Go to the documentation of this file.
16 #ifndef __SPECKLEY_WAVE_ASSEMBLER_3D_H__
17 #define __SPECKLEY_WAVE_ASSEMBLER_3D_H__
19 #include <speckley/Speckley.h>
20 #include <speckley/AbstractAssembler.h>
21 #include <speckley/Brick.h>
22 #include <speckley/SpeckleyException.h>
24 #include <escript/Data.h>
31 class WaveAssembler3D :
public AbstractAssembler
36 : AbstractAssembler(),
41 domain = REFCOUNTNS::static_pointer_cast<const Brick>(dom);
43 DataMap::const_iterator a = c.find(
"c12"), b = c.find(
"c23");
44 if (c.find(
"c11") == c.end()
45 || c.find(
"c13") == c.end() || c.find(
"c33") == c.end()
46 || c.find(
"c44") == c.end() || c.find(
"c66") == c.end()
47 || (a == c.end() && b == c.end()))
50 if (a != c.end() && b != c.end()) {
52 }
else if (a == c.end()) {
61 }
else if (b == c.end()) {
65 throw SpeckleyException(
"C tensor elements must be reduced");
68 throw SpeckleyException(
"C tensor elements must not be empty");
71 c11 = c.find(
"c11")->second;
72 c13 = c.find(
"c13")->second;
73 c33 = c.find(
"c33")->second;
74 c44 = c.find(
"c44")->second;
75 c66 = c.find(
"c66")->second;
81 throw SpeckleyException(
"C tensor elements must be reduced");
88 throw SpeckleyException(
"C tensor elements must not be empty");
187 #endif // __SPECKLEY_WAVE_ASSEMBLER_3D_H__
escript::Data c11
Definition: speckley/src/WaveAssembler3D.h:205
escript::Data c33
Definition: speckley/src/WaveAssembler3D.h:205
escript::Data c44
Definition: speckley/src/WaveAssembler3D.h:205
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:514
WaveAssembler3D(escript::const_Domain_ptr dom, const double *dx, const dim_t *NE, const dim_t *NN, const DataMap &c)
Definition: speckley/src/WaveAssembler3D.h:57
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:535
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:461
bool isHTI
Definition: speckley/src/WaveAssembler3D.h:204
Definition: AbstractAssembler.cpp:17
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:521
Data represents a collection of datapoints.
Definition: Data.h:62
void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/WaveAssembler3D.cpp:114
index_t dim_t
Definition: DataTypes.h:87
escript::Data c23
Definition: speckley/src/WaveAssembler3D.h:205
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:53
Brick is the 3-dimensional implementation of a SpeckleyDomain.
Definition: speckley/src/Brick.h:43
escript::Data c66
Definition: speckley/src/WaveAssembler3D.h:205
Definition: Speckley.h:75
bool isEmpty() const
Definition: Data.cpp:1135
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:528
POINTER_WRAPPER_CLASS(const Brick) domain
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/WaveAssembler3D.cpp:543
const dim_t * m_NE
Definition: speckley/src/WaveAssembler3D.h:202
escript::Data c13
Definition: speckley/src/WaveAssembler3D.h:205
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/WaveAssembler3D.cpp:209
const dim_t * m_NN
Definition: speckley/src/WaveAssembler3D.h:203
bool isVTI
Definition: speckley/src/WaveAssembler3D.h:204
int getTypeCode() const
Returns the function space type code.
Definition: FunctionSpace.cpp:93
~WaveAssembler3D()
Definition: speckley/src/WaveAssembler3D.h:115
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
const double * m_dx
Definition: speckley/src/WaveAssembler3D.h:201
escript::Data c12
Definition: speckley/src/WaveAssembler3D.h:205
SpeckleyException exception class.
Definition: SpeckleyException.h:40
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/WaveAssembler3D.cpp:504
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/WaveAssembler3D.cpp:552