ModNewtonRootFinding.cppGo to the documentation of this file.00001 /* ****************************************************************** ** 00002 ** OpenSees - Open System for Earthquake Engineering Simulation ** 00003 ** Pacific Earthquake Engineering Research Center ** 00004 ** ** 00005 ** ** 00006 ** (C) Copyright 2001, The Regents of the University of California ** 00007 ** All Rights Reserved. ** 00008 ** ** 00009 ** Commercial use of this program without express permission of the ** 00010 ** University of California, Berkeley, is strictly prohibited. See ** 00011 ** file 'COPYRIGHT' in main directory for information on usage and ** 00012 ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. ** 00013 ** ** 00014 ** Developed by: ** 00015 ** Frank McKenna (fmckenna@ce.berkeley.edu) ** 00016 ** Gregory L. Fenves (fenves@ce.berkeley.edu) ** 00017 ** Filip C. Filippou (filippou@ce.berkeley.edu) ** 00018 ** ** 00019 ** Reliability module developed by: ** 00020 ** Terje Haukaas (haukaas@ce.berkeley.edu) ** 00021 ** Armen Der Kiureghian (adk@ce.berkeley.edu) ** 00022 ** ** 00023 ** ****************************************************************** */ 00024 00025 // $Revision: 1.1 $ 00026 // $Date: 2003/10/27 23:45:44 $ 00027 // $Source: /usr/local/cvs/OpenSees/SRC/reliability/analysis/rootFinding/ModNewtonRootFinding.cpp,v $ 00028 00029 00030 // 00031 // Written by Terje Haukaas (haukaas@ce.berkeley.edu) 00032 // 00033 00034 #include <ModNewtonRootFinding.h> 00035 #include <RootFinding.h> 00036 #include <GFunEvaluator.h> 00037 #include <ProbabilityTransformation.h> 00038 #include <ReliabilityDomain.h> 00039 #include <RandomVariable.h> 00040 #include <math.h> 00041 #include <Vector.h> 00042 00043 00044 ModNewtonRootFinding::ModNewtonRootFinding( 00045 ReliabilityDomain *passedReliabilityDomain, 00046 ProbabilityTransformation *passedProbabilityTransformation, 00047 GFunEvaluator *passedGFunEvaluator, 00048 int passedMaxIter, 00049 double ptol, 00050 double pmaxStepLength) 00051 :RootFinding() 00052 { 00053 theReliabilityDomain = passedReliabilityDomain; 00054 theProbabilityTransformation = passedProbabilityTransformation; 00055 theGFunEvaluator = passedGFunEvaluator; 00056 maxIter = passedMaxIter; 00057 tol = ptol; 00058 maxStepLength = pmaxStepLength; 00059 } 00060 00061 ModNewtonRootFinding::~ModNewtonRootFinding() 00062 { 00063 } 00064 00065 00066 Vector 00067 ModNewtonRootFinding::findLimitStateSurface(int space, double g, Vector pDirection, Vector thePoint) 00068 { 00069 opserr << "Currently, the Modified Newton root-finding algorithm is not " << endln 00070 << "implemented. This algorithm requires the directional gradient in " << endln 00071 << "the direction of the root-finding search. For now; it is " << endln 00072 << "recommended to use the Secant Root-Finding Algorithm." << endln; 00073 return thePoint; 00074 } 00075 00076 |