Homec4science
Diffusion akantu (bugfix/69-segmentation-fault-caused-by-model-applybc-aka-fromtraction-trac-traction) Browse

rAKA/src/mesh_utils/mesh_partition.hh
bf637b29f909bugfix/69-segmentation…

Last Change
Raw File

mesh_partition.hh
View Options

/**
* Copyright (©) 2010-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_csr.hh"
#include "mesh.hh"
/* -------------------------------------------------------------------------- */
#ifndef AKANTU_MESH_PARTITION_HH_
#define AKANTU_MESH_PARTITION_HH_
namespace akantu {
class MeshPartition {
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public:
MeshPartition(Mesh & mesh, Int spatial_dimension,
const ID & id = "MeshPartitioner");
virtual ~MeshPartition();
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
public:
/// define a partition of the mesh
virtual void partitionate(
Int nb_part,
const std::function<Int(const Element &, const Element &)> &
edge_load_func =
[](auto && /*unused*/, auto && /*unused*/) { return 1; },
const std::function<Int(const Element &)> & vertex_load_func =
[](auto && /*unused*/) { return 1; }) = 0;
/// reorder the nodes to reduce the filling during the factorization of a
/// matrix that has a profil based on the connectivity of the mesh
virtual void reorder() = 0;
/// fill the partitions array with a given linearized partition information
void fillPartitionInformation(const Mesh & mesh,
const Int * linearized_partitions);
virtual void printself(std::ostream & stream, int indent = 0) const;
protected:
/// build the dual graph of the mesh, for all element of spatial_dimension
void
buildDualGraph(Array<Int> & dxadj, Array<Int> & dadjncy,
Array<Int> & edge_loads,
const std::function<Int(const Element &, const Element &)> &
edge_load_func,
Array<Int> & vertex_loads,
const std::function<Int(const Element &)> & vertex_load_func);
/// tweak the mesh to handle the PBC pairs
void tweakConnectivity();
/// restore the mesh that has been tweaked
void restoreConnectivity();
/* ------------------------------------------------------------------------ */
/* Accessors */
/* ------------------------------------------------------------------------ */
public:
bool hasPartitions(ElementType type, GhostType ghost_type);
AKANTU_GET_MACRO_AUTO(Partitions, partitions);
AKANTU_GET_MACRO_BY_ELEMENT_TYPE_CONST(Partition, partitions, Idx);
AKANTU_GET_MACRO_AUTO(GhostPartitionCSR, ghost_partitions_csr);
AKANTU_GET_MACRO_AUTO(NbPartition, nb_partitions);
AKANTU_SET_MACRO(NbPartition, nb_partitions, Int);
protected:
Idx linearized(const Element & element);
Element unlinearized(Idx lin_element);
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected:
/// id
ID id;
/// the mesh to partition
Mesh & mesh;
/// dimension of the elements to consider in the mesh
Int spatial_dimension;
/// number of partitions
Int nb_partitions;
/// partition numbers
ElementTypeMapArray<Idx> partitions;
ElementTypeMap<CSR<Idx>> ghost_partitions_csr;
ElementTypeMapArray<Idx> ghost_partitions;
ElementTypeMapArray<Idx> ghost_partitions_offset;
Array<Int> * permutation;
ElementTypeMapArray<Idx> saved_connectivity;
// vector of pair to ensure the iteration order
std::vector<std::pair<ElementType, Idx>> linearized_offsets;
};
/// standard output stream operator
inline std::ostream & operator<<(std::ostream & stream,
const MeshPartition & _this) {
_this.printself(stream);
return stream;
}
} // namespace akantu
#ifdef AKANTU_USE_SCOTCH
#include "mesh_partition_scotch.hh"
#endif
#endif /* AKANTU_MESH_PARTITION_HH_ */
c4science · Help