17#ifndef __SPECKLEY_WAVE_ASSEMBLER_3D_H__
18#define __SPECKLEY_WAVE_ASSEMBLER_3D_H__
20#include <speckley/Speckley.h>
21#include <speckley/AbstractAssembler.h>
22#include <speckley/Brick.h>
23#include <speckley/SpeckleyException.h>
25#include <escript/Data.h>
42 domain = REFCOUNTNS::static_pointer_cast<const Brick>(dom);
44 DataMap::const_iterator a = c.find(
"c12"), b = c.find(
"c23");
45 if (c.find(
"c11") == c.end()
46 || c.find(
"c13") == c.end() || c.find(
"c33") == c.end()
47 || c.find(
"c44") == c.end() || c.find(
"c66") == c.end()
48 || (a == c.end() && b == c.end()))
51 if (a != c.end() && b != c.end()) {
53 }
else if (a == c.end()) {
62 }
else if (b == c.end()) {
72 c11 = c.find(
"c11")->second;
73 c13 = c.find(
"c13")->second;
74 c33 = c.find(
"c33")->second;
75 c44 = c.find(
"c44")->second;
76 c66 = c.find(
"c66")->second;
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:44
Data represents a collection of datapoints.
Definition: Data.h:64
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:463
bool isEmpty() const
Definition: Data.cpp:1132
int getTypeCode() const
Returns the function space type code.
Definition: FunctionSpace.cpp:95
Definition: speckley/src/AbstractAssembler.h:35
Brick is the 3-dimensional implementation of a SpeckleyDomain.
Definition: speckley/src/Brick.h:34
SpeckleyException exception class.
Definition: SpeckleyException.h:31
Definition: speckley/src/WaveAssembler3D.h:33
escript::Data c33
Definition: speckley/src/WaveAssembler3D.h:183
escript::Data c23
Definition: speckley/src/WaveAssembler3D.h:183
escript::Data c12
Definition: speckley/src/WaveAssembler3D.h:183
bool isHTI
Definition: speckley/src/WaveAssembler3D.h:182
~WaveAssembler3D()
Definition: speckley/src/WaveAssembler3D.h:93
const dim_t * m_NE
Definition: speckley/src/WaveAssembler3D.h:180
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:554
bool isVTI
Definition: speckley/src/WaveAssembler3D.h:182
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:211
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:545
POINTER_WRAPPER_CLASS(const Brick) domain
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:537
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:530
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:523
escript::Data c13
Definition: speckley/src/WaveAssembler3D.h:183
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:35
escript::Data c11
Definition: speckley/src/WaveAssembler3D.h:183
escript::Data c66
Definition: speckley/src/WaveAssembler3D.h:183
escript::Data c44
Definition: speckley/src/WaveAssembler3D.h:183
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:516
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:506
const double * m_dx
Definition: speckley/src/WaveAssembler3D.h:179
const dim_t * m_NN
Definition: speckley/src/WaveAssembler3D.h:181
void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/WaveAssembler3D.cpp:116
index_t dim_t
Definition: DataTypes.h:66
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:44
Definition: AbstractAssembler.cpp:19
@ ReducedElements
Definition: Speckley.h:64
std::map< std::string, escript::Data > DataMap
Definition: speckley/src/domainhelpers.h:25