- Timestamp:
- 11/04/16 13:48:43 (8 years ago)
- Location:
- issm/trunk
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
-
issm/trunk/src
- Property svn:mergeinfo changed
-
issm/trunk/src/c
- Property svn:ignore
-
issm/trunk/src/c/modules/ModelProcessorx/ElementsAndVerticesPartitioning.cpp
r17806 r21341 44 44 45 45 /*Fetch parameters: */ 46 iomodel-> Constant(&numrifts,RiftsNumriftsEnum);46 iomodel->FindConstant(&numrifts,"md.rifts.numrifts"); 47 47 48 48 /*First, check that partitioning has not yet been carryed out. Just check whether my_elements pointers is not already assigned a value: */ … … 66 66 case PentaEnum: 67 67 elements_width=6; 68 iomodel->FetchData(&elements2d,NULL,NULL, MeshElements2dEnum);69 iomodel-> Constant(&numberofelements2d,MeshNumberofelements2dEnum);70 iomodel-> Constant(&numberofvertices2d,MeshNumberofvertices2dEnum);71 iomodel-> Constant(&numlayers,MeshNumberoflayersEnum);68 iomodel->FetchData(&elements2d,NULL,NULL,"md.mesh.elements2d"); 69 iomodel->FindConstant(&numberofelements2d,"md.mesh.numberofelements2d"); 70 iomodel->FindConstant(&numberofvertices2d,"md.mesh.numberofvertices2d"); 71 iomodel->FindConstant(&numlayers,"md.mesh.numberoflayers"); 72 72 break; 73 73 default: 74 _error_("mesh not supported yet");74 _error_("mesh elements "<< EnumToStringx(iomodel->meshelementtype) <<" not supported yet"); 75 75 } 76 76 … … 82 82 /*Deal with rifts, they have to be included into one partition only, not several: */ 83 83 if(numrifts){ 84 iomodel->FetchData(&riftinfo,&numrifts,NULL, RiftsRiftstructEnum);84 iomodel->FetchData(&riftinfo,&numrifts,NULL,"md.rifts.riftstruct"); 85 85 for(i=0;i<numrifts;i++){ 86 86 el1=reCast<int>(*(riftinfo+RIFTINFOSIZE*i+2))-1; //matlab indexing to c indexing … … 88 88 epart[el2]=epart[el1]; //ensures that this pair of elements will be in the same partition, as well as the corresponding vertices; 89 89 } 90 iomodel->DeleteData(riftinfo, RiftsRiftstructEnum);90 iomodel->DeleteData(riftinfo,"md.rifts.riftstruct"); 91 91 } 92 92 … … 114 114 * penpair has 2 nodes that are poointing toward 2 vertices. 115 115 * The 2 vertices must be in the same cpu as the penpair*/ 116 iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL, StressbalanceVertexPairingEnum);116 iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL,"md.stressbalance.vertex_pairing"); 117 117 for(i=0;i<numvertex_pairing;i++){ 118 118 if(my_vertices[vertex_pairing[2*i+0]-1] && !my_vertices[vertex_pairing[2*i+1]-1]){ … … 121 121 } 122 122 xDelete<int>(vertex_pairing); 123 iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL, MasstransportVertexPairingEnum);123 iomodel->FetchData(&vertex_pairing,&numvertex_pairing,NULL,"md.masstransport.vertex_pairing"); 124 124 for(i=0;i<numvertex_pairing;i++){ 125 125 if(my_vertices[vertex_pairing[2*i+0]-1] && !my_vertices[vertex_pairing[2*i+1]-1]){
Note:
See TracChangeset
for help on using the changeset viewer.