Function getBestQueryGenerator

src/Initializer/ParameterDB.cpp:651–682

}

QueryGenerator* getBestQueryGenerator(bool anelasticity,
                                      bool plasticity,
                                      bool anisotropy,
                                      bool poroelasticity,
                                      bool useCellHomogenizedMaterial,
                                      const CellToVertexArray& cellToVertex) {
  QueryGenerator* queryGen = nullptr;
  if (!useCellHomogenizedMaterial) {
    queryGen = new ElementBarycentreGenerator(cellToVertex);
  } else {
    const auto rank = MPI::mpi.rank();
    if (anisotropy) {
      logWarning(rank)
          << "Material Averaging is not implemented for anisotropic materials. Falling back to "
             "material properties sampled from the element barycenters instead.";
      queryGen = new ElementBarycentreGenerator(cellToVertex);
    } else if (plasticity) {
      logWarning(rank)
          << "Material Averaging is not implemented for plastic materials. Falling back to "
             "material properties sampled from the element barycenters instead.";
      queryGen = new ElementBarycentreGenerator(cellToVertex);
    } else if (poroelasticity) {
      logWarning(rank)
          << "Material Averaging is not implemented for poroelastic materials. Falling back to "
             "material properties sampled from the element barycenters instead.";
      queryGen = new ElementBarycentreGenerator(cellToVertex);
    } else {
      queryGen = new ElementAverageGenerator(cellToVertex);
    }
  }
  return queryGen;
}

template class MaterialParameterDB<seissol::model::AnisotropicMaterial>;
template class MaterialParameterDB<seissol::model::ElasticMaterial>;