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

Diff of /trunk/finley/src/CPPAdapter/TransportProblemAdapter.cpp

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

revision 1512 by gross, Mon Feb 25 04:45:48 2008 UTC revision 1513 by gross, Tue Apr 15 08:47:57 2008 UTC
# Line 118  void TransportProblemAdapter::copyConstr Line 118  void TransportProblemAdapter::copyConstr
118      /* insert 0 rows and cols into transport matrix */      /* insert 0 rows and cols into transport matrix */
119      q.expand();      q.expand();
120      double* q_dp=q.getSampleData(0);      double* q_dp=q.getSampleData(0);
121      Paso_SystemMatrix_nullifyRowsAndCols(transp->transport_matrix,q_dp,q_dp,0.);      Paso_SystemMatrix_nullifyRows(transp->transport_matrix,q_dp, 0.);
122      checkPasoError();      checkPasoError();
123    
124      /* insert 0 rows amd 1 in main diagonal into mass matrix */      /* insert 0 rows amd 1 in main diagonal into mass matrix */
125      Paso_SystemMatrix_nullifyRows(transp->mass_matrix,q_dp, 1.);      Paso_SystemMatrix_nullifyRowsAndCols(transp->mass_matrix,q_dp,q_dp,1.);
126      checkPasoError();      checkPasoError();
127    
128      source.copyWithMask(escript::Data(0.,q.getDataPointShape(),q.getFunctionSpace()),q);      source.copyWithMask(escript::Data(0.,q.getDataPointShape(),q.getFunctionSpace()),q);

Legend:
Removed from v.1512  
changed lines
  Added in v.1513

  ViewVC Help
Powered by ViewVC 1.1.26