v0.11.0
Public Member Functions | Public Attributes | List of all members
CommonData Struct Reference

#include <users_modules/topology_optimization/src/Header.hpp>

Collaboration diagram for CommonData:
[legend]

Public Member Functions

 CommonData (MoFEM::Interface &m_field)
 
MoFEMErrorCode getParameters ()
 
MoFEMErrorCode getSensitivityFromTag (const EntityHandle fe_ent, const int nb_gauss_pts)
 
MoFEMErrorCode setSensitivityToTag (const EntityHandle fe_ent, const double sens)
 
MoFEMErrorCode getDensityFromTag (const EntityHandle fe_ent, const int nb_gauss_pts)
 
MoFEMErrorCode setDensityToTag (const EntityHandle fe_ent, const double rho)
 
MoFEMErrorCode createTags (moab::Interface &moab)
 
MoFEMErrorCode build_main_dm (DM &dm, DM &sub_dm_elastic, DM &sub_dm_rho, MoFEM::Interface &m_field, BitRefLevel &bit_level0, PetscBool is_partitioned)
 
MoFEMErrorCode build_sub_dms (DM &dm, DM &sub_dm_elastic, DM &sub_dm_rho, MoFEM::Interface &m_field, BitRefLevel &bit_level0, PetscBool is_partitioned)
 
MoFEMErrorCode getBoundBox (const double material_coords[3], std::vector< double > &coords, double &boxSize)
 

Public Attributes

MatrixDouble dotEdge
 
MatrixDouble dotEleLeft
 
MatrixDouble dotEleRight
 
MoFEM::Interface & mField
 
Tag thSensitivity
 
Tag thDens
 
Vec massVec
 
Vec changeVec
 
boost::shared_ptr< VectorDouble > dataRhotGaussPtr
 
boost::shared_ptr< MatrixDouble > dispGradAtGaussPts
 
MatrixDouble D
 
VectorDouble sTrain
 
vector< VectorDouble > sTress
 
double volumeFrac
 
MatrixDouble cOmplianceMat
 
VectorDouble sEnsitivityVec
 
VectorDouble rhoVecAtTag
 
MatrixDouble anisLambdaMat
 
double edgeLength
 
double lAmbda
 
double mU
 
double lAmbda0
 
double rHo
 
double yOung
 
double pOisson
 
int penalty
 
double mOve
 
double fIlter_radius
 
double cHange
 
double cUrrent_mass
 
double L_mid
 
bool uPdate_density
 

Detailed Description

Examples
ContactOps.hpp, NavierStokesElement.hpp, PlasticOps.hpp, Remodeling.hpp, SmallStrainPlasticity.hpp, continuity_check_on_contact_prism_side_ele.cpp, continuity_check_on_skeleton_3d.cpp, continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mixed_poisson.cpp, and reaction_diffusion.cpp.

Definition at line 40 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

Constructor & Destructor Documentation

◆ CommonData()

CommonData::CommonData ( MoFEM::Interface &  m_field)

Definition at line 56 of file Header.hpp.

56  : mField(m_field) {
57 
58  dataRhotGaussPtr = boost::shared_ptr<VectorDouble>(new VectorDouble());
59  dispGradAtGaussPts = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
60  volumeFrac = 0.5;
61  yOung = 10;
62  fIlter_radius = 0;
63  penalty = 3;
64  cHange = 0;
65  edgeLength = 1.;
66  lAmbda0 = 0.;
67  cUrrent_mass = 0;
68  L_mid = 0; // langange multiplier
69  uPdate_density = false;
70 
71  anisLambdaMat.resize(3,3,false);
72  anisLambdaMat.clear();
73 
74  ierr = createTags(m_field.get_moab());
75  }
MoFEM::Exceptions::ierr
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
Definition: Exceptions.hpp:87
MoFEM::Types::VectorDouble
ublas::vector< double, DoubleAllocator > VectorDouble
Definition: Types.hpp:74
MoFEM::Types::MatrixDouble
ublas::matrix< double, ublas::row_major, DoubleAllocator > MatrixDouble
Definition: Types.hpp:76
CommonData::dataRhotGaussPtr
boost::shared_ptr< VectorDouble > dataRhotGaussPtr
Definition: Header.hpp:28
CommonData::mField
MoFEM::Interface & mField
Definition: Header.hpp:21
CommonData::volumeFrac
double volumeFrac
Definition: Header.hpp:34
CommonData::L_mid
double L_mid
Definition: Header.hpp:53
CommonData::anisLambdaMat
MatrixDouble anisLambdaMat
Definition: Header.hpp:38
CommonData::yOung
double yOung
Definition: Header.hpp:45
CommonData::cHange
double cHange
Definition: Header.hpp:50
CommonData::uPdate_density
bool uPdate_density
Definition: Header.hpp:54
CommonData::dispGradAtGaussPts
boost::shared_ptr< MatrixDouble > dispGradAtGaussPts
Definition: Header.hpp:29
CommonData::cUrrent_mass
double cUrrent_mass
Definition: Header.hpp:52
CommonData::lAmbda0
double lAmbda0
Definition: Header.hpp:43
CommonData::fIlter_radius
double fIlter_radius
Definition: Header.hpp:49
CommonData::penalty
int penalty
Definition: Header.hpp:47
CommonData::edgeLength
double edgeLength
Definition: Header.hpp:40
CommonData::createTags
MoFEMErrorCode createTags(moab::Interface &moab)
Definition: Header.hpp:232
MoFEM::CoreInterface::get_moab
virtual moab::Interface & get_moab()=0

Member Function Documentation

◆ build_main_dm()

MoFEMErrorCode CommonData::build_main_dm ( DM &  dm,
DM &  sub_dm_elastic,
DM &  sub_dm_rho,
MoFEM::Interface &  m_field,
BitRefLevel &  bit_level0,
PetscBool  is_partitioned 
)

Definition at line 244 of file Header.hpp.

246  {
247  MoFEMFunctionBegin;
248  DMType dm_name = "DM_TOPO";
249  // Register DM problem
250  CHKERR DMRegister_MoFEM(dm_name);
251  CHKERR DMCreate(PETSC_COMM_WORLD, &dm);
252  CHKERR DMSetType(dm, dm_name);
253  CHKERR DMMoFEMCreateMoFEM(dm, &m_field, dm_name, bit_level0);
254  CHKERR DMSetFromOptions(dm);
255  CHKERR DMMoFEMAddElement(dm, "ELASTIC");
256  CHKERR DMMoFEMAddElement(dm, "RHO_APPROX");
257  if (m_field.check_finite_element("PRESSURE_FE"))
258  CHKERR DMMoFEMAddElement(dm, "PRESSURE_FE");
259  if (m_field.check_finite_element("FORCE_FE"))
260  CHKERR DMMoFEMAddElement(dm, "FORCE_FE");
261  CHKERR DMMoFEMSetIsPartitioned(dm, is_partitioned);
262  m_field.getInterface<ProblemsManager>()->buildProblemFromFields = PETSC_TRUE;
263  CHKERR DMSetUp(dm);
264  m_field.getInterface<ProblemsManager>()->buildProblemFromFields = PETSC_FALSE;
265 
266 
267  MoFEMFunctionReturn(0);
268  }
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:415
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:485
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:604
MoFEM::DMMoFEMSetIsPartitioned
PetscErrorCode DMMoFEMSetIsPartitioned(DM dm, PetscBool is_partitioned)
Definition: DMMMoFEM.cpp:982
MoFEM::DMMoFEMCreateMoFEM
PetscErrorCode DMMoFEMCreateMoFEM(DM dm, MoFEM::Interface *m_field_ptr, const char problem_name[], const MoFEM::BitRefLevel bit_level, const MoFEM::BitRefLevel bit_mask=MoFEM::BitRefLevel().set())
Must be called by user to set MoFEM data structures.
Definition: DMMMoFEM.cpp:105
MoFEM::DMMoFEMAddElement
PetscErrorCode DMMoFEMAddElement(DM dm, const char fe_name[])
add element to dm
Definition: DMMMoFEM.cpp:425
MoFEM::DMRegister_MoFEM
PetscErrorCode DMRegister_MoFEM(const char sname[])
Register MoFEM problem.
Definition: DMMMoFEM.cpp:48
MoFEM::CoreInterface::check_finite_element
virtual bool check_finite_element(const std::string &name) const =0
Check if finite element is in database.
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(const MOFEMuuid &uuid, IFACE *&iface) const
Get interface by uuid and return reference to pointer of interface.
Definition: UnknownInterface.hpp:130

◆ build_sub_dms()

MoFEMErrorCode CommonData::build_sub_dms ( DM &  dm,
DM &  sub_dm_elastic,
DM &  sub_dm_rho,
MoFEM::Interface &  m_field,
BitRefLevel &  bit_level0,
PetscBool  is_partitioned 
)

Definition at line 270 of file Header.hpp.

273  {
274 
275  MoFEMFunctionBegin;
276  DMType dm_name_sub_elast = "DM_ELASTIC";
277  CHKERR DMRegister_MoFEM(dm_name_sub_elast);
278  CHKERR DMCreate(PETSC_COMM_WORLD, &sub_dm_elastic);
279  CHKERR DMSetType(sub_dm_elastic, dm_name_sub_elast);
280  // m_field.getInterface<ProblemsManager>()->buildProblemFromFields =
281  // PETSC_TRUE;
282  CHKERR DMMoFEMCreateSubDM(sub_dm_elastic, dm, "ELASTIC_PROB");
283  CHKERR DMMoFEMSetDestroyProblem(sub_dm_elastic, PETSC_TRUE);
284  CHKERR DMMoFEMAddSubFieldRow(sub_dm_elastic, "DISPLACEMENTS");
285  CHKERR DMMoFEMAddSubFieldCol(sub_dm_elastic, "DISPLACEMENTS");
286  CHKERR DMMoFEMAddElement(sub_dm_elastic, "ELASTIC");
287  if (m_field.check_finite_element("PRESSURE_FE"))
288  CHKERR DMMoFEMAddElement(sub_dm_elastic, "PRESSURE_FE");
289  if (m_field.check_finite_element("FORCE_FE"))
290  CHKERR DMMoFEMAddElement(sub_dm_elastic, "FORCE_FE");
291 
292  CHKERR DMMoFEMSetSquareProblem(sub_dm_elastic, PETSC_TRUE);
293  CHKERR DMMoFEMSetIsPartitioned(sub_dm_elastic, is_partitioned);
294  CHKERR DMSetUp(sub_dm_elastic);
295 
296  DMType dm_name_sub_rho = "DM_DENSITY";
297  CHKERR DMRegister_MoFEM(dm_name_sub_rho);
298  CHKERR DMCreate(PETSC_COMM_WORLD, &sub_dm_rho);
299  CHKERR DMSetType(sub_dm_rho, dm_name_sub_rho);
300  // m_field.getInterface<ProblemsManager>()->buildProblemFromFields =
301  // PETSC_TRUE;
302  CHKERR DMMoFEMCreateSubDM(sub_dm_rho, dm, "DENSITY_PROB");
303  CHKERR DMMoFEMSetDestroyProblem(sub_dm_rho, PETSC_TRUE);
304  CHKERR DMMoFEMAddSubFieldRow(sub_dm_rho, "RHO");
305  CHKERR DMMoFEMAddSubFieldCol(sub_dm_rho, "RHO");
306  CHKERR DMMoFEMAddElement(sub_dm_rho, "RHO_APPROX");
307  CHKERR DMMoFEMSetSquareProblem(sub_dm_rho, PETSC_TRUE);
308  CHKERR DMMoFEMSetIsPartitioned(sub_dm_rho, is_partitioned);
309  CHKERR DMSetUp(sub_dm_rho);
310  MoFEMFunctionReturn(0);
311  }
MoFEM::DMMoFEMCreateSubDM
PetscErrorCode DMMoFEMCreateSubDM(DM subdm, DM dm, const char problem_name[])
Must be called by user to set Sub DM MoFEM data structures.
Definition: DMMMoFEM.cpp:177
MoFEM::DMMoFEMSetSquareProblem
PetscErrorCode DMMoFEMSetSquareProblem(DM dm, PetscBool square_problem)
set squared problem
Definition: DMMMoFEM.cpp:378
MoFEM::DMMoFEMAddSubFieldRow
PetscErrorCode DMMoFEMAddSubFieldRow(DM dm, const char field_name[], EntityType lo_type=MBVERTEX, EntityType hi_type=MBMAXTYPE)
Definition: DMMMoFEM.cpp:198
MoFEM::DMMoFEMAddSubFieldCol
PetscErrorCode DMMoFEMAddSubFieldCol(DM dm, const char field_name[], EntityType lo_type=MBVERTEX, EntityType hi_type=MBMAXTYPE)
Definition: DMMMoFEM.cpp:221
MoFEM::DMMoFEMSetDestroyProblem
PetscErrorCode DMMoFEMSetDestroyProblem(DM dm, PetscBool destroy_problem)
Definition: DMMMoFEM.cpp:360

◆ createTags()

MoFEMErrorCode CommonData::createTags ( moab::Interface &  moab)

Definition at line 232 of file Header.hpp.

232  {
233  MoFEMFunctionBegin;
234  std::vector<double> def_val(9, 0);
235  CHKERR mField.get_moab().tag_get_handle(
236  "SENSITIVITY", 1, MB_TYPE_DOUBLE, thSensitivity,
237  MB_TAG_CREAT | MB_TAG_VARLEN | MB_TAG_SPARSE, &def_val[0]);
238  CHKERR mField.get_moab().tag_get_handle(
239  "RHO", 1, MB_TYPE_DOUBLE, thDens,
240  MB_TAG_CREAT | MB_TAG_VARLEN | MB_TAG_SPARSE, &def_val[0]);
241  MoFEMFunctionReturn(0);
242  }
CommonData::thDens
Tag thDens
Definition: Header.hpp:23
CommonData::thSensitivity
Tag thSensitivity
Definition: Header.hpp:22

◆ getBoundBox()

MoFEMErrorCode CommonData::getBoundBox ( const double  material_coords[3],
std::vector< double > &  coords,
double &  boxSize 
)

Definition at line 313 of file Header.hpp.

314  {
315  MoFEMFunctionBeginHot;
316  FTensor::Tensor1<double, 3> t_material_coords(
317  material_coords[0], material_coords[1], material_coords[2]);
318  FTensor::Tensor1<double *, 3> t_coords(&coords[0], &coords[1], &coords[2],
319  3);
320  FTensor::Index<'i', 3> i;
321  double max_xyz[] = {0, 0, 0};
322  for (int bb = 0; bb != coords.size() / 3; bb++) {
323  t_coords(i) -= t_material_coords(i);
324  for (int dd = 0; dd != 3; dd++) {
325  max_xyz[dd] = (max_xyz[dd] < fabs(t_coords(dd))) ? fabs(t_coords(dd))
326  : max_xyz[dd];
327  }
328  ++t_coords;
329  }
330  //take the max of the box
331  boxSize = (max_xyz[0] > max_xyz[1]) ? max_xyz[0] : max_xyz[1];
332  boxSize = (boxSize > max_xyz[2]) ? boxSize : max_xyz[2];
333  MoFEMFunctionReturnHot(0);
334  }
MoFEMFunctionReturnHot
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:516
MoFEMFunctionBeginHot
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:509
i
FTensor::Index< 'i', 3 > i
Definition: matrix_function.cpp:18
FTensor::dd
const Tensor2_symmetric_Expr< const ddTensor0< T, Dim, i, j >, typename promote< T, double >::V, Dim, i, j > dd(const Tensor0< T * > &a, const Index< i, Dim > index1, const Index< j, Dim > index2, const Tensor1< int, Dim > &d_ijk, const Tensor1< double, Dim > &d_xyz)
Definition: ddTensor0.hpp:33

◆ getDensityFromTag()

MoFEMErrorCode CommonData::getDensityFromTag ( const EntityHandle  fe_ent,
const int  nb_gauss_pts 
)

Definition at line 198 of file Header.hpp.

198  {
199  MoFEMFunctionBegin;
200  double *tag_data;
201  int tag_size;
202  CHKERR mField.get_moab().tag_get_by_ptr(
203  thDens, &fe_ent, 1, (const void **)&tag_data, &tag_size);
204 
205  if (tag_size == 1) {
206  rhoVecAtTag.resize(nb_gauss_pts , false);
207  rhoVecAtTag.clear();
208  void const *tag_data[] = {&*rhoVecAtTag.data().begin()};
209  const int tag_size = rhoVecAtTag.data().size();
210  CHKERR mField.get_moab().tag_set_by_ptr(thDens, &fe_ent, 1,
211  tag_data, &tag_size);
212  } else if (tag_size != nb_gauss_pts) {
213  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
214  "Wrong size of the tag, data inconsistency");
215  } else {
216  VectorAdaptor tag_vec = VectorAdaptor(
217  tag_size, ublas::shallow_array_adaptor<double>(tag_size, tag_data));
218  rhoVecAtTag = tag_vec;
219  }
220  MoFEMFunctionReturn(0);
221  }
MOFEM_DATA_INCONSISTENCY
@ MOFEM_DATA_INCONSISTENCY
Definition: definitions.h:123
MoFEM::Types::VectorAdaptor
VectorShallowArrayAdaptor< double > VectorAdaptor
Definition: Types.hpp:109
CommonData::rhoVecAtTag
VectorDouble rhoVecAtTag
Definition: Header.hpp:37

◆ getParameters()

MoFEMErrorCode CommonData::getParameters ( )

Definition at line 76 of file Header.hpp.

76  {
77  MoFEMFunctionBegin;
78  ierr = PetscOptionsBegin(PETSC_COMM_WORLD, "", "Topo Config", "none");
79  CHKERR(ierr);
80 
81  CHKERR PetscOptionsInt("-penalty", "default penalty ", "", penalty, &penalty,
82  PETSC_NULL);
83 
84  CHKERR PetscOptionsScalar("-young", "default yOung modulus", "", yOung, &yOung,
85  PETSC_NULL);
86 
87  CHKERR PetscOptionsScalar("-poisson", "poisson ratio", "", 0.1, &pOisson,
88  PETSC_NULL);
89 
90  CHKERR PetscOptionsScalar("-move", "max move", "", 0.1, &mOve, PETSC_NULL);
91  // CHKERR PetscOptionsScalar("-lambda", "lambda smooth", "", lAmbda0, &lAmbda0, PETSC_NULL);
92 
93  double lambda_coef[3] = {0.,0.,0.};
94  int nmax = 3;
95  PetscBool is_lambda_defined = PETSC_FALSE;
96  PetscOptionsGetRealArray(NULL, "", "-lambda", lambda_coef,
97  &nmax, &is_lambda_defined);
98 
99  switch (nmax) {
100 
101  case 0:
102  case 1:
103  anisLambdaMat(0, 0) = anisLambdaMat(1, 1) = anisLambdaMat(2, 2) =
104  lambda_coef[0];
105  break;
106  case 3:
107  anisLambdaMat(0, 0) = lambda_coef[0];
108  anisLambdaMat(1, 1) = lambda_coef[1];
109  anisLambdaMat(2, 2) = lambda_coef[2];
110  break;
111  default:
112 
113  SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
114  "Wrong number of lambda parameters, should be 1 or 3, data "
115  "inconsistency");
116  }
117 
118  CHKERR PetscOptionsScalar("-volume_frac", "volume frac", "", volumeFrac,
119  &volumeFrac, PETSC_NULL);
120 
121  // CHKERR PetscOptionsScalar("-filter_r", "default filter radius", "", 0,
122  // &fIlter_radius, PETSC_NULL);
123  ierr = PetscOptionsEnd();
124  CHKERRG(ierr);
125 
126  lAmbda = (yOung * pOisson) / ((1. + pOisson) * (1. - 2. * pOisson));
127  mU = yOung / (2. * (1. + pOisson));
128 
129  const double coefficient = yOung / ((1 + pOisson) * (1 - 2 * pOisson));
130  D.resize(36, 1, false);
131  D.clear();
132 
133  FTensor::Ddg<FTensor::PackPtr<double *, 1>, 3, 3> t_D(MAT_TO_DDG(&D));
134 
135  FTensor::Index<'i', 3> i;
136  FTensor::Index<'j', 3> j;
137  FTensor::Index<'k', 3> k;
138  FTensor::Index<'l', 3> l;
139 
140  t_D(i, j, k, l) = 0.;
141 
142  t_D(0, 0, 0, 0) = 1 - pOisson;
143  t_D(1, 1, 1, 1) = 1 - pOisson;
144  t_D(2, 2, 2, 2) = 1 - pOisson;
145 
146  t_D(0, 1, 0, 1) = 0.5 * (1 - 2 * pOisson);
147  t_D(0, 2, 0, 2) = 0.5 * (1 - 2 * pOisson);
148  t_D(1, 2, 1, 2) = 0.5 * (1 - 2 * pOisson);
149 
150  t_D(0, 0, 1, 1) = pOisson;
151  t_D(1, 1, 0, 0) = pOisson;
152  t_D(0, 0, 2, 2) = pOisson;
153  t_D(2, 2, 0, 0) = pOisson;
154  t_D(1, 1, 2, 2) = pOisson;
155  t_D(2, 2, 1, 1) = pOisson;
156 
157  t_D(i, j, k, l) *= coefficient;
158 
159  MoFEMFunctionReturn(0);
160  }
MAT_TO_DDG
#define MAT_TO_DDG(SM)
Definition: HookeElement.hpp:156
CHKERRG
#define CHKERRG(n)
Check error code of MoFEM/MOAB/PETSc function.
Definition: definitions.h:552
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
MoFEM::PetscOptionsGetRealArray
PetscErrorCode PetscOptionsGetRealArray(PetscOptions *, const char pre[], const char name[], PetscReal dval[], PetscInt *nmax, PetscBool *set)
Definition: DeprecatedPetsc.hpp:203
CommonData::lAmbda
double lAmbda
Definition: Header.hpp:41
CommonData::mOve
double mOve
Definition: Header.hpp:48
CommonData::mU
double mU
Definition: Header.hpp:42
CommonData::pOisson
double pOisson
Definition: Header.hpp:46
CommonData::D
MatrixDouble D
Definition: Header.hpp:31

◆ getSensitivityFromTag()

MoFEMErrorCode CommonData::getSensitivityFromTag ( const EntityHandle  fe_ent,
const int  nb_gauss_pts 
)

Definition at line 163 of file Header.hpp.

163  {
164  MoFEMFunctionBegin;
165  double *tag_data;
166  int tag_size;
167  CHKERR mField.get_moab().tag_get_by_ptr(
168  thSensitivity, &fe_ent, 1, (const void **)&tag_data, &tag_size);
169 
170  if (tag_size == 1) {
171  sEnsitivityVec.resize(nb_gauss_pts , false);
172  sEnsitivityVec.clear();
173  void const *tag_data[] = {&*sEnsitivityVec.data().begin()};
174  const int tag_size = sEnsitivityVec.data().size();
175  CHKERR mField.get_moab().tag_set_by_ptr(thSensitivity, &fe_ent, 1,
176  tag_data, &tag_size);
177  } else if (tag_size != nb_gauss_pts) {
178  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
179  "Wrong size of the tag, data inconsistency");
180  } else {
181  VectorAdaptor tag_vec = VectorAdaptor(
182  tag_size, ublas::shallow_array_adaptor<double>(tag_size, tag_data));
183  sEnsitivityVec = tag_vec;
184  }
185  MoFEMFunctionReturn(0);
186  }
CommonData::sEnsitivityVec
VectorDouble sEnsitivityVec
Definition: Header.hpp:36

◆ setDensityToTag()

MoFEMErrorCode CommonData::setDensityToTag ( const EntityHandle  fe_ent,
const double  rho 
)

Definition at line 223 of file Header.hpp.

223  {
224  MoFEMFunctionBegin;
225  void const *tag_data[] = {&*rhoVecAtTag.data().begin()};
226  const int tag_size = rhoVecAtTag.data().size();
227  CHKERR mField.get_moab().tag_set_by_ptr(thDens, &fe_ent, 1,
228  tag_data, &tag_size);
229  MoFEMFunctionReturn(0);
230  }

◆ setSensitivityToTag()

MoFEMErrorCode CommonData::setSensitivityToTag ( const EntityHandle  fe_ent,
const double  sens 
)

Definition at line 189 of file Header.hpp.

189  {
190  MoFEMFunctionBegin;
191  void const *tag_data[] = {&*sEnsitivityVec.data().begin()};
192  const int tag_size = sEnsitivityVec.data().size();
193  CHKERR mField.get_moab().tag_set_by_ptr(thSensitivity, &fe_ent, 1,
194  tag_data, &tag_size);
195  MoFEMFunctionReturn(0);
196  }

Member Data Documentation

◆ anisLambdaMat

MatrixDouble CommonData::anisLambdaMat

Definition at line 38 of file Header.hpp.

◆ cHange

double CommonData::cHange

Definition at line 50 of file Header.hpp.

◆ changeVec

Vec CommonData::changeVec

Definition at line 26 of file Header.hpp.

◆ cOmplianceMat

MatrixDouble CommonData::cOmplianceMat

Definition at line 35 of file Header.hpp.

◆ cUrrent_mass

double CommonData::cUrrent_mass

Definition at line 52 of file Header.hpp.

◆ D

MatrixDouble CommonData::D

Definition at line 31 of file Header.hpp.

◆ dataRhotGaussPtr

boost::shared_ptr<VectorDouble> CommonData::dataRhotGaussPtr

Definition at line 28 of file Header.hpp.

◆ dispGradAtGaussPts

boost::shared_ptr<MatrixDouble> CommonData::dispGradAtGaussPts

Definition at line 29 of file Header.hpp.

◆ dotEdge

MatrixDouble CommonData::dotEdge
Examples
continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, and continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

Definition at line 41 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

◆ dotEleLeft

MatrixDouble CommonData::dotEleLeft
Examples
continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, and continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

Definition at line 42 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

◆ dotEleRight

MatrixDouble CommonData::dotEleRight
Examples
continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, and continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

Definition at line 43 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

◆ edgeLength

double CommonData::edgeLength

Definition at line 40 of file Header.hpp.

◆ fIlter_radius

double CommonData::fIlter_radius

Definition at line 49 of file Header.hpp.

◆ L_mid

double CommonData::L_mid

Definition at line 53 of file Header.hpp.

◆ lAmbda

double CommonData::lAmbda

Definition at line 41 of file Header.hpp.

◆ lAmbda0

double CommonData::lAmbda0

Definition at line 43 of file Header.hpp.

◆ massVec

Vec CommonData::massVec

Definition at line 25 of file Header.hpp.

◆ mField

MoFEM::Interface& CommonData::mField

Definition at line 21 of file Header.hpp.

◆ mOve

double CommonData::mOve

Definition at line 48 of file Header.hpp.

◆ mU

double CommonData::mU

Definition at line 42 of file Header.hpp.

◆ penalty

int CommonData::penalty

Definition at line 47 of file Header.hpp.

◆ pOisson

double CommonData::pOisson

Definition at line 46 of file Header.hpp.

◆ rHo

double CommonData::rHo

Definition at line 44 of file Header.hpp.

◆ rhoVecAtTag

VectorDouble CommonData::rhoVecAtTag

Definition at line 37 of file Header.hpp.

◆ sEnsitivityVec

VectorDouble CommonData::sEnsitivityVec

Definition at line 36 of file Header.hpp.

◆ sTrain

VectorDouble CommonData::sTrain

Definition at line 32 of file Header.hpp.

◆ sTress

vector<VectorDouble> CommonData::sTress

Definition at line 33 of file Header.hpp.

◆ thDens

Tag CommonData::thDens

Definition at line 23 of file Header.hpp.

◆ thSensitivity

Tag CommonData::thSensitivity

Definition at line 22 of file Header.hpp.

◆ uPdate_density

bool CommonData::uPdate_density

Definition at line 54 of file Header.hpp.

◆ volumeFrac

double CommonData::volumeFrac

Definition at line 34 of file Header.hpp.

◆ yOung

double CommonData::yOung

Definition at line 45 of file Header.hpp.

The documentation for this struct was generated from the following files:
Generated on Sun Apr 11 2021 23:15:53 for MoFEM by Doxygen 1.9.1 and hosted at University of Glasgow