MumpsParallelSOE.hGo 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.2 $ 00022 // $Date: 2006/04/13 20:58:07 $ 00023 // $Source: /usr/local/cvs/OpenSees/SRC/system_of_eqn/linearSOE/mumps/MumpsParallelSOE.h,v $ 00024 00025 #ifndef MumpsParallelSOE_h 00026 #define MumpsParallelSOE_h 00027 00028 // Written: fmk 00029 // Description: This file contains the class definition for MumpsParallelSOE 00030 // MumpsParallelSOE is a subclass of LinearSOE. It uses the sparse column 00031 // storage scheme. The matrix A is kept distributed, X and B kept on all processors. 00032 // 00033 // matrix types (matType): 0 Unsymmetrc 00034 // 1 Symmetrix positive definite 00035 // 2 General Symmetric 00036 00037 // What: "@(#) MumpsParallelSOE.h, revA" 00038 00039 00040 #include <MumpsSOE.h> 00041 #include <Vector.h> 00042 00043 class MumpsParallelSolver; 00044 00045 class MumpsParallelSOE : public MumpsSOE 00046 { 00047 public: 00048 MumpsParallelSOE(MumpsParallelSolver &theSolver, int matType=2); 00049 00050 ~MumpsParallelSOE(); 00051 00052 // these methods need to be rewritten 00053 int setSize(Graph &theGraph); 00054 00055 int addB(const Vector &, const ID &, double fact = 1.0); 00056 int setB(const Vector &, double fact = 1.0); 00057 const Vector &getB(void); 00058 void zeroB(void); 00059 int solve(void); 00060 00061 int sendSelf(int commitTag, Channel &theChannel); 00062 int recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker); 00063 friend class MumpsParallelSolver; 00064 00065 int setProcessID(int processTag); 00066 int setChannels(int numChannels, Channel **theChannels); 00067 00068 protected: 00069 00070 private: 00071 int processID; 00072 int numChannels; 00073 Channel **theChannels; 00074 ID **localCol; 00075 00076 double *workArea; 00077 int sizeWork; 00078 double *myB; 00079 Vector *myVectB; 00080 }; 00081 00082 00083 #endif 00084 |