Changeset 24109
- Timestamp:
 - 07/29/19 05:27:37 (6 years ago)
 - File:
 - 
      
- 1 edited
 
 
Legend:
- Unmodified
 - Added
 - Removed
 
- 
      
issm/trunk-jpl/src/c/modules/SystemMatricesx/SystemMatricesx.cpp
r24108 r24109 120 120 //Kff->AllocationInfo(); 121 121 //Kfs->AllocationInfo(); 122 if(analysisenum==GLheightadvectionAnalysisEnum){123 PetscViewer viewer; PetscViewerASCIIOpen(PETSC_COMM_WORLD,"KFF.m",&viewer); PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_MATLAB);124 MatView(Kff->pmatrix->matrix,viewer);125 PetscViewerPopFormat(viewer);126 PetscViewerDestroy(&viewer);127 for(int e=0;e<femmodel->elements->Size();e++){128 Element* element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(e));129 int numdofs;130 for(int i=0;i<element->GetNumberOfNodes();i++){131 numdofs=element->nodes[i]->GetNumberOfDofs(NoneApproximationEnum,FsetEnum);132 int *doflist = xNew<int>(numdofs);133 element->nodes[i]->GetDofList(doflist,NoneApproximationEnum,FsetEnum);134 for(int j=0;j<numdofs;j++){135 if(doflist[j]==11796-1){136 printf("node: %i, element: %i, element->nodes[i]->Sid: %i \n", i, element->id, element->nodes[i]->Sid());137 }138 }139 xDelete<int>(doflist);140 }141 }142 }143 122 144 123 /*cleanu up and assign output pointers: */  
  Note:
 See   TracChangeset
 for help on using the changeset viewer.
  ![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)