escript
Revision_
|
00001 00002 /***************************************************************************** 00003 * 00004 * Copyright (c) 2003-2014 by University of Queensland 00005 * http://www.uq.edu.au 00006 * 00007 * Primary Business: Queensland, Australia 00008 * Licensed under the Open Software License version 3.0 00009 * http://www.opensource.org/licenses/osl-3.0.php 00010 * 00011 * Development until 2012 by Earth Systems Science Computational Center (ESSCC) 00012 * Development 2012-2013 by School of Earth Sciences 00013 * Development from 2014 by Centre for Geoscience Computing (GeoComp) 00014 * 00015 *****************************************************************************/ 00016 #ifndef __RIPLEY_DEFAULTASSEMBLER3D_H__ 00017 #define __RIPLEY_DEFAULTASSEMBLER3D_H__ 00018 00019 #include <map> 00020 #include <escript/Data.h> 00021 #include <ripley/Ripley.h> 00022 #include <ripley/RipleyException.h> 00023 #include <ripley/AbstractAssembler.h> 00024 #include <ripley/Brick.h> 00025 00026 namespace ripley { 00027 00028 00029 class DefaultAssembler3D : public AbstractAssembler { 00030 public: 00031 DefaultAssembler3D(Brick *dom, double *m_dx, dim_t *m_NX, dim_t *m_NE, 00032 dim_t *m_NN) : AbstractAssembler() { 00033 domain = dom; 00034 this->m_dx = m_dx; 00035 this->m_NX = m_NX; 00036 this->m_NE = m_NE; 00037 this->m_NN = m_NN; 00038 } 00039 ~DefaultAssembler3D(){}; 00040 /* The default RipleyDomain assemblers, with original signatures */ 00041 00044 virtual void assemblePDESingle(paso::SystemMatrix_ptr mat, 00045 escript::Data& rhs, 00046 const escript::Data& A, const escript::Data& B, 00047 const escript::Data& C, const escript::Data& D, 00048 const escript::Data& X, const escript::Data& Y) const; 00049 00052 virtual void assemblePDEBoundarySingle(paso::SystemMatrix_ptr mat, 00053 escript::Data& rhs, const escript::Data& d, 00054 const escript::Data& y) const; 00055 00058 virtual void assemblePDESingleReduced(paso::SystemMatrix_ptr mat, 00059 escript::Data& rhs, const escript::Data& A, const escript::Data& B, 00060 const escript::Data& C, const escript::Data& D, 00061 const escript::Data& X, const escript::Data& Y) const; 00062 00065 virtual void assemblePDEBoundarySingleReduced(paso::SystemMatrix_ptr mat, 00066 escript::Data& rhs, const escript::Data& d, 00067 const escript::Data& y) const; 00068 00071 virtual void assemblePDESystem(paso::SystemMatrix_ptr mat, escript::Data& rhs, 00072 const escript::Data& A, const escript::Data& B, 00073 const escript::Data& C, const escript::Data& D, 00074 const escript::Data& X, const escript::Data& Y) const; 00075 00078 virtual void assemblePDEBoundarySystem(paso::SystemMatrix_ptr mat, 00079 escript::Data& rhs, const escript::Data& d, 00080 const escript::Data& y) const; 00081 00084 virtual void assemblePDESystemReduced(paso::SystemMatrix_ptr mat, 00085 escript::Data& rhs, const escript::Data& A, const escript::Data& B, 00086 const escript::Data& C, const escript::Data& D, 00087 const escript::Data& X, const escript::Data& Y) const; 00088 00091 virtual void assemblePDEBoundarySystemReduced(paso::SystemMatrix_ptr mat, 00092 escript::Data& rhs, const escript::Data& d, 00093 const escript::Data& y) const; 00094 00095 /* The new interface for assemblers */ 00096 virtual void assemblePDESingle(paso::SystemMatrix_ptr mat, escript::Data& rhs, 00097 std::map<std::string, escript::Data> coefs) const; 00098 virtual void assemblePDEBoundarySingle(paso::SystemMatrix_ptr mat, 00099 escript::Data& rhs, 00100 std::map<std::string, escript::Data> coefs) const; 00101 virtual void assemblePDESingleReduced(paso::SystemMatrix_ptr mat, 00102 escript::Data& rhs, 00103 std::map<std::string, escript::Data> coefs) const; 00104 virtual void assemblePDEBoundarySingleReduced(paso::SystemMatrix_ptr mat, 00105 escript::Data& rhs, 00106 std::map<std::string, escript::Data> coefs) const; 00107 virtual void assemblePDESystem(paso::SystemMatrix_ptr mat, escript::Data& rhs, 00108 std::map<std::string, escript::Data> coefs) const; 00109 virtual void assemblePDEBoundarySystem(paso::SystemMatrix_ptr mat, 00110 escript::Data& rhs, 00111 std::map<std::string, escript::Data> coefs) const; 00112 virtual void assemblePDESystemReduced(paso::SystemMatrix_ptr mat, 00113 escript::Data& rhs, 00114 std::map<std::string, escript::Data> coefs) const; 00115 virtual void assemblePDEBoundarySystemReduced(paso::SystemMatrix_ptr mat, 00116 escript::Data& rhs, 00117 std::map<std::string, escript::Data> coefs) const; 00118 00119 void collateFunctionSpaceTypes(std::vector<int>& fsTypes, 00120 std::map<std::string, escript::Data> coefs) const; 00121 protected: 00122 const Brick *domain; 00123 const double *m_dx; 00124 const dim_t *m_NX; 00125 const dim_t *m_NE; 00126 const dim_t *m_NN; 00127 }; 00128 00129 } // namespace ripley 00130 00131 #endif // __RIPLEY_DEFAULTASSEMBLER3D_H__ 00132