-// Copyright (C) 2007-2013 CEA/DEN, EDF R&D
+// Copyright (C) 2007-2016 CEA/DEN, EDF R&D
//
// This library 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 2.1 of the License.
+// version 2.1 of the License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
#include <fvm_locator.h>
}
-namespace ParaMEDMEM
+namespace MEDCoupling
{
/*!
- \defgroup noncoincidentdec NonCoincidentDEC
+ \anchor NonCoincidentDEC-det
+ \class NonCoincidentDEC
- \section overview Overview
-
- \c NonCoincidentDEC enables nonconservative remapping of fields
+ \c NonCoincidentDEC enables non-conservative remapping of fields
between two parallel codes.
The computation is possible for 3D meshes and 2D meshes.
- It is not available for 3D surfaces. The computation enables fast parallel localization, and is based on a point in element search, followed
+ It is not available for 3D surfaces.
+
+ The computation enables fast parallel localization, and is based on a point in element search, followed
by a field evaluation at the point location. Thus, it is typically
- faster than the \ref interpkerneldec which gives a \ref conservativeremapping.
- It is particularly true for the initialisation phase (synchronize)
- which is very computationnaly intensive in \ref interpkerneldec.
+ faster than the \ref InterpKernelDEC-det "InterpKernelDEC" which uses a
+ \ref InterpKerRemapGlobal "conservative remapping" (i.e. the same algorithms of volume
+ intersection as in the \ref remapper "sequential remapper")
+ It is particularly true for the initialisation phase (synchronize() method)
+ which has a significant computation cost in \ref InterpKernelDEC-det.
- In the present version, only fields lying on elements are considered.
+ In the present version, only fields lying on elements ("P0") are considered.
The value is estimated by locating the barycenter of the target
side cell in a source cell and sending the value of this source cell
as the value of the target cell.
fvm_nodal_t* medmemMeshToFVMMesh(const MEDMEM::MESH* mesh)
{
// create an FVM structure from the paramesh structure
- fvm_nodal_t * fvm_nodal = fvm_nodal_create(mesh->getName().c_str(),mesh->getMeshDimension());
+ std::string meshName(mesh->getName());//this line avoid that mesh->getName() object killed before fvm_nodal_create read the const char *.
+ fvm_nodal_t * fvm_nodal = fvm_nodal_create(meshName.c_str(),mesh->getMeshDimension());
//loop on cell types
int nbtypes = mesh->getNumberOfTypes(MED_EN::MED_CELL);
{
// create an FVM structure from the paramesh structure
- fvm_nodal_t * fvm_nodal = fvm_nodal_create(support->getName().c_str(),1);
+ std::string supportName(support->getName());//this line avoid that support->getName() object killed before fvm_nodal_create read the const char *.
+ fvm_nodal_t * fvm_nodal = fvm_nodal_create(supportName.c_str(),1);
const MEDMEM::MESH* mesh= support->getMesh();
{
}
- /*!
- \addtogroup noncoincidentdec
- @{
- */
-
- /*! Constructor of a non coincident \ref dec with
+ /*! Constructor of a non coincident \ref para-dec "DEC" with
* a source group on which lies a field lying on a mesh and a
* target group on which lies a mesh.
*
if (_source_group->containsMyRank())
{
MEDMEM::MESH* mesh = _local_field->getField()->getSupport()->getMesh();
- fvm_nodal_t* source_nodal = ParaMEDMEM::medmemMeshToFVMMesh(mesh);
+ fvm_nodal_t* source_nodal = MEDCoupling::medmemMeshToFVMMesh(mesh);
int target_size = _target_group->size() ;
int start_rank= _source_group->size();
{
MEDMEM::MESH* mesh = _local_field->getField()->getSupport()->getMesh();
- fvm_nodal_t* target_nodal = ParaMEDMEM::medmemMeshToFVMMesh(mesh);
+ fvm_nodal_t* target_nodal = MEDCoupling::medmemMeshToFVMMesh(mesh);
int source_size = _source_group->size();
int start_rank= 0 ;
const MPI_Comm* comm = (dynamic_cast<const MPIProcessorGroup*> (_union_group))->getComm();
renormalizeTargetField();
}
- /*!
- @}
- */
}