escript  Revision_
WaveAssembler2D.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/Rectangle.h>
00025 
00026 namespace ripley {
00027 
00028 
00029 class WaveAssembler2D : public AbstractAssembler {
00030 public:
00031     WaveAssembler2D(Rectangle *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     ~WaveAssembler2D(){};
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 {throw RipleyException("This assembly not supported by this assembler");}
00042     void assemblePDEBoundarySingle(paso::SystemMatrix_ptr mat, escript::Data& rhs,
00043             std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");}
00044     void assemblePDESingleReduced(paso::SystemMatrix_ptr mat, escript::Data& rhs,
00045             std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");}
00046     void assemblePDEBoundarySingleReduced(paso::SystemMatrix_ptr mat,
00047             escript::Data& rhs,
00048             std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");}
00049     void assemblePDEBoundarySystem(paso::SystemMatrix_ptr mat, escript::Data& rhs,
00050             std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");}
00051     void assemblePDESystemReduced(paso::SystemMatrix_ptr mat, escript::Data& rhs,
00052             std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");}
00053     void assemblePDEBoundarySystemReduced(paso::SystemMatrix_ptr mat,
00054             escript::Data& rhs,
00055             std::map<std::string, escript::Data> coefs) const {throw RipleyException("This assembly not supported by this assembler");}
00056 
00057     void collateFunctionSpaceTypes(std::vector<int>& fsTypes, 
00058             std::map<std::string, escript::Data> coefs) const;
00059 
00060 private:
00061     std::map<std::string, escript::Data> c;
00062     Rectangle *domain;
00063     double *m_dx;
00064     dim_t *m_NX;
00065     dim_t *m_NE;
00066     dim_t *m_NN;
00067     escript::Data c11, c12, c13, c23, c33, c44, c66;
00068     bool isVTI, isHTI;
00069 };
00070 
00071 } // namespace ripley
00072 
00073 #endif // __RIPLEY_WAVEASSEMBLER2D_H__
00074