66 {
68
70 if (
bit->getName().compare(0, 3,
"ROD") == 0) {
71
72 const int id =
bit->getMeshsetId();
76
77 std::vector<double> attributes;
78 bit->getAttributes(attributes);
79 if (attributes.size() < 3) {
80 SETERRQ(
82 "Input mesh for ROD should have 3 attributes but there is %ld",
83 attributes.size());
84 }
89
90
91 CHKERR PetscPrintf(PETSC_COMM_WORLD,
"\nSimple rod block %d\n",
id);
92 CHKERR PetscPrintf(PETSC_COMM_WORLD,
"\tYoung's modulus %3.4g\n",
93 attributes[0]);
94 CHKERR PetscPrintf(PETSC_COMM_WORLD,
"\tCross-section area %3.4g\n",
95 attributes[1]);
96 CHKERR PetscPrintf(PETSC_COMM_WORLD,
"\tPrestress %3.4g\n",
97 attributes[2]);
98 }
99 }
100
102 }
@ MOFEM_ATOM_TEST_INVALID
#define CHKERR
Inline error check.
#define _IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(MESHSET_MANAGER, CUBITBCTYPE, IT)
Iterator that loops over a specific Cubit MeshSet having a particular BC meshset in a moFEM field.
std::map< int, BlockOptionDataSimpleRods > mapSimpleRod
virtual moab::Interface & get_moab()=0