/[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 4818 by caltinay, Mon Mar 31 00:16:20 2014 UTC revision 4836 by caltinay, Mon Apr 7 05:51:55 2014 UTC
# Line 752  void MeshAdapter::addPDEToSystem( Line 752  void MeshAdapter::addPDEToSystem(
752          throw FinleyAdapterException("finley only supports Paso system matrices.");          throw FinleyAdapterException("finley only supports Paso system matrices.");
753    
754      Mesh* mesh=m_finleyMesh.get();      Mesh* mesh=m_finleyMesh.get();
755      Paso_SystemMatrix* S = smat->getPaso_SystemMatrix();      paso::SystemMatrix_ptr S(smat->getPaso_SystemMatrix());
756      Assemble_PDE(mesh->Nodes, mesh->Elements, S, rhs, A, B, C, D, X, Y);      Assemble_PDE(mesh->Nodes, mesh->Elements, S, rhs, A, B, C, D, X, Y);
757      checkFinleyError();      checkFinleyError();
758    
# Line 795  void MeshAdapter::addPDEToRHS(escript::D Line 795  void MeshAdapter::addPDEToRHS(escript::D
795          const escript::Data& y_contact, const escript::Data& y_dirac) const          const escript::Data& y_contact, const escript::Data& y_dirac) const
796  {  {
797      Mesh* mesh=m_finleyMesh.get();      Mesh* mesh=m_finleyMesh.get();
798      Assemble_PDE(mesh->Nodes, mesh->Elements, 0, rhs,      Assemble_PDE(mesh->Nodes, mesh->Elements, paso::SystemMatrix_ptr(), rhs,
799              escript::Data(), escript::Data(), escript::Data(), escript::Data(),              escript::Data(), escript::Data(), escript::Data(), escript::Data(),
800              X, Y);              X, Y);
801      checkFinleyError();      checkFinleyError();
802    
803      Assemble_PDE(mesh->Nodes, mesh->FaceElements, 0, rhs,      Assemble_PDE(mesh->Nodes, mesh->FaceElements, paso::SystemMatrix_ptr(),
804              escript::Data(), escript::Data(), escript::Data(), escript::Data(),              rhs, escript::Data(), escript::Data(), escript::Data(),
805              escript::Data(), y);              escript::Data(), escript::Data(), y);
806      checkFinleyError();      checkFinleyError();
807    
808      Assemble_PDE(mesh->Nodes, mesh->ContactElements, 0, rhs,      Assemble_PDE(mesh->Nodes, mesh->ContactElements, paso::SystemMatrix_ptr(),
809              escript::Data(), escript::Data(), escript::Data(),              rhs, escript::Data(), escript::Data(), escript::Data(),
810              escript::Data(), escript::Data(), y_contact);              escript::Data(), escript::Data(), y_contact);
811      checkFinleyError();      checkFinleyError();
812    
813      Assemble_PDE(mesh->Nodes, mesh->Points, 0, rhs,      Assemble_PDE(mesh->Nodes, mesh->Points, paso::SystemMatrix_ptr(), rhs,
814              escript::Data(), escript::Data(), escript::Data(), escript::Data(),              escript::Data(), escript::Data(), escript::Data(), escript::Data(),
815              escript::Data(), y_dirac);              escript::Data(), y_dirac);
816      checkFinleyError();      checkFinleyError();
# Line 1458  escript::ASM_ptr MeshAdapter::newSystemM Line 1458  escript::ASM_ptr MeshAdapter::newSystemM
1458      paso::SystemMatrixPattern_ptr pattern = getFinley_Mesh()->getPattern(      paso::SystemMatrixPattern_ptr pattern = getFinley_Mesh()->getPattern(
1459              reduceRowOrder, reduceColOrder);              reduceRowOrder, reduceColOrder);
1460      checkFinleyError();      checkFinleyError();
1461      Paso_SystemMatrix* fsystemMatrix;      paso::SystemMatrix_ptr fsystemMatrix;
1462      const int trilinos = 0;      const int trilinos = 0;
1463      if (trilinos) {      if (trilinos) {
1464  #ifdef TRILINOS  #ifdef TRILINOS
1465          // FIXME: Allocation Epetra_VrbMatrix here...          // FIXME: Allocation Epetra_VrbMatrix here...
1466  #endif  #endif
1467      } else {      } else {
1468          fsystemMatrix=Paso_SystemMatrix_alloc(type, pattern, row_blocksize,          fsystemMatrix.reset(new paso::SystemMatrix(type, pattern,
1469                  column_blocksize, FALSE);                              row_blocksize, column_blocksize, false));
1470      }      }
1471      checkPasoError();      checkPasoError();
1472      SystemMatrixAdapter* sma=new SystemMatrixAdapter(fsystemMatrix, row_blocksize, row_functionspace, column_blocksize, column_functionspace);      SystemMatrixAdapter* sma=new SystemMatrixAdapter(fsystemMatrix, row_blocksize, row_functionspace, column_blocksize, column_functionspace);

Legend:
Removed from v.4818  
changed lines
  Added in v.4836

  ViewVC Help
Powered by ViewVC 1.1.26