16 #ifndef __SPECKLEY_WAVE_ASSEMBLER_3D_H__ 17 #define __SPECKLEY_WAVE_ASSEMBLER_3D_H__ 20 #include <escript/Data.h> 21 #include <speckley/Speckley.h> 22 #include <speckley/SpeckleyException.h> 23 #include <speckley/AbstractAssembler.h> 24 #include <speckley/Brick.h> 39 domain = boost::static_pointer_cast<
const Brick>(dom);
41 DataMap::const_iterator a = c.find(
"c12"), b = c.find(
"c23");
42 if (c.find(
"c11") == c.end()
43 || c.find(
"c13") == c.end() || c.find(
"c33") == c.end()
44 || c.find(
"c44") == c.end() || c.find(
"c66") == c.end()
45 || (a == c.end() && b == c.end()))
48 if (a != c.end() && b != c.end()) {
50 }
else if (a == c.end()) {
59 }
else if (b == c.end()) {
69 c11 = c.find(
"c11")->second;
70 c13 = c.find(
"c13")->second;
71 c33 = c.find(
"c33")->second;
72 c44 = c.find(
"c44")->second;
73 c66 = c.find(
"c66")->second;
185 #endif // __SPECKLEY_WAVE_ASSEMBLER_3D_H__ escript::Data c33
Definition: speckley/src/WaveAssembler3D.h:180
escript::Data c44
Definition: speckley/src/WaveAssembler3D.h:180
Definition: AbstractAssembler.cpp:22
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:532
bool isVTI
Definition: speckley/src/WaveAssembler3D.h:179
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:213
boost::shared_ptr< const Brick > domain
Definition: speckley/src/WaveAssembler3D.h:175
void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/WaveAssembler3D.cpp:118
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:518
escript::Data c23
Definition: speckley/src/WaveAssembler3D.h:180
std::map< std::string, escript::Data > DataMap
Definition: speckley/src/domainhelpers.h:24
const dim_t * m_NN
Definition: speckley/src/WaveAssembler3D.h:178
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:508
escript::Data c11
Definition: speckley/src/WaveAssembler3D.h:180
~WaveAssembler3D()
Definition: speckley/src/WaveAssembler3D.h:90
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:556
Definition: Speckley.h:49
const double * m_dx
Definition: speckley/src/WaveAssembler3D.h:176
escript::Data c13
Definition: speckley/src/WaveAssembler3D.h:180
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:455
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:547
Data represents a collection of datapoints.
Definition: Data.h:68
escript::Data c12
Definition: speckley/src/WaveAssembler3D.h:180
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:539
SpeckleyException exception class.
Definition: SpeckleyException.h:29
int getTypeCode() const
Returns the function space type code.
Definition: FunctionSpace.cpp:99
bool isEmpty() const
Return true if this Data holds an instance of DataEmpty. This is _not_ the same as asking if the obj...
Definition: Data.cpp:853
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:525
const dim_t * m_NE
Definition: speckley/src/WaveAssembler3D.h:177
Definition: speckley/src/AbstractAssembler.h:33
index_t dim_t
Definition: types.h:27
bool isHTI
Definition: speckley/src/WaveAssembler3D.h:179
Brick is the 3-dimensional implementation of a SpeckleyDomain.
Definition: speckley/src/Brick.h:32
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:37
Definition: speckley/src/WaveAssembler3D.h:29
escript::Data c66
Definition: speckley/src/WaveAssembler3D.h:180
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:32
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:39