SuperLU.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.3 $ 00022 // $Date: 2005/03/17 20:47:09 $ 00023 // $Source: /usr/local/cvs/OpenSees/SRC/system_of_eqn/linearSOE/sparseGEN/SuperLU.h,v $ 00024 00025 00026 #ifndef SuperLU_h 00027 #define SuperLU_h 00028 00029 // File: ~/system_of_eqn/linearSOE/sparseGEN/SuperLU.h 00030 // 00031 // Written: fmk 00032 // Created: 11/96 00033 // 00034 // Description: This file contains the class definition for SuperLU. 00035 // A SuperLU object can be constructed to solve a SparseGenColLinSOE 00036 // object. It obtains the solution by making calls on the 00037 // the SuperLU library developed at UC Berkeley by Prof. James Demmel, 00038 // Xiaoye S. Li and John R. Gilbert. 00039 // The SuperLU library contains a set of subroutines to solve a sparse 00040 // linear system $AX=B$. It uses Gaussian elimination with partial 00041 // pivoting (GEPP). The columns of A may be preordered before 00042 // factorization; the preordering for sparsity is completely separate 00043 // from the factorization and a number of ordering schemes are provided. 00044 // 00045 // What: "@(#) SuperLU.h, revA" 00046 00047 #include <SparseGenColLinSolver.h> 00048 #include <dsp_defs.h> 00049 #include <supermatrix.h> 00050 00051 class SuperLU : public SparseGenColLinSolver 00052 { 00053 public: 00054 SuperLU(int permSpec = 0, 00055 double drop_tol = 0.0, 00056 int panelSize = 6, 00057 int relax = 6, 00058 char symmetric = 'N'); 00059 ~SuperLU(); 00060 00061 int solve(void); 00062 int setSize(void); 00063 00064 int sendSelf(int commitTag, Channel &theChannel); 00065 int recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker); 00066 00067 protected: 00068 00069 private: 00070 SuperMatrix A,L,U,B,AC; 00071 int *perm_r; 00072 int *perm_c; 00073 int *etree; 00074 int sizePerm; 00075 int relax, permSpec, panelSize; 00076 double drop_tol; 00077 char symmetric; 00078 superlu_options_t options; 00079 SuperLUStat_t stat; 00080 }; 00081 00082 #endif 00083 |