ripley/src/WaveAssembler3D.h
Go to the documentation of this file.
/*****************************************************************************
*
* Copyright (c) 2003-2020 by The University of Queensland
* http://www.uq.edu.au
*
* Primary Business: Queensland, Australia
* Licensed under the Apache License, version 2.0
* http://www.apache.org/licenses/LICENSE-2.0
*
* Development until 2012 by Earth Systems Science Computational Center (ESSCC)
* Development 2012-2013 by School of Earth Sciences
* Development from 2014-2017 by Centre for Geoscience Computing (GeoComp)
* Development from 2019 by School of Earth and Environmental Sciences
**
*****************************************************************************/
#ifndef __RIPLEY_WAVEASSEMBLER2D_H__
#define __RIPLEY_WAVEASSEMBLER2D_H__
#include <ripley/Brick.h>
namespace ripley {
class WaveAssembler3D : public AbstractAssembler
{
public:
WaveAssembler3D(escript::const_Domain_ptr dom, const double *dx,
const dim_t *NE, const dim_t *NN,
const DataMap& c);
~WaveAssembler3D() {}
/* The only assembly function we care about right now*/
void assemblePDESystem(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const;
void assemblePDESingle(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const {
throw escript::NotImplementedError("assemblePDESingle() not supported by this assembler");
}
void assemblePDEBoundarySingle(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const {
throw escript::NotImplementedError("assemblePDEBoundarySingle() not supported by this assembler");
}
void assemblePDESingleReduced(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const {
throw escript::NotImplementedError("assemblePDESingleReduced() not supported by this assembler");
}
void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const {
throw escript::NotImplementedError("assemblePDEBoundarySingleReduced() not supported by this assembler");
}
void assemblePDEBoundarySystem(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const {
throw escript::NotImplementedError("assemblePDEBoundarySystem() not supported by this assembler");
}
void assemblePDESystemReduced(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const {
throw escript::NotImplementedError("assemblePDESystemReduced() not supported by this assembler");
}
void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix* mat,
escript::Data& rhs, const DataMap& coefs) const {
throw escript::NotImplementedError("assemblePDEBoundarySystemReduced() not supported by this assembler");
}
void collateFunctionSpaceTypes(std::vector<int>& fsTypes,
const DataMap& coefs) const;
private:
DataMap c;
POINTER_WRAPPER_CLASS(const Brick) domain;
const double *m_dx;
const dim_t *m_NE;
const dim_t *m_NN;
escript::Data c11, c12, c13, c23, c33, c44, c66;
bool isVTI, isHTI;
};
} // namespace ripley
#endif // __RIPLEY_WAVEASSEMBLER2D_H__
ripley::WaveAssembler3D::c23
escript::Data c23
Definition: ripley/src/WaveAssembler3D.h:105
speckley::WaveAssembler3D::c11
escript::Data c11
Definition: speckley/src/WaveAssembler3D.h:211
ripley::isNotEmpty
bool isNotEmpty(const std::string target, const DataMap &mapping)
Definition: ripley/src/domainhelpers.h:54
ripley
Definition: ripley/src/AbstractAssembler.h:26
ripley::DataMap
std::map< std::string, escript::Data > DataMap
Definition: ripley/src/domainhelpers.h:39
escript::Data::requireWrite
void requireWrite()
Ensures data is ready for write access. This means that the data will be resolved if lazy and will be...
Definition: Data.cpp:1239
speckley::WaveAssembler3D::c33
escript::Data c33
Definition: speckley/src/WaveAssembler3D.h:211
INDEX4
#define INDEX4(_X1_, _X2_, _X3_, _X4_, _N1_, _N2_, _N3_)
Definition: index.h:27
ripley::AbstractAssembler
Definition: ripley/src/AbstractAssembler.h:47
INDEX2
#define INDEX2(_X1_, _X2_, _N1_)
Definition: index.h:23
ripley::RipleyException
RipleyException exception class.
Definition: RipleyException.h:44
ripley::WaveAssembler3D::~WaveAssembler3D
~WaveAssembler3D()
Definition: ripley/src/WaveAssembler3D.h:60
speckley::WaveAssembler3D::c44
escript::Data c44
Definition: speckley/src/WaveAssembler3D.h:211
escript::Data::getSampleDataRO
const DataTypes::real_t * getSampleDataRO(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0) const
Return the sample data for the given sample no. Please do not use this unless you NEED to access samp...
Definition: Data.h:1989
ripley::WaveAssembler3D::assemblePDESingleReduced
void assemblePDESingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.h:75
ripley::WaveAssembler3D::assemblePDESystem
void assemblePDESystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.cpp:88
speckley::WaveAssembler3D::assemblePDEBoundarySingle
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:516
escript::NotImplementedError
An exception class for features which are not (yet) implemented.
Definition: EsysException.h:93
ripley::WaveAssembler3D::POINTER_WRAPPER_CLASS
POINTER_WRAPPER_CLASS(const Brick) domain
speckley::WaveAssembler3D::assemblePDEBoundarySystemReduced
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:537
INDEX5
#define INDEX5(_X1_, _X2_, _X3_, _X4_, _X5_, _N1_, _N2_, _N3_, _N4_)
Definition: index.h:29
escript::Data::getDataPointSize
int getDataPointSize() const
Return the size of the data point. It is the product of the data point shape dimensions.
Definition: Data.cpp:1360
ripley::WaveAssembler3D::isVTI
bool isVTI
Definition: ripley/src/WaveAssembler3D.h:106
ripley::unpackData
const escript::Data unpackData(const std::string target, const DataMap &mapping)
Definition: ripley/src/domainhelpers.h:43
ripley::WaveAssembler3D::assemblePDEBoundarySystemReduced
void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.h:91
all_weights
const double all_weights[][11]
Definition: speckley/src/WaveAssembler3D.cpp:23
escript::Data::getFunctionSpace
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:463
ripley::WaveAssembler3D::assemblePDEBoundarySingleReduced
void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.h:79
speckley::WaveAssembler3D::isHTI
bool isHTI
Definition: speckley/src/WaveAssembler3D.h:210
ripley::WaveAssembler3D::assemblePDEBoundarySystem
void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.h:83
speckley
Definition: AbstractAssembler.cpp:19
ripley::WaveAssembler3D::WaveAssembler3D
WaveAssembler3D(escript::const_Domain_ptr dom, const double *dx, const dim_t *NE, const dim_t *NN, const DataMap &c)
Definition: ripley/src/WaveAssembler3D.cpp:27
ripley::WaveAssembler3D::m_NN
const dim_t * m_NN
Definition: ripley/src/WaveAssembler3D.h:104
speckley::WaveAssembler3D::assemblePDEBoundarySingleReduced
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
Data represents a collection of datapoints.
Definition: Data.h:64
speckley::WaveAssembler3D::collateFunctionSpaceTypes
void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/WaveAssembler3D.cpp:116
INDEX3
#define INDEX3(_X1_, _X2_, _X3_, _N1_, _N2_)
Definition: index.h:25
escript::DataTypes::dim_t
index_t dim_t
Definition: DataTypes.h:93
speckley::WaveAssembler3D::c23
escript::Data c23
Definition: speckley/src/WaveAssembler3D.h:211
escript::AbstractSystemMatrix
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:57
ripley::WaveAssembler3D::c12
escript::Data c12
Definition: ripley/src/WaveAssembler3D.h:105
ripley::WaveAssembler3D::c11
escript::Data c11
Definition: ripley/src/WaveAssembler3D.h:105
escript::Data::getSampleDataRW
DataTypes::real_t * getSampleDataRW(DataTypes::RealVectorType::size_type sampleNo, DataTypes::real_t dummy=0)
Return the sample data for the given sample no. Please do not use this unless you NEED to access samp...
Definition: Data.h:1954
ripley::WaveAssembler3D::m_NE
const dim_t * m_NE
Definition: ripley/src/WaveAssembler3D.h:103
speckley::unpackData
escript::Data unpackData(std::string target, std::map< std::string, escript::Data > mapping)
Definition: AbstractAssembler.cpp:35
ripley::WaveAssembler3D::c66
escript::Data c66
Definition: ripley/src/WaveAssembler3D.h:105
ripley::WaveAssembler3D::c13
escript::Data c13
Definition: ripley/src/WaveAssembler3D.h:105
SI
#define SI(_x_, _y_, _qx_, _qy_, _qz_)
ripley::WaveAssembler3D::c33
escript::Data c33
Definition: ripley/src/WaveAssembler3D.h:105
speckley::WaveAssembler3D::c66
escript::Data c66
Definition: speckley/src/WaveAssembler3D.h:211
escript::Data::isEmpty
bool isEmpty() const
Definition: Data.cpp:1132
speckley::WaveAssembler3D::assemblePDEBoundarySystem
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler3D.cpp:530
ripley::WaveAssembler3D::assemblePDESingle
void assemblePDESingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.h:67
speckley::WaveAssembler3D::assemblePDESingle
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
ripley::WaveAssembler3D::c
DataMap c
Definition: ripley/src/WaveAssembler3D.h:100
speckley::WaveAssembler3D::m_NE
const dim_t * m_NE
Definition: speckley/src/WaveAssembler3D.h:208
ripley::Brick
Brick is the 3-dimensional implementation of a RipleyDomain.
Definition: ripley/src/Brick.h:43
escript::DataTypes::index_t
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:88
speckley::WaveAssembler3D::c13
escript::Data c13
Definition: speckley/src/WaveAssembler3D.h:211
speckley::WaveAssembler3D::assemblePDESystem
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
speckley::WaveAssembler3D::m_NN
const dim_t * m_NN
Definition: speckley/src/WaveAssembler3D.h:209
ripley::WaveAssembler3D::isHTI
bool isHTI
Definition: ripley/src/WaveAssembler3D.h:106
escript::FunctionSpace::getTypeCode
int getTypeCode() const
Returns the function space type code.
Definition: FunctionSpace.cpp:95
escript::Data::actsExpanded
bool actsExpanded() const
Return true if this Data is expanded or resolves to expanded. That is, if it has a separate value for...
Definition: Data.cpp:1118
all_lagrange_derivs
const double all_lagrange_derivs[][11][11]
Definition: speckley/src/WaveAssembler3D.cpp:35
escript::AbstractSystemMatrix::getRowBlockSize
int getRowBlockSize() const
returns the row block size
Definition: AbstractSystemMatrix.h:128
speckley::DataMap
std::map< std::string, escript::Data > DataMap
Definition: speckley/src/domainhelpers.h:39
escript::AbstractSystemMatrix::getColumnBlockSize
int getColumnBlockSize() const
returns the column block size
Definition: AbstractSystemMatrix.h:139
escript::const_Domain_ptr
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:58
speckley::isNotEmpty
bool isNotEmpty(const std::string target, const DataMap &mapping)
Definition: speckley/src/domainhelpers.h:54
speckley::WaveAssembler3D::m_dx
const double * m_dx
Definition: speckley/src/WaveAssembler3D.h:207
escript::ValueError
An exception class that signals an invalid argument value.
Definition: EsysException.h:103
ripley::WaveAssembler3D::c44
escript::Data c44
Definition: ripley/src/WaveAssembler3D.h:105
speckley::WaveAssembler3D::c12
escript::Data c12
Definition: speckley/src/WaveAssembler3D.h:211
speckley::SpeckleyException
SpeckleyException exception class.
Definition: SpeckleyException.h:44
ripley::WaveAssembler3D::collateFunctionSpaceTypes
void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.cpp:70
speckley::WaveAssembler3D::assemblePDESystemReduced
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
speckley::WaveAssembler3D::assemblePDESingleReduced
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
ripley::WaveAssembler3D::m_dx
const double * m_dx
Definition: ripley/src/WaveAssembler3D.h:102
ripley::WaveAssembler3D::assemblePDEBoundarySingle
void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.h:71
ripley::WaveAssembler3D::assemblePDESystemReduced
void assemblePDESystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const DataMap &coefs) const
Definition: ripley/src/WaveAssembler3D.h:87