From 1e304e99364236ebca83fde4648cd41bfbc0c2b5 Mon Sep 17 00:00:00 2001 From: geay Date: Thu, 12 Jun 2014 15:23:30 +0200 Subject: [PATCH] Waiting for Gauthier feedback. --- src/ParaMEDMEM/CMakeLists.txt | 2 - src/ParaMEDMEM/ICoCoTrioField.cxx | 395 ------------------------------ src/ParaMEDMEM/ICoCoTrioField.hxx | 74 ------ src/ParaMEDMEM/TrioDEC.cxx | 44 ---- src/ParaMEDMEM/TrioDEC.hxx | 33 --- 5 files changed, 548 deletions(-) delete mode 100644 src/ParaMEDMEM/ICoCoTrioField.cxx delete mode 100644 src/ParaMEDMEM/ICoCoTrioField.hxx delete mode 100644 src/ParaMEDMEM/TrioDEC.cxx delete mode 100644 src/ParaMEDMEM/TrioDEC.hxx diff --git a/src/ParaMEDMEM/CMakeLists.txt b/src/ParaMEDMEM/CMakeLists.txt index d5bed9bfe..295313064 100644 --- a/src/ParaMEDMEM/CMakeLists.txt +++ b/src/ParaMEDMEM/CMakeLists.txt @@ -58,8 +58,6 @@ SET(paramedmem_SOURCES ParaFIELD.cxx ParaGRID.cxx BlockTopology.cxx - ICoCoTrioField.cxx - TrioDEC.cxx ) ADD_LIBRARY(paramedmem SHARED ${paramedmem_SOURCES}) diff --git a/src/ParaMEDMEM/ICoCoTrioField.cxx b/src/ParaMEDMEM/ICoCoTrioField.cxx deleted file mode 100644 index ea0b81e4a..000000000 --- a/src/ParaMEDMEM/ICoCoTrioField.cxx +++ /dev/null @@ -1,395 +0,0 @@ -// Copyright (C) 2007-2014 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, 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 -// 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 this library; if not, write to the Free Software -// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA -// -// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com -// - -// ICoCo file common to several codes -// ICoCoTrioField.cxx -// version 1.2 10/05/2010 - -#include - -#include "ICoCoMEDField.hxx" -#include "MEDCouplingUMesh.hxx" -#include "MEDCouplingFieldDouble.hxx" -#include "MEDCouplingAutoRefCountObjectPtr.hxx" - -#include -#include -#include - -using namespace ICoCo; -using namespace std; - -// Default constructor -TrioField::TrioField() : - _type(0), - _mesh_dim(0), - _space_dim(0), - _nbnodes(0), - _nodes_per_elem(0), - _nb_elems(0), - _itnumber(0), - _connectivity(0), - _coords(0), - _time1(0.), - _time2(0.), - _nb_field_components(0), - _field(0), - _has_field_ownership(false) { } - -// Copy constructor -TrioField::TrioField(const TrioField& OtherField):_connectivity(0),_coords(0),_field(0) { - (*this)=OtherField; -} - -// Destructor -TrioField::~TrioField() { - clear(); -} - -// After the call to clear(), all pointers are null and field ownership is false. -// Arrays are deleted if necessary -void TrioField::clear() { - if (_connectivity) - delete[] _connectivity; - if (_coords) - delete[] _coords; - if (_field && _has_field_ownership) - delete[] _field; - _connectivity=0; - _coords=0; - _field=0; - _has_field_ownership=false; -} - -// Returns the number of value locations -// The size of field is nb_values()*_nb_field_components -int TrioField::nb_values() const { - if (_type==0) - return _nb_elems; - else if (_type==1) - return _nbnodes; - throw 0; - return -1; -} - -// Save field to a .field file (loadable by visit!) -void TrioField::save(ostream& os) const{ - - os << setprecision(12); - os << getName() << endl; - os << _type << endl; - os << _mesh_dim << endl; - os << _space_dim << endl; - os << _nbnodes << endl; - os << _nodes_per_elem << endl; - os << _nb_elems << endl; - - os<< _itnumber<> name; - setName(name); - in >> _type; - in >> _mesh_dim; - in >> _space_dim; - in >> _nbnodes; - in >> _nodes_per_elem; - in >> _nb_elems; - - in >> _itnumber; - if (_connectivity) - delete [] _connectivity; - _connectivity=new int[_nodes_per_elem*_nb_elems]; - for (int i=0;i<_nb_elems;i++) { - for (int j=0;j<_nodes_per_elem;j++) - in >> _connectivity[i*_nodes_per_elem+j]; - } - if (_coords) - delete [] _coords; - _coords=new double[_nbnodes*_space_dim]; - for (int i=0;i<_nbnodes;i++) { - for (int j=0;j<_space_dim;j++) - in >> _coords[i*_space_dim+j]; - } - - in >> _time1; - in >> _time2; - in >> _nb_field_components; - int test; - in >> test; - if (test) { - if (_field) - delete [] _field; - _field=new double[_nb_field_components*nb_values()]; - for (int i=0;i> _field[i*_nb_field_components+j]; - } - } - else - _field=0; - - in >> _has_field_ownership; -} - - -// After the call to set_standalone(), field ownership is true and field is allocated -// to the size _nb_field_components*nb_values(). -// The values of the field have been copied if necessary. -void TrioField::set_standalone() { - if (!_field) { - _field=new double[_nb_field_components*nb_values()]; - _has_field_ownership=true; - - } - else if (!_has_field_ownership) { - double *tmp_field=new double[_nb_field_components*nb_values()]; - memcpy(tmp_field,_field,_nb_field_components*nb_values()*sizeof(double)); - _field=tmp_field; - _has_field_ownership=true; - } -} - -// Used to simulate a 0D geometry (Cathare/Trio for example). -void TrioField::dummy_geom() { - _type=0; - _mesh_dim=2; - _space_dim=2; - _nbnodes=3; - _nodes_per_elem=3; - _nb_elems=1; - _itnumber=0; - if (_connectivity) - delete[] _connectivity; - _connectivity=new int[3]; - _connectivity[0]=0; - _connectivity[1]=1; - _connectivity[2]=2; - if (_coords) - delete[] _coords; - _coords=new double[6]; - _coords[0]=0; - _coords[1]=0; - _coords[2]=1; - _coords[3]=0; - _coords[4]=0; - _coords[5]=1; - _time1=0; - _time2=1; - _nb_field_components=1; - if (_field && _has_field_ownership) - delete[] _field; - _has_field_ownership=false; - _field=0; -} - -// Overloading operator = for TrioField -// This becomes an exact copy of NewField. -// If NewField._has_field_ownership is false, they point to the same values. -// Otherwise the values are copied. -TrioField& TrioField::operator=(const TrioField& NewField){ - - clear(); - - _type=NewField._type; - _mesh_dim=NewField._mesh_dim; - _space_dim=NewField._space_dim; - _nbnodes=NewField._nbnodes; - _nodes_per_elem=NewField._nodes_per_elem; - _nb_elems=NewField._nb_elems; - _itnumber=NewField._itnumber; - _time1=NewField._time1; - _time2=NewField._time2; - _nb_field_components=NewField._nb_field_components; - - if (!NewField._connectivity) - _connectivity=0; - else { - _connectivity=new int[_nodes_per_elem*_nb_elems]; - memcpy( _connectivity,NewField._connectivity,_nodes_per_elem*_nb_elems*sizeof(int)); - } - - if (!NewField._coords) - _coords=0; - else { - _coords=new double[_nbnodes*_space_dim]; - memcpy( _coords,NewField._coords,_nbnodes*_space_dim*sizeof(double)); - } - - //Copie des valeurs du champ - _has_field_ownership=NewField._has_field_ownership; - if (_has_field_ownership) { - _field=new double[nb_values()*_nb_field_components]; - memcpy(_field,NewField._field,nb_values()*_nb_field_components*sizeof(double)); - } - else - _field=NewField._field; - - return(*this); - -} - -/*! - * This method is non const only due to this->_field that can be modified (to point to the same zone than returned object). - * So \b warning, to access to \a this->_field only when the returned object is alive. - */ -MEDField *TrioField::build_medfield() -{ - ParaMEDMEM::MEDCouplingAutoRefCountObjectPtr mesh(ParaMEDMEM::MEDCouplingUMesh::New("",_mesh_dim)); - ParaMEDMEM::MEDCouplingAutoRefCountObjectPtr coo(ParaMEDMEM::DataArrayDouble::New()); coo->alloc(_nbnodes,_space_dim); - mesh->setCoords(coo); - double *ptr(coo->getPointer()); - std::copy(_coords,_coords+_space_dim*_nbnodes,ptr); - mesh->allocateCells(_nb_elems); - INTERP_KERNEL::NormalizedCellType elemtype; - switch(_mesh_dim) - { - case 1: - { - switch (_nodes_per_elem) - { - case 2: - elemtype=INTERP_KERNEL::NORM_SEG2; - break; - default: - throw INTERP_KERNEL::Exception("incompatible Trio field - wrong nb of nodes per elem"); - } - break; - } - case 2: - { - switch (_nodes_per_elem) - { - case 3: - elemtype=INTERP_KERNEL::NORM_TRI3; - break; - case 4 : - elemtype=INTERP_KERNEL::NORM_QUAD4; - break; - default: - throw INTERP_KERNEL::Exception("incompatible Trio field - wrong nb of nodes per elem"); - } - break; - } - case 3: - { - switch (_nodes_per_elem) - { - case 4: - elemtype=INTERP_KERNEL::NORM_TETRA4; - break; - case 8 : - elemtype=INTERP_KERNEL::NORM_HEXA8; - break; - default: - throw INTERP_KERNEL::Exception("incompatible Trio field - wrong nb of nodes per elem"); - } - break; - default: - throw INTERP_KERNEL::Exception("incompatible Trio field - wrong mesh dimension"); - } - } - //creating a connectivity table that complies to MED (1 indexing) - //and passing it to _mesh - ParaMEDMEM::MEDCouplingAutoRefCountObjectPtr field; - int *conn(new int[_nodes_per_elem]); - for (int i=0; i<_nb_elems;i++) - { - for(int j=0;j<_nodes_per_elem;j++) - { - conn[j]=_connectivity[i*_nodes_per_elem+j]; - } - mesh->insertNextCell(elemtype,_nodes_per_elem,conn); - } - delete [] conn; - mesh->finishInsertingCells(); - // - //field on the sending end - int nb_case=nb_values(); - if (_type==0) - { - field = ParaMEDMEM::MEDCouplingFieldDouble::New(ParaMEDMEM::ON_CELLS,ParaMEDMEM::ONE_TIME); - } - else - { - field = ParaMEDMEM::MEDCouplingFieldDouble::New(ParaMEDMEM::ON_NODES,ParaMEDMEM::ONE_TIME ); - } - field->setMesh(mesh); - field->setNature(ParaMEDMEM::ConservativeVolumic); - ParaMEDMEM::MEDCouplingAutoRefCountObjectPtr fieldArr(ParaMEDMEM::DataArrayDouble::New()); - fieldArr->alloc(field->getNumberOfTuplesExpected(),_nb_field_components); - field->setName(getName().c_str()); - std::string meshName("SupportOf_"); meshName+=getName(); - mesh->setName(meshName.c_str()); - field->setTime(_time1,0,_itnumber); - if (_field!=0) - { - for (int i =0; isetIJ(i,j,_field[i*_nb_field_components+j]); - } - } - //field on the receiving end - else - { - // the trio field points to the pointer inside the MED field - _field=fieldArr->getPointer(); - for (int i=0; i<_nb_field_components*nb_case;i++) - _field[i]=0.0; - } - field->setArray(fieldArr); - return new MEDField(field); -} - - - diff --git a/src/ParaMEDMEM/ICoCoTrioField.hxx b/src/ParaMEDMEM/ICoCoTrioField.hxx deleted file mode 100644 index c0b81b128..000000000 --- a/src/ParaMEDMEM/ICoCoTrioField.hxx +++ /dev/null @@ -1,74 +0,0 @@ -// Copyright (C) 2007-2014 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, 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 -// 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 this library; if not, write to the Free Software -// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA -// -// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com -// - -// ICoCo file common to several codes -// ICoCoTrioField.h -// version 1.2 10/05/2010 - -#ifndef _ICoCoTrioField_included_ -#define _ICoCoTrioField_included_ - -#include -namespace ICoCo { - class MEDField; - ////////////////////////////////////////////////////////////////////////////// - // - // .DESCRIPTION - // class TrioField, used for coupling Trio codes via the ICoCo interface - // This structure contains all the necessary information - // for constructing a ParaMEDMEM::ParaFIELD (with the addition of the MPI - // communicator). - // This structure can either own or not _field values (_has_field_ownership) - // For _coords, _connectivity and _field, a null pointer means no data allocated. - // _coords and _connectivity tables, when allocated, are always owned by the TrioField. - // - ////////////////////////////////////////////////////////////////////////////// - class TrioField:public Field { - public: - - TrioField(); - TrioField(const TrioField& OtherField); - ~TrioField(); - void clear(); - void set_standalone(); - void dummy_geom(); - TrioField& operator=(const TrioField& NewField); - void save(std::ostream& os) const; - void restore(std::istream& in); - int nb_values() const ; - MEDField *build_medfield(); - public: - int _type ; // 0 elem 1 nodes - int _mesh_dim; - int _space_dim; - int _nbnodes; - int _nodes_per_elem; - int _nb_elems; - int _itnumber; - int* _connectivity; - double* _coords; - - double _time1,_time2; - int _nb_field_components; - double* _field; - bool _has_field_ownership; - }; -} - -#endif diff --git a/src/ParaMEDMEM/TrioDEC.cxx b/src/ParaMEDMEM/TrioDEC.cxx deleted file mode 100644 index 6c9451c87..000000000 --- a/src/ParaMEDMEM/TrioDEC.cxx +++ /dev/null @@ -1,44 +0,0 @@ -// Data exchange channel for trio -// TrioDEC.cxx -// version 0.0 06/06/2014 - -#include "TrioDEC.hxx" - -#include "ICoCoMEDField.hxx" -#include "ICoCoTrioField.hxx" - -using namespace ParaMEDMEM; -using namespace ICoCo; - -TrioDEC::TrioDEC():_my_traduced_field(0) -{ -} - -TrioDEC::TrioDEC(ProcessorGroup& source_group, ProcessorGroup& target_group):InterpKernelDEC(source_group,target_group),_my_traduced_field(0) -{ -} - -TrioDEC::TrioDEC(const std::set& src_ids, const std::set& trg_ids, const MPI_Comm& world_comm):InterpKernelDEC(src_ids,trg_ids,world_comm),_my_traduced_field(0) -{ -} - -void TrioDEC::attachLocalField(ICoCo::TrioField *field) -{ - if(!field) - throw INTERP_KERNEL::Exception("TrioDEC::attachLocalField : The input trio Field is NULL !"); - releaseInternalPointer(); - _my_traduced_field=field->build_medfield(); - DisjointDEC::attachLocalField(_my_traduced_field); -} - -void TrioDEC::releaseInternalPointer() -{ - if(_my_traduced_field) - delete _my_traduced_field; - _my_traduced_field=0; -} - -TrioDEC::~TrioDEC() -{ - releaseInternalPointer(); -} diff --git a/src/ParaMEDMEM/TrioDEC.hxx b/src/ParaMEDMEM/TrioDEC.hxx deleted file mode 100644 index 711bfe661..000000000 --- a/src/ParaMEDMEM/TrioDEC.hxx +++ /dev/null @@ -1,33 +0,0 @@ -// Data exchange channel for trio -// TrioDEC.h -// version 0.0 06/06/2014 - -#ifndef _TrioDEC_included_ -#define _TrioDEC_included_ - -#include "InterpKernelDEC.hxx" - -namespace ICoCo -{ - class MEDField; - class TrioField; -} - -namespace ParaMEDMEM -{ - class TrioDEC : public InterpKernelDEC - { - public: - TrioDEC(); - TrioDEC(ProcessorGroup& source_group, ProcessorGroup& target_group); - TrioDEC(const std::set& src_ids, const std::set& trg_ids, const MPI_Comm& world_comm=MPI_COMM_WORLD); - void attachLocalField(ICoCo::TrioField *field); - virtual ~TrioDEC(); - private: - void releaseInternalPointer(); - private : - ICoCo::MEDField *_my_traduced_field; - }; -} - -#endif -- 2.39.2