Homec4science
Diffusion tamaas (feature/mpi_fft_engine) Browse

tamaas/src/model/model_template.cpp
0670ef204b92feature/mpi_fft_engine

Last Change
Raw File

model_template.cpp
View Options

/**
* @file
* @section LICENSE
*
* Copyright (©) 2016-2020 EPFL (École Polytechnique Fédérale de Lausanne),
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* This program is free software: you can redistribute it and/or modify
* it under the terms of the GNU Affero General Public License as published
* by the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Affero General Public License for more details.
*
* You should have received a copy of the GNU Affero General Public License
* along with this program. If not, see <https://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#include "model_template.hh"
#include "influence.hh"
#include "partitioner.hh"
#include "westergaard.hh"
/* -------------------------------------------------------------------------- */
namespace tamaas {
/* -------------------------------------------------------------------------- */
template <model_type type>
ModelTemplate<type>::ModelTemplate(std::vector<Real> system_size,
std::vector<UInt> discretization)
: Model(std::move(system_size), std::move(discretization)) {
constexpr UInt dim = trait::dimension;
constexpr UInt dim_b = trait::boundary_dimension;
constexpr UInt nb_components = trait::components;
if (this->system_size.size() != dim)
TAMAAS_EXCEPTION("System size does not match model type");
if (this->discretization.size() != dim)
TAMAAS_EXCEPTION("Discretization size does not match model type");
// Copying sizes for traction grid
std::array<UInt, dim_b> traction_size;
auto disc_it = this->discretization.begin() + (dim > dim_b);
std::copy(disc_it, this->discretization.end(), traction_size.begin());
// Adjust MPI sizes
traction_size = Partitioner<dim_b>::local_size(traction_size);
*disc_it = traction_size.front();
// Allocating
auto traction =
std::make_unique<Grid<Real, dim_b>>(traction_size, nb_components);
auto displacement =
std::make_unique<Grid<Real, dim>>(this->discretization, nb_components);
registerField("traction", std::move(traction));
registerField("displacement", std::move(displacement));
this->initializeBEEngine();
}
/* -------------------------------------------------------------------------- */
template <model_type type>
void ModelTemplate<type>::initializeBEEngine() {
engine = std::make_unique<BEEngineTmpl<type>>(this);
}
namespace {
template <typename T>
std::vector<T> extract(const std::vector<T>& v, UInt n) {
std::vector<T> ex(v.size() - n);
std::copy(v.begin() + n, v.end(), ex.begin());
return ex;
}
} // namespace
template <model_type type>
std::vector<UInt> ModelTemplate<type>::getGlobalDiscretization() const {
auto global_bdisc = Partitioner<trait::boundary_dimension>::global_size(
getBoundaryDiscretization());
if (trait::dimension != trait::boundary_dimension)
global_bdisc.insert(global_bdisc.begin(), getDiscretization().front());
return global_bdisc;
}
template <model_type type>
std::vector<UInt> ModelTemplate<type>::getBoundaryDiscretization() const {
return getDiscretization();
}
template <>
std::vector<UInt>
ModelTemplate<model_type::volume_1d>::getBoundaryDiscretization() const {
return extract(getDiscretization(), 1);
}
template <>
std::vector<UInt>
ModelTemplate<model_type::volume_2d>::getBoundaryDiscretization() const {
return extract(getDiscretization(), 1);
}
template <model_type type>
std::vector<Real> ModelTemplate<type>::getBoundarySystemSize() const {
return getSystemSize();
}
template <>
std::vector<Real>
ModelTemplate<model_type::volume_1d>::getBoundarySystemSize() const {
return extract(getSystemSize(), 1);
}
template <>
std::vector<Real>
ModelTemplate<model_type::volume_2d>::getBoundarySystemSize() const {
return extract(getSystemSize(), 1);
}
/* -------------------------------------------------------------------------- */
template <model_type type>
void ModelTemplate<type>::applyElasticity(GridBase<Real>& stress,
const GridBase<Real>& strain) const {
const influence::ElasticHelper<dim> elasticity(getShearModulus(),
getPoissonRatio());
// Checking number of components
if (strain.getNbComponents() == dim * dim) {
Loop::loop([&elasticity](auto sigma,
auto epsilon) { sigma = elasticity(epsilon); },
range<MatrixProxy<Real, dim, dim>>(stress),
range<MatrixProxy<const Real, dim, dim>>(strain));
}
// in case we have voigt notation
else if (strain.getNbComponents() == voigt_size<dim>::value) {
Loop::loop([&elasticity](auto sigma,
auto epsilon) { sigma = elasticity(epsilon); },
range<SymMatrixProxy<Real, dim>>(stress),
range<SymMatrixProxy<const Real, dim>>(strain));
}
else
TAMAAS_EXCEPTION("Strain components do not match dimension");
}
/* -------------------------------------------------------------------------- */
/* Template instanciation */
/* -------------------------------------------------------------------------- */
template class ModelTemplate<model_type::basic_1d>;
template class ModelTemplate<model_type::basic_2d>;
template class ModelTemplate<model_type::surface_1d>;
template class ModelTemplate<model_type::surface_2d>;
template class ModelTemplate<model_type::volume_1d>;
template class ModelTemplate<model_type::volume_2d>;
} // namespace tamaas
c4science · Help