escript  Revision_
LameAssembler3D.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_LAMEASSEMBLER3D_H__
00017 #define __RIPLEY_LAMEASSEMBLER3D_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 LameAssembler3D : public AbstractAssembler {
00030 public:
00031     LameAssembler3D(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     ~LameAssembler3D(){};
00040     
00041     /* The new interface for assemblers */
00042     virtual void assemblePDESingle(paso::SystemMatrix_ptr mat,
00043             escript::Data& rhs,
00044             std::map<std::string, escript::Data> coefs) const;
00045     virtual void assemblePDEBoundarySingle(paso::SystemMatrix_ptr mat,
00046             escript::Data& rhs,
00047             std::map<std::string, escript::Data> coefs) const;
00048     virtual void assemblePDESingleReduced(paso::SystemMatrix_ptr mat,
00049             escript::Data& rhs,
00050             std::map<std::string, escript::Data> coefs) const;
00051     virtual void assemblePDEBoundarySingleReduced(paso::SystemMatrix_ptr mat,
00052             escript::Data& rhs,
00053             std::map<std::string, escript::Data> coefs) const;
00054     virtual void assemblePDESystem(paso::SystemMatrix_ptr mat,
00055             escript::Data& rhs,
00056             std::map<std::string, escript::Data> coefs) const;
00057     virtual void assemblePDEBoundarySystem(paso::SystemMatrix_ptr mat,
00058             escript::Data& rhs,
00059             std::map<std::string, escript::Data> coefs) const;
00060     virtual void assemblePDESystemReduced(paso::SystemMatrix_ptr mat,
00061             escript::Data& rhs,
00062             std::map<std::string, escript::Data> coefs) const;
00063     virtual void assemblePDEBoundarySystemReduced(paso::SystemMatrix_ptr mat,
00064             escript::Data& rhs,
00065             std::map<std::string, escript::Data> coefs) const;
00066             
00067     void collateFunctionSpaceTypes(std::vector<int>& fsTypes, 
00068             std::map<std::string, escript::Data> coefs) const;
00069 protected:
00070     const Brick *domain;
00071     const double *m_dx;
00072     const dim_t *m_NX;
00073     const dim_t *m_NE;
00074     const dim_t *m_NN;
00075 };
00076 
00077 } // namespace ripley
00078 
00079 
00080 #endif // __RIPLEY_LAMEASSEMBLER3D_H__
00081