Homec4science
Diffusion akantu (features/72-czm-stability) Browse

rAKA/src/mesh_utils/cohesive_element_inserter_helper.hh
a67b29746942features/72-czm…

Last Change
Raw File

cohesive_element_inserter_helper.hh
View Options

/**
* Copyright (©) 2020-2023 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* This file is part of Akantu
*
* 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