escript  Revision_
AbstractContinuousDomain.h
Go to the documentation of this file.
1 
2 /*****************************************************************************
3 *
4 * Copyright (c) 2003-2016 by The University of Queensland
5 * http://www.uq.edu.au
6 *
7 * Primary Business: Queensland, Australia
8 * Licensed under the Apache License, version 2.0
9 * http://www.apache.org/licenses/LICENSE-2.0
10 *
11 * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12 * Development 2012-2013 by School of Earth Sciences
13 * Development from 2014 by Centre for Geoscience Computing (GeoComp)
14 *
15 *****************************************************************************/
16 
17 
18 #if !defined escript_AbstractContinuousDomain_20040528_H
19 #define escript_AbstractContinuousDomain_20040528_H
20 
21 #include "system_dep.h"
22 #include "AbstractDomain.h"
23 
24 #include <string>
25 #include <vector>
26 
27 #include "AbstractSystemMatrix.h"
29 
30 namespace escript {
31 
32 //
33 // Forward declaration
34 class Data;
35 
47 {
48 
49  public:
50 
65 
79  virtual ~AbstractContinuousDomain();
80 
86  virtual bool isValidFunctionSpaceType(int functionSpaceType) const;
87 
92  virtual std::string getDescription() const;
93 
98  virtual int getContinuousFunctionCode() const;
99 
104  virtual int getReducedContinuousFunctionCode() const;
105 
110  virtual int getFunctionCode() const;
111 
116  virtual int getReducedFunctionCode() const;
117 
122  virtual int getFunctionOnBoundaryCode() const;
123 
128  virtual int getReducedFunctionOnBoundaryCode() const;
129 
130 
135  virtual int getFunctionOnContactZeroCode() const;
136 
141  virtual int getReducedFunctionOnContactZeroCode() const;
142 
147  virtual int getFunctionOnContactOneCode() const;
148 
153  virtual int getReducedFunctionOnContactOneCode() const;
154 
159  virtual int getSolutionCode() const;
160 
165  virtual int getReducedSolutionCode() const;
166 
171  virtual int getDiracDeltaFunctionsCode() const;
172 
178  virtual int getSystemMatrixTypeId(const boost::python::object& options) const;
179 
186  virtual int getTransportTypeId(int solver, int preconditioner, int package, bool symmetry) const;
187 
194  virtual void setToIntegrals(std::vector<double>& integrals,const escript::Data& arg) const;
195 
196 // /**
197 // \brief
198 // Return the domain as const AbstractContinuousDomain&
199 // */
200 // static const AbstractContinuousDomain& asAbstractContinuousDomain(const AbstractDomain& domain);
201 
202 
203 
204 
209  virtual void addPDEToSystem(
211  const escript::Data& A, const escript::Data& B, const escript::Data& C,
212  const escript::Data& D, const escript::Data& X, const escript::Data& Y,
213  const escript::Data& d, const escript::Data& y,
214  const escript::Data& d_contact, const escript::Data& y_contact,
215  const escript::Data& d_dirac, const escript::Data& y_dirac) const;
216 
217 // We do not require this method at this level since the python side checks to ensure it exists
218 // before calling it.
219 
220 // /**
221 // \brief
222 // adds a PDE onto the lumped stiffness matrix matrix
223 // */
224 // virtual void addPDEToLumpedSystem(
225 // escript::Data& mat,
226 // const escript::Data& D,
227 // const escript::Data& d) const;
228 
233  virtual void addPDEToRHS(escript::Data& rhs,
234  const escript::Data& X, const escript::Data& Y,
235  const escript::Data& y, const escript::Data& y_contact, const escript::Data& y_dirac) const;
241  virtual void addPDEToTransportProblem(
243  const escript::Data& M,
244  const escript::Data& A, const escript::Data& B, const escript::Data& C,const escript::Data& D,
245  const escript::Data& X,const escript::Data& Y,
246  const escript::Data& d, const escript::Data& y,
247  const escript::Data& d_contact,const escript::Data& y_contact,
248  const escript::Data& d_dirac,const escript::Data& y_dirac) const;
249 
254  virtual ASM_ptr newSystemMatrix(
255  const int row_blocksize,
256  const escript::FunctionSpace& row_functionspace,
257  const int column_blocksize,
258  const escript::FunctionSpace& column_functionspace,
259  const int type) const;
266  virtual ATP_ptr newTransportProblem(
267  const int blocksize,
268  const escript::FunctionSpace& functionspace,
269  const int type) const;
270 
275  virtual dim_t getNumDataPointsGlobal() const;
276 
282  virtual std::pair<int,dim_t> getDataShape(int functionSpaceCode) const;
283 
288  virtual void setNewX(const escript::Data& arg);
289 
294  virtual void Print_Mesh_Info(const bool full=false) const;
295 
296 
297  protected:
298 
299  private:
300 
301 };
302 
303 } // end of namespace
304 
305 #endif
AbstractContinuousDomain, base class for continuous domains.
Definition: AbstractContinuousDomain.h:46
Definition: FunctionSpace.h:34
Definition: AbstractContinuousDomain.cpp:24
static dim_t M
Definition: SparseMatrix_saveHB.cpp:37
boost::shared_ptr< AbstractTransportProblem > ATP_ptr
Definition: AbstractTransportProblem.h:162
Data represents a collection of datapoints.
Definition: Data.h:68
Give a short description of what AbstractTransportProblem does.
Definition: AbstractTransportProblem.h:45
#define ESCRIPT_DLL_API
Definition: escriptcore/src/system_dep.h:54
index_t dim_t
Definition: types.h:27
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:37
boost::shared_ptr< AbstractSystemMatrix > ASM_ptr
Definition: AbstractSystemMatrix.h:170
Base class for all escript domains.
Definition: AbstractDomain.h:45