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