Nektar++
GlobalLinSysStaticCond.h
Go to the documentation of this file.
1 ///////////////////////////////////////////////////////////////////////////////
2 //
3 // File: GlobalLinSysStaticCond.h
4 //
5 // For more information, please see: http://www.nektar.info
6 //
7 // The MIT License
8 //
9 // Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
10 // Department of Aeronautics, Imperial College London (UK), and Scientific
11 // Computing and Imaging Institute, University of Utah (USA).
12 //
13 // License for the specific language governing rights and limitations under
14 // Permission is hereby granted, free of charge, to any person obtaining a
15 // copy of this software and associated documentation files (the "Software"),
16 // to deal in the Software without restriction, including without limitation
17 // the rights to use, copy, modify, merge, publish, distribute, sublicense,
18 // and/or sell copies of the Software, and to permit persons to whom the
19 // Software is furnished to do so, subject to the following conditions:
20 //
21 // The above copyright notice and this permission notice shall be included
22 // in all copies or substantial portions of the Software.
23 //
24 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
25 // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
26 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
27 // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
28 // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
29 // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
30 // DEALINGS IN THE SOFTWARE.
31 //
32 // Description: A collection of routines common to statically condensed systems.
33 //
34 ///////////////////////////////////////////////////////////////////////////////
35 
36 #ifndef NEKTAR_LIB_MULTIREGIONS_GLOBALLINSYSSTATICCOND_H
37 #define NEKTAR_LIB_MULTIREGIONS_GLOBALLINSYSSTATICCOND_H
38 
42 
43 namespace Nektar
44 {
45  namespace MultiRegions
46  {
47  // Forward declarations
48  class ExpList;
50 
51  typedef boost::shared_ptr<GlobalLinSysStaticCond>
53 
54  /// A global linear system.
55  class GlobalLinSysStaticCond : virtual public GlobalLinSys
56  {
57  public:
58  /// Constructor for full direct matrix solve.
60  const GlobalLinSysKey &mkey,
61  const boost::weak_ptr<ExpList> &pExpList,
62  const boost::shared_ptr<AssemblyMap> &locToGloMap);
63 
64  virtual ~GlobalLinSysStaticCond();
65 
66  protected:
68  int scLevel,
69  NekVector<NekDouble> &F_GlobBnd)
70  {
71  return m_schurCompl;
72  }
73 
74  virtual void v_BasisTransform(
75  Array<OneD, NekDouble>& pInOut,
76  int offset)
77  {
78 
79  }
80 
81  virtual void v_BasisInvTransform(
82  Array<OneD, NekDouble>& pInOut)
83  {
84 
85  }
86 
88  boost::shared_ptr<AssemblyMap> pLoctoGloMap)
89  {
90 
91  }
92 
93  virtual int v_GetNumBlocks();
94 
95  virtual GlobalLinSysStaticCondSharedPtr v_Recurse(
96  const GlobalLinSysKey &mkey,
97  const boost::weak_ptr<ExpList> &pExpList,
98  const DNekScalBlkMatSharedPtr pSchurCompl,
99  const DNekScalBlkMatSharedPtr pBinvD,
100  const DNekScalBlkMatSharedPtr pC,
101  const DNekScalBlkMatSharedPtr pInvD,
102  const boost::shared_ptr<AssemblyMap> &locToGloMap) = 0;
103 
104  /// Schur complement for Direct Static Condensation.
105  GlobalLinSysStaticCondSharedPtr m_recursiveSchurCompl;
106  /// Block Schur complement matrix.
108  /// Block \f$ BD^{-1} \f$ matrix.
110  /// Block \f$ C \f$ matrix.
112  /// Block \f$ D^{-1} \f$ matrix.
114  /// Local to global map.
115  boost::shared_ptr<AssemblyMap> m_locToGloMap;
116  /// Workspace array for matrix multiplication
118 
119  /// Solve the linear system for given input and output vectors
120  /// using a specified local to global map.
121  virtual void v_Solve(
124  const AssemblyMapSharedPtr &locToGloMap,
125  const Array<OneD, const NekDouble> &dirForcing
127 
128  virtual void v_InitObject();
129 
130  /// Initialise this object
131  virtual void v_Initialise(
132  const boost::shared_ptr<AssemblyMap>& locToGloMap);
133 
134  /// Set up the storage for the Schur complement or the top level
135  /// of the multi-level Schur complement.
136  void SetupTopLevel(
137  const boost::shared_ptr<AssemblyMap>& locToGloMap);
138 
139  ///
141  const boost::shared_ptr<AssemblyMap>& locToGloMap);
142  };
143  }
144 }
145 
146 #endif
Array< OneD, NekDouble > m_wsp
Workspace array for matrix multiplication.
virtual void v_Solve(const Array< OneD, const NekDouble > &in, Array< OneD, NekDouble > &out, const AssemblyMapSharedPtr &locToGloMap, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
Solve the linear system for given input and output vectors using a specified local to global map...
virtual GlobalLinSysStaticCondSharedPtr v_Recurse(const GlobalLinSysKey &mkey, const boost::weak_ptr< ExpList > &pExpList, const DNekScalBlkMatSharedPtr pSchurCompl, const DNekScalBlkMatSharedPtr pBinvD, const DNekScalBlkMatSharedPtr pC, const DNekScalBlkMatSharedPtr pInvD, const boost::shared_ptr< AssemblyMap > &locToGloMap)=0
static Array< OneD, NekDouble > NullNekDouble1DArray
boost::shared_ptr< AssemblyMap > AssemblyMapSharedPtr
Definition: AssemblyMap.h:53
DNekScalBlkMatSharedPtr m_schurCompl
Block Schur complement matrix.
GlobalLinSysStaticCondSharedPtr m_recursiveSchurCompl
Schur complement for Direct Static Condensation.
DNekScalBlkMatSharedPtr m_C
Block matrix.
virtual void v_BasisTransform(Array< OneD, NekDouble > &pInOut, int offset)
virtual void v_Initialise(const boost::shared_ptr< AssemblyMap > &locToGloMap)
Initialise this object.
DNekScalBlkMatSharedPtr m_invD
Block matrix.
virtual void v_AssembleSchurComplement(boost::shared_ptr< AssemblyMap > pLoctoGloMap)
virtual DNekScalBlkMatSharedPtr v_PreSolve(int scLevel, NekVector< NekDouble > &F_GlobBnd)
boost::shared_ptr< DNekScalBlkMat > DNekScalBlkMatSharedPtr
Definition: NekTypeDefs.hpp:74
void SetupTopLevel(const boost::shared_ptr< AssemblyMap > &locToGloMap)
Set up the storage for the Schur complement or the top level of the multi-level Schur complement...
Describe a linear system.
virtual int v_GetNumBlocks()
Get the number of blocks in this system.
A global linear system.
Definition: GlobalLinSys.h:70
GlobalLinSysStaticCond(const GlobalLinSysKey &mkey, const boost::weak_ptr< ExpList > &pExpList, const boost::shared_ptr< AssemblyMap > &locToGloMap)
Constructor for full direct matrix solve.
boost::shared_ptr< AssemblyMap > m_locToGloMap
Local to global map.
virtual void v_BasisInvTransform(Array< OneD, NekDouble > &pInOut)
DNekScalBlkMatSharedPtr m_BinvD
Block matrix.
void ConstructNextLevelCondensedSystem(const boost::shared_ptr< AssemblyMap > &locToGloMap)
boost::shared_ptr< GlobalLinSysStaticCond > GlobalLinSysStaticCondSharedPtr