Go to the documentation of this file.
16 #ifndef __SPECKLEY_WAVE_ASSEMBLER_2D_H__
17 #define __SPECKLEY_WAVE_ASSEMBLER_2D_H__
19 #include <speckley/Rectangle.h>
24 class WaveAssembler2D :
public AbstractAssembler
29 : AbstractAssembler(),
34 domain = REFCOUNTNS::static_pointer_cast<const Rectangle>(dom);
36 DataMap::const_iterator a = c.find(
"c12"), b = c.find(
"c23");
37 if (c.find(
"c11") == c.end()
38 || c.find(
"c13") == c.end() || c.find(
"c33") == c.end()
39 || c.find(
"c44") == c.end() || c.find(
"c66") == c.end()
40 || (a == c.end() && b == c.end()))
43 if (a != c.end() && b != c.end()) {
45 }
else if (a == c.end()) {
54 }
else if (b == c.end()) {
58 throw SpeckleyException(
"C tensor elements must be reduced");
61 throw SpeckleyException(
"C tensor elements must not be empty");
64 c11 = c.find(
"c11")->second;
65 c13 = c.find(
"c13")->second;
66 c33 = c.find(
"c33")->second;
67 c44 = c.find(
"c44")->second;
68 c66 = c.find(
"c66")->second;
74 throw SpeckleyException(
"C tensor elements must be reduced");
81 throw SpeckleyException(
"C tensor elements must not be empty");
184 #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:395
POINTER_WRAPPER_CLASS(const Rectangle) domain
escript::Data c33
Definition: speckley/src/WaveAssembler2D.h:202
escript::Data c12
Definition: speckley/src/WaveAssembler2D.h:202
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:376
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:50
escript::Data c11
Definition: speckley/src/WaveAssembler2D.h:202
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:461
~WaveAssembler2D()
Definition: speckley/src/WaveAssembler2D.h:108
escript::Data c44
Definition: speckley/src/WaveAssembler2D.h:202
Definition: AbstractAssembler.cpp:17
Data represents a collection of datapoints.
Definition: Data.h:62
index_t dim_t
Definition: DataTypes.h:87
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:53
const dim_t * m_NE
Definition: speckley/src/WaveAssembler2D.h:199
Definition: Speckley.h:75
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:389
bool isEmpty() const
Definition: Data.cpp:1135
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:412
Rectangle is the 2-dimensional implementation of a SpeckleyDomain.
Definition: speckley/src/Rectangle.h:43
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:382
bool isHTI
Definition: speckley/src/WaveAssembler2D.h:201
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:403
escript::Data c13
Definition: speckley/src/WaveAssembler2D.h:202
const double * m_dx
Definition: speckley/src/WaveAssembler2D.h:198
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:367
int getTypeCode() const
Returns the function space type code.
Definition: FunctionSpace.cpp:93
bool isVTI
Definition: speckley/src/WaveAssembler2D.h:201
virtual void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/WaveAssembler2D.cpp:114
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:210
escript::Data c23
Definition: speckley/src/WaveAssembler2D.h:202
escript::Data c66
Definition: speckley/src/WaveAssembler2D.h:202
std::map< std::string, escript::Data > DataMap
Definition: speckley/src/domainhelpers.h:35
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:50
SpeckleyException exception class.
Definition: SpeckleyException.h:40
const dim_t * m_NN
Definition: speckley/src/WaveAssembler2D.h:200