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_WAVEASSEMBLER2D_H__ 00017 #define __RIPLEY_WAVEASSEMBLER2D_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 WaveAssembler3D : public AbstractAssembler { 00030 public: 00031 WaveAssembler3D(Brick *dom, double *m_dx, dim_t *m_NX, dim_t *m_NE, 00032 dim_t *m_NN, std::map<std::string, escript::Data> c); 00033 ~WaveAssembler3D(){}; 00034 00035 /* The only assembly function we care about right now*/ 00036 void assemblePDESystem(paso::SystemMatrix_ptr mat, escript::Data& rhs, 00037 std::map<std::string, escript::Data> coefs) const; 00038 00039 00040 void assemblePDESingle(paso::SystemMatrix_ptr mat, escript::Data& rhs, 00041 std::map<std::string, escript::Data> coefs) const { 00042 throw RipleyException("This assembly not supported by this assembler"); 00043 } 00044 void assemblePDEBoundarySingle(paso::SystemMatrix_ptr mat, 00045 escript::Data& rhs, 00046 std::map<std::string, escript::Data> coefs) const { 00047 throw RipleyException("This assembly not supported by this assembler"); 00048 } 00049 void assemblePDESingleReduced(paso::SystemMatrix_ptr mat, 00050 escript::Data& rhs, 00051 std::map<std::string, escript::Data> coefs) const { 00052 throw RipleyException("This assembly not supported by this assembler"); 00053 } 00054 void assemblePDEBoundarySingleReduced(paso::SystemMatrix_ptr mat, 00055 escript::Data& rhs, 00056 std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");} 00057 void assemblePDEBoundarySystem(paso::SystemMatrix_ptr mat, escript::Data& rhs, 00058 std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");} 00059 void assemblePDESystemReduced(paso::SystemMatrix_ptr mat, escript::Data& rhs, 00060 std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");} 00061 void assemblePDEBoundarySystemReduced(paso::SystemMatrix_ptr mat, 00062 escript::Data& rhs, 00063 std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");} 00064 00065 void collateFunctionSpaceTypes(std::vector<int>& fsTypes, 00066 std::map<std::string, escript::Data> coefs) const; 00067 00068 private: 00069 std::map<std::string, escript::Data> c; 00070 Brick *domain; 00071 double *m_dx; 00072 dim_t *m_NX; 00073 dim_t *m_NE; 00074 dim_t *m_NN; 00075 escript::Data c11, c12, c13, c23, c33, c44, c66; 00076 bool isVTI, isHTI; 00077 }; 00078 00079 } // namespace ripley 00080 00081 00082 #endif // __RIPLEY_WAVEASSEMBLER2D_H__ 00083