/[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 964 by gross, Tue Feb 13 05:10:26 2007 UTC revision 969 by ksteube, Tue Feb 13 23:02:23 2007 UTC
# Line 839  SystemMatrixAdapter MeshAdapter::newSyst Line 839  SystemMatrixAdapter MeshAdapter::newSyst
839            
840      Paso_SystemMatrixPattern* fsystemMatrixPattern=Finley_getPattern(getFinley_Mesh(),reduceRowOrder,reduceColOrder);      Paso_SystemMatrixPattern* fsystemMatrixPattern=Finley_getPattern(getFinley_Mesh(),reduceRowOrder,reduceColOrder);
841      checkFinleyError();      checkFinleyError();
842      Paso_SystemMatrix* fsystemMatrix=Paso_SystemMatrix_alloc(type,fsystemMatrixPattern,row_blocksize,column_blocksize);      printf("ksteube in finley/src/CPPAdapter/MeshAdapter.cpp newSystemMatrix\n");
843        Paso_SystemMatrix* fsystemMatrix;
844        int trilinos = 0;
845        if (trilinos) {
846    #ifdef TRILINOS
847          /* printf("ksteube newSystemMatrix package=%d\n", row_functionspace.domain.package); */
848          /* Allocation Epetra_VrbMatrix here */
849    #endif
850        }
851        else {
852          fsystemMatrix=Paso_SystemMatrix_alloc(type,fsystemMatrixPattern,row_blocksize,column_blocksize);
853        }
854      checkPasoError();      checkPasoError();
855      Paso_SystemMatrixPattern_dealloc(fsystemMatrixPattern);      Paso_SystemMatrixPattern_dealloc(fsystemMatrixPattern);
856      return SystemMatrixAdapter(fsystemMatrix,row_blocksize,row_functionspace,column_blocksize,column_functionspace);      return SystemMatrixAdapter(fsystemMatrix,row_blocksize,row_functionspace,column_blocksize,column_functionspace);
# Line 1071  int* MeshAdapter::borrowSampleReferenceI Line 1082  int* MeshAdapter::borrowSampleReferenceI
1082      break;      break;
1083    case(DegreesOfFreedom):    case(DegreesOfFreedom):
1084      out=mesh->Nodes->degreeOfFreedomId;      out=mesh->Nodes->degreeOfFreedomId;
 /*  
     for (i=0;i<mesh->Nodes->numNodes; ++i) {  
        if (mesh->Nodes->degreeOfFreedom[i]==sampleNo) {  
           out=mesh->Nodes->degreeOfFreedomIds[i];  
           break;  
        }  
     }  
 */  
1085      break;      break;
1086    case(ReducedDegreesOfFreedom):    case(ReducedDegreesOfFreedom):
1087      out=mesh->Nodes->reducedDegreeOfFreedomId;      out=mesh->Nodes->reducedDegreeOfFreedomId;
 /*  
     for (i=0;i<mesh->Nodes->numNodes; ++i) {  
        if (mesh->Nodes->reducedDegreeOfFreedom[i]==sampleNo) {  
           out=mesh->Nodes->Id[i];  
           break;  
        }  
     }  
 */  
1088      break;      break;
1089    default:    default:
1090      stringstream temp;      stringstream temp;

Legend:
Removed from v.964  
changed lines
  Added in v.969

  ViewVC Help
Powered by ViewVC 1.1.26