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.