DistributedSuperLU.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.1 $ 00022 // $Date: 2005/12/06 22:21:03 $ 00023 // $Source: /usr/local/cvs/OpenSees/SRC/system_of_eqn/linearSOE/sparseGEN/DistributedSuperLU.h,v $ 00024 00025 00026 #ifndef DistributedSuperLU_h 00027 #define DistributedSuperLU_h 00028 00029 // Written: fmk 00030 // 00031 // Description: This file contains the class definition for DistributedSuperLU. 00032 // A DistributedSuperLU object can be constructed to solve a SparseGenColLinSOE 00033 // object. It obtains the solution by making calls on the 00034 // the SuperLU library developed at UC Berkeley by Prof. James Demmel, 00035 // Xiaoye S. Li and John R. Gilbert. 00036 // The SuperLU library contains a set of subroutines to solve a sparse 00037 // linear system $AX=B$. It uses Gaussian elimination with partial 00038 // pivoting (GEPP). The columns of A may be preordered before 00039 // factorization; the preordering for sparsity is completely separate 00040 // from the factorization and a number of ordering schemes are provided. 00041 // 00042 // What: "@(#) DistributedSuperLU.h, revA" 00043 00044 #include <superlu_ddefs.h> 00045 #include <SparseGenColLinSolver.h> 00046 00047 00048 class DistributedSuperLU : public SparseGenColLinSolver 00049 { 00050 public: 00051 DistributedSuperLU(int npRow, int npCol); 00052 DistributedSuperLU(); 00053 ~DistributedSuperLU(); 00054 00055 int solve(void); 00056 int setSize(void); 00057 00058 virtual int setProcessID(int domainTag); 00059 virtual int setChannels(int numChannels, Channel **theChannels); 00060 int sendSelf(int commitTag, Channel &theChannel); 00061 int recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker); 00062 00063 protected: 00064 00065 private: 00066 superlu_options_t_Distributed options; 00067 SuperLUStat_t stat; 00068 SuperMatrix A; 00069 ScalePermstruct_t ScalePermstruct; 00070 LUstruct_t LUstruct; 00071 gridinfo_t grid; 00072 00073 bool gridInit; 00074 int npRow, npCol; 00075 00076 int processID; 00077 int numChannels; 00078 Channel **theChannels; 00079 00080 MPI_Comm comm_SuperLU; 00081 double *b; 00082 int *rowA; 00083 }; 00084 00085 #endif 00086 |