/[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 797 by dhawcroft, Sat Jul 29 19:40:22 2006 UTC revision 798 by gross, Fri Aug 4 01:05:36 2006 UTC
# Line 252  void MeshAdapter::addPDEToSystem( Line 252  void MeshAdapter::addPDEToSystem(
252                       const escript::Data& d, const escript::Data& y,                       const escript::Data& d, const escript::Data& y,
253                       const escript::Data& d_contact,const escript::Data& y_contact) const                       const escript::Data& d_contact,const escript::Data& y_contact) const
254  {  {
255    escriptDataC _rhs=rhs.getDataC();     escriptDataC _rhs=rhs.getDataC();
256    escriptDataC _A  =A.getDataC();     escriptDataC _A  =A.getDataC();
257    escriptDataC _B=B.getDataC();     escriptDataC _B=B.getDataC();
258    escriptDataC _C=C.getDataC();     escriptDataC _C=C.getDataC();
259    escriptDataC _D=D.getDataC();     escriptDataC _D=D.getDataC();
260    escriptDataC _X=X.getDataC();     escriptDataC _X=X.getDataC();
261    escriptDataC _Y=Y.getDataC();     escriptDataC _Y=Y.getDataC();
262    escriptDataC _d=d.getDataC();     escriptDataC _d=d.getDataC();
263    escriptDataC _y=y.getDataC();     escriptDataC _y=y.getDataC();
264    escriptDataC _d_contact=d_contact.getDataC();     escriptDataC _d_contact=d_contact.getDataC();
265    escriptDataC _y_contact=y_contact.getDataC();     escriptDataC _y_contact=y_contact.getDataC();
266    
267     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
    Finley_Assemble_PDE(mesh->Nodes,mesh->Elements,mat.getPaso_SystemMatrix(), &_rhs, &_A, &_B, &_C, &_D, &_X, &_Y );  
268    
269       Finley_Assemble_PDE(mesh->Nodes,mesh->Elements,mat.getPaso_SystemMatrix(), &_rhs, &_A, &_B, &_C, &_D, &_X, &_Y );
270     checkFinleyError();     checkFinleyError();
271    
272     Finley_Assemble_RobinCondition(mesh->Nodes,mesh->FaceElements, mat.getPaso_SystemMatrix(), &_rhs, &_d, &_y, Finley_Assemble_handelShapeMissMatch_Mean_out);     Finley_Assemble_PDE(mesh->Nodes,mesh->FaceElements, mat.getPaso_SystemMatrix(), &_rhs, 0, 0, 0, &_d, 0, &_y );
273     checkFinleyError();     checkFinleyError();
274    
275     Finley_Assemble_RobinCondition(mesh->Nodes,mesh->ContactElements, mat.getPaso_SystemMatrix(), &_rhs , &_d_contact, &_y_contact ,             Finley_Assemble_handelShapeMissMatch_Step_out);     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, mat.getPaso_SystemMatrix(), &_rhs , 0, 0, 0, &_d_contact, 0, &_y_contact );
276     checkFinleyError();     checkFinleyError();
277  }  }
278    
# Line 283  void MeshAdapter::addPDEToRHS( escript:: Line 283  void MeshAdapter::addPDEToRHS( escript::
283  {  {
284     Finley_Mesh* mesh=m_finleyMesh.get();     Finley_Mesh* mesh=m_finleyMesh.get();
285    
286     // Finley_Assemble_PDE_RHS(mesh->Nodes,mesh->Elements,&(rhs.getDataC()),&(X.getDataC()),&(Y.getDataC()));     escriptDataC _rhs=rhs.getDataC();
287     Finley_Assemble_PDE(mesh->Nodes,mesh->Elements,0,&(rhs.getDataC()),0,0,0,0,&(X.getDataC()),&(Y.getDataC()));     escriptDataC _X=X.getDataC();
288     checkFinleyError();     escriptDataC _Y=Y.getDataC();
289       escriptDataC _y=y.getDataC();
290       escriptDataC _y_contact=y_contact.getDataC();
291    
292     // Finley_Assemble_RobinCondition_RHS(mesh->Nodes,mesh->FaceElements,&(rhs.getDataC()),&(y.getDataC()),Finley_Assemble_handelShapeMissMatch_Mean_out);     Finley_Assemble_PDE(mesh->Nodes,mesh->Elements, 0, &_rhs, 0, 0, 0, 0, &_X, &_Y );
293     Finley_Assemble_RobinCondition(mesh->Nodes,mesh->FaceElements,0,&(rhs.getDataC()),0,&(y.getDataC()),Finley_Assemble_handelShapeMissMatch_Mean_out);     checkFinleyError();
294    
295       Finley_Assemble_PDE(mesh->Nodes,mesh->FaceElements, 0, &_rhs, 0, 0, 0, 0, 0, &_y );
296     checkFinleyError();     checkFinleyError();
297     Finley_Assemble_RobinCondition(mesh->Nodes,mesh->ContactElements,0,&(rhs.getDataC()),0,&(y_contact.getDataC()),Finley_Assemble_handelShapeMissMatch_Step_out);  
298     // Finley_Assemble_RobinCondition_RHS(mesh->Nodes,mesh->ContactElements,&(rhs.getDataC()),&(y_contact.getDataC()),Finley_Assemble_handelShapeMissMatch_Step_out);     Finley_Assemble_PDE(mesh->Nodes,mesh->ContactElements, 0, &_rhs , 0, 0, 0, 0, 0, &_y_contact );
299     checkFinleyError();     checkFinleyError();
300  }  }
301    

Legend:
Removed from v.797  
changed lines
  Added in v.798

  ViewVC Help
Powered by ViewVC 1.1.26