escript  Revision_
WaveAssembler3D.h
Go to the documentation of this file.
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