escript  Revision_
DefaultAssembler3D.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_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