Homec4science
Diffusion akantu (features/38-damage-energy-computation-missing-in-phasefield-model) Browse

rAKA/src/synchronizer/facet_synchronizer_inline_impl.cc
ba6837e25b9bfeatures/38-damage…

Last Change
Raw File

facet_synchronizer_inline_impl.cc
View Options

/**
* @file facet_synchronizer_inline_impl.cc
*
* @author Nicolas Richart <nicolas.richart@epfl.ch>
* @author Marco Vocialta <marco.vocialta@epfl.ch>
*
* @date creation: Wed Nov 05 2014
* @date last modification: Tue Nov 07 2017
*
* @brief facet synchronizer inline implementation
*
* @section LICENSE
*
* Copyright (©) 2015-2018 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/>.
*
*/
/* -------------------------------------------------------------------------- */
// template<GhostType ghost_facets>
// inline void FacetSynchronizer::getFacetGlobalConnectivity(const
// DistributedSynchronizer & distributed_synchronizer,
// const
// ElementTypeMapArray<UInt>
// & rank_to_facet,
// const Array<Element> *
// elements,
// Array<ElementTypeMapArray<UInt>
// *> & connectivity,
// Array<ElementTypeMapArray<UInt>
// *> & facets) {
// AKANTU_DEBUG_IN();
// UInt spatial_dimension = mesh.getSpatialDimension();
// /// init facet check tracking
// ElementTypeMapArray<bool> facet_checked("facet_checked", id);
// mesh.initElementTypeMapArray(facet_checked, 1, spatial_dimension - 1);
// Mesh::type_iterator first = mesh.firstType(spatial_dimension - 1,
// ghost_facets);
// Mesh::type_iterator last = mesh.lastType(spatial_dimension - 1,
// ghost_facets);
// for (; first != last; ++first) {
// ElementType type = *first;
// Array<bool> & f_checked = facet_checked(type, ghost_facets);
// UInt nb_facet = mesh.getNbElement(type, ghost_facets);
// f_checked.resize(nb_facet);
// }
// const Array<UInt> & nodes_global_ids =
// distributed_synchronizer.mesh.getGlobalNodesIds();
// /// loop on every processor
// for (UInt p = 0; p < nb_proc; ++p) {
// if (p == rank) continue;
// /// reset facet check
// Mesh::type_iterator first = mesh.firstType(spatial_dimension - 1,
// ghost_facets);
// Mesh::type_iterator last = mesh.lastType(spatial_dimension - 1,
// ghost_facets);
// for (; first != last; ++first) {
// ElementType type = *first;
// Array<bool> & f_checked = facet_checked(type, ghost_facets);
// f_checked.clear();
// }
// ElementTypeMapArray<UInt> & global_conn = (*connectivity(p));
// const Array<Element> & elem = elements[p];
// ElementTypeMapArray<UInt> & facet_list = (*facets(p));
// UInt nb_element = elem.getSize();
// /// loop on every send/recv element
// for (UInt el = 0; el < nb_element; ++el) {
// ElementType type = elem(el).type;
// GhostType gt = elem(el).ghost_type;
// UInt el_index = elem(el).element;
// const Array<Element> & facet_to_element =
// mesh.getSubelementToElement(type, gt);
// UInt nb_facets_per_element = Mesh::getNbFacetsPerElement(type);
// ElementType facet_type = Mesh::getFacetType(type);
// UInt nb_nodes_per_facet = Mesh::getNbNodesPerElement(facet_type);
// Vector<UInt> conn_tmp(nb_nodes_per_facet);
// /// loop on every facet of the element
// for (UInt f = 0; f < nb_facets_per_element; ++f) {
// const Element & facet = facet_to_element(el_index, f);
// if (facet == ElementNull) continue;
// UInt facet_index = facet.element;
// GhostType facet_gt = facet.ghost_type;
// const Array<UInt> & t_to_f = rank_to_facet(facet_type, facet_gt);
// /// exclude not ghost facets, facets assigned to other
// /// processors
// if (facet_gt != ghost_facets) continue;
// if ((facet_gt == _ghost) && (t_to_f(facet_index) != p)) continue;
// /// exclude facets that have already been added
// Array<bool> & f_checked = facet_checked(facet_type, facet_gt);
// if (f_checked(facet_index)) continue;
// else f_checked(facet_index) = true;
// /// add facet index
// Array<UInt> & f_list = facet_list(facet_type, facet_gt);
// f_list.push_back(facet_index);
// /// add sorted facet global connectivity
// const Array<UInt> & conn = mesh.getConnectivity(facet_type, facet_gt);
// Array<UInt> & g_conn = global_conn(facet_type, facet_gt);
// for (UInt n = 0; n < nb_nodes_per_facet; ++n)
// conn_tmp(n) = nodes_global_ids(conn(facet_index, n));
// std::sort(conn_tmp.storage(), conn_tmp.storage() + nb_nodes_per_facet);
// g_conn.push_back(conn_tmp);
// }
// }
// }
// AKANTU_DEBUG_OUT();
// }
c4science · Help