ThreadedSuperLU.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: 2002/01/25 20:27:11 $ 00023 // $Source: /usr/local/cvs/OpenSees/SRC/system_of_eqn/linearSOE/sparseGEN/ThreadedSuperLU.h,v $ 00024 00025 00026 #ifndef ThreadedSuperLU_h 00027 #define ThreadedSuperLU_h 00028 00029 // File: ~/system_of_eqn/linearSOE/sparseGEN/ThreadedSuperLU.h 00030 // 00031 // Written: fmk 00032 // Created: 11/96 00033 // Revision: A 00034 // 00035 // Description: This file contains the class interface for ThreadedSuperLU. 00036 // This is a class that uses the threads version of SuperLU 00037 // 00038 // What: "@(#) ThreadedSuperLU.h, revA" 00039 00040 #include <SparseGenColLinSolver.h> 00041 #include <supermatrix.h> 00042 #include <pdsp_defs.h> 00043 00044 class ThreadedSuperLU : public SparseGenColLinSolver 00045 { 00046 public: 00047 ThreadedSuperLU(int numThreads = 2, 00048 int permSpec = 0, int panelSize = 6, 00049 int relax = 6, double thresh = 0.0); 00050 00051 ~ThreadedSuperLU(); 00052 00053 int solve(void); 00054 int setSize(void); 00055 00056 int sendSelf(int commitTag, Channel &theChannel); 00057 int recvSelf(int commitTag, Channel &theChannel, 00058 FEM_ObjectBroker &theBroker); 00059 00060 protected: 00061 00062 private: 00063 SuperMatrix A,L,U,B,AC; 00064 int *perm_r; 00065 int *perm_c; 00066 int *etree; 00067 int sizePerm; 00068 int relax, permSpec, panelSize; 00069 float thresh; 00070 00071 int numThreads; 00072 pdgstrf_options_t pdgstrf_options; 00073 yes_no_t refact, usepr; 00074 fact_t fact; 00075 trans_t trans; 00076 void *work; 00077 int lwork; 00078 Gstat_t gStat; 00079 }; 00080 00081 #endif 00082 00083 |