Go to the documentation of this file.
17 #ifndef __SPECKLEY_WAVE_ASSEMBLER_2D_H__
18 #define __SPECKLEY_WAVE_ASSEMBLER_2D_H__
20 #include <speckley/Rectangle.h>
25 class WaveAssembler2D :
public AbstractAssembler
30 : AbstractAssembler(),
35 domain = REFCOUNTNS::static_pointer_cast<const Rectangle>(dom);
37 DataMap::const_iterator a = c.find(
"c12"), b = c.find(
"c23");
38 if (c.find(
"c11") == c.end()
39 || c.find(
"c13") == c.end() || c.find(
"c33") == c.end()
40 || c.find(
"c44") == c.end() || c.find(
"c66") == c.end()
41 || (a == c.end() && b == c.end()))
44 if (a != c.end() && b != c.end()) {
46 }
else if (a == c.end()) {
55 }
else if (b == c.end()) {
59 throw SpeckleyException(
"C tensor elements must be reduced");
62 throw SpeckleyException(
"C tensor elements must not be empty");
65 c11 = c.find(
"c11")->second;
66 c13 = c.find(
"c13")->second;
67 c33 = c.find(
"c33")->second;
68 c44 = c.find(
"c44")->second;
69 c66 = c.find(
"c66")->second;
75 throw SpeckleyException(
"C tensor elements must be reduced");
82 throw SpeckleyException(
"C tensor elements must not be empty");
185 #endif // __SPECKLEY_DEFAULTASSEMBLER2D_H__
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:396
POINTER_WRAPPER_CLASS(const Rectangle) domain
escript::Data c33
Definition: speckley/src/WaveAssembler2D.h:205
escript::Data c12
Definition: speckley/src/WaveAssembler2D.h:205
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:377
WaveAssembler2D(escript::const_Domain_ptr dom, const double *dx, const dim_t *NE, const dim_t *NN, const DataMap &c)
Definition: speckley/src/WaveAssembler2D.h:53
escript::Data c11
Definition: speckley/src/WaveAssembler2D.h:205
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:462
~WaveAssembler2D()
Definition: speckley/src/WaveAssembler2D.h:111
escript::Data c44
Definition: speckley/src/WaveAssembler2D.h:205
Definition: AbstractAssembler.cpp:18
Data represents a collection of datapoints.
Definition: Data.h:63
index_t dim_t
Definition: DataTypes.h:90
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:55
const dim_t * m_NE
Definition: speckley/src/WaveAssembler2D.h:202
Definition: Speckley.h:77
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:390
bool isEmpty() const
Definition: Data.cpp:1136
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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:413
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:383
bool isHTI
Definition: speckley/src/WaveAssembler2D.h:204
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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:404
escript::Data c13
Definition: speckley/src/WaveAssembler2D.h:205
const double * m_dx
Definition: speckley/src/WaveAssembler2D.h:201
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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:368
int getTypeCode() const
Returns the function space type code.
Definition: FunctionSpace.cpp:94
bool isVTI
Definition: speckley/src/WaveAssembler2D.h:204
virtual void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/WaveAssembler2D.cpp:115
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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:211
escript::Data c23
Definition: speckley/src/WaveAssembler2D.h:205
escript::Data c66
Definition: speckley/src/WaveAssembler2D.h:205
std::map< std::string, escript::Data > DataMap
Definition: speckley/src/domainhelpers.h:37
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:56
SpeckleyException exception class.
Definition: SpeckleyException.h:42
const dim_t * m_NN
Definition: speckley/src/WaveAssembler2D.h:203