Homec4science
Diffusion akantu (features/changes-to-address-joss-review) Browse

rAKA/src/mesh_utils/cohesive_element_inserter_helper.hh
cb6de5f8a3bcfeatures/changes-to…

Last Change
Raw File

cohesive_element_inserter_helper.hh
View Options

/**
* @file cohesive_element_inserter_helper.hh
*
* @author Nicolas Richart <nicolas.richart@epfl.ch>
*
* @date creation: Tue Sep 08 2020
* @date last modification: Wed Nov 11 2020
*
* @brief An helper class to handle cohesive element insertion
*
*
* @section LICENSE
*
* Copyright (©) 2018-2021 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* Akantu is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Akantu 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 Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#include "aka_array.hh"
#include "mesh.hh"
/* -------------------------------------------------------------------------- */
#ifndef __AKANTU_COHESIVE_ELEMENT_INSERTER_HELPER_HH__
#define __AKANTU_COHESIVE_ELEMENT_INSERTER_HELPER_HH__
namespace akantu {
/* -------------------------------------------------------------------------- */
class CohesiveElementInserterHelper {
public:
CohesiveElementInserterHelper(
Mesh & mesh, const ElementTypeMapArray<bool> & facet_insertion);
Int insertCohesiveElement();
Int insertFacetsOnly();
private:
template <Int dim> Int insertFacetsOnlyImpl();
template <Int dim> void doubleFacets();
template <Int dim> void findSubfacetToDouble();
void doubleNodes(const std::vector<Idx> & old_nodes);
bool findElementsAroundSubfacet(
const Element & starting_element, const Element & end_facet,
const Vector<Idx> & subfacet_connectivity,
std::vector<Element> & element_list, std::vector<Element> & facet_list,
std::vector<Element> * subfacet_list = nullptr);
static inline bool hasElement(const Vector<Idx> & nodes_element,
const Vector<Idx> & nodes);
static inline bool
removeElementsInVector(const std::vector<Element> & elem_to_remove,
std::vector<Element> & elem_list);
void updateElementalConnectivity(
Mesh & mesh, Idx old_node, Idx new_node,
const std::vector<Element> & element_list,
const std::vector<Element> * facet_list = nullptr);
// update functions
void updateElementToSubelement(Int dim, bool facet_mode);
void updateSubelementToElement(Int dim, bool facet_mode);
void updateQuadraticSegments(Int dim);
void updateCohesiveData();
void doublePointFacet();
template <Int spatial_dimension> void doubleSubfacet();
decltype(auto) elementsOfDimToElementsOfDim(Int dim1, Int dim2) {
AKANTU_DEBUG_ASSERT(dim1 >= 0 and dim1 <= 3,
"dimension of target element out of range");
AKANTU_DEBUG_ASSERT(dim2 >= 0 and dim2 <= 3,
"dimension of source element out of range");
auto & array = dimelements_to_dimelements[dim1][dim2];
if (not array) {
array = std::make_unique<Array<std::vector<Element>>>();
}
return (*array);
}
public:
decltype(auto) getNewElements() const { return (new_elements); }
decltype(auto) getDoubledNodes() const { return (doubled_nodes); }
private:
std::array<std::unique_ptr<Array<Element>>, 3> facets_to_double_by_dim;
std::array<std::array<std::unique_ptr<Array<std::vector<Element>>>, 2>, 4>
dimelements_to_dimelements;
Array<Idx> doubled_nodes;
Array<Element> new_elements;
Mesh & mesh;
Mesh & mesh_facets;
ElementTypeMap<Int> nb_new_facets;
};
} // namespace akantu
#endif /* __AKANTU_COHESIVE_ELEMENT_INSERTER_HELPER_HH__ */
c4science · Help