xc
ProfileSPDLinSOE.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.3 $
49 // $Date: 2005/12/06 22:06:04 $
50 // $Source: /usr/local/cvs/OpenSees/SRC/system_of_eqn/linearSOE/profileSPD/ProfileSPDLinSOE.h,v $
51 
52 
53 // Written: fmk
54 // Revision: A
55 //
56 // Description: This file contains the class definition for ProfileSPDLinSOE
57 // ProfileSPDLinSOE is a subclass of LinearSOE. It uses the LAPACK Upper storage
58 // scheme to store the components of the A matrix.
59 
60 // What: "@(#) ProfileSPDLinSOE.h, revA"
61 
62 #ifndef ProfileSPDLinSOE_h
63 #define ProfileSPDLinSOE_h
64 
65 #include <solution/system_of_eqn/linearSOE/FactoredSOEBase.h>
66 #include "utility/matrix/Vector.h"
67 #include "utility/matrix/ID.h"
68 
69 namespace XC {
70 class ProfileSPDLinSolver;
71 
121  {
122  protected:
123  int profileSize;
124  Vector A;
125  ID iDiagLoc;
126  bool isAcondensed;
127  int numInt;
128  protected:
129  virtual bool setSolver(LinearSOESolver *);
130 
131  friend class SolutionStrategy;
132  friend class FEM_ObjectBroker;
134  ProfileSPDLinSOE(SolutionStrategy *,int classTag);
135  ProfileSPDLinSOE(SolutionStrategy *,int N, int *iLoc, ProfileSPDLinSolver *theSolver= nullptr);
136  SystemOfEqn *getCopy(void) const;
137  public:
138  virtual int setSize(Graph &theGraph);
139  virtual int addA(const Matrix &, const ID &, double fact = 1.0);
140 
141  virtual void zeroA(void);
142 
143  virtual int sendSelf(Communicator &);
144  virtual int recvSelf(const Communicator &);
145 
146  friend class ProfileSPDLinSolver;
147  friend class ProfileSPDLinDirectSolver;
148  friend class ProfileSPDLinDirectBlockSolver;
149  friend class ProfileSPDLinDirectThreadSolver;
150  friend class ProfileSPDLinDirectSkypackSolver;
151  friend class ProfileSPDLinSubstrSolver;
152  friend class ProfileSPDLinSubstrThreadSolver;
153  };
155  { return new ProfileSPDLinSOE(*this); }
156 } // end of XC namespace
157 
158 
159 #endif
160 
161 
162 
Skypakc based linear SOE solver.
Definition: ProfileSPDLinDirectSkypackSolver.h:90
virtual int sendSelf(Communicator &)
Send the object.
Definition: ProfileSPDLinSOE.cpp:338
Float vector abstraction.
Definition: Vector.h:94
It solves a ProfileSPDLinSOE object using the LDL^t factorization.
Definition: ProfileSPDLinDirectSolver.h:84
Communication parameters between processes.
Definition: Communicator.h:66
SystemOfEqn * getCopy(void) const
Virtual constructor.
Definition: ProfileSPDLinSOE.h:154
FEM_ObjectBroker is is an object broker class for the finite element method.
Definition: FEM_ObjectBroker.h:151
Vector of integers.
Definition: ID.h:95
solves a ProfileSPDLinSOE object using the LDL^t factorization.
Definition: ProfileSPDLinDirectBlockSolver.h:85
System of equations base class.
Definition: SystemOfEqn.h:90
virtual int addA(const Matrix &, const ID &, double fact=1.0)
Assembles the product of m by fact into A.
Definition: ProfileSPDLinSOE.cpp:258
virtual bool setSolver(LinearSOESolver *)
Set the solver to use.
Definition: ProfileSPDLinSOE.cpp:128
The Graph class provides the abstraction of a graph.
Definition: Graph.h:94
Base class for factored systems of equations.
Definition: FactoredSOEBase.h:41
A ProfileSPDLinSubstrSolver object will perform the numerical substructuring operations on a ProfileS...
Definition: ProfileSPDLinSubstrSolver.h:81
virtual int recvSelf(const Communicator &)
Receive the object.
Definition: ProfileSPDLinSOE.cpp:342
Solution strategy for the finite element problem.
Definition: SolutionStrategy.h:94
Profile matrix system of equations.
Definition: ProfileSPDLinSOE.h:120
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
Base class for profile matrix linear SOE solvers.
Definition: ProfileSPDLinSolver.h:77
Matrix of floats.
Definition: Matrix.h:111
solves a ProfileSPDLinSOE object using the LDL^t factorization (threaded version).
Definition: ProfileSPDLinDirectThreadSolver.h:86
virtual void zeroA(void)
Zeros the entries in the 1d array for and marks the system as not having been factored.
Definition: ProfileSPDLinSOE.cpp:331
ProfileSPDLinSOE(SolutionStrategy *)
Constructor.
Definition: ProfileSPDLinSOE.cpp:69
virtual int setSize(Graph &theGraph)
Sets the system size.
Definition: ProfileSPDLinSOE.cpp:161