escript  Revision_
speckley/src/WaveAssembler2D.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_2D_H__
17 #define __SPECKLEY_WAVE_ASSEMBLER_2D_H__
18 
19 #include <speckley/Rectangle.h>
20 
21 namespace speckley {
22 
23 
25 {
26 public:
28  const dim_t *NE, const dim_t *NN, const DataMap& c)
30  m_dx(dx),
31  m_NE(NE),
32  m_NN(NN)
33  {
34  domain = REFCOUNTNS::static_pointer_cast<const Rectangle>(dom);
35  isHTI = isVTI = false;
36  DataMap::const_iterator a = c.find("c12"), b = c.find("c23");
37  if (c.find("c11") == c.end()
38  || c.find("c13") == c.end() || c.find("c33") == c.end()
39  || c.find("c44") == c.end() || c.find("c66") == c.end()
40  || (a == c.end() && b == c.end()))
41  throw SpeckleyException("required constants missing for WaveAssembler");
42 
43  if (a != c.end() && b != c.end()) {
44  throw SpeckleyException("WaveAssembler3D() doesn't support general form waves");
45  } else if (a == c.end()) {
46  c23 = b->second;
47  isHTI = true;
49  throw SpeckleyException("C tensor elements must be reduced");
50  }
51  if (c23.isEmpty()) {
52  throw SpeckleyException("C tensor elements must not be empty");
53  }
54  } else if (b == c.end()) {
55  c12 = a->second;
56  isVTI = true;
58  throw SpeckleyException("C tensor elements must be reduced");
59  }
60  if (c12.isEmpty()) {
61  throw SpeckleyException("C tensor elements must not be empty");
62  }
63  } // final else case taken care of with the missing constants above
64  c11 = c.find("c11")->second;
65  c13 = c.find("c13")->second;
66  c33 = c.find("c33")->second;
67  c44 = c.find("c44")->second;
68  c66 = c.find("c66")->second;
74  throw SpeckleyException("C tensor elements must be reduced");
75  }
76  if (c11.isEmpty()
77  || c13.isEmpty()
78  || c33.isEmpty()
79  || c44.isEmpty()
80  || c66.isEmpty()) {
81  throw SpeckleyException("C tensor elements must not be empty");
82  }
83  }
84 
86 
87  /* The default SpeckleyDomain assemblers, with original signatures */
88 
92  const escript::Data& A, const escript::Data& B,
93  const escript::Data& C, const escript::Data& D,
94  const escript::Data& du, const escript::Data& Y) const;
95 
99  escript::Data& rhs, const escript::Data& d,
100  const escript::Data& y) const;
101 
105  escript::Data& rhs, const escript::Data& A, const escript::Data& B,
106  const escript::Data& C, const escript::Data& D,
107  const escript::Data& du, const escript::Data& Y) const;
108 
112  escript::Data& rhs, const escript::Data& d,
113  const escript::Data& y) const;
114 
118  const escript::Data& A, const escript::Data& B,
119  const escript::Data& C, const escript::Data& D,
120  const escript::Data& du, const escript::Data& Y) const;
121 
125  escript::Data& rhs, const escript::Data& d,
126  const escript::Data& y) const;
127 
131  escript::Data& rhs, const escript::Data& A, const escript::Data& B,
132  const escript::Data& C, const escript::Data& D,
133  const escript::Data& du, const escript::Data& Y) const;
134 
138  escript::Data& rhs, const escript::Data& d,
139  const escript::Data& y) const;
140 
141  /* The new interface for assemblers */
142 
144  escript::Data& rhs,
145  const DataMap& coefs) const;
147  escript::Data& rhs,
148  const DataMap& coefs) const;
150  escript::Data& rhs,
151  const DataMap& coefs) const;
154  escript::Data& rhs,
155  const DataMap& coefs) const;
157  escript::Data& rhs,
158  const DataMap& coefs) const;
160  escript::Data& rhs,
161  const DataMap& coefs) const;
163  escript::Data& rhs,
164  const DataMap& coefs) const;
167  escript::Data& rhs,
168  const DataMap& coefs) const;
169 
170  virtual void collateFunctionSpaceTypes(std::vector<int>& fsTypes,
171  const DataMap& coefs) const;
172 
173 protected:
174  POINTER_WRAPPER_CLASS(const Rectangle) domain;
175  const double *m_dx;
176  const dim_t *m_NE;
177  const dim_t *m_NN;
178  bool isHTI, isVTI;
180 };
181 
182 } // namespace speckley
183 
184 #endif // __SPECKLEY_DEFAULTASSEMBLER2D_H__
185 
Definition: AbstractAssembler.cpp:18
Definition: Speckley.h:52
virtual void assemblePDEBoundarySingleReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:383
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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:404
escript::Data c23
Definition: speckley/src/WaveAssembler2D.h:179
escript::Data c44
Definition: speckley/src/WaveAssembler2D.h:179
~WaveAssembler2D()
Definition: speckley/src/WaveAssembler2D.h:85
bool isVTI
Definition: speckley/src/WaveAssembler2D.h:178
std::map< std::string, escript::Data > DataMap
Definition: speckley/src/domainhelpers.h:24
escript::Data c12
Definition: speckley/src/WaveAssembler2D.h:179
const dim_t * m_NE
Definition: speckley/src/WaveAssembler2D.h:176
bool isHTI
Definition: speckley/src/WaveAssembler2D.h:178
escript::Data c66
Definition: speckley/src/WaveAssembler2D.h:179
virtual void assemblePDEBoundarySingle(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:377
POINTER_WRAPPER_CLASS(const Rectangle) domain
Definition: speckley/src/WaveAssembler2D.h:24
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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:413
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:1015
Rectangle is the 2-dimensional implementation of a SpeckleyDomain.
Definition: speckley/src/Rectangle.h:32
virtual void assemblePDEBoundarySystem(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:390
virtual void collateFunctionSpaceTypes(std::vector< int > &fsTypes, const DataMap &coefs) const
Definition: speckley/src/WaveAssembler2D.cpp:115
const double * m_dx
Definition: speckley/src/WaveAssembler2D.h:175
const FunctionSpace & getFunctionSpace() const
Return the function space.
Definition: Data.h:452
escript::Data c11
Definition: speckley/src/WaveAssembler2D.h:179
Data represents a collection of datapoints.
Definition: Data.h:63
int getTypeCode() const
Returns the function space type code.
Definition: FunctionSpace.cpp:94
WaveAssembler2D(escript::const_Domain_ptr dom, const double *dx, const dim_t *NE, const dim_t *NN, const DataMap &c)
Definition: speckley/src/WaveAssembler2D.h:27
SpeckleyException exception class.
Definition: SpeckleyException.h:29
Definition: speckley/src/AbstractAssembler.h:33
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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:368
escript::Data c13
Definition: speckley/src/WaveAssembler2D.h:179
Base class for escript system matrices.
Definition: AbstractSystemMatrix.h:42
const dim_t * m_NN
Definition: speckley/src/WaveAssembler2D.h:177
escript::Data c33
Definition: speckley/src/WaveAssembler2D.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 &du, const escript::Data &Y) const
Definition: speckley/src/WaveAssembler2D.cpp:211
boost::shared_ptr< const AbstractDomain > const_Domain_ptr
Definition: AbstractDomain.h:39
virtual void assemblePDEBoundarySystemReduced(escript::AbstractSystemMatrix *mat, escript::Data &rhs, const escript::Data &d, const escript::Data &y) const
Definition: speckley/src/WaveAssembler2D.cpp:396
index_t dim_t
Definition: DataTypes.h:64