xc
ItpackLinSolver.h
1 // -*-c++-*-
2 //----------------------------------------------------------------------------
3 // XC program; finite element analysis code
4 // for structural analysis and design.
5 //
6 // Copyright (C) Luis C. Pérez Tato
7 //
8 // This program derives from OpenSees <http://opensees.berkeley.edu>
9 // developed by the «Pacific earthquake engineering research center».
10 //
11 // Except for the restrictions that may arise from the copyright
12 // of the original program (see copyright_opensees.txt)
13 // XC is free software: you can redistribute it and/or modify
14 // it under the terms of the GNU General Public License as published by
15 // the Free Software Foundation, either version 3 of the License, or
16 // (at your option) any later version.
17 //
18 // This software is distributed in the hope that it will be useful, but
19 // WITHOUT ANY WARRANTY; without even the implied warranty of
20 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
21 // GNU General Public License for more details.
22 //
23 //
24 // You should have received a copy of the GNU General Public License
25 // along with this program.
26 // If not, see <http://www.gnu.org/licenses/>.
27 //----------------------------------------------------------------------------
28 /* ****************************************************************** **
29 ** OpenSees - Open System for Earthquake Engineering Simulation **
30 ** Pacific Earthquake Engineering Research Center **
31 ** **
32 ** **
33 ** (C) Copyright 1999, The Regents of the University of California **
34 ** All Rights Reserved. **
35 ** **
36 ** Commercial use of this program without express permission of the **
37 ** University of California, Berkeley, is strictly prohibited. See **
38 ** file 'COPYRIGHT' in main directory for information on usage and **
39 ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. **
40 ** **
41 ** Developed by: **
42 ** Frank McKenna (fmckenna@ce.berkeley.edu) **
43 ** Gregory L. Fenves (fenves@ce.berkeley.edu) **
44 ** Filip C. Filippou (filippou@ce.berkeley.edu) **
45 ** **
46 ** ****************************************************************** */
47 
48 // $Revision: 1.2 $
49 // $Date: 2002/06/08 16:17:26 $
50 // $Source: /usr/local/cvs/OpenSees/SRC/system_of_eqn/linearSOE/itpack/ItpackLinSolver.h,v $
51 
52 #ifndef ItpackLinSolver_h
53 #define ItpackLinSolver_h
54 
55 // Written: MHS
56 // Created: Sept 2001
57 //
58 // Description: This file contains the class definition for ItpackLinSolver.
59 // ItpackLinSolver is a concrete subclass of LinearSOE. It stores full
60 // unsymmetric linear system of equations using 1d arrays in Fortran style
61 
62 #include <solution/system_of_eqn/linearSOE/LinearSOESolver.h>
63 
64 // Adaptive methods
65 #define ItpackJCG 1
66 #define ItpackJSI 2
67 #define ItpackSOR 3
68 #define ItpackSSORCG 4
69 #define ItpackSSORSI 5
70 #define ItpackRSCG 6
71 #define ItpackRSSI 7
72 
73 // Textbook methods
74 #define ItpackJ 8
75 #define ItpackGS 9
76 #define ItpackSORFixed 10
77 #define ItpackSSORFixed 11
78 #define ItpackRS 12
79 
80 namespace XC {
81 class ItpackLinSOE;
82 
84 //
87  {
88  private:
89 
90  ItpackLinSOE *theSOE;
91 
92  ID IA;
93  ID JA;
94 
95  int n; // Size of system, i.e., the number of equations
96  // Parameter arrays sent to ITPACK subroutines
97  int iparm[12];
98  double rparm[12];
99  // Workspace arrays sent to ITPACK subroutines
100  ID iwksp;
101  Vector wksp;
102 
103  int nwksp;
104  int maxIter;
105  int method;
106  double omega;
107 
108  ItpackLinSolver(int method, int maxIter = 100, double omega = 1.0);
109  ItpackLinSolver(void);
110  virtual LinearSOESolver *getCopy(void) const;
111  virtual bool setLinearSOE(LinearSOE *theSOE);
112  public:
113 
114  int solve(void);
115  int setSize(void);
116  bool setLinearSOE(ItpackLinSOE &theSOE);
117 
118  int sendSelf(Communicator &);
119  int recvSelf(const Communicator &);
120  };
121 
122 inline LinearSOESolver *ItpackLinSolver::getCopy(void) const
123  { return new ItpackLinSolver(*this); }
124 } // end of XC namespace
125 
126 #endif
Float vector abstraction.
Definition: Vector.h:94
Communication parameters between processes.
Definition: Communicator.h:66
Linear system of equations.
Definition: LinearSOE.h:92
Vector of integers.
Definition: ID.h:95
int sendSelf(Communicator &)
Send the object.
Definition: ItpackLinSolver.cpp:168
ITPACK based systems of equations.
Definition: ItpackLinSOE.h:78
int solve(void)
Causes the solver to compute the solution of the system of equations.
Definition: ItpackLinSolver.cpp:206
int setSize(void)
Sets the number of equations.
Definition: ItpackLinSolver.cpp:93
Open source finite element program for structural analysis.
Definition: ContinuaReprComponent.h:35
Base class of the solvers for linear system of equations.
Definition: LinearSOESolver.h:87
ITPACK based linear SOE solver.
Definition: ItpackLinSolver.h:86
int recvSelf(const Communicator &)
Receive the object.
Definition: ItpackLinSolver.cpp:171