+++ /dev/null
-// Copyright (C) 2007-2013 CEA/DEN, EDF R&D, OPEN CASCADE
-//
-// Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
-// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
-//
-// 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.
-//
-// 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
-//
-
-// File : CalculatorEngine.cxx
-// Author : Laurent DADA, CEA
-// Module : CalculatorComponent
-// $Header$
-//
-#include <unistd.h>
-
-#include <iostream>
-#include <sstream>
-
-#include "COMPONENT_version.h"
-#include "CalculatorEngine.hxx"
-#include "MEDMEM_Support_i.hxx"
-#include "SUPPORTClient.hxx"
-
-#include "MEDMEM_define.hxx"
-#include "MEDMEM_STRING.hxx"
-#include "MEDMEM_Exception.hxx"
-#include "MEDMEM_Unit.hxx"
-
-#include "utilities.h"
-
-using namespace MEDMEM;
-using namespace std;
-
-#ifndef MED_TAILLE_PNOM
-#define MED_TAILLE_PNOM MED_SNAME_SIZE
-#endif
-
-typedef FIELD<double,MEDMEM::FullInterlace> TFieldDouble;
-typedef FIELDTEMPLATE_I<double,MEDMEM::FullInterlace> TFieldDouble_i;
-
-//================================================================================
-// static PrintFieldValues - shows field contents
-//================================================================================
-
-static void PrintFieldValues (TFieldDouble * field, int until_index)
-{
- const double * values = field -> getValue();
- int nb_comp = field -> getNumberOfComponents();
- MESSAGE( "Field : " << field -> getName() );
- MESSAGE( " Description : " << field -> getDescription() );
- MESSAGE( " Number of components : " << nb_comp );
-
- for (int i=0 ; i<nb_comp ; i++) {
- MESSAGE ( " Component " << field -> getComponentName(i+1) << endl \
- << " Unit : " << field -> getMEDComponentUnit(i+1) << endl \
- << " Description : " << field -> getComponentDescription(i+1) );
- }
- MESSAGE( " First " << until_index << " values : " );
-
- for (int j=0 ; j < until_index ; j++) {
- for (int i=0;i<nb_comp;i++) MESSAGE( values[i+j*nb_comp] << " " );
- }
-}
-
-
-//================================================================================
-// Constructors
-//================================================================================
-
-CalculatorEngine::CalculatorEngine(CORBA::ORB_ptr orb,
- PortableServer::POA_ptr poa,
- PortableServer::ObjectId * contId,
- const char *instanceName,
- const char *interfaceName) :
- Engines_Component_i(orb, poa, contId, instanceName, interfaceName,true)
-{
- MESSAGE("CalculatorEngine::CalculatorEngine activate object instanceName("
- << instanceName << ") interfaceName(" << interfaceName << ")" )
- _thisObj = this ;
- _id = _poa->activate_object(_thisObj);
- SALOME_NamingService *NS = new SALOME_NamingService(orb);
- _NS=NS;
-}
-
-CalculatorEngine::CalculatorEngine()
-{
-}
-
-//================================================================================
-// Destructors
-//================================================================================
-
-CalculatorEngine::~CalculatorEngine()
-{
- delete _NS;
-}
-
-char* CalculatorEngine::getVersion()
-{
-#if COMPONENT_DEVELOPMENT
- return CORBA::string_dup(COMPONENT_VERSION_STR"dev");
-#else
- return CORBA::string_dup(COMPONENT_VERSION_STR);
-#endif
-}
-
-static omni_mutex aPutToStudyMutex;
-SALOME_MED::FIELDDOUBLE_ptr CalculatorEngine::PutToStudy(SALOME_MED::FIELDDOUBLE_ptr theField1,
- CORBA::Long theStudyId) {
- omni_mutex_lock aLock(aPutToStudyMutex);
-
- CORBA::Object_var anObj = _NS->Resolve("/myStudyManager");
- SALOMEDS::StudyManager_var aManager = SALOMEDS::StudyManager::_narrow(anObj);
-// SALOMEDS::ListOfOpenStudies_var aList = aManager->GetOpenStudies();
-// SALOMEDS::Study_var aStudy = aManager->GetStudyByName(aList[0]);
- SALOMEDS::Study_var aStudy = aManager->GetStudyByID(theStudyId);
- SALOMEDS::StudyBuilder_var aBuilder = aStudy->NewBuilder();
-
- SALOMEDS::SComponent_var aComponent = aStudy->FindComponent("Calculator");
- if (aComponent->_is_nil()) aComponent = aBuilder->NewComponent("Calculator");
- SALOMEDS::GenericAttribute_var anAttr = aBuilder->FindOrCreateAttribute(aComponent, "AttributeName");
- SALOMEDS::AttributeName_var aName = SALOMEDS::AttributeName::_narrow(anAttr);
- aName->SetValue("Calculator");
-
- SALOMEDS::SObject_var aSO = aBuilder->NewObject(aComponent);
- anAttr = aBuilder->FindOrCreateAttribute(aSO, "AttributeName");
- aName = SALOMEDS::AttributeName::_narrow(anAttr);
- char aFieldName[10];
- sprintf(aFieldName, "Field_%d", (int)aSO->Tag());
- aName->SetValue(strdup(aFieldName));
-
- ORB_INIT &init = *SINGLETON_<ORB_INIT>::Instance() ;
- ASSERT(SINGLETON_<ORB_INIT>::IsAlreadyExisting()) ;
- CORBA::ORB_var &orb = init(0,0);
-
- string iorStr = orb->object_to_string(theField1);
- anAttr = aBuilder->FindOrCreateAttribute(aSO, "AttributeIOR");
- SALOMEDS::AttributeIOR_var anIOR = SALOMEDS::AttributeIOR::_narrow(anAttr);
- anIOR->SetValue(iorStr.c_str());
-
- return SALOME_MED::FIELDDOUBLE::_duplicate(theField1);
-}
-
-
-//=======================================================================================
-// Addition of 2 FieldNodeDouble (based on Med library)
-//
-// We assume that :
-// - support
-// - number of components
-// - name, type, unit and description of components
-// are the same for the 2 fields.
-//
-// Otherwise,
-// the result will be a copy of the first field.
-//
-//=======================================================================================
-
-static omni_mutex anAddMutex;
-SALOME_MED::FIELDDOUBLE_ptr CalculatorEngine::Add(SALOME_MED::FIELDDOUBLE_ptr FirstField,
- SALOME_MED::FIELDDOUBLE_ptr SecondField)
-{
- omni_mutex_lock aLock(anAddMutex);
-
- ostringstream tmp;
- bool same_support = true;
-
- beginService("CalculatorEngine::Add");
-
- sendMessage(NOTIF_TRACE, "CalculatorEngine::Add is Computing");
-
- const char* LOC = "CalculatorEngine::Add(SALOME_MED::FIELDDOUBLE_ptr FirstField,SALOME_MED::FIELDDOUBLE_ptr SecondField)";
- BEGIN_OF(LOC);
-
- SCRUTE(FirstField);
- SCRUTE(SecondField);
-
- CORBA::String_var field1_name = FirstField -> getName();
- CORBA::String_var field2_name = SecondField -> getName();
- CORBA::String_var field_description = FirstField -> getDescription();
-
- string message = string("\n CalculatorEngine::Add : ") +
- string("\n first field name --> ") + string(field1_name.in()) +
- string("\n second field name --> ") + string(field2_name.in());
-
- sendMessage(NOTIF_TRACE, message.c_str());
-
- // Number of components
- int nb_comp1 = FirstField -> getNumberOfComponents();
- int nb_comp2 = SecondField -> getNumberOfComponents();
- if (nb_comp1 != nb_comp2) same_support = false;
- MESSAGE( " nb_comp1 --> " << nb_comp1 <<endl );
- MESSAGE( " nb_comp2 --> " << nb_comp2 <<endl );
-
-
- /* Name Description and Unit of each component of the field */
-
- string * component_name = new string[nb_comp1];
- string * component_unit = new string[nb_comp1];
-
- for (int i=0 ; i<nb_comp1 ; i++) {
- CORBA::String_var name = FirstField -> getComponentName(i+1);
- CORBA::String_var unit = FirstField -> getComponentUnit(i+1);
- component_name[i].assign(name.in());
- component_unit[i].assign(unit.in());
-
- if (same_support) {
- CORBA::String_var name2 = SecondField -> getComponentName(i+1);
- CORBA::String_var unit2 = SecondField -> getComponentUnit(i+1);
- if ( (component_name[i] != name2.in()) || (component_unit[i] != unit2.in()))
- same_support = false;
- }
- }
-
- MESSAGE("CalculatorEngine::Add Number of entities in the Support of each fields");
-
- // Number of entities in the Support of the field
-
- int len_value1 = 0;
- int len_value2 = 0;
- SALOME_MED::SUPPORT_var support1 = FirstField -> getSupport();
- SALOME_MED::SUPPORT_var support2 = SecondField -> getSupport();
-
- SCRUTE(support1);
-
- SCRUTE(support2);
-
- SALOME_MED::GMESH_var mesh1 = support1 -> getMesh();
- SALOME_MED::GMESH_var mesh2 = support2 -> getMesh();
-
- SCRUTE(mesh1);
- SCRUTE(mesh2);
-
- if (same_support && support1 != support2) same_support = false;
-
- if (support1 -> isOnAllElements())
- len_value1 = mesh1-> getNumberOfElements(support1 -> getEntity(),SALOME_MED::MED_ALL_ELEMENTS);
- else
- len_value1 = support1 -> getNumberOfElements(SALOME_MED::MED_ALL_ELEMENTS);
-
- SCRUTE(len_value1);
-
- if (support2 -> isOnAllElements())
- len_value2 = mesh2-> getNumberOfElements(support2 -> getEntity(),SALOME_MED::MED_ALL_ELEMENTS);
- else
- len_value2 = support2 -> getNumberOfElements(SALOME_MED::MED_ALL_ELEMENTS);
-
- if (same_support && len_value1 != len_value2) same_support = false;
-
- SCRUTE(len_value2);
-
- // Values of new field
- SALOME_TYPES::ListOfDouble * first_value = FirstField -> getValue(SALOME_MED::MED_FULL_INTERLACE);
- SALOME_TYPES::ListOfDouble * second_value;
- if (same_support) {
- second_value = SecondField -> getValue(SALOME_MED::MED_FULL_INTERLACE);
- tmp << "\n CalculatorEngine::Add : " ;
- tmp << "\n Number of entities in the Support of first field = " << len_value1 ;
- tmp << "\n Number of entities in the Support of second field = " << len_value2 ;
- message = tmp.str();
-
- sendMessage(NOTIF_TRACE, message.c_str());
-
- }
-
- int totalLength = nb_comp1*len_value1;
- double * new_value = new double[totalLength];
- if (same_support)
- {
- sendMessage(NOTIF_TRACE, "CalculatorEngine::Add - field1 and field2 have the same support");
-
- for(int i=0 ; i<totalLength ; i++)
- new_value[i] = ((*first_value)[i]) + ((*second_value)[i]);
- }
- else {
- for(int i=0 ; i<totalLength ; i++)
- new_value[i] = ((*first_value)[i]) ;
- }
-
- // Creation of a Local Field
-
- MESSAGE("CalculatorEngine::Add Creation of the local field "<< totalLength);
-
- TFieldDouble * fieldloc = new TFieldDouble();
- fieldloc -> allocValue(nb_comp1,len_value1);
- fieldloc -> setValue(new_value);
- fieldloc -> setName("-new_Add-");
- fieldloc -> setDescription( field_description.in() );
- fieldloc -> setComponentsNames(component_name);
- fieldloc -> setMEDComponentsUnits(component_unit);
-
- delete [] new_value;
- delete first_value;
- if (same_support)
- delete second_value;
- delete [] component_name;
- delete [] component_unit;
- mesh1->UnRegister();
- mesh2->UnRegister();
- support2->UnRegister();
-
- // Control
- int until_index = ( 20 > len_value1)? len_value1 : 20;
- PrintFieldValues (fieldloc,until_index);
-
- sendMessage(NOTIF_TRACE, "CalculatorEngine::Add - new field created");
-
- // Export result
- // Creation of the new CORBA field
-
- MESSAGE("CalculatorEngine::Add Creation of the CORBA field");
-
- SUPPORT *support1Clt=new SUPPORTClient(support1);
- fieldloc->setSupport(support1Clt);
- support1Clt->removeReference();
-
- TFieldDouble_i * NewField = new TFieldDouble_i(fieldloc,true) ;
- SALOME_MED::FIELDDOUBLE_ptr myFieldIOR = NewField->_this() ;
-
- //sleep(4);
-
- END_OF(LOC);
-
- endService("CalculatorEngine::Add");
- return myFieldIOR;
-}
-
-
-//===============================================================================
-// Multiplication of a FieldNodeDouble by a double value (based on Med library)
-//===============================================================================
-
-static omni_mutex aMulMutex;
-SALOME_MED::FIELDDOUBLE_ptr CalculatorEngine::Mul(SALOME_MED::FIELDDOUBLE_ptr OldField,
- CORBA::Double x1)
-{
- omni_mutex_lock aLock(aMulMutex);
-
- beginService("CalculatorEngine::Mul");
- const char* LOC = "CalculatorEngine::Mul(SALOME_MED::FIELDDOUBLE_ptr OldField,CORBA::Double x1)";
- BEGIN_OF(LOC);
-
- SCRUTE(OldField);
- SCRUTE(x1);
-
- // Name and description of field
- CORBA::String_var field_name = OldField -> getName();
- CORBA::String_var field_description = OldField -> getDescription();
-
- // Number of components
- int nb_comp = OldField -> getNumberOfComponents();
- string * component_name = new string[nb_comp];
- string * component_unit = new string[nb_comp];
-
- for (int i=0 ; i<nb_comp ; i++) {
- CORBA::String_var name = OldField -> getComponentName(i+1);
- component_name[i] = name.in();
- CORBA::String_var unit = OldField -> getComponentUnit(i+1);
- component_unit[i] = unit.in();
- }
-
- MESSAGE("CalculatorEngine::Mul Number of entities in the Support of the field");
-
- // Number of entities in the Support of the field
- int len_value = 0;
-
- SALOME_MED::SUPPORT_var support = OldField -> getSupport();
-
- SCRUTE(support);
-
- SALOME_MED::GMESH_var mesh = support -> getMesh();
-
- SCRUTE(mesh);
-
- if (support -> isOnAllElements())
- len_value = mesh-> getNumberOfElements(support -> getEntity(),SALOME_MED::MED_ALL_ELEMENTS);
- else
- len_value = support -> getNumberOfElements(SALOME_MED::MED_ALL_ELEMENTS);
-
- SCRUTE(len_value);
-
- // Values of new field
- SALOME_TYPES::ListOfDouble * old_value = OldField -> getValue(SALOME_MED::MED_FULL_INTERLACE);
-
- int totalLength = nb_comp*len_value;
- double * new_value = new double[totalLength];
- for(int i=0 ; i<totalLength ; i++) {
- new_value[i] = x1 * ((*old_value)[i]);
- }
-
- // Creation of a Local Field
-
- MESSAGE("CalculatorEngine::Mul Creation of the local field "<< totalLength);
-
- TFieldDouble * fieldloc = new TFieldDouble();
- fieldloc -> allocValue(nb_comp,len_value);
- fieldloc -> setValue(new_value);
- fieldloc -> setName("-new_Mul-");
- fieldloc -> setDescription(field_description.in());
- fieldloc -> setComponentsNames(component_name);
- fieldloc -> setMEDComponentsUnits(component_unit);
-
- delete old_value;
- delete [] new_value;
- delete [] component_name;
- delete [] component_unit;
- mesh->UnRegister();
-
- // Creation of the new CORBA field
-
- MESSAGE("CalculatorEngine::Mul Creation of the CORBA field");
-
- SUPPORT *supportClt=new SUPPORTClient(support);
- fieldloc->setSupport(supportClt);
- supportClt->removeReference();
- TFieldDouble_i * NewField = new TFieldDouble_i(fieldloc,true) ;
- SALOME_MED::FIELDDOUBLE_ptr myFieldIOR = NewField->_this() ;
-
- //sleep(4);
-
- END_OF(LOC);
-
- endService("CalculatorEngine::Mul");
- return myFieldIOR;
-}
-
-
-//================================================================================
-// Build a constant field of double based on first field support (Med library)
-//================================================================================
-
-static omni_mutex aConstantMutex;
-SALOME_MED::FIELDDOUBLE_ptr CalculatorEngine::Constant(SALOME_MED::FIELDDOUBLE_ptr FirstField,
- CORBA::Double x1)
-{
- omni_mutex_lock aLock(aConstantMutex);
-
- beginService("CalculatorEngine::Const");
- const char* LOC = "CalculatorEngine::Constant(SALOME_MED::FIELDDOUBLE_ptr FirstField,CORBA::Double x1)";
- BEGIN_OF(LOC);
-
- // Name and description of field
- CORBA::String_var field_name = FirstField -> getName();
- CORBA::String_var field_description = FirstField -> getDescription();
-
- SALOME_MED::SUPPORT_var FirstSupport = FirstField -> getSupport();
-
- // Number of components
- int nb_comp = FirstField -> getNumberOfComponents();
-
- // Type, name, unit and description of components
- string * component_name = new string[nb_comp];
- string * component_unit = new string[nb_comp];
-
- for (int i = 0 ; i<nb_comp ; i++) {
- CORBA::String_var name = FirstField -> getComponentName(i+1);
- component_name[i] = name.in();
- CORBA::String_var unit = FirstField -> getComponentUnit(i+1);
- component_unit[i] = unit.in();
- }
-
- MESSAGE("CalculatorEngine::Constant Number of entities in the Support of the field");
-
- // Number of entities in the Support of the field
- int len_value = 0;
-
- SCRUTE(FirstSupport);
-
- SALOME_MED::GMESH_var mesh = FirstSupport -> getMesh();
-
- SCRUTE(mesh);
-
- if ( FirstSupport -> isOnAllElements() )
- len_value = mesh -> getNumberOfElements(FirstSupport -> getEntity(),SALOME_MED::MED_ALL_ELEMENTS);
- else
- len_value = FirstSupport -> getNumberOfElements(SALOME_MED::MED_ALL_ELEMENTS);
-
- // Values of new field
- int totalLength = nb_comp*len_value;
- double * new_value = new double[totalLength];
-
- for(int i=0 ; i<totalLength ; i++)
- new_value[i] = x1 ;
-
-
- // Creation of a Local Field
-
- MESSAGE("CalculatorEngine::Constant Creation of the local field "<< totalLength);
-
- TFieldDouble * fieldloc = new TFieldDouble();
- fieldloc -> allocValue(nb_comp,len_value);
- fieldloc -> setValue(new_value);
- fieldloc -> setName("-new_Const_Field-");
- fieldloc -> setDescription(field_description.in());
- fieldloc -> setComponentsNames(component_name);
- fieldloc -> setMEDComponentsUnits(component_unit);
-
- delete [] new_value;
- delete [] component_name;
- delete [] component_unit;
- mesh->UnRegister();
-
- // Control
- int until_index = ( 20 > len_value)? len_value : 20 ;
- PrintFieldValues (fieldloc,until_index);
-
- // Creation of the new CORBA field
-
- MESSAGE("CalculatorEngine::Constant Creation of the CORBA field");
-
- SUPPORT *supportClt=new SUPPORTClient(FirstSupport);
- fieldloc->setSupport(supportClt);
- supportClt->removeReference();
- TFieldDouble_i * NewField = new TFieldDouble_i(fieldloc,true) ;
-
- SALOME_MED::FIELDDOUBLE_ptr myFieldIOR = NewField->_this() ;
-
- //sleep(4);
-
- endService("CalculatorEngine::Const");
-
- END_OF(LOC);
-
- return myFieldIOR;
-}
-
-//================================================================================
-// write a field in a MED file
-//================================================================================
-
-static omni_mutex aWriteMEDfileMutex;
-void CalculatorEngine::writeMEDfile(SALOME_MED::FIELDDOUBLE_ptr field, const char *filename)
-{
- omni_mutex_lock aLock(aWriteMEDfileMutex);
-
- const char * LOC = "CalculatorEngine::writeMEDfile ";
-
- beginService("CalculatorEngine::writeMEDfile");
- BEGIN_OF(LOC);
- MESSAGE("fichier d'impression du champ resultat:"<<filename);
-
- // get support of the field
-
- SALOME_MED::SUPPORT_ptr support = field -> getSupport();
-
- // get mesh from this support
-
- SALOME_MED::GMESH_ptr mesh = support -> getMesh();
-
- // write the mesh
-
- int index_m = mesh->addDriver(SALOME_MED::MED_DRIVER,filename,mesh->getName());
-
- mesh -> write(index_m,"");
-
- //write the field
-
- MESSAGE("fichier :"<<filename);
-
- med_2_3::med_idt _medIdt = med_2_3::MEDfileOpen(const_cast <char *> (filename) , med_2_3::MED_ACC_RDWR);
- SCRUTE(_medIdt);
-
- if (_medIdt<0) return;
-
- int err ;
-
- int component_count= field->getNumberOfComponents();
- string component_name(component_count*MED_TAILLE_PNOM,' ') ;
- string component_unit(component_count*MED_TAILLE_PNOM,' ') ;
-
- for (int i=0; i < component_count ; i++) {
- component_name.replace(i*MED_TAILLE_PNOM,MED_TAILLE_PNOM,(field -> getComponentName(i+1)),0,MED_TAILLE_PNOM);
- component_unit.replace(i*MED_TAILLE_PNOM,MED_TAILLE_PNOM,(field -> getComponentUnit(i+1)),0,MED_TAILLE_PNOM);
- }
-
- MESSAGE("component_name=|"<<component_name<<"|");
- MESSAGE("component_unit=|"<<component_unit<<"|");
-
- // already existing ?
-
- char * champName = new char[MED_NAME_SIZE+1] ;
- char * meshName = new char[MED_NAME_SIZE+1] ;
- char * compDtunit = new char[MED_SNAME_SIZE+1];
- char * compName, * compUnit;
- bool Find = false ;
-
- int n = med_2_3::MEDnField(_medIdt);
-
- med_2_3::med_int nbComp ;
- med_2_3::med_field_type type ;
- med_2_3::med_bool _local;
-
- for (int i = 1; i <= n; i++)
- {
- nbComp = med_2_3::MEDnField(_medIdt);
- compName = new char[MED_TAILLE_PNOM*nbComp+1];
- compUnit = new char[MED_TAILLE_PNOM*nbComp+1];
-
- err = med_2_3::MEDfieldInfo(_medIdt,i,champName,meshName,&_local,&type,compName,compUnit,compDtunit,&nbComp);
- if (err == 0) {
- if (strcmp(champName, field->getName())==0) { // Found !
- Find = true ;
- break ;
- }
- }
- delete[] compName ;
- delete[] compUnit ;
- }
-
-
- if (Find)
- {
- // the same ?
- if (nbComp != component_count)
- throw MEDEXCEPTION ( LOCALIZED (STRING(LOC)
- <<": Field exist in file, but number of component are different : " \
- <<nbComp<<" in file and "<<component_count<<" in CORBA object."
- )
- );
-
- // component name and unit
- MESSAGE(LOC<<" Component name in file : "<<compName);
- MESSAGE(LOC<<" Component name in CORBA object : "<<component_name);
- MESSAGE(LOC<<" Component unit in file : "<<compUnit);
- MESSAGE(LOC<<" Component unit in CORBA object : "<<component_unit);
- delete[] compName ;
- delete[] compUnit ;
-
- }
- else
- {
- // Verify the field doesn't exist
- string dataGroupName = "/CHA/"+ string(field->getName());
-
- MESSAGE(LOC << "|" << dataGroupName << "|" );
-
- med_2_3::med_idt gid = H5Gopen(_medIdt, dataGroupName.c_str() );
-
- if ( gid < 0 )
- {
- // create field :
- err=med_2_3::MEDfieldCr(_medIdt,
- const_cast <char*> (field->getName()),
- med_2_3::MED_FLOAT64,
- component_count,
- const_cast <char*> ( component_name.c_str() ),
- const_cast <char*> ( component_unit.c_str() ),
- compDtunit,
- meshName);
- if ( err < 0 )
- throw MEDEXCEPTION( LOCALIZED (STRING(LOC) << ": Error MEDfieldCr : "<<err ));
- }
- else
- H5Gclose(gid);
- }
-
- SALOME_MED::SUPPORT_ptr mySupport = field -> getSupport();
- SCRUTE(mySupport);
-
- if (! mySupport->isOnAllElements())
- throw MEDEXCEPTION ( LOCALIZED (STRING(LOC)
- <<": Field must be on all entity"
- )
- );
- SALOME_MED::MESH_ptr myMesh = mySupport -> getMesh()->convertInMESH();
- SCRUTE(mesh);
-
- SALOME_MED::medGeometryElement_array* Types = mySupport->getTypes() ;
-
- int NumberOfType = Types->length() ;
-
- for ( int i = 0; i < NumberOfType; i++ )
- {
- int NumberOfElements = mySupport->getNumberOfElements ( (*Types)[i] ) ;
- //Temporarily commented to avoid compilation errors:int NumberOfGaussPoint = mySupport->getNumberOfGaussPoint ( (*Types)[i] ) ;
- int NumberOfGaussPoint = 1; //Temporary line instead of the above one
-
- MESSAGE( " " << field->getName() );
- MESSAGE( " " << NumberOfElements );
- //!!!tmp commented:MESSAGE( " " << NumberOfGaussPoint );
- MESSAGE( "Attention! NumberOfGaussPoint is invalid: " << NumberOfGaussPoint ); //tmp line!!!
- MESSAGE( " " << (int) (convertIdlEntToMedEnt(mySupport->getEntity())) );
- MESSAGE( " " << (int)(convertIdlEltToMedElt((*Types)[i])) );
- MESSAGE( " " << field->getIterationNumber() );
- MESSAGE( " " << field->getTime() );
- MESSAGE( " " << field->getOrderNumber() );
- MESSAGE( "MEDchampEcr :" << myMesh->getName() );
-
- SALOME_TYPES::ListOfDouble * value = field->getValue( SALOME_MED::MED_FULL_INTERLACE ) ;
- double *locvalue = new double[NumberOfElements];
- for (int k = 0; k < NumberOfElements; k++)
- locvalue[k] = (*value) [k];
-
- err=med_2_3::MEDfieldValueWithProfileWr(_medIdt,
- const_cast <char*> (field->getName()),
- MED_NO_DT,
- MED_NO_IT,
- field->getTime(),
- (med_2_3::med_entity_type)convertIdlEntToMedEnt(mySupport->getEntity()),
- (med_2_3::med_geometry_type)(convertIdlEltToMedElt((*Types)[i])),
- med_2_3::MED_COMPACT_PFLMODE,
- MED_ALLENTITIES_PROFILE,
- " ",
- med_2_3::MED_FULL_INTERLACE,
- SALOME_MED::MED_ALL_ELEMENTS,
- 1,
- (unsigned char*)locvalue
- );
-
- delete locvalue;
-
- if (err < MED_VALID )
- throw MEDEXCEPTION(LOCALIZED( STRING(LOC)
- <<": Error in writing Field "<< field->getName() <<", type "<<(*Types)[i]
- )
- );
- }
- END_OF(LOC);
-
- SCRUTE( err );
- if (err < 0 ) return;
-
- med_2_3::MEDfileClose(_medIdt);
-
- endService("CalculatorEngine::writeMEDfile");
- return;
-}
-
-extern "C"
-{
- //===========================================================================
- // CalculatorEngine_factory
- //===========================================================================
- PortableServer::ObjectId * CalculatorEngine_factory
- (CORBA::ORB_ptr orb,
- PortableServer::POA_ptr poa,
- PortableServer::ObjectId * contId,
- const char *instanceName,
- const char *interfaceName)
- {
- MESSAGE("CalculatorEngine_factory CalculatorEngine ("
- << instanceName << "," << interfaceName << ")");
- CalculatorEngine * myCalculator
- = new CalculatorEngine(orb, poa, contId, instanceName, interfaceName);
- return myCalculator->getId() ;
- }
-}
