OpenSees » SRC » analysis » algorithm » equiSolnAlgo

KrylovNewton.h

Go 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 1999, 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 ** ****************************************************************** */
00020                                                                         
00021 // $Revision: 1.8 $
00022 // $Date: 2005/11/29 22:42:41 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/algorithm/equiSolnAlgo/KrylovNewton.h,v $
00024                                                                         
00025 #ifndef KrylovNewton_h
00026 #define KrylovNewton_h
00027 
00028 // Written: MHS
00029 // Created: June 2001
00030 //
00031 // Description: This file contains the class definition for 
00032 // KrylovNewton.  KrylovNewton is a class which uses a Krylov
00033 // subspace accelerator on the modified Newton method.
00034 // The accelerator is described by Carlson and Miller in
00035 // "Design and Application of a 1D GWMFE Code"
00036 // from SIAM Journal of Scientific Computing (Vol. 19, No. 3,
00037 // pp. 728-765, May 1998)
00038 
00039 #include <EquiSolnAlgo.h>
00040 #include <Vector.h>
00041 
00042 class KrylovNewton: public EquiSolnAlgo
00043 {
00044   public:
00045     KrylovNewton(int tangent = CURRENT_TANGENT, int maxDim = 3);    
00046     KrylovNewton(ConvergenceTest &theTest, int tangent = CURRENT_TANGENT, int maxDim = 3);
00047     ~KrylovNewton();
00048 
00049     int solveCurrentStep(void);    
00050     int setConvergenceTest(ConvergenceTest *theNewTest);
00051     ConvergenceTest *getConvergenceTest(void);     
00052     
00053     virtual int sendSelf(int commitTag, Channel &theChannel);
00054     virtual int recvSelf(int commitTag, Channel &theChannel, 
00055                          FEM_ObjectBroker &theBroker);
00056     void Print(OPS_Stream &s, int flag =0);    
00057     
00058   protected:
00059     
00060   private:
00061     ConvergenceTest *theTest;
00062     int tangent;
00063 
00064     // Storage for update vectors
00065     Vector **v;
00066 
00067     // Storage for subspace vectors
00068     Vector **Av;
00069 
00070     // Array data sent to LAPACK subroutine
00071     double *AvData;
00072     double *rData;
00073     double *work;
00074 
00075     // Length of work array
00076     int lwork;
00077 
00078     // Size information
00079     int numEqns;
00080     int maxDimension;
00081 
00082     // Private lsq routine to do Krylov updates
00083     // dimension is the current dimension of the subspace
00084     int leastSquares(int dimension);
00085 };
00086 
00087 #endif
00088 
00089

Generated on Mon Oct 23 15:04:57 2006 for OpenSees by doxygen 1.5.0