/[escript]/branches/domexper/dudley/src/CPPAdapter/MeshAdapter.cpp
ViewVC logotype

Diff of /branches/domexper/dudley/src/CPPAdapter/MeshAdapter.cpp

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

revision 3090 by jfenwick, Wed Aug 11 00:51:28 2010 UTC revision 3114 by jfenwick, Fri Aug 27 05:26:25 2010 UTC
# Line 684  pair<int,int> MeshAdapter::getDataShape( Line 684  pair<int,int> MeshAdapter::getDataShape(
684  void MeshAdapter::addPDEToSystem(  void MeshAdapter::addPDEToSystem(
685                                   SystemMatrixAdapter& mat, escript::Data& rhs,                                   SystemMatrixAdapter& mat, escript::Data& rhs,
686                                   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,
687                                   const escript::Data& d, const escript::Data& y,                                   const escript::Data& d, const escript::Data& y) const
                                  const escript::Data& d_contact,const escript::Data& y_contact) const  
688  {  {
689     escriptDataC _rhs=rhs.getDataC();     escriptDataC _rhs=rhs.getDataC();
690     escriptDataC _A  =A.getDataC();     escriptDataC _A  =A.getDataC();
# Line 696  void MeshAdapter::addPDEToSystem( Line 695  void MeshAdapter::addPDEToSystem(
695     escriptDataC _Y=Y.getDataC();     escriptDataC _Y=Y.getDataC();
696     escriptDataC _d=d.getDataC();     escriptDataC _d=d.getDataC();
697     escriptDataC _y=y.getDataC();     escriptDataC _y=y.getDataC();
    escriptDataC _d_contact=d_contact.getDataC();  
    escriptDataC _y_contact=y_contact.getDataC();  
698    
699     Dudley_Mesh* mesh=m_dudleyMesh.get();     Dudley_Mesh* mesh=m_dudleyMesh.get();
700    
701     Dudley_Assemble_PDE(mesh->Nodes,mesh->Elements,mat.getPaso_SystemMatrix(), &_rhs, &_A, &_B, &_C, &_D, &_X, &_Y );     Dudley_Assemble_PDE(mesh->Nodes,mesh->Elements,mat.getPaso_SystemMatrix(), &_rhs, &_A, &_B, &_C, &_D, &_X, &_Y );
702     checkDudleyError();     checkDudleyError();
703    
704    
705     Dudley_Assemble_PDE(mesh->Nodes,mesh->FaceElements, mat.getPaso_SystemMatrix(), &_rhs, 0, 0, 0, &_d, 0, &_y );     Dudley_Assemble_PDE(mesh->Nodes,mesh->FaceElements, mat.getPaso_SystemMatrix(), &_rhs, 0, 0, 0, &_d, 0, &_y );
706     checkDudleyError();     checkDudleyError();
707    
708    
709     checkDudleyError();     checkDudleyError();
710  }  }
711    
# Line 741  void MeshAdapter::addPDEToRHS( escript:: Line 740  void MeshAdapter::addPDEToRHS( escript::
740     escriptDataC _X=X.getDataC();     escriptDataC _X=X.getDataC();
741     escriptDataC _Y=Y.getDataC();     escriptDataC _Y=Y.getDataC();
742     escriptDataC _y=y.getDataC();     escriptDataC _y=y.getDataC();
743     escriptDataC _y_contact=y_contact.getDataC();  //   escriptDataC _y_contact=y_contact.getDataC();
744    
745     Dudley_Assemble_PDE(mesh->Nodes,mesh->Elements, 0, &_rhs, 0, 0, 0, 0, &_X, &_Y );     Dudley_Assemble_PDE(mesh->Nodes,mesh->Elements, 0, &_rhs, 0, 0, 0, 0, &_X, &_Y );
746     checkDudleyError();     checkDudleyError();

Legend:
Removed from v.3090  
changed lines
  Added in v.3114

  ViewVC Help
Powered by ViewVC 1.1.26