escript  Revision_
ripley/src/DefaultAssembler2D.h
Go to the documentation of this file.
1 
2 /*****************************************************************************
3 *
4 * Copyright (c) 2003-2020 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-2017 by Centre for Geoscience Computing (GeoComp)
14 * Development from 2019 by School of Earth and Environmental Sciences
15 **
16 *****************************************************************************/
17 #ifndef __RIPLEY_DEFAULTASSEMBLER2D_H__
18 #define __RIPLEY_DEFAULTASSEMBLER2D_H__
19 
20 #include <ripley/Rectangle.h>
21 
22 namespace ripley {
23 
24 template<class Scalar = double>
26 {
27 public:
29  const dim_t* NE, const dim_t* NN)
31  m_dx(dx),
32  m_NE(NE),
33  m_NN(NN)
34  {
35  domain = REFCOUNTNS::static_pointer_cast<const Rectangle>(dom);
36  }
37 
39 
40  /* The default RipleyDomain assemblers, with original signatures */
41 
45  const escript::Data& A, const escript::Data& B,
46  const escript::Data& C, const escript::Data& D,
47  const escript::Data& X, const escript::Data& Y) const;
48 
52  escript::Data& rhs, const escript::Data& d,
53  const escript::Data& y) const;
54 
58  escript::Data& rhs, const escript::Data& A, const escript::Data& B,
59  const escript::Data& C, const escript::Data& D,
60  const escript::Data& X, const escript::Data& Y) const;
61 
65  escript::Data& rhs, const escript::Data& d,
66  const escript::Data& y) const;
67 
71  const escript::Data& A, const escript::Data& B,
72  const escript::Data& C, const escript::Data& D,
73  const escript::Data& X, const escript::Data& Y) const;
74 
78  escript::Data& rhs, const escript::Data& d,
79  const escript::Data& y) const;
80 
84  escript::Data& rhs, const escript::Data& A, const escript::Data& B,
85  const escript::Data& C, const escript::Data& D,
86  const escript::Data& X, const escript::Data& Y) const;
87 
91  escript::Data& rhs, const escript::Data& d,
92  const escript::Data& y) const;
93 
94  /* The new interface for assemblers */
95 
97  escript::Data& rhs,
98  const DataMap& coefs) const;
100  escript::Data& rhs,
101  const DataMap& coefs) const;
103  escript::Data& rhs,
104  const DataMap& coefs) const;
107  escript::Data& rhs,
108  const DataMap& coefs) const;
110  escript::Data& rhs,
111  const DataMap& coefs) const;
113  escript::Data& rhs,
114  const DataMap& coefs) const;
116  escript::Data& rhs,
117  const DataMap& coefs) const;
120  escript::Data& rhs,
121  const DataMap& coefs) const;
122 
123  virtual void collateFunctionSpaceTypes(std::vector<int>& fsTypes,
124  const DataMap& coefs) const;
125 
126 protected:
128  const double *m_dx;
129  const dim_t *m_NE;
130  const dim_t *m_NN;
131 };
132 
133 } // namespace ripley
134 
135 #endif // __RIPLEY_DEFAULTASSEMBLER2D_H__
136 
ripley::isNotEmpty
bool isNotEmpty(const std::string target, const DataMap &mapping)
Definition: ripley/src/domainhelpers.h:40
escript::Data::isComplex
bool isComplex() const
True if components of this data are stored as complex.
Definition: Data.cpp:1160
ripley::DefaultAssembler2D::assemblePDESingleReduced
virtual void assemblePDESingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: ripley/src/DefaultAssembler2D.cpp:842
ripley
Definition: ripley/src/AbstractAssembler.h:26
ripley::DataMap
std::map< std::string, escript::Data > DataMap
Definition: ripley/src/domainhelpers.h:25
speckley::DefaultAssembler2D::assemblePDEBoundarySingleReduced
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:580
escript::Data::requireWrite
void requireWrite()
Ensures data is ready for write access. This means that the data will be resolved if lazy and will be...
Definition: Data.cpp:1239
INDEX4
#define INDEX4(_X1_, _X2_, _X3_, _X4_, _N1_, _N2_, _N3_)
Definition: index.h:27
ripley::DefaultAssembler2D::m_NN
const dim_t * m_NN
Definition: ripley/src/DefaultAssembler2D.h:130
ripley::AbstractAssembler
Definition: ripley/src/AbstractAssembler.h:34
ripley::DefaultAssembler2D::assemblePDESingle
virtual void assemblePDESingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: ripley/src/DefaultAssembler2D.cpp:152
INDEX2
#define INDEX2(_X1_, _X2_, _N1_)
Definition: index.h:23
escript::Data::getSampleDataRO
const DataTypes::real_t * getSampleDataRO(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0) const
Return the sample data for the given sample no. Please do not use this unless you NEED to access samp...
Definition: Data.h:1989
speckley::DefaultAssembler2D::assembleComplexPDESystemReduced
virtual void assembleComplexPDESystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:559
speckley::DefaultAssembler2D::assembleComplexPDESystem
virtual void assembleComplexPDESystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:407
speckley::DefaultAssembler2D::assemblePDEBoundarySingle
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:568
ripley::DefaultAssembler2D::assemblePDEBoundarySingleReduced
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: ripley/src/DefaultAssembler2D.cpp:1014
speckley::DefaultAssembler2D::assembleComplexPDEBoundarySingle
virtual void assembleComplexPDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:574
INDEX5
#define INDEX5(_X1_, _X2_, _X3_, _X4_, _X5_, _N1_, _N2_, _N3_, _N4_)
Definition: index.h:29
escript::Data::getDataPointSize
int getDataPointSize() const
Return the size of the data point. It is the product of the data point shape dimensions.
Definition: Data.cpp:1360
ripley::unpackData
const escript::Data unpackData(const std::string target, const DataMap &mapping)
Definition: ripley/src/domainhelpers.h:29
escript::Data::complicate
void complicate()
Definition: Data.cpp:6057
speckley::DefaultAssembler2D::m_NE
const dim_t * m_NE
Definition: speckley/src/DefaultAssembler2D.h:157
speckley::DefaultAssembler2D::m_NN
const dim_t * m_NN
Definition: speckley/src/DefaultAssembler2D.h:158
speckley
Definition: AbstractAssembler.cpp:19
speckley::DefaultAssembler2D::assemblePDESystem
virtual void assemblePDESystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:271
ripley::Rectangle
Rectangle is the 2-dimensional implementation of a RipleyDomain.
Definition: ripley/src/Rectangle.h:30
escript::Data
Data represents a collection of datapoints.
Definition: Data.h:64
speckley::DefaultAssembler2D::assemblePDESingle
virtual void assemblePDESingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:621
speckley::DefaultAssembler2D::assembleComplexPDESingleReduced
virtual void assembleComplexPDESingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:851
INDEX3
#define INDEX3(_X1_, _X2_, _X3_, _N1_, _N2_)
Definition: index.h:25
escript::DataTypes::dim_t
index_t dim_t
Definition: DataTypes.h:65
speckley::DefaultAssembler2D::assembleComplexPDEBoundarySystem
virtual void assembleComplexPDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:600
escript::AbstractSystemMatrix
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:44
ripley::DefaultAssembler2D::assemblePDEBoundarySingle
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: ripley/src/DefaultAssembler2D.cpp:601
all_weights
const double all_weights[][11]
Definition: speckley/src/DefaultAssembler2D.cpp:23
escript::Data::getSampleDataRW
DataTypes::real_t * getSampleDataRW(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0)
Return the sample data for the given sample no. Please do not use this unless you NEED to access samp...
Definition: Data.h:1954
ripley::DefaultAssembler2D::assemblePDESystemReduced
virtual void assemblePDESystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: ripley/src/DefaultAssembler2D.cpp:1983
speckley::DefaultAssembler2D::assemblePDESystemReduced
virtual void assemblePDESystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:550
speckley::DefaultAssembler2D::m_dx
const double * m_dx
Definition: speckley/src/DefaultAssembler2D.h:156
speckley::unpackData
escript::Data unpackData(std::string target, std::map< std::string, escript::Data > mapping)
Definition: AbstractAssembler.cpp:21
ripley::DefaultAssembler2D::~DefaultAssembler2D
~DefaultAssembler2D()
Definition: ripley/src/DefaultAssembler2D.h:38
speckley::DefaultAssembler2D::assembleComplexPDESingle
virtual void assembleComplexPDESingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:727
escript::Data::isEmpty
bool isEmpty() const
Definition: Data.cpp:1132
all_lagrange_derivs
const double all_lagrange_derivs[][11][11]
Definition: speckley/src/DefaultAssembler2D.cpp:35
speckley::DefaultAssembler2D::assemblePDEBoundarySystemReduced
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:606
speckley::DefaultAssembler2D::assemblePDEBoundarySystem
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:594
ripley::DefaultAssembler2D
Definition: ripley/src/DefaultAssembler2D.h:26
ripley::DefaultAssembler2D::m_dx
const double * m_dx
Definition: ripley/src/DefaultAssembler2D.h:128
speckley::DefaultAssembler2D::assembleComplexPDEBoundarySingleReduced
virtual void assembleComplexPDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:587
ripley::DefaultAssembler2D::POINTER_WRAPPER_CLASS
POINTER_WRAPPER_CLASS(const Rectangle) domain
ripley::DefaultAssembler2D::collateFunctionSpaceTypes
virtual void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: ripley/src/DefaultAssembler2D.cpp:31
escript::DataTypes::index_t
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:60
speckley::DefaultAssembler2D::assembleComplexPDEBoundarySystemReduced
virtual void assembleComplexPDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/DefaultAssembler2D.cpp:613
ripley::DefaultAssembler2D::assemblePDEBoundarySystemReduced
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: ripley/src/DefaultAssembler2D.cpp:2179
ripley::DefaultAssembler2D::assemblePDESystem
virtual void assemblePDESystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: ripley/src/DefaultAssembler2D.cpp:1186
speckley::DefaultAssembler2D::collateFunctionSpaceTypes
virtual void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/DefaultAssembler2D.cpp:116
escript::Data::actsExpanded
bool actsExpanded() const
Return true if this Data is expanded or resolves to expanded. That is, if it has a separate value for...
Definition: Data.cpp:1118
escript::DataTypes::cplx_t
std::complex< real_t > cplx_t
complex data type
Definition: DataTypes.h:54
escript::AbstractSystemMatrix::getRowBlockSize
int getRowBlockSize() const
returns the row block size
Definition: AbstractSystemMatrix.h:114
escript::Scalar
Data Scalar(double value, const FunctionSpace &what, bool expanded)
A collection of factory functions for creating Data objects which contain data points of various shap...
Definition: DataFactory.cpp:49
speckley::DataMap
std::map< std::string, escript::Data > DataMap
Definition: speckley/src/domainhelpers.h:25
escript::AbstractSystemMatrix::getColumnBlockSize
int getColumnBlockSize() const
returns the column block size
Definition: AbstractSystemMatrix.h:125
escript::const_Domain_ptr
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:44
ripley::DefaultAssembler2D::assemblePDEBoundarySystem
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: ripley/src/DefaultAssembler2D.cpp:1682
speckley::isNotEmpty
bool isNotEmpty(const std::string target, const DataMap &mapping)
Definition: speckley/src/domainhelpers.h:40
speckley::SpeckleyException
SpeckleyException exception class.
Definition: SpeckleyException.h:31
ripley::DefaultAssembler2D::m_NE
const dim_t * m_NE
Definition: ripley/src/DefaultAssembler2D.h:129
speckley::DefaultAssembler2D::assemblePDESingleReduced
virtual void assemblePDESingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &A, const escript::Data &B, const escript::Data &C, const escript::Data &D, const escript::Data &X, const escript::Data &Y) const
Definition: speckley/src/DefaultAssembler2D.cpp:843