Changeset 22647
- Timestamp:
- 03/28/18 14:18:01 (7 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
r22515 r22647 24 24 } 25 25 /*}}}*/ 26 Penta::Penta(int penta_id, int penta_sid, int index,IoModel* iomodel,int nummodels)/*{{{*/27 :ElementHook(nummodels, index+1,NUMVERTICES,iomodel){26 Penta::Penta(int penta_id, int penta_sid,IoModel* iomodel,int nummodels)/*{{{*/ 27 :ElementHook(nummodels,penta_id,NUMVERTICES,iomodel){ 28 28 29 29 int penta_elements_ids[2]; … … 38 38 39 39 /*Build neighbors list*/ 40 if (xIsNan<IssmDouble>(iomodel->Data("md.mesh.upperelements")[ index]) || iomodel->Data("md.mesh.upperelements")[index]==-1.) penta_elements_ids[1]=this->id; //upper penta is the same penta41 else penta_elements_ids[1]=reCast<int,IssmDouble>((iomodel->Data("md.mesh.upperelements")[ index]));42 if (xIsNan<IssmDouble>(iomodel->Data("md.mesh.lowerelements")[ index]) || iomodel->Data("md.mesh.lowerelements")[index]==-1.) penta_elements_ids[0]=this->id; //lower penta is the same penta43 else penta_elements_ids[0]=reCast<int,IssmDouble>((iomodel->Data("md.mesh.lowerelements")[ index]));40 if (xIsNan<IssmDouble>(iomodel->Data("md.mesh.upperelements")[penta_sid]) || iomodel->Data("md.mesh.upperelements")[penta_sid]==-1.) penta_elements_ids[1]=this->id; //upper penta is the same penta 41 else penta_elements_ids[1]=reCast<int,IssmDouble>((iomodel->Data("md.mesh.upperelements")[penta_sid])); 42 if (xIsNan<IssmDouble>(iomodel->Data("md.mesh.lowerelements")[penta_sid]) || iomodel->Data("md.mesh.lowerelements")[penta_sid]==-1.) penta_elements_ids[0]=this->id; //lower penta is the same penta 43 else penta_elements_ids[0]=reCast<int,IssmDouble>((iomodel->Data("md.mesh.lowerelements")[penta_sid])); 44 44 this->InitHookNeighbors(penta_elements_ids); 45 45 … … 2801 2801 /*Checks if debuging*/ 2802 2802 _assert_(iomodel->elements); 2803 _assert_(index==this->sid); 2803 2804 2804 2805 /*Recover element type*/ -
issm/trunk-jpl/src/c/classes/Elements/Penta.h
r22625 r22647 36 36 /*Penta constructors and destructor: {{{*/ 37 37 Penta(){}; 38 Penta(int penta_id,int penta_sid, int i,IoModel* iomodel,int nummodels);38 Penta(int penta_id,int penta_sid,IoModel* iomodel,int nummodels); 39 39 ~Penta(); 40 40 /*}}}*/ -
issm/trunk-jpl/src/c/classes/Elements/Seg.cpp
r21899 r22647 19 19 #define NUMVERTICES 2 20 20 /*Constructors/destructor/copy*/ 21 Seg::Seg(int seg_id, int seg_sid, int index,IoModel* iomodel,int nummodels)/*{{{*/22 :ElementHook(nummodels, index+1,NUMVERTICES,iomodel){21 Seg::Seg(int seg_id, int seg_sid,IoModel* iomodel,int nummodels)/*{{{*/ 22 :ElementHook(nummodels,seg_id,NUMVERTICES,iomodel){ 23 23 24 24 /*id: */ -
issm/trunk-jpl/src/c/classes/Elements/Seg.h
r22625 r22647 32 32 /*Seg constructors, destructors {{{*/ 33 33 Seg(){}; 34 Seg(int seg_id,int seg_sid, int i,IoModel* iomodel,int nummodels);34 Seg(int seg_id,int seg_sid,IoModel* iomodel,int nummodels); 35 35 ~Seg(); 36 36 /*}}}*/ -
issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp
r22515 r22647 20 20 21 21 /*Constructors/destructor/copy*/ 22 Tetra::Tetra(int seg_id, int seg_sid, int index,IoModel* iomodel,int nummodels)/*{{{*/23 :ElementHook(nummodels, index+1,NUMVERTICES,iomodel){22 Tetra::Tetra(int tet_id, int tet_sid,IoModel* iomodel,int nummodels)/*{{{*/ 23 :ElementHook(nummodels,tet_id,NUMVERTICES,iomodel){ 24 24 25 25 /*id: */ 26 this->id = seg_id;27 this->sid = seg_sid;26 this->id = tet_id; 27 this->sid = tet_sid; 28 28 29 29 //this->parameters: we still can't point to it, it may not even exist. Configure will handle this. … … 904 904 /*Checks if debuging*/ 905 905 _assert_(iomodel->elements); 906 _assert_(index==this->sid); 906 907 907 908 /*Recover element type*/ -
issm/trunk-jpl/src/c/classes/Elements/Tetra.h
r22625 r22647 32 32 /*Tetra constructors, destructors {{{*/ 33 33 Tetra(){}; 34 Tetra(int seg_id,int seg_sid,int i,IoModel* iomodel,int nummodels);34 Tetra(int tet_id,int tet_sid,IoModel* iomodel,int nummodels); 35 35 ~Tetra(); 36 36 /*}}}*/ -
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r22642 r22647 25 25 26 26 /*Constructors/destructor/copy*/ 27 Tria::Tria(int tria_id, int tria_sid, int index,IoModel* iomodel,int nummodels)/*{{{*/28 :ElementHook(nummodels, index+1,NUMVERTICES,iomodel){27 Tria::Tria(int tria_id, int tria_sid, IoModel* iomodel,int nummodels)/*{{{*/ 28 :ElementHook(nummodels,tria_id,NUMVERTICES,iomodel){ 29 29 30 30 /*id: */ … … 3604 3604 /*Checks if debuging*/ 3605 3605 _assert_(iomodel->elements); 3606 _assert_(index==this->sid); 3606 3607 3607 3608 /*Recover element type*/ -
issm/trunk-jpl/src/c/classes/Elements/Tria.h
r22625 r22647 34 34 /*Tria constructors, destructors {{{*/ 35 35 Tria(){}; 36 Tria(int tria_id,int tria_sid, int i,IoModel* iomodel,int nummodels);36 Tria(int tria_id,int tria_sid,IoModel* iomodel,int nummodels); 37 37 ~Tria(); 38 38 /*}}}*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp
r22515 r22647 37 37 case TriaEnum: 38 38 for(i=0;i<iomodel->numberofelements;i++){ 39 if(iomodel->my_elements[i]) elements->AddObject(new Tria(i+1,i,i ,iomodel,nummodels));39 if(iomodel->my_elements[i]) elements->AddObject(new Tria(i+1,i,iomodel,nummodels)); 40 40 } 41 41 break; 42 42 case TetraEnum: 43 43 for(i=0;i<iomodel->numberofelements;i++){ 44 if(iomodel->my_elements[i]) elements->AddObject(new Tetra(i+1,i,i ,iomodel,nummodels));44 if(iomodel->my_elements[i]) elements->AddObject(new Tetra(i+1,i,iomodel,nummodels)); 45 45 } 46 46 break; … … 48 48 iomodel->FetchData(2,"md.mesh.upperelements","md.mesh.lowerelements"); 49 49 for(i=0;i<iomodel->numberofelements;i++){ 50 if(iomodel->my_elements[i]) elements->AddObject(new Penta(i+1,i,i ,iomodel,nummodels));50 if(iomodel->my_elements[i]) elements->AddObject(new Penta(i+1,i,iomodel,nummodels)); 51 51 } 52 52 break;
Note:
See TracChangeset
for help on using the changeset viewer.