escript  Revision_
speckley/src/WaveAssembler3D.h
Go to the documentation of this file.
1 
2 /*****************************************************************************
3 *
4 * Copyright (c) 2003-2016 by The University of Queensland
5 * http://www.uq.edu.au
6 *
7 * Primary Business: Queensland, Australia
8 * Licensed under the Apache License, version 2.0
9 * http://www.apache.org/licenses/LICENSE-2.0
10 *
11 * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12 * Development 2012-2013 by School of Earth Sciences
13 * Development from 2014 by Centre for Geoscience Computing (GeoComp)
14 *
15 *****************************************************************************/
16 #ifndef __SPECKLEY_WAVE_ASSEMBLER_3D_H__
17 #define __SPECKLEY_WAVE_ASSEMBLER_3D_H__
18 
19 #include <map>
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>
25 
26 namespace speckley {
27 
28 
30 {
31 public:
33  const dim_t *NE, const dim_t *NN, const DataMap& c)
35  m_dx(dx),
36  m_NE(NE),
37  m_NN(NN)
38  {
39  domain = boost::static_pointer_cast<const Brick>(dom);
40  isHTI = isVTI = false;
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()))
46  throw SpeckleyException("required constants missing for WaveAssembler");
47 
48  if (a != c.end() && b != c.end()) {
49  throw SpeckleyException("WaveAssembler3D() doesn't support general form waves");
50  } else if (a == c.end()) {
51  c23 = b->second;
52  isHTI = true;
54  throw SpeckleyException("C tensor elements must be reduced");
55  }
56  if (c23.isEmpty()) {
57  throw SpeckleyException("C tensor elements must not be empty");
58  }
59  } else if (b == c.end()) {
60  c12 = a->second;
61  isVTI = true;
63  throw SpeckleyException("C tensor elements must be reduced");
64  }
65  if (c12.isEmpty()) {
66  throw SpeckleyException("C tensor elements must not be empty");
67  }
68  } // final else case taken care of with the missing constants above
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;
79  throw SpeckleyException("C tensor elements must be reduced");
80  }
81  if (c11.isEmpty()
82  || c13.isEmpty()
83  || c33.isEmpty()
84  || c44.isEmpty()
85  || c66.isEmpty()) {
86  throw SpeckleyException("C tensor elements must not be empty");
87  }
88  }
89 
91 
92  /* The default SpeckleyDomain assemblers, with original signatures */
93 
97  escript::Data& rhs, const escript::Data& A,
98  const escript::Data& B, const escript::Data& C,
99  const escript::Data& D, const escript::Data& X,
100  const escript::Data& Y) const;
101 
105  escript::Data& rhs, const escript::Data& d,
106  const escript::Data& y) const;
107 
111  escript::Data& rhs, const escript::Data& A, const escript::Data& B,
112  const escript::Data& C, const escript::Data& D,
113  const escript::Data& X, const escript::Data& Y) const;
114 
119  escript::Data& rhs, const escript::Data& d,
120  const escript::Data& y) const;
121 
125  escript::Data& rhs, const escript::Data& A,
126  const escript::Data& B, const escript::Data& C,
127  const escript::Data& D, const escript::Data& X,
128  const escript::Data& Y) const;
129 
133  escript::Data& rhs, const escript::Data& d,
134  const escript::Data& y) const;
135 
139  escript::Data& rhs, const escript::Data& A, const escript::Data& B,
140  const escript::Data& C, const escript::Data& D,
141  const escript::Data& X, const escript::Data& Y) const;
142 
147  escript::Data& rhs, const escript::Data& d,
148  const escript::Data& y) const;
149 
150  /* The new interface for assemblers */
151 
153  escript::Data& rhs, const DataMap& coefs) const;
155  escript::Data& rhs, const DataMap& coefs) const;
157  escript::Data& rhs, const DataMap& coefs) const;
160  escript::Data& rhs, const DataMap& coefs) const;
162  escript::Data& rhs, const DataMap& coefs) const;
164  escript::Data& rhs, const DataMap& coefs) const;
166  escript::Data& rhs, const DataMap& coefs) const;
169  escript::Data& rhs, const DataMap& coefs) const;
170 
171  void collateFunctionSpaceTypes(std::vector<int>& fsTypes,
172  const DataMap& coefs) const;
173 
174 protected:
175  boost::shared_ptr<const Brick> domain;
176  const double *m_dx;
177  const dim_t *m_NE;
178  const dim_t *m_NN;
179  bool isHTI, isVTI;
181 };
182 
183 } // namespace speckley
184 
185 #endif // __SPECKLEY_WAVE_ASSEMBLER_3D_H__
186 
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