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

Diff of /trunk/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 3981 by jfenwick, Fri Sep 21 02:47:54 2012 UTC
# Line 1  Line 1 
1    
2  /*******************************************************  /*****************************************************************************
3  *  *
4  * Copyright (c) 2003-2010 by University of Queensland  * Copyright (c) 2003-2012 by University of Queensland
5  * Earth Systems Science Computational Center (ESSCC)  * http://www.uq.edu.au
 * http://www.uq.edu.au/esscc  
6  *  *
7  * Primary Business: Queensland, Australia  * Primary Business: Queensland, Australia
8  * Licensed under the Open Software License version 3.0  * Licensed under the Open Software License version 3.0
9  * http://www.opensource.org/licenses/osl-3.0.php  * http://www.opensource.org/licenses/osl-3.0.php
10  *  *
11  *******************************************************/  * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12    * Development since 2012 by School of Earth Sciences
13    *
14    *****************************************************************************/
15    
16    #include <pasowrap/PasoException.h>
17    #include <pasowrap/TransportProblemAdapter.h>
18  #include "MeshAdapter.h"  #include "MeshAdapter.h"
19  #include "escript/Data.h"  #include "escript/Data.h"
20  #include "escript/DataFactory.h"  #include "escript/DataFactory.h"
# Line 23  extern "C" { Line 26  extern "C" {
26  }  }
27    
28  #include <boost/python/import.hpp>  #include <boost/python/import.hpp>
 #include <boost/python/tuple.hpp>  
29    
30  using namespace std;  using namespace std;
31  using namespace escript;  using namespace escript;
32    using namespace paso;
33    
34  namespace finley {  namespace finley {
35    
# Line 648  int MeshAdapter::getReducedSolutionCode( Line 651  int MeshAdapter::getReducedSolutionCode(
651     return ReducedDegreesOfFreedom;     return ReducedDegreesOfFreedom;
652  }  }
653    
654  int MeshAdapter::getDiracDeltaFunctionCode() const  int MeshAdapter::getDiracDeltaFunctionsCode() const
655  {  {
656     return Points;     return Points;
657  }  }
# Line 772  void MeshAdapter::addPDEToSystem( Line 775  void MeshAdapter::addPDEToSystem(
775                                   AbstractSystemMatrix& mat, escript::Data& rhs,                                   AbstractSystemMatrix& mat, escript::Data& rhs,
776                                   const escript::Data& A, const escript::Data& B, const escript::Data& C,const  escript::Data& D,const  escript::Data& X,const  escript::Data& Y,                                   const escript::Data& A, const escript::Data& B, const escript::Data& C,const  escript::Data& D,const  escript::Data& X,const  escript::Data& Y,
777                                   const escript::Data& d, const escript::Data& y,                                   const escript::Data& d, const escript::Data& y,
778                                   const escript::Data& d_contact,const escript::Data& y_contact) const                                   const escript::Data& d_contact,const escript::Data& y_contact,
779                                     const escript::Data& d_dirac,const escript::Data& y_dirac) const
780  {  {
781     SystemMatrixAdapter* smat=dynamic_cast<SystemMatrixAdapter*>(&mat);     SystemMatrixAdapter* smat=dynamic_cast<SystemMatrixAdapter*>(&mat);
782     if (smat==0)     if (smat==0)
783     {     {
784      throw FinleyAdapterException("finley only supports its own system matrices.");      throw FinleyAdapterException("finley only supports Paso system matrices.");
785     }     }
786     escriptDataC _rhs=rhs.getDataC();     escriptDataC _rhs=rhs.getDataC();
787     escriptDataC _A  =A.getDataC();     escriptDataC _A  =A.getDataC();
# Line 790  void MeshAdapter::addPDEToSystem( Line 794  void MeshAdapter::addPDEToSystem(
794     escriptDataC _y=y.getDataC();     escriptDataC _y=y.getDataC();
795     escriptDataC _d_contact=d_contact.getDataC();     escriptDataC _d_contact=d_contact.getDataC();
796     escriptDataC _y_contact=y_contact.getDataC();     escriptDataC _y_contact=y_contact.getDataC();
797       escriptDataC _d_dirac=d_dirac.getDataC();
798       escriptDataC _y_dirac=y_dirac.getDataC();
799    
800     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
801    
# Line 801  void MeshAdapter::addPDEToSystem( Line 807  void MeshAdapter::addPDEToSystem(
807    
808     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, smat->getPaso_SystemMatrix(), &_rhs , 0, 0, 0, &_d_contact, 0, &_y_contact );     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, smat->getPaso_SystemMatrix(), &_rhs , 0, 0, 0, &_d_contact, 0, &_y_contact );
809     checkFinleyError();     checkFinleyError();
810    
811        Finley_Assemble_PDE(mesh->Nodes,mesh->Points, smat->getPaso_SystemMatrix(), &_rhs , 0, 0, 0, &_d_dirac, 0, &_y_dirac );
812       checkFinleyError();
813  }  }
814    
815  void  MeshAdapter::addPDEToLumpedSystem(  void  MeshAdapter::addPDEToLumpedSystem(
816                                          escript::Data& mat,                                          escript::Data& mat,
817                                          const escript::Data& D,                                          const escript::Data& D,
818                                          const escript::Data& d) const                                          const escript::Data& d,
819                                            const escript::Data& d_dirac,
820                                            const bool useHRZ) const
821  {  {
822     escriptDataC _mat=mat.getDataC();     escriptDataC _mat=mat.getDataC();
823     escriptDataC _D=D.getDataC();     escriptDataC _D=D.getDataC();
824     escriptDataC _d=d.getDataC();     escriptDataC _d=d.getDataC();
825       escriptDataC _d_dirac=d_dirac.getDataC();
826    
827     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
828    
829     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->Elements,&_mat, &_D);     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->Elements,&_mat, &_D, useHRZ);
830     checkFinleyError();     checkFinleyError();
831        
832     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->FaceElements,&_mat, &_d);     Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->FaceElements,&_mat, &_d, useHRZ);
833       checkFinleyError();
834    
835       Finley_Assemble_LumpedSystem(mesh->Nodes,mesh->Points,&_mat, &_d_dirac, useHRZ);
836     checkFinleyError();     checkFinleyError();
837    
838  }  }
# Line 826  void  MeshAdapter::addPDEToLumpedSystem( Line 841  void  MeshAdapter::addPDEToLumpedSystem(
841  //  //
842  // adds linear PDE of second order into the right hand side only  // adds linear PDE of second order into the right hand side only
843  //  //
844  void MeshAdapter::addPDEToRHS( escript::Data& rhs, const  escript::Data& X,const  escript::Data& Y, const escript::Data& y, const escript::Data& y_contact) const  void MeshAdapter::addPDEToRHS( escript::Data& rhs, const  escript::Data& X,const  escript::Data& Y, const escript::Data& y, const escript::Data& y_contact,  const escript::Data& y_dirac) const
845  {  {
846     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
847    
# Line 835  void MeshAdapter::addPDEToRHS( escript:: Line 850  void MeshAdapter::addPDEToRHS( escript::
850     escriptDataC _Y=Y.getDataC();     escriptDataC _Y=Y.getDataC();
851     escriptDataC _y=y.getDataC();     escriptDataC _y=y.getDataC();
852     escriptDataC _y_contact=y_contact.getDataC();     escriptDataC _y_contact=y_contact.getDataC();
853       escriptDataC _y_dirac=y_dirac.getDataC();
854    
855     Finley_Assemble_PDE(mesh->Nodes,mesh->Elements, 0, &_rhs, 0, 0, 0, 0, &_X, &_Y );     Finley_Assemble_PDE(mesh->Nodes,mesh->Elements, 0, &_rhs, 0, 0, 0, 0, &_X, &_Y );
856     checkFinleyError();     checkFinleyError();
# Line 844  void MeshAdapter::addPDEToRHS( escript:: Line 860  void MeshAdapter::addPDEToRHS( escript::
860    
861     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, 0, &_rhs , 0, 0, 0, 0, 0, &_y_contact );     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, 0, &_rhs , 0, 0, 0, 0, 0, &_y_contact );
862     checkFinleyError();     checkFinleyError();
863    
864       Finley_Assemble_PDE(mesh->Nodes,mesh->Points, 0, &_rhs , 0, 0, 0, 0, 0, &_y_dirac );
865       checkFinleyError();
866    
867  }  }
868  //  //
869  // adds PDE of second order into a transport problem  // adds PDE of second order into a transport problem
# Line 853  void MeshAdapter::addPDEToTransportProbl Line 873  void MeshAdapter::addPDEToTransportProbl
873                                             const escript::Data& A, const escript::Data& B, const escript::Data& C,                                             const escript::Data& A, const escript::Data& B, const escript::Data& C,
874                                             const  escript::Data& D,const  escript::Data& X,const  escript::Data& Y,                                             const  escript::Data& D,const  escript::Data& X,const  escript::Data& Y,
875                                             const escript::Data& d, const escript::Data& y,                                             const escript::Data& d, const escript::Data& y,
876                                             const escript::Data& d_contact,const escript::Data& y_contact) const                                             const escript::Data& d_contact,const escript::Data& y_contact,
877                                               const escript::Data& d_dirac, const escript::Data& y_dirac) const
878  {  {
879     TransportProblemAdapter* tpa=dynamic_cast<TransportProblemAdapter*>(&tp);     TransportProblemAdapter* tpa=dynamic_cast<TransportProblemAdapter*>(&tp);
880     if (tpa==0)     if (tpa==0)
881     {     {
882      throw FinleyAdapterException("finley only supports its own transport problems.");      throw FinleyAdapterException("finley only supports Paso transport problems.");
883     }     }
884    
885    
# Line 876  void MeshAdapter::addPDEToTransportProbl Line 897  void MeshAdapter::addPDEToTransportProbl
897     escriptDataC _y=y.getDataC();     escriptDataC _y=y.getDataC();
898     escriptDataC _d_contact=d_contact.getDataC();     escriptDataC _d_contact=d_contact.getDataC();
899     escriptDataC _y_contact=y_contact.getDataC();     escriptDataC _y_contact=y_contact.getDataC();
900       escriptDataC _d_dirac=d_dirac.getDataC();
901       escriptDataC _y_dirac=y_dirac.getDataC();
902    
903    
904     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
905     Paso_TransportProblem* _tp = tpa->getPaso_TransportProblem();     Paso_TransportProblem* _tp = tpa->getPaso_TransportProblem();
# Line 891  void MeshAdapter::addPDEToTransportProbl Line 915  void MeshAdapter::addPDEToTransportProbl
915    
916     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, _tp->transport_matrix, &_source , 0, 0, 0, &_d_contact, 0, &_y_contact );     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, _tp->transport_matrix, &_source , 0, 0, 0, &_d_contact, 0, &_y_contact );
917     checkFinleyError();     checkFinleyError();
918    
919       Finley_Assemble_PDE(mesh->Nodes,mesh->Points, _tp->transport_matrix, &_source , 0, 0, 0, &_d_dirac, 0, &_y_dirac );
920       checkFinleyError();
921    
922  }  }
923    
924  //  //
# Line 1075  void MeshAdapter::interpolateOnDomain(es Line 1103  void MeshAdapter::interpolateOnDomain(es
1103        break;        break;
1104        case(Points):        case(Points):
1105        if (getMPISize()>1) {        if (getMPISize()>1) {
1106           escript::Data temp=escript::Data( in,  continuousFunction(*this) );           //escript::Data temp=escript::Data( in,  continuousFunction(*this) );
1107           escriptDataC _in2 = temp.getDataC();           //escriptDataC _in2 = temp.getDataC();
1108        } else {        } else {
1109           Finley_Assemble_interpolate(mesh->Nodes,mesh->Points,&_in,&_target);           Finley_Assemble_interpolate(mesh->Nodes,mesh->Points,&_in,&_target);
1110        }        }
# Line 1486  void MeshAdapter::setNewX(const escript: Line 1514  void MeshAdapter::setNewX(const escript:
1514         tmp = new_x.getDataC();         tmp = new_x.getDataC();
1515         Finley_Mesh_setCoordinates(mesh,&tmp);         Finley_Mesh_setCoordinates(mesh,&tmp);
1516     } else {     } else {
1517         escript::Data new_x_inter=escript::Data( new_x,  continuousFunction(*this) );         throw FinleyAdapterException("As of version escript3.3 - SetNewX only accepts ContinuousFunction arguments please interpolate.");    
1518    /*       escript::Data new_x_inter=escript::Data( new_x,  continuousFunction(*this) );
1519         tmp = new_x_inter.getDataC();         tmp = new_x_inter.getDataC();
1520         Finley_Mesh_setCoordinates(mesh,&tmp);         Finley_Mesh_setCoordinates(mesh,&tmp);*/
1521     }     }
1522     checkFinleyError();     checkFinleyError();
1523  }  }
# Line 1549  void MeshAdapter::saveDX(const string& f Line 1578  void MeshAdapter::saveDX(const string& f
1578  //  //
1579  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
1580  {  {
1581      boost::python::object pySaveVTK = boost::python::import("esys.weipa").attr("saveVTK");      boost::python::object pySaveVTK = boost::python::import("esys.weipa").attr("_saveVTK");
1582      pySaveVTK(*boost::python::make_tuple(filename,      pySaveVTK(filename, const_cast<MeshAdapter*>(this)->getPtr(),
1583                 const_cast<MeshAdapter*>(this)->getPtr(),                metadata, metadata_schema, arg);
                metadata, metadata_schema), **arg);  
1584  }  }
1585    
1586  bool MeshAdapter::ownSample(int fs_code, index_t id) const  bool MeshAdapter::ownSample(int fs_code, index_t id) const
# Line 1640  ASM_ptr MeshAdapter::newSystemMatrix( Line 1668  ASM_ptr MeshAdapter::newSystemMatrix(
1668  // creates a TransportProblemAdapter  // creates a TransportProblemAdapter
1669  //  //
1670  ATP_ptr MeshAdapter::newTransportProblem(  ATP_ptr MeshAdapter::newTransportProblem(
                                                          const bool useBackwardEuler,  
1671                                                           const int blocksize,                                                           const int blocksize,
1672                                                           const escript::FunctionSpace& functionspace,                                                           const escript::FunctionSpace& functionspace,
1673                                                           const int type) const                                                           const int type) const
# Line 1663  ATP_ptr MeshAdapter::newTransportProblem Line 1690  ATP_ptr MeshAdapter::newTransportProblem
1690     Paso_SystemMatrixPattern* fsystemMatrixPattern=Finley_getPattern(getFinley_Mesh(),reduceOrder,reduceOrder);     Paso_SystemMatrixPattern* fsystemMatrixPattern=Finley_getPattern(getFinley_Mesh(),reduceOrder,reduceOrder);
1691     checkFinleyError();     checkFinleyError();
1692     Paso_TransportProblem* transportProblem;     Paso_TransportProblem* transportProblem;
1693     transportProblem=Paso_TransportProblem_alloc(useBackwardEuler,fsystemMatrixPattern,blocksize);     transportProblem=Paso_TransportProblem_alloc(fsystemMatrixPattern,blocksize);
1694     checkPasoError();     checkPasoError();
1695     Paso_SystemMatrixPattern_free(fsystemMatrixPattern);     Paso_SystemMatrixPattern_free(fsystemMatrixPattern);
1696     TransportProblemAdapter* tpa=new TransportProblemAdapter(transportProblem,useBackwardEuler,blocksize,functionspace);     TransportProblemAdapter* tpa=new TransportProblemAdapter(transportProblem,blocksize,functionspace);
1697     return ATP_ptr(tpa);     return ATP_ptr(tpa);
1698  //   return TransportProblemAdapter(transportProblem,useBackwardEuler,blocksize,functionspace);  //   return TransportProblemAdapter(transportProblem, blocksize,functionspace);
1699  }  }
1700    
1701  //  //
# Line 2027  bool MeshAdapter::operator!=(const Abstr Line 2054  bool MeshAdapter::operator!=(const Abstr
2054  int MeshAdapter::getSystemMatrixTypeId(const int solver, const int preconditioner, const int package, const bool symmetry) const  int MeshAdapter::getSystemMatrixTypeId(const int solver, const int preconditioner, const int package, const bool symmetry) const
2055  {  {
2056     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
2057     int out=Paso_SystemMatrix_getSystemMatrixTypeId(SystemMatrixAdapter::mapOptionToPaso(solver),SystemMatrixAdapter::mapOptionToPaso(preconditioner), SystemMatrixAdapter::mapOptionToPaso(package),symmetry?1:0, mesh->MPIInfo);     return SystemMatrixAdapter::getSystemMatrixTypeId(solver, preconditioner,
2058     checkPasoError();             package, symmetry, mesh->MPIInfo);
    return out;  
2059  }  }
2060    
2061  int MeshAdapter::getTransportTypeId(const int solver, const int preconditioner, const int package, const bool symmetry) const  int MeshAdapter::getTransportTypeId(const int solver, const int preconditioner, const int package, const bool symmetry) const
2062  {  {
2063     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
2064     int out=Paso_TransportProblem_getTypeId(SystemMatrixAdapter::mapOptionToPaso(solver),SystemMatrixAdapter::mapOptionToPaso(preconditioner), SystemMatrixAdapter::mapOptionToPaso(package),symmetry?1:0, mesh->MPIInfo);     return TransportProblemAdapter::getTransportTypeId(solver, preconditioner,
2065     checkPasoError();             package, symmetry, mesh->MPIInfo);
    return out;  
2066  }  }
2067    
2068  escript::Data MeshAdapter::getX() const  escript::Data MeshAdapter::getX() const
# Line 2217  void MeshAdapter::setTagMap(const string Line 2243  void MeshAdapter::setTagMap(const string
2243  {  {
2244     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
2245     Finley_Mesh_addTagMap(mesh, name.c_str(),tag);     Finley_Mesh_addTagMap(mesh, name.c_str(),tag);
2246     checkPasoError();     checkFinleyError();
2247     // throwStandardException("MeshAdapter::set TagMap is not implemented.");     // throwStandardException("MeshAdapter::set TagMap is not implemented.");
2248  }  }
2249    
# Line 2226  int MeshAdapter::getTag(const string& na Line 2252  int MeshAdapter::getTag(const string& na
2252     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
2253     int tag=0;     int tag=0;
2254     tag=Finley_Mesh_getTag(mesh, name.c_str());     tag=Finley_Mesh_getTag(mesh, name.c_str());
2255     checkPasoError();     checkFinleyError();
2256     // throwStandardException("MeshAdapter::getTag is not implemented.");     // throwStandardException("MeshAdapter::getTag is not implemented.");
2257     return tag;     return tag;
2258  }  }
# Line 2414  ReferenceElementSetWrapper::~ReferenceEl Line 2440  ReferenceElementSetWrapper::~ReferenceEl
2440    Finley_ReferenceElementSet_dealloc(m_refSet);    Finley_ReferenceElementSet_dealloc(m_refSet);
2441  }  }
2442    
2443    // points will be flattened
2444    void MeshAdapter:: addDiracPoints(const std::vector<double>& points, const std::vector<int>& tags) const
2445    {
2446          const int dim = getDim();
2447          int numPoints=points.size()/dim;
2448          int numTags=tags.size();
2449          Finley_Mesh* mesh=m_finleyMesh.get();
2450          
2451          if ( points.size() % dim != 0 )
2452          {
2453        throw FinleyAdapterException("Error - number of coords does not appear to be a multiple of dimension.");
2454          }
2455          
2456          if  ( (numTags > 0) && ( numPoints !=  numTags ) )
2457         throw FinleyAdapterException("Error - if tags are given number of tags and points must match.");
2458          
2459          double* points_ptr=TMPMEMALLOC(numPoints * dim, double);
2460          int*    tags_ptr= TMPMEMALLOC(numPoints, int);
2461          
2462          for (int i=0;i<numPoints;++i) {
2463           points_ptr[ i * dim     ] = points[i * dim ];
2464           if ( dim > 1 ) points_ptr[ i * dim + 1 ] = points[i * dim + 1];
2465           if ( dim > 2 ) points_ptr[ i * dim + 2 ] = points[i * dim + 2];    
2466               tags_ptr[i]=tags[i];
2467          }
2468          
2469          Finley_Mesh_addPoints(mesh, numPoints, points_ptr, tags_ptr);
2470          checkFinleyError();
2471          
2472          TMPMEMFREE(points_ptr);
2473          TMPMEMFREE(tags_ptr);
2474    }
2475    
2476    
2477    // void MeshAdapter:: addDiracPoints(const boost::python::list& points, const boost::python::list& tags) const
2478    // {
2479    //       const int dim = getDim();
2480    //       int numPoints=boost::python::extract<int>(points.attr("__len__")());
2481    //       int numTags=boost::python::extract<int>(tags.attr("__len__")());
2482    //       Finley_Mesh* mesh=m_finleyMesh.get();
2483    //      
2484    //       if  ( (numTags > 0) && ( numPoints !=  numTags ) )
2485    //   throw FinleyAdapterException("Error - if tags are given number of tags and points must match.");
2486    //      
2487    //       double* points_ptr=TMPMEMALLOC(numPoints * dim, double);
2488    //       int*    tags_ptr= TMPMEMALLOC(numPoints, int);
2489    //      
2490    //       for (int i=0;i<numPoints;++i) {
2491    //     int tag_id=-1;
2492    //     int numComps=boost::python::extract<int>(points[i].attr("__len__")());
2493    //     if  ( numComps !=   dim ) {
2494    //                stringstream temp;            
2495    //                temp << "Error - illegal number of components " << numComps << " for point " << i;              
2496    //                throw FinleyAdapterException(temp.str());
2497    //     }
2498    //     points_ptr[ i * dim     ] = boost::python::extract<double>(points[i][0]);
2499    //     if ( dim > 1 ) points_ptr[ i * dim + 1 ] = boost::python::extract<double>(points[i][1]);
2500    //     if ( dim > 2 ) points_ptr[ i * dim + 2 ] = boost::python::extract<double>(points[i][2]);
2501    //    
2502    //     if ( numTags > 0) {
2503    //            boost::python::extract<string> ex_str(tags[i]);
2504    //        if  ( ex_str.check() ) {
2505    //            tag_id=getTag( ex_str());
2506    //        } else {
2507    //             boost::python::extract<int> ex_int(tags[i]);
2508    //             if ( ex_int.check() ) {
2509    //                 tag_id=ex_int();
2510    //             } else {
2511    //              stringstream temp;          
2512    //                  temp << "Error - unable to extract tag for point " << i;
2513    //              throw FinleyAdapterException(temp.str());
2514    //            }
2515    //        }
2516    //     }      
2517    //            tags_ptr[i]=tag_id;
2518    //       }
2519    //      
2520    //       Finley_Mesh_addPoints(mesh, numPoints, points_ptr, tags_ptr);
2521    //       checkPasoError();
2522    //      
2523    //       TMPMEMFREE(points_ptr);
2524    //       TMPMEMFREE(tags_ptr);
2525    // }
2526    
2527    /*
2528    void MeshAdapter:: addDiracPoint( const boost::python::list& point, const int tag) const
2529    {  
2530        boost::python::list points =  boost::python::list();
2531        boost::python::list tags =  boost::python::list();
2532        points.append(point);
2533        tags.append(tag);
2534        addDiracPoints(points, tags);
2535    }
2536    */
2537    
2538    /*
2539    void MeshAdapter:: addDiracPointWithTagName( const boost::python::list& point, const std::string& tag) const
2540    {
2541            boost::python::list points =   boost::python::list();
2542            boost::python::list tags =   boost::python::list();
2543            points.append(point);
2544            tags.append(tag);
2545            addDiracPoints(points, tags);
2546    }
2547    */
2548  }  // end of namespace  }  // end of namespace

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

  ViewVC Help
Powered by ViewVC 1.1.26