/[escript]/branches/diaplayground/finley/src/CPPAdapter/MeshAdapter.cpp
ViewVC logotype

Diff of /branches/diaplayground/finley/src/CPPAdapter/MeshAdapter.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3344 by caltinay, Thu Nov 11 23:26:52 2010 UTC revision 3490 by caltinay, Wed Mar 30 02:24:33 2011 UTC
# Line 23  extern "C" { Line 23  extern "C" {
23  }  }
24    
25  #include <boost/python/import.hpp>  #include <boost/python/import.hpp>
 #include <boost/python/tuple.hpp>  
26    
27  using namespace std;  using namespace std;
28  using namespace escript;  using namespace escript;
# Line 806  void MeshAdapter::addPDEToSystem( Line 805  void MeshAdapter::addPDEToSystem(
805  void  MeshAdapter::addPDEToLumpedSystem(  void  MeshAdapter::addPDEToLumpedSystem(
806                                          escript::Data& mat,                                          escript::Data& mat,
807                                          const escript::Data& D,                                          const escript::Data& D,
808                                          const escript::Data& d) const                                          const escript::Data& d,
809                                            const bool useHRZ) const
810  {  {
811     escriptDataC _mat=mat.getDataC();     escriptDataC _mat=mat.getDataC();
812     escriptDataC _D=D.getDataC();     escriptDataC _D=D.getDataC();
# Line 814  void  MeshAdapter::addPDEToLumpedSystem( Line 814  void  MeshAdapter::addPDEToLumpedSystem(
814    
815     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
816    
817     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->Elements,&_mat, &_D);     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->Elements,&_mat, &_D, useHRZ);
818     checkFinleyError();     checkFinleyError();
819        
820     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->FaceElements,&_mat, &_d);     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->FaceElements,&_mat, &_d, useHRZ);
821     checkFinleyError();     checkFinleyError();
822    
823  }  }
# Line 1075  void MeshAdapter::interpolateOnDomain(es Line 1075  void MeshAdapter::interpolateOnDomain(es
1075        break;        break;
1076        case(Points):        case(Points):
1077        if (getMPISize()>1) {        if (getMPISize()>1) {
1078           escript::Data temp=escript::Data( in,  continuousFunction(*this) );           //escript::Data temp=escript::Data( in,  continuousFunction(*this) );
1079           escriptDataC _in2 = temp.getDataC();           //escriptDataC _in2 = temp.getDataC();
1080        } else {        } else {
1081           Finley_Assemble_interpolate(mesh->Nodes,mesh->Points,&_in,&_target);           Finley_Assemble_interpolate(mesh->Nodes,mesh->Points,&_in,&_target);
1082        }        }
# Line 1549  void MeshAdapter::saveDX(const string& f Line 1549  void MeshAdapter::saveDX(const string& f
1549  //  //
1550  void MeshAdapter::saveVTK(const string& filename,const boost::python::dict& arg,  const string& metadata, const string& metadata_schema) const  void MeshAdapter::saveVTK(const string& filename,const boost::python::dict& arg,  const string& metadata, const string& metadata_schema) const
1551  {  {
1552      boost::python::object pySaveVTK = boost::python::import("esys.weipa").attr("saveVTK");      boost::python::object pySaveVTK = boost::python::import("esys.weipa").attr("_saveVTK");
1553      pySaveVTK(*boost::python::make_tuple(filename,      pySaveVTK(filename, const_cast<MeshAdapter*>(this)->getPtr(),
1554                 const_cast<MeshAdapter*>(this)->getPtr(),                metadata, metadata_schema, arg);
                metadata, metadata_schema), **arg);  
1555  }  }
1556    
1557  bool MeshAdapter::ownSample(int fs_code, index_t id) const  bool MeshAdapter::ownSample(int fs_code, index_t id) const

Legend:
Removed from v.3344  
changed lines
  Added in v.3490

  ViewVC Help
Powered by ViewVC 1.1.26