Changeset 3567 for issm/trunk
- Timestamp:
- 04/19/10 09:49:52 (15 years ago)
- Location:
- issm/trunk/src
- Files:
-
- 178 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/DataSet/DataSet.cpp
r3556 r3567 188 188 memcpy(&enum_type,marshalled_dataset,sizeof(int)); marshalled_dataset+=sizeof(int); 189 189 190 if(enum_type==NodeEnum ()){190 if(enum_type==NodeEnum){ 191 191 Node* node=NULL; 192 192 node=new Node(); … … 194 194 dataset->AddObject(node); 195 195 } 196 else if(enum_type==VertexEnum ()){196 else if(enum_type==VertexEnum){ 197 197 Vertex* vertex=NULL; 198 198 vertex=new Vertex(); … … 200 200 dataset->AddObject(vertex); 201 201 } 202 else if(enum_type==ParamEnum ()){202 else if(enum_type==ParamEnum){ 203 203 Param* param=NULL; 204 204 param=new Param(); … … 206 206 dataset->AddObject(param); 207 207 } 208 else if(enum_type==NumparEnum ()){208 else if(enum_type==NumparEnum){ 209 209 Numpar* numpar=NULL; 210 210 numpar=new Numpar(); … … 212 212 dataset->AddObject(numpar); 213 213 } 214 else if(enum_type==TriaEnum ()){214 else if(enum_type==TriaEnum){ 215 215 Tria* tria=NULL; 216 216 tria=new Tria(); … … 218 218 dataset->AddObject(tria); 219 219 } 220 else if(enum_type==SingEnum ()){220 else if(enum_type==SingEnum){ 221 221 Sing* sing=NULL; 222 222 sing=new Sing(); … … 224 224 dataset->AddObject(sing); 225 225 } 226 else if(enum_type==BeamEnum ()){226 else if(enum_type==BeamEnum){ 227 227 Beam* beam=NULL; 228 228 beam=new Beam(); … … 230 230 dataset->AddObject(beam); 231 231 } 232 else if(enum_type==PentaEnum ()){232 else if(enum_type==PentaEnum){ 233 233 Penta* penta=NULL; 234 234 penta=new Penta(); … … 236 236 dataset->AddObject(penta); 237 237 } 238 else if(enum_type==MaticeEnum ()){238 else if(enum_type==MaticeEnum){ 239 239 Matice* matice=NULL; 240 240 matice=new Matice(); … … 242 242 dataset->AddObject(matice); 243 243 } 244 else if(enum_type==MatparEnum ()){244 else if(enum_type==MatparEnum){ 245 245 Matpar* matpar=NULL; 246 246 matpar=new Matpar(); … … 248 248 dataset->AddObject(matpar); 249 249 } 250 else if(enum_type==SpcEnum ()){250 else if(enum_type==SpcEnum){ 251 251 Spc* spc=NULL; 252 252 spc=new Spc(); … … 254 254 dataset->AddObject(spc); 255 255 } 256 else if(enum_type==PengridEnum ()){256 else if(enum_type==PengridEnum){ 257 257 Pengrid* pengrid=NULL; 258 258 pengrid=new Pengrid(); … … 260 260 dataset->AddObject(pengrid); 261 261 } 262 else if(enum_type==PenpairEnum ()){262 else if(enum_type==PenpairEnum){ 263 263 Penpair* penpair=NULL; 264 264 penpair=new Penpair(); … … 266 266 dataset->AddObject(penpair); 267 267 } 268 else if(enum_type==IcefrontEnum ()){268 else if(enum_type==IcefrontEnum){ 269 269 Icefront* icefront=NULL; 270 270 icefront=new Icefront(); … … 272 272 dataset->AddObject(icefront); 273 273 } 274 else if(enum_type==NumericalfluxEnum ()){274 else if(enum_type==NumericalfluxEnum){ 275 275 Numericalflux* numericalflux=NULL; 276 276 numericalflux=new Numericalflux(); … … 278 278 dataset->AddObject(numericalflux); 279 279 } 280 else if(enum_type==RgbEnum ()){280 else if(enum_type==RgbEnum){ 281 281 Rgb* rgb=NULL; 282 282 rgb=new Rgb(); … … 284 284 dataset->AddObject(rgb); 285 285 } 286 else if(enum_type==RiftfrontEnum ()){286 else if(enum_type==RiftfrontEnum){ 287 287 Riftfront* riftfront=NULL; 288 288 riftfront=new Riftfront(); … … 385 385 386 386 /*Find param type objects: */ 387 if((*object)->Enum()==ParamEnum ()){387 if((*object)->Enum()==ParamEnum){ 388 388 389 389 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 416 416 417 417 /*Find param type objects: */ 418 if((*object)->Enum()==ParamEnum ()){418 if((*object)->Enum()==ParamEnum){ 419 419 420 420 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 446 446 447 447 /*Find param type objects: */ 448 if((*object)->Enum()==ParamEnum ()){448 if((*object)->Enum()==ParamEnum){ 449 449 450 450 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 477 477 478 478 /*Find param type objects: */ 479 if((*object)->Enum()==ParamEnum ()){479 if((*object)->Enum()==ParamEnum){ 480 480 481 481 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 509 509 510 510 /*Find param type objects: */ 511 if((*object)->Enum()==ParamEnum ()){511 if((*object)->Enum()==ParamEnum){ 512 512 513 513 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 542 542 543 543 /*Find param type objects: */ 544 if((*object)->Enum()==ParamEnum ()){544 if((*object)->Enum()==ParamEnum){ 545 545 546 546 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 573 573 574 574 /*Find param type objects: */ 575 if((*object)->Enum()==ParamEnum ()){575 if((*object)->Enum()==ParamEnum){ 576 576 577 577 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 602 602 603 603 /*Find param type objects: */ 604 if((*object)->Enum()==ParamEnum ()){604 if((*object)->Enum()==ParamEnum){ 605 605 606 606 /*Ok, this object is a parameter, recover it and ask which name it has: */ … … 630 630 631 631 /*Find param type objects: */ 632 if((*object)->Enum()==ResultEnum ()){632 if((*object)->Enum()==ResultEnum){ 633 633 634 634 /*Ok, this object is a result,recover it and ask which name it has: */ … … 660 660 661 661 /*Find param type objects: */ 662 if((*object)->Enum()==ResultEnum ()){662 if((*object)->Enum()==ResultEnum){ 663 663 664 664 /*Ok, this object is a result, recover it and ask which name it has: */ … … 678 678 } 679 679 /*}}}*/ 680 /*FUNCTION DataSet::GetEnum (){{{1*/680 /*FUNCTION DataSet::GetEnum{{{1*/ 681 681 int DataSet::GetEnum(){ 682 682 return enum_type; … … 861 861 } 862 862 863 if((*object)->Enum()==NodeEnum ()){863 if((*object)->Enum()==NodeEnum){ 864 864 node=(Node*)(*object); 865 865 node->Configure(nodes,vertices); 866 866 } 867 if((*object)->Enum()==NumparEnum ()){867 if((*object)->Enum()==NumparEnum){ 868 868 numpar=(Numpar*)(*object); 869 869 numpar->Configure(parameters); … … 1080 1080 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1081 1081 1082 if((*object)->Enum()==NodeEnum ()){1082 if((*object)->Enum()==NodeEnum){ 1083 1083 node=(Node*)(*object); 1084 1084 node->FieldAverageOntoVertices(fieldsum,connectivity,field); … … 1097 1097 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1098 1098 1099 if((*object)->Enum()==NodeEnum ()){1099 if((*object)->Enum()==NodeEnum){ 1100 1100 1101 1101 node=(Node*)(*object); … … 1115 1115 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1116 1116 1117 if((*object)->Enum()==PentaEnum ()){1117 if((*object)->Enum()==PentaEnum){ 1118 1118 1119 1119 penta=(Penta*)(*object); … … 1121 1121 1122 1122 } 1123 if((*object)->Enum()==NodeEnum ()){1123 if((*object)->Enum()==NodeEnum){ 1124 1124 1125 1125 node=(Node*)(*object); … … 1193 1193 1194 1194 /*Check this is a single point constraint (spc): */ 1195 if((*object)->Enum()==NodeEnum ()){1195 if((*object)->Enum()==NodeEnum){ 1196 1196 1197 1197 node=(Node*)(*object); … … 1250 1250 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1251 1251 1252 if((*object)->Enum()==PengridEnum ()){1252 if((*object)->Enum()==PengridEnum){ 1253 1253 found=1; 1254 1254 break; … … 1282 1282 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1283 1283 1284 if((*object)->Enum()==PengridEnum ()){1284 if((*object)->Enum()==PengridEnum){ 1285 1285 1286 1286 pengrid=(Pengrid*)(*object); … … 1343 1343 1344 1344 /*Check this is a node: */ 1345 if((*object)->Enum()==NodeEnum ()){1345 if((*object)->Enum()==NodeEnum){ 1346 1346 1347 1347 node=(Node*)(*object); … … 1372 1372 1373 1373 /*Check this is a single point constraint (spc): */ 1374 if((*object)->Enum()==RgbEnum ()){ //we assume uniqueness of all Rgbs, no error checking here.1374 if((*object)->Enum()==RgbEnum){ //we assume uniqueness of all Rgbs, no error checking here. 1375 1375 1376 1376 count++; … … 1390 1390 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1391 1391 1392 if((*object)->Enum()==RiftfrontEnum ()){1392 if((*object)->Enum()==RiftfrontEnum){ 1393 1393 1394 1394 riftfront=(Riftfront*)(*object); … … 1448 1448 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1449 1449 1450 if((*object)->Enum()==PenpairEnum ()){1450 if((*object)->Enum()==PenpairEnum){ 1451 1451 1452 1452 penpair=(Penpair*)(*object); … … 1484 1484 1485 1485 /*Check this is a single point constraint (spc): */ 1486 if((*object)->Enum()==RgbEnum ()){ //we assume uniqueness of all Rgbs, no error checking here.1486 if((*object)->Enum()==RgbEnum){ //we assume uniqueness of all Rgbs, no error checking here. 1487 1487 1488 1488 /*we found an rgb, increment counter, so that row index for Rmg is up to date: */ … … 1545 1545 1546 1546 /*Check this is a single point constraint (spc): */ 1547 if((*object)->Enum()==SpcEnum ()){1547 if((*object)->Enum()==SpcEnum){ 1548 1548 1549 1549 spc=(Spc*)(*object); … … 1619 1619 material->UpdateFromInputs(inputs); 1620 1620 } 1621 if((*object)->Enum()==NodeEnum ()){1621 if((*object)->Enum()==NodeEnum){ 1622 1622 node=(Node*)(*object); 1623 1623 node->UpdateFromInputs(inputs); … … 1635 1635 for ( object=objects.begin() ; object < objects.end(); object++ ){ 1636 1636 1637 if((*object)->Enum()==VertexEnum ()){1637 if((*object)->Enum()==VertexEnum){ 1638 1638 1639 1639 vertex=(Vertex*)(*object); -
issm/trunk/src/c/EnumDefinitions/EnumDefinitions.cpp
r3469 r3567 6 6 #include "./EnumDefinitions.h" 7 7 8 /*Datasets: */9 int DatasetsEnum(void){ return 100; }10 int ElementsEnum(void){ return 101; }11 int NodesEnum(void){ return 102; }12 int ConstraintsEnum(void){ return 103; }13 int LoadsEnum(void){ return 104; }14 int MaterialsEnum(void){ return 105; }15 int ParametersEnum(void){ return 106; }16 int ResultsEnum(void){ return 107; }17 int VerticesEnum(void){ return 108; }18 19 /*ANALYSIS TYPES: */20 int AnalysisEnum(void){ return 200; }21 //diagnostic22 int DiagnosticAnalysisEnum(void){ return 210; }23 int HorizAnalysisEnum(void){ return 211; }24 int StokesAnalysisEnum(void){ return 212; }25 int HutterAnalysisEnum(void){ return 213; }26 int VertAnalysisEnum(void){ return 214; }27 //control28 int ControlAnalysisEnum(void){ return 220; }29 int AdjointAnalysisEnum(void){ return 221; }30 int InverseAnalysisEnum(void){ return 222; }31 int GradientAnalysisEnum(void){ return 223; }32 //thermal33 int ThermalAnalysisEnum(void){ return 230; }34 //transient35 int TransientAnalysisEnum(void){ return 231; }36 int SteadyAnalysisEnum(void){ return 232; } //NEEDED by nightly runs37 //slope38 int SlopecomputeAnalysisEnum(void){ return 240; }39 int SurfaceXAnalysisEnum(void){ return 241; }40 int SurfaceYAnalysisEnum(void){ return 242; }41 int BedXAnalysisEnum(void){ return 243; }42 int BedYAnalysisEnum(void){ return 244; }43 //prognostic44 int PrognosticAnalysisEnum(void){ return 250; }45 int Prognostic2AnalysisEnum(void){ return 251; }46 int BalancedthicknessAnalysisEnum(void){ return 252; }47 int BalancedvelocitiesAnalysisEnum(void){ return 253; }48 //melting49 int MeltingAnalysisEnum(void){ return 260; }50 //mesh2grid51 int Mesh2gridAnalysisEnum(void){ return 270; }52 //parameters53 int ParametersAnalysisEnum(void){ return 280; }54 //steadystate55 int SteadystateAnalysisEnum(void){ return 281; }56 //none57 int NoneAnalysisEnum(void){ return 290; }58 59 60 /*Formulations: */61 int FormulationEnum(void){ return 300; }62 int NoneFormulationEnum(void){ return 301; }63 int HutterFormulationEnum(void){ return 302; }64 int MacAyealFormulationEnum(void){ return 303; }65 int PattynFormulationEnum(void){ return 304; }66 int StokesFormulationEnum(void){ return 305; }67 68 /*OBJECTS: */69 int ObjectEnum(void){ return 400; }70 /*Elements: */71 int HookEnum(void){ return 409; }72 int ElementEnum(void){ return 410; }73 int TriaEnum(void){ return 411; }74 int ElementPropertiesEnum(void){ return 412; }75 int NodePropertiesEnum(void){ return 413; }76 int PentaEnum(void){ return 414; }77 int SingEnum(void){ return 415; }78 int BeamEnum(void){ return 416; }79 int DofIndexingEnum(void){ return 417; }80 81 /*Grids: */82 int VertexEnum(void){ return 420; }83 int NodeEnum(void){ return 421; }84 /*Loads: */85 int LoadEnum(void){ return 430; }86 int IcefrontEnum(void){ return 431; }87 int RiftfrontEnum(void){ return 432; }88 int PenpairEnum(void){ return 433; }89 int PengridEnum(void){ return 434; }90 int NumericalfluxEnum(void){ return 435; }91 /*Materials: */92 int MaterialEnum(void){ return 440; }93 int MaticeEnum(void){ return 441; }94 int MatparEnum(void){ return 442; }95 int NumparEnum(void){ return 443; }96 /*Inputs: */97 int InputEnum(void){ return 450; }98 /*Params: */99 int ParamEnum(void){ return 460; }100 /*Results: */101 int ResultEnum(void){ return 470; }102 /*Rgb: */103 int RgbEnum(void){ return 480; }104 /*Spc: */105 int SpcEnum(void){ return 490; }106 /*DofVec: */107 int DofVecEnum(void){ return 495; }108 109 110 /*GEOGRAPHY:*/111 int GeographyEnum(void){ return 500; }112 int IceSheetEnum(void){ return 502; }113 int IceShelfEnum(void){ return 502; }114 115 /*FILL:*/116 int WaterEnum(void){ return 601; }117 int IceEnum(void){ return 602; }118 int AirEnum(void){ return 603; }119 int MelangeEnum(void){ return 604; }120 121 8 /*functions on enums: */ 122 9 int EnumIsElement(int en){ 123 if ((en>= 410) && (en<=419))return 1;10 if ((en>=ElementEnum) && (en<=SingEnum))return 1; 124 11 else return 0; 125 12 } 126 13 127 14 int EnumIsLoad(int en){ 128 if ((en>= 430) && (en<=439))return 1;15 if ((en>=LoadEnum) && (en<=NumericalfluxEnum))return 1; 129 16 else return 0; 130 17 } 131 18 int EnumIsMaterial(int en){ 132 if ((en>= 440) && (en<=449))return 1;19 if ((en>=MaterialEnum) && (en<=NumparEnum))return 1; 133 20 else return 0; 134 21 } -
issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h
r3417 r3567 7 7 #define _ENUM_DEFINITIONS_ 8 8 9 /*Datasets: */ 10 int DatasetsEnum(void); 11 int ElementsEnum(void); 12 int NodesEnum(void); 13 int ConstraintsEnum(void); 14 int LoadsEnum(void); 15 int MaterialsEnum(void); 16 int ParametersEnum(void); 17 int ResultsEnum(void); 18 int VerticesEnum(void); 9 enum definitions{ 19 10 20 /*ANALYSIS TYPES: */ 21 int AnalysisEnum(void); 22 //diagnostic 23 int DiagnosticAnalysisEnum(void); 24 int HorizAnalysisEnum(void); 25 int StokesAnalysisEnum(void); 26 int HutterAnalysisEnum(void); 27 int VertAnalysisEnum(void); 28 //control 29 int ControlAnalysisEnum(void); 30 int AdjointAnalysisEnum(void); 31 int InverseAnalysisEnum(void); 32 int GradientAnalysisEnum(void); 33 //thermal 34 int ThermalAnalysisEnum(void); 35 //transient 36 int TransientAnalysisEnum(void); 37 int SteadyAnalysisEnum(void); 38 //slope 39 int SlopecomputeAnalysisEnum(void); 40 int SurfaceXAnalysisEnum(void); 41 int SurfaceYAnalysisEnum(void); 42 int BedXAnalysisEnum(void); 43 int BedYAnalysisEnum(void); 44 //prognostic 45 int PrognosticAnalysisEnum(void); 46 int Prognostic2AnalysisEnum(void); 47 int BalancedthicknessAnalysisEnum(void); 48 int BalancedvelocitiesAnalysisEnum(void); 49 //melting 50 int MeltingAnalysisEnum(void); 51 //mesh2grid 52 int Mesh2gridAnalysisEnum(void); 53 //parameters 54 int ParametersAnalysisEnum(void); 55 //steadystate 56 int SteadystateAnalysisEnum(void); 57 //none 58 int NoneAnalysisEnum(void); 11 /*Datasets {{{1*/ 12 DatasetsEnum, 13 ElementsEnum, 14 NodesEnum, 15 ConstraintsEnum, 16 LoadsEnum, 17 MaterialsEnum, 18 ParametersEnum, 19 ResultsEnum, 20 VerticesEnum, 21 /*}}}*/ 22 /*Analysis types {{{1 */ 23 AnalysisEnum, 24 //diagnostic 25 DiagnosticAnalysisEnum, 26 HorizAnalysisEnum, 27 StokesAnalysisEnum, 28 HutterAnalysisEnum, 29 VertAnalysisEnum, 30 //control 31 ControlAnalysisEnum, 32 AdjointAnalysisEnum, 33 InverseAnalysisEnum, 34 GradientAnalysisEnum, 35 //thermal 36 ThermalAnalysisEnum, 37 //transient 38 TransientAnalysisEnum, 39 SteadyAnalysisEnum, 40 //slope 41 SlopecomputeAnalysisEnum, 42 SurfaceXAnalysisEnum, 43 SurfaceYAnalysisEnum, 44 BedXAnalysisEnum, 45 BedYAnalysisEnum, 46 //prognostic 47 PrognosticAnalysisEnum, 48 Prognostic2AnalysisEnum, 49 BalancedthicknessAnalysisEnum, 50 BalancedvelocitiesAnalysisEnum, 51 //melting 52 MeltingAnalysisEnum, 53 //mesh2grid 54 Mesh2gridAnalysisEnum, 55 //parameters 56 ParametersAnalysisEnum, 57 //steadystate 58 SteadystateAnalysisEnum, 59 //none 60 NoneAnalysisEnum, 61 /*}}}*/ 62 /*Formulations {{{1*/ 63 FormulationEnum, 64 NoneFormulationEnum, 65 HutterFormulationEnum, 66 MacAyealFormulationEnum, 67 PattynFormulationEnum, 68 StokesFormulationEnum, 69 /*}}}*/ 70 /*Objects {{{1*/ 71 ObjectEnum, 72 HookEnum, 73 /*Elements: */ 74 ElementEnum, 75 TriaEnum, 76 ElementPropertiesEnum, 77 NodePropertiesEnum, 78 PentaEnum, 79 BeamEnum, 80 SingEnum, 81 DofIndexingEnum, 82 /*Grids: */ 83 NodeEnum, 84 VertexEnum, 85 /*Loads: */ 86 LoadEnum, 87 IcefrontEnum, 88 RiftfrontEnum, 89 PenpairEnum, 90 PengridEnum, 91 NumericalfluxEnum, 92 /*Materials: */ 93 MaterialEnum, 94 MaticeEnum, 95 MatparEnum, 96 NumparEnum, 97 /*Inputs: */ 98 InputEnum, 99 /*Params: */ 100 ParamEnum, 101 /*Results: */ 102 ResultEnum, 103 /*Rgb: */ 104 RgbEnum, 105 /*Spc: */ 106 SpcEnum, 107 /*DofVec: */ 108 DofVecEnum, 109 /*}}}*/ 110 /*Geography {{{1*/ 111 GeographyEnum, 112 IceSheetEnum, 113 IceShelfEnum, 114 /*}}}*/ 115 /*Fill {{{1*/ 116 WaterEnum, 117 IceEnum, 118 AirEnum, 119 MelangeEnum, 120 /*}}}*/ 59 121 60 61 /*Formulations: */ 62 int FormulationEnum(void); 63 int NoneFormulationEnum(void); 64 int HutterFormulationEnum(void); 65 int MacAyealFormulationEnum(void); 66 int PattynFormulationEnum(void); 67 int StokesFormulationEnum(void); 68 69 /*OBJECTS: */ 70 int ObjectEnum(void); 71 /*Elements: */ 72 int HookEnum(void); 73 int ElementEnum(void); 74 int TriaEnum(void); 75 int ElementPropertiesEnum(void); 76 int NodePropertiesEnum(void); 77 int PentaEnum(void); 78 int SingEnum(void); 79 int BeamEnum(void); 80 int DofIndexingEnum(void); 81 /*Grids: */ 82 int NodeEnum(void); 83 int VertexEnum(void); 84 /*Loads: */ 85 int LoadEnum(void); 86 int IcefrontEnum(void); 87 int RiftfrontEnum(void); 88 int PenpairEnum(void); 89 int PengridEnum(void); 90 int NumericalfluxEnum(void); 91 /*Materials: */ 92 int MaterialEnum(void); 93 int MaticeEnum(void); 94 int MatparEnum(void); 95 int NumparEnum(void); 96 /*Inputs: */ 97 int InputEnum(void); 98 /*Params: */ 99 int ParamEnum(void); 100 /*Results: */ 101 int ResultEnum(void); 102 /*Rgb: */ 103 int RgbEnum(void); 104 /*Spc: */ 105 int SpcEnum(void); 106 /*DofVec: */ 107 int DofVecEnum(void); 108 109 /*GEOGRAPHY:*/ 110 int GeographyEnum(void); 111 int IceSheetEnum(void); 112 int IceShelfEnum(void); 113 114 /*FILL: */ 115 int WaterEnum(void); 116 int IceEnum(void); 117 int AirEnum(void); 118 int MelangeEnum(void); 122 }; 119 123 120 124 /*Functions on enums: */ … … 122 126 int EnumIsLoad(int en); 123 127 int EnumIsMaterial(int en); 128 124 129 #endif -
issm/trunk/src/c/EnumDefinitions/SynchronizeMatlabEnum.sh
r1901 r3567 7 7 rm $ISSM_DIR/src/m/enum/*.m 8 8 9 #Get the list of interesting lines of EnumDefinitions.cpp and introduce line number10 9 echo "Synchronizong with Matlab's Enums..." 11 cat EnumDefinitions.cpp | grep "return" | grep "int" | awk '{ printf "line%so %s\n", NR, $0 }' | sed -e "s/(void){//" > temp 10 #Get all lines of EnumDefinitions2.h which hold Enum, | remove all comas | add line number in the first column > put everything in file temp 11 cat EnumDefinitions.h | grep "Enum," |sed -e "s/,//g" | awk '{ printf "%s %s\n", NR, $0 }' > temp 12 13 #Build File AnalysisTypeFromEnum.m 12 14 cat <<END > AnalysisTypeFromEnum.m 13 15 function string=AnalysisTypeFromEnum(enum) … … 23 25 %initialize output 24 26 string='not found'; 25 26 27 END 27 28 28 29 #get number of lines in temp 29 NUMBEROFLINES= `wc -l temp | awk '{printf("%s",$1);}'`;30 NUMBEROFLINES=$(wc -l temp | awk '{printf("%s",$1);}'); 30 31 31 32 # go through the lines of temp … … 33 34 34 35 #Get name and enum of the line i 35 NAME=`cat temp | awk '/line'$i'o/ {printf("%s\n",$3);}';`; 36 ENUM=`cat temp | awk '/line'$i'o/ {printf("%s\n",$5);}';`; 37 echo writing file $NAME \($i/$NUMBEROFLINES\) 36 NAME=$(cat temp | grep "^[ ]*$i " | awk '{printf("%s",$2);}'); 37 ENUM=$i; 38 #offset Enum by one (Enum begins with 0 and not 1!) 39 let ENUM=$ENUM-1 38 40 39 41 #write corresponding Matlab file 40 cat <<END > $(echo $NAME.m) 42 echo "writing file "$NAME".m ($i/$NUMBEROFLINES)" 43 cat <<END > $(echo $NAME".m") 41 44 function macro=$(echo $NAME)() 42 45 %$(echo `echo $NAME | sed -e "s/[a-z]/\U&/g"`) - Enum of $(echo `echo $NAME | sed -e "s/Enum//"`) … … 47 50 % macro=$NAME() 48 51 49 macro=$ENUM 52 macro=$ENUM; 50 53 END 51 54 … … 62 65 63 66 #move file to m/enum 64 mv $NAME .m $ISSM_DIR/src/m/enum/67 mv $NAME".m" $ISSM_DIR"/src/m/enum/" 65 68 66 69 done -
issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp
r3452 r3567 17 17 18 18 /*Create constraints: */ 19 constraints = new DataSet(ConstraintsEnum ());19 constraints = new DataSet(ConstraintsEnum); 20 20 21 21 /*Fetch data: */ -
issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateElementsNodesAndMaterialsBalancedthickness.cpp
r3476 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Partition elements and vertices and nodes: */ -
issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateLoadsBalancedthickness.cpp
r3332 r3567 16 16 17 17 /*Create loads: */ 18 loads = new DataSet(LoadsEnum ());18 loads = new DataSet(LoadsEnum); 19 19 20 20 -
issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateConstraintsBalancedvelocities.cpp
r3456 r3567 19 19 20 20 /*Create constraints: */ 21 constraints = new DataSet(ConstraintsEnum ());21 constraints = new DataSet(ConstraintsEnum); 22 22 23 23 /*Fetch data: */ -
issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateElementsNodesAndMaterialsBalancedvelocities.cpp
r3434 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Partition elements and vertices and nodes: */ -
issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateLoadsBalancedvelocities.cpp
r3332 r3567 16 16 17 17 /*Create loads: */ 18 loads = new DataSet(LoadsEnum ());18 loads = new DataSet(LoadsEnum); 19 19 20 20 -
issm/trunk/src/c/ModelProcessorx/CreateDataSets.cpp
r3417 r3567 22 22 23 23 /*This is just a high level driver: */ 24 if (iomodel->analysis_type==DiagnosticAnalysisEnum ()){24 if (iomodel->analysis_type==DiagnosticAnalysisEnum){ 25 25 26 if (iomodel->sub_analysis_type==HorizAnalysisEnum ()){26 if (iomodel->sub_analysis_type==HorizAnalysisEnum){ 27 27 28 28 CreateElementsNodesAndMaterialsDiagnosticHoriz(pelements,pnodes, pvertices, pmaterials, iomodel,iomodel_handle); … … 32 32 33 33 } 34 else if (iomodel->sub_analysis_type==VertAnalysisEnum ()){34 else if (iomodel->sub_analysis_type==VertAnalysisEnum){ 35 35 36 36 CreateElementsNodesAndMaterialsDiagnosticVert(pelements,pnodes, pvertices, pmaterials, iomodel,iomodel_handle); … … 39 39 40 40 } 41 else if (iomodel->sub_analysis_type==StokesAnalysisEnum ()){41 else if (iomodel->sub_analysis_type==StokesAnalysisEnum){ 42 42 43 43 CreateElementsNodesAndMaterialsDiagnosticStokes(pelements,pnodes, pvertices, pmaterials, iomodel,iomodel_handle); … … 46 46 47 47 } 48 else if (iomodel->sub_analysis_type==HutterAnalysisEnum ()){48 else if (iomodel->sub_analysis_type==HutterAnalysisEnum){ 49 49 50 50 CreateElementsNodesAndMaterialsDiagnosticHutter(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle); … … 54 54 } 55 55 } 56 else if (iomodel->analysis_type==SlopecomputeAnalysisEnum ()){56 else if (iomodel->analysis_type==SlopecomputeAnalysisEnum){ 57 57 58 58 CreateElementsNodesAndMaterialsSlopeCompute(pelements,pnodes, pvertices,pmaterials, iomodel,iomodel_handle); … … 61 61 62 62 } 63 else if (iomodel->analysis_type==ThermalAnalysisEnum ()){63 else if (iomodel->analysis_type==ThermalAnalysisEnum){ 64 64 65 65 CreateElementsNodesAndMaterialsThermal(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle); … … 69 69 70 70 } 71 else if (iomodel->analysis_type==MeltingAnalysisEnum ()){71 else if (iomodel->analysis_type==MeltingAnalysisEnum){ 72 72 73 73 CreateElementsNodesAndMaterialsMelting(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle); … … 76 76 CreateParametersMelting(pparameters,iomodel,iomodel_handle); 77 77 } 78 else if (iomodel->analysis_type==PrognosticAnalysisEnum ()){78 else if (iomodel->analysis_type==PrognosticAnalysisEnum){ 79 79 80 80 CreateElementsNodesAndMaterialsPrognostic(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle); … … 84 84 85 85 } 86 else if (iomodel->analysis_type==Prognostic2AnalysisEnum ()){86 else if (iomodel->analysis_type==Prognostic2AnalysisEnum){ 87 87 88 88 CreateElementsNodesAndMaterialsPrognostic2(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle); … … 92 92 93 93 } 94 else if (iomodel->analysis_type==BalancedthicknessAnalysisEnum ()){94 else if (iomodel->analysis_type==BalancedthicknessAnalysisEnum){ 95 95 96 96 CreateElementsNodesAndMaterialsBalancedthickness(pelements,pnodes,pvertices, pmaterials, iomodel,iomodel_handle); … … 100 100 101 101 } 102 else if (iomodel->analysis_type==BalancedvelocitiesAnalysisEnum ()){102 else if (iomodel->analysis_type==BalancedvelocitiesAnalysisEnum){ 103 103 104 104 CreateElementsNodesAndMaterialsBalancedvelocities(pelements,pnodes,pvertices,pmaterials, iomodel,iomodel_handle); -
issm/trunk/src/c/ModelProcessorx/CreateParameters.cpp
r3538 r3567 30 30 31 31 /*Initialize dataset: */ 32 parameters = new DataSet(ParametersEnum ());32 parameters = new DataSet(ParametersEnum); 33 33 34 34 //solution parameters: always 1'st, to speed-up its lookup by elements.. … … 245 245 count++; 246 246 param= new Param(count,"numberofnodes",DOUBLE); 247 if (iomodel->analysis_type==Prognostic2AnalysisEnum ()) param->SetDouble(3*iomodel->numberofelements);247 if (iomodel->analysis_type==Prognostic2AnalysisEnum) param->SetDouble(3*iomodel->numberofelements); 248 248 else param->SetDouble(iomodel->numberofvertices); 249 249 parameters->AddObject(param); -
issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp
r3507 r3567 23 23 24 24 /*Create constraints: */ 25 constraints = new DataSet(ConstraintsEnum ());25 constraints = new DataSet(ConstraintsEnum); 26 26 27 27 /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp
r3476 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */ … … 58 58 if(iomodel->my_elements[i]){ 59 59 60 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum ()){ //elements of type 1 are Hutter type Tria. Don't create this elements.60 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum){ //elements of type 1 are Hutter type Tria. Don't create this elements. 61 61 62 62 /*Create and add tria element to elements dataset: */ … … 106 106 for (i=0;i<iomodel->numberofelements;i++){ 107 107 if(iomodel->my_elements[i]){ 108 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum () | *(iomodel->elements_type+2*i+0)==PattynFormulationEnum()){ //elements of type 1 are Hutter type Tria. Don't create this elements.108 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum | *(iomodel->elements_type+2*i+0)==PattynFormulationEnum){ //elements of type 1 are Hutter type Tria. Don't create this elements. 109 109 /*Create and add penta element to elements dataset: */ 110 110 elements->AddObject(new Penta(i,iomodel)); -
issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateLoadsDiagnosticHoriz.cpp
r3488 r3567 21 21 22 22 /*Create loads: */ 23 loads = new DataSet(LoadsEnum ());23 loads = new DataSet(LoadsEnum); 24 24 25 25 /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */ … … 45 45 46 46 /*Do not create ice front if Hutter or Stokes elements*/ 47 if ((int)*(iomodel->elements_type+2*element+0)==(HutterFormulationEnum () || StokesFormulationEnum())) continue;47 if ((int)*(iomodel->elements_type+2*element+0)==(HutterFormulationEnum || StokesFormulationEnum)) continue; 48 48 49 49 /*Create and add load: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp
r3477 r3567 19 19 20 20 /*Create constraints: */ 21 constraints = new DataSet(ConstraintsEnum ());21 constraints = new DataSet(ConstraintsEnum); 22 22 23 23 /*Now, is the flag ishutter on? otherwise, do nothing: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp
r3508 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Now, is the flag ishutter on? otherwise, do nothing: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateLoadsDiagnosticHutter.cpp
r3456 r3567 15 15 16 16 /*Create loads: */ 17 loads = new DataSet(LoadsEnum ());17 loads = new DataSet(LoadsEnum); 18 18 19 19 /*Now, is the flag ishutter on? otherwise, do nothing: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateConstraintsDiagnosticStokes.cpp
r3456 r3567 18 18 19 19 /*Create constraints: */ 20 constraints = new DataSet(ConstraintsEnum ());20 constraints = new DataSet(ConstraintsEnum); 21 21 22 22 /*Now, is the flag ishutter on? otherwise, do nothing: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateElementsNodesAndMaterialsDiagnosticStokes.cpp
r3550 r3567 25 25 26 26 /*First create the elements, nodes and material properties: */ 27 elements = new DataSet(ElementsEnum ());28 nodes = new DataSet(NodesEnum ());29 vertices = new DataSet(VerticesEnum ());30 materials = new DataSet(MaterialsEnum ());27 elements = new DataSet(ElementsEnum); 28 nodes = new DataSet(NodesEnum); 29 vertices = new DataSet(VerticesEnum); 30 materials = new DataSet(MaterialsEnum); 31 31 32 32 /*Now, do we have Stokes elements?*/ … … 59 59 60 60 if(iomodel->my_elements[i]){ 61 if (*(iomodel->elements_type+2*i+1)==StokesFormulationEnum ()){61 if (*(iomodel->elements_type+2*i+1)==StokesFormulationEnum){ 62 62 63 63 /*Create and add penta element to elements dataset: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp
r3488 r3567 25 25 26 26 /*Create loads: */ 27 loads = new DataSet(LoadsEnum ());27 loads = new DataSet(LoadsEnum); 28 28 29 29 /*Now, is the flag ishutter on? otherwise, do nothing: */ … … 49 49 50 50 /*Do not create ice front if it is not a Stokes element*/ 51 if (!(int)*(iomodel->elements_type+2*element+1)==StokesFormulationEnum ()){51 if (!(int)*(iomodel->elements_type+2*element+1)==StokesFormulationEnum){ 52 52 continue; 53 53 } -
issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateConstraintsDiagnosticVert.cpp
r3444 r3567 17 17 18 18 /*Create constraints: */ 19 constraints = new DataSet(ConstraintsEnum ());19 constraints = new DataSet(ConstraintsEnum); 20 20 21 21 /*return if 2d mesh*/ -
issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateElementsNodesAndMaterialsDiagnosticVert.cpp
r3476 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */ -
issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateLoadsDiagnosticVert.cpp
r3443 r3567 15 15 16 16 /*Create loads: */ 17 loads = new DataSet(LoadsEnum ());17 loads = new DataSet(LoadsEnum); 18 18 19 19 /*Assign output pointer: */ -
issm/trunk/src/c/ModelProcessorx/Melting/CreateConstraintsMelting.cpp
r3456 r3567 18 18 19 19 /*Create constraints: */ 20 constraints = new DataSet(ConstraintsEnum ());20 constraints = new DataSet(ConstraintsEnum); 21 21 22 22 /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these -
issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp
r3494 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Partition elements and vertices and nodes: */ -
issm/trunk/src/c/ModelProcessorx/Melting/CreateLoadsMelting.cpp
r3456 r3567 21 21 22 22 /*Create loads: */ 23 loads = new DataSet(LoadsEnum ());23 loads = new DataSet(LoadsEnum); 24 24 25 25 //create penalties for grids: no grid can have a temperature over the melting point -
issm/trunk/src/c/ModelProcessorx/Partitioning.cpp
r3559 r3567 22 22 void Partitioning(bool** pmy_elements, bool** pmy_vertices, bool** pmy_nodes, bool** pmy_bordervertices, IoModel* iomodel, ConstDataHandle iomodel_handle){ 23 23 24 if (iomodel->analysis_type==Prognostic2AnalysisEnum ())24 if (iomodel->analysis_type==Prognostic2AnalysisEnum) 25 25 DiscontinuousGalerkinPartitioning(pmy_elements, pmy_vertices, pmy_nodes, pmy_bordervertices, iomodel, iomodel_handle); 26 26 else -
issm/trunk/src/c/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp
r3456 r3567 19 19 20 20 /*Create constraints: */ 21 constraints = new DataSet(ConstraintsEnum ());21 constraints = new DataSet(ConstraintsEnum); 22 22 23 23 /*Fetch data: */ -
issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp
r3434 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Partition elements and vertices and nodes: */ -
issm/trunk/src/c/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp
r3456 r3567 15 15 16 16 /*Create loads: */ 17 loads = new DataSet(LoadsEnum ());17 loads = new DataSet(LoadsEnum); 18 18 19 19 /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these -
issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateConstraintsPrognostic2.cpp
r3359 r3567 15 15 16 16 /*Create constraints: */ 17 constraints = new DataSet(ConstraintsEnum ());17 constraints = new DataSet(ConstraintsEnum); 18 18 19 19 /*Assign output pointer: */ -
issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateElementsNodesAndMaterialsPrognostic2.cpp
r3534 r3567 27 27 28 28 /*First create the elements, nodes and material properties: */ 29 elements = new DataSet(ElementsEnum ());30 nodes = new DataSet(NodesEnum ());31 vertices = new DataSet(VerticesEnum ());32 materials = new DataSet(MaterialsEnum ());29 elements = new DataSet(ElementsEnum); 30 nodes = new DataSet(NodesEnum); 31 vertices = new DataSet(VerticesEnum); 32 materials = new DataSet(MaterialsEnum); 33 33 34 34 /*Partition elements and vertices and nodes: */ -
issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateLoadsPrognostic2.cpp
r3549 r3567 29 29 30 30 /*Create loads: */ 31 loads = new DataSet(LoadsEnum ());31 loads = new DataSet(LoadsEnum); 32 32 33 33 /*Get edges and elements*/ -
issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateConstraintsSlopeCompute.cpp
r3456 r3567 15 15 16 16 /*Create constraints: */ 17 constraints = new DataSet(ConstraintsEnum ());17 constraints = new DataSet(ConstraintsEnum); 18 18 19 19 /*Assign output pointer: */ -
issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateElementsNodesAndMaterialsSlopeCompute.cpp
r3425 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Partition elements and vertices and nodes: */ -
issm/trunk/src/c/ModelProcessorx/SlopeCompute/CreateLoadsSlopeCompute.cpp
r3456 r3567 15 15 16 16 /*Create loads: */ 17 loads = new DataSet(LoadsEnum ());17 loads = new DataSet(LoadsEnum); 18 18 19 19 /*Now, is the flag isstokes on? otherwise, do nothing: */ -
issm/trunk/src/c/ModelProcessorx/Thermal/CreateConstraintsThermal.cpp
r3455 r3567 19 19 20 20 /*Create constraints: */ 21 constraints = new DataSet(ConstraintsEnum ());21 constraints = new DataSet(ConstraintsEnum); 22 22 23 23 /*return if 2d mesh*/ -
issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp
r3491 r3567 24 24 25 25 /*First create the elements, nodes and material properties: */ 26 elements = new DataSet(ElementsEnum ());27 nodes = new DataSet(NodesEnum ());28 vertices = new DataSet(VerticesEnum ());29 materials = new DataSet(MaterialsEnum ());26 elements = new DataSet(ElementsEnum); 27 nodes = new DataSet(NodesEnum); 28 vertices = new DataSet(VerticesEnum); 29 materials = new DataSet(MaterialsEnum); 30 30 31 31 /*Partition elements and vertices and nodes: */ -
issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp
r3456 r3567 22 22 23 23 /*Create loads: */ 24 loads = new DataSet(LoadsEnum ());24 loads = new DataSet(LoadsEnum); 25 25 26 26 //create penalties for grids: no grid can have a temperature over the melting point -
issm/trunk/src/c/PenaltyConstraintsx/RiftConstraints.cpp
r3332 r3567 46 46 for (i=0;i<loads->Size();i++){ 47 47 48 if(RiftfrontEnum ()==loads->GetEnum(i)){48 if(RiftfrontEnum==loads->GetEnum(i)){ 49 49 50 50 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 77 77 for (i=0;i<loads->Size();i++){ 78 78 79 if(RiftfrontEnum ()==loads->GetEnum(i)){79 if(RiftfrontEnum==loads->GetEnum(i)){ 80 80 found=1; 81 81 break; … … 104 104 for (i=0;i<loads->Size();i++){ 105 105 106 if(RiftfrontEnum ()==loads->GetEnum(i)){106 if(RiftfrontEnum==loads->GetEnum(i)){ 107 107 108 108 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 142 142 for (i=0;i<loads->Size();i++){ 143 143 144 if(RiftfrontEnum ()==loads->GetEnum(i)){144 if(RiftfrontEnum==loads->GetEnum(i)){ 145 145 146 146 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 164 164 for (i=0;i<loads->Size();i++){ 165 165 166 if (RiftfrontEnum ()==loads->GetEnum(i)){166 if (RiftfrontEnum==loads->GetEnum(i)){ 167 167 168 168 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 199 199 for (i=0;i<loads->Size();i++){ 200 200 201 if (RiftfrontEnum ()==loads->GetEnum(i)){201 if (RiftfrontEnum==loads->GetEnum(i)){ 202 202 203 203 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 230 230 for (i=0;i<loads->Size();i++){ 231 231 232 if (RiftfrontEnum ()==loads->GetEnum(i)){232 if (RiftfrontEnum==loads->GetEnum(i)){ 233 233 234 234 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 253 253 for (i=0;i<loads->Size();i++){ 254 254 255 if (RiftfrontEnum ()==loads->GetEnum(i)){255 if (RiftfrontEnum==loads->GetEnum(i)){ 256 256 257 257 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 291 291 for (i=0;i<loads->Size();i++){ 292 292 293 if (RiftfrontEnum ()==loads->GetEnum(i)){293 if (RiftfrontEnum==loads->GetEnum(i)){ 294 294 295 295 riftfront=(Riftfront*)loads->GetObjectByOffset(i); … … 330 330 for (i=0;i<loads->Size();i++){ 331 331 332 if (RiftfrontEnum ()==loads->GetEnum(i)){332 if (RiftfrontEnum==loads->GetEnum(i)){ 333 333 334 334 riftfront=(Riftfront*)loads->GetObjectByOffset(i); -
issm/trunk/src/c/ProcessParamsx/ProcessParamsx.cpp
r3534 r3567 41 41 param_type=parameters->GetEnum(i); 42 42 43 if(param_type==ParamEnum ()){43 if(param_type==ParamEnum){ 44 44 45 45 Param* param=(Param*)parameters->GetObjectByOffset(i); -
issm/trunk/src/c/Qmux/DakotaResponses.cpp
r3446 r3567 238 238 239 239 if(isstokes){ 240 femmodel=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());240 femmodel=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 241 241 } 242 242 if(ismacayealpattyn){ 243 femmodel=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());243 femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 244 244 } 245 245 if(ishutter){ 246 femmodel=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());246 femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 247 247 } 248 248 … … 300 300 301 301 if(isstokes){ 302 femmodel=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());302 femmodel=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 303 303 } 304 304 if(ismacayealpattyn){ 305 femmodel=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());305 femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 306 306 } 307 307 if(ishutter){ 308 femmodel=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());308 femmodel=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 309 309 } 310 310 -
issm/trunk/src/c/Qmux/SpawnCoreParallel.cpp
r3332 r3567 60 60 61 61 /*First off, recover the response descriptors for the response functions: */ 62 param=(Param*)model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum())->parameters->FindParamObject("responsedescriptors");62 param=(Param*)model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum)->parameters->FindParamObject("responsedescriptors"); 63 63 if(!param)ISSMERROR(" could not find response descriptors!"); 64 64 … … 123 123 124 124 //initialize results: 125 results=new DataSet(ResultsEnum ());125 results=new DataSet(ResultsEnum); 126 126 127 127 /*Modify core inputs to reflect the dakota variables inputs: */ 128 inputs->UpdateFromDakota(variables,variables_descriptors,numvariables,model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum())->parameters,qmu_part,qmu_npart); //diagnostic horiz model is the one holding the parameters for Dakota.128 inputs->UpdateFromDakota(variables,variables_descriptors,numvariables,model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum)->parameters,qmu_part,qmu_npart); //diagnostic horiz model is the one holding the parameters for Dakota. 129 129 130 130 /*Run the analysis core solution sequence, with the updated inputs: */ 131 if(analysis_type==DiagnosticAnalysisEnum ()){131 if(analysis_type==DiagnosticAnalysisEnum){ 132 132 133 133 if(verbose)_printf_("Starting diagnostic core\n"); … … 136 136 137 137 } 138 else if(analysis_type==ThermalAnalysisEnum ()){138 else if(analysis_type==ThermalAnalysisEnum){ 139 139 140 140 if(verbose)_printf_("Starting thermal core\n"); … … 142 142 143 143 } 144 else if(analysis_type==PrognosticAnalysisEnum ()){144 else if(analysis_type==PrognosticAnalysisEnum){ 145 145 146 146 if(verbose)_printf_("Starting prognostic core\n"); … … 148 148 149 149 } 150 else if(analysis_type==TransientAnalysisEnum ()){150 else if(analysis_type==TransientAnalysisEnum){ 151 151 152 152 if(verbose)_printf_("Starting transient core\n"); -
issm/trunk/src/c/UpdateGeometryx/UpdateGeometryx.cpp
r3500 r3567 57 57 for(i=0;i<materials->Size();i++){ 58 58 Object* object=materials->GetObjectByOffset(i); 59 if (object->Enum()==MatparEnum ()){59 if (object->Enum()==MatparEnum){ 60 60 matpar=(Matpar*)object; 61 61 break; -
issm/trunk/src/c/objects/Beam.cpp
r3543 r3567 205 205 206 206 /*get enum type of Beam: */ 207 enum_type=BeamEnum ();207 enum_type=BeamEnum; 208 208 209 209 /*marshall enum: */ … … 377 377 378 378 /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */ 379 if (analysis_type==DiagnosticAnalysisEnum ()) {380 381 if (sub_analysis_type==HutterAnalysisEnum ()) {379 if (analysis_type==DiagnosticAnalysisEnum) { 380 381 if (sub_analysis_type==HutterAnalysisEnum) { 382 382 383 383 CreateKMatrixDiagnosticHutter( Kgg,inputs,analysis_type,sub_analysis_type); … … 430 430 431 431 /*Just branch to the correct load generator, according to the type of analysis we are carrying out: */ 432 if (analysis_type==DiagnosticAnalysisEnum ()) {433 if (sub_analysis_type==HutterAnalysisEnum ()) {432 if (analysis_type==DiagnosticAnalysisEnum) { 433 if (sub_analysis_type==HutterAnalysisEnum) { 434 434 CreatePVectorDiagnosticHutter( pg,inputs,analysis_type,sub_analysis_type); 435 435 } … … 577 577 int Beam::Enum(void){ 578 578 579 return BeamEnum ();579 return BeamEnum; 580 580 581 581 } -
issm/trunk/src/c/objects/DofIndexing.cpp
r3483 r3567 150 150 151 151 /*get enum type of DofIndexing: */ 152 enum_type=DofIndexingEnum ();152 enum_type=DofIndexingEnum; 153 153 154 154 /*marshall enum: */ -
issm/trunk/src/c/objects/DofVec.cpp
r3470 r3567 191 191 192 192 /*get enum type of DofVec: */ 193 enum_type=DofVecEnum ();193 enum_type=DofVecEnum; 194 194 195 195 /*marshall enum: */ … … 227 227 int DofVec::Enum(void){ 228 228 229 return DofVecEnum ();229 return DofVecEnum; 230 230 231 231 } -
issm/trunk/src/c/objects/Element.cpp
r3470 r3567 9 9 int Element::Enum(void){ 10 10 11 return ElementEnum ();11 return ElementEnum; 12 12 } -
issm/trunk/src/c/objects/ElementProperties.cpp
r3483 r3567 257 257 258 258 /*get enum type of ElementProperties: */ 259 enum_type=ElementPropertiesEnum ();259 enum_type=ElementPropertiesEnum; 260 260 261 261 /*marshall enum: */ -
issm/trunk/src/c/objects/Hook.cpp
r3483 r3567 117 117 118 118 /*get enum type of Hook: */ 119 enum_type=HookEnum ();119 enum_type=HookEnum; 120 120 121 121 /*marshall enum: */ -
issm/trunk/src/c/objects/Icefront.cpp
r3511 r3567 85 85 if (strcmp(iomodel->meshtype,"2d")==0){ 86 86 segment_width=4; 87 element_type=TriaEnum ();87 element_type=TriaEnum; 88 88 } 89 89 else{ 90 90 segment_width=6; 91 element_type=PentaEnum ();91 element_type=PentaEnum; 92 92 } 93 93 element=(int)(*(iomodel->pressureload+segment_width*i+segment_width-2)-1); //element is in the penultimate column (grid1 grid2 ... elem fill) … … 112 112 icefront_b[1]=iomodel->bed[i2-1]; 113 113 114 if (iomodel->sub_analysis_type==HorizAnalysisEnum ()){115 if ((int)*(iomodel->elements_type+2*element+0)==MacAyealFormulationEnum ()){ //this is a collapsed 3d element, icefront will be 2d114 if (iomodel->sub_analysis_type==HorizAnalysisEnum){ 115 if ((int)*(iomodel->elements_type+2*element+0)==MacAyealFormulationEnum){ //this is a collapsed 3d element, icefront will be 2d 116 116 strcpy(icefront_type,"segment"); 117 117 } 118 else if ((int)*(iomodel->elements_type+2*element+0)==PattynFormulationEnum ()){ //this is a real 3d element, icefront will be 3d.118 else if ((int)*(iomodel->elements_type+2*element+0)==PattynFormulationEnum){ //this is a real 3d element, icefront will be 3d. 119 119 strcpy(icefront_type,"quad"); 120 120 i3=(int)*(iomodel->pressureload+segment_width*i+2); … … 133 133 } 134 134 } 135 else if (iomodel->sub_analysis_type==StokesAnalysisEnum ()){135 else if (iomodel->sub_analysis_type==StokesAnalysisEnum){ 136 136 //We have a Stokes element, so we need a 3d Icefront 137 137 strcpy(icefront_type,"quad"); … … 209 209 210 210 /*get enum type of Icefront: */ 211 enum_type=IcefrontEnum ();211 enum_type=IcefrontEnum; 212 212 213 213 /*marshall enum: */ … … 306 306 307 307 /*Just branch to the correct element icefront vector generator, according to the type of analysis we are carrying out: */ 308 if (analysis_type==ControlAnalysisEnum ()){308 if (analysis_type==ControlAnalysisEnum){ 309 309 310 310 CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type); 311 311 312 312 } 313 else if (analysis_type==DiagnosticAnalysisEnum ()){314 315 if (sub_analysis_type==HorizAnalysisEnum ()){313 else if (analysis_type==DiagnosticAnalysisEnum){ 314 315 if (sub_analysis_type==HorizAnalysisEnum){ 316 316 317 317 CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type); 318 318 319 319 } 320 else if (sub_analysis_type==StokesAnalysisEnum ()){320 else if (sub_analysis_type==StokesAnalysisEnum){ 321 321 322 322 CreatePVectorDiagnosticStokes( pg,inputs,analysis_type,sub_analysis_type); … … 380 380 381 381 //check that the element is onbed (collapsed formulation) otherwise:pe=0 382 if(element_type==PentaEnum ()){382 if(element_type==PentaEnum){ 383 383 if (!element->GetOnBed()){ 384 384 return; … … 389 389 390 390 /*Identify which grids are comprised in the segment. */ 391 if(element_type==TriaEnum ())element_nodes=(Node**)xmalloc(3*sizeof(Node*));392 if(element_type==PentaEnum ())element_nodes=(Node**)xmalloc(6*sizeof(Node*));391 if(element_type==TriaEnum)element_nodes=(Node**)xmalloc(3*sizeof(Node*)); 392 if(element_type==PentaEnum)element_nodes=(Node**)xmalloc(6*sizeof(Node*)); 393 393 element->GetNodes((void**)element_nodes); 394 394 … … 496 496 497 497 /*check icefront is associated to a pentaelem: */ 498 if(element_type!=PentaEnum ()){498 if(element_type!=PentaEnum){ 499 499 ISSMERROR(" quad icefront is associated to a TriaElem!"); 500 500 } … … 510 510 511 511 //Identify which grids are comprised in the quad: 512 if(element_type==PentaEnum ())element_nodes=(Node**)xmalloc(6*sizeof(Node*));512 if(element_type==PentaEnum)element_nodes=(Node**)xmalloc(6*sizeof(Node*)); 513 513 element->GetNodes((void**)element_nodes); 514 514 … … 653 653 654 654 /*check icefront is associated to a pentaelem: */ 655 if(element_type!=PentaEnum ()){655 if(element_type!=PentaEnum){ 656 656 ISSMERROR(" quad icefront is associated to a TriaElem!"); 657 657 } … … 667 667 668 668 //Identify which grids are comprised in the quad: 669 if(element_type==PentaEnum ())element_nodes=(Node**)xmalloc(6*sizeof(Node*));669 if(element_type==PentaEnum)element_nodes=(Node**)xmalloc(6*sizeof(Node*)); 670 670 element->GetNodes((void**)element_nodes); 671 671 … … 855 855 int Icefront::Enum(void){ 856 856 857 return IcefrontEnum ();857 return IcefrontEnum; 858 858 859 859 } … … 1128 1128 1129 1129 //Now deal with water pressure: 1130 if(fill==WaterEnum ()){ //icefront ends in water1130 if(fill==WaterEnum){ //icefront ends in water 1131 1131 water_level_above_g_tria=min(0,z_g[i]);//0 if the gaussian point is above water level 1132 1132 water_pressure_tria=rho_water*gravity*water_level_above_g_tria; 1133 1133 } 1134 else if(fill==AirEnum ()){1134 else if(fill==AirEnum){ 1135 1135 water_pressure_tria=0; 1136 1136 } … … 1382 1382 1383 1383 //Now deal with water pressure: 1384 if(fill==WaterEnum ()){ //icefront ends in water1384 if(fill==WaterEnum){ //icefront ends in water 1385 1385 water_level_above_g_tria=min(0,z_g[i]);//0 if the gaussian point is above water level 1386 1386 water_pressure_tria=rho_water*gravity*water_level_above_g_tria; 1387 1387 } 1388 else if(fill==AirEnum ()){1388 else if(fill==AirEnum){ 1389 1389 water_pressure_tria=0; 1390 1390 } … … 1463 1463 bed=b1*(1+segment_gauss_coord[ig])/2+b2*(1-segment_gauss_coord[ig])/2; 1464 1464 1465 if (fill==WaterEnum ()){1465 if (fill==WaterEnum){ 1466 1466 //icefront ends in water: 1467 1467 ice_pressure=1.0/2.0*gravity*rho_ice*pow(thickness,2); … … 1473 1473 water_pressure=1.0/2.0*gravity*rho_water*(pow(surface_under_water,2) - pow(base_under_water,2)); 1474 1474 } 1475 else if (fill==AirEnum ()){1475 else if (fill==AirEnum){ 1476 1476 ice_pressure=1.0/2.0*gravity*rho_ice*pow(thickness,2); 1477 1477 air_pressure=0; -
issm/trunk/src/c/objects/Input.cpp
r3463 r3567 153 153 int Input::Enum(void){ 154 154 155 return InputEnum ();155 return InputEnum; 156 156 157 157 } -
issm/trunk/src/c/objects/Load.cpp
r246 r3567 10 10 int Load::Enum(void){ 11 11 12 return LoadEnum ();12 return LoadEnum; 13 13 } 14 14 -
issm/trunk/src/c/objects/Material.cpp
r246 r3567 10 10 int Material::Enum(void){ 11 11 12 return MaterialEnum ();12 return MaterialEnum; 13 13 } 14 14 -
issm/trunk/src/c/objects/Matice.cpp
r3541 r3567 102 102 103 103 /*get enum type of Matice: */ 104 enum_type=MaticeEnum ();104 enum_type=MaticeEnum; 105 105 106 106 /*marshall enum: */ … … 175 175 int Matice::Enum(void){ 176 176 177 return MaticeEnum ();177 return MaticeEnum; 178 178 179 179 } -
issm/trunk/src/c/objects/Matpar.cpp
r3509 r3567 46 46 double matpar_g; 47 47 48 if (iomodel->analysis_type==DiagnosticAnalysisEnum () && iomodel->sub_analysis_type==HutterAnalysisEnum()){48 if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==HutterAnalysisEnum){ 49 49 if (strcmp(iomodel->meshtype,"2d")==0){ 50 50 matpar_mid=iomodel->numberofvertices+1; //put it at the end of the materials … … 133 133 134 134 /*get enum type of Matpar: */ 135 enum_type=MatparEnum ();135 enum_type=MatparEnum; 136 136 137 137 /*marshall enum: */ … … 221 221 int Matpar::Enum(void){ 222 222 223 return MatparEnum ();223 return MatparEnum; 224 224 225 225 } -
issm/trunk/src/c/objects/Model.cpp
r3417 r3567 143 143 } 144 144 /*sub_analysis_type should default to none, as it was not specified: */ 145 sub_analysis_type=NoneAnalysisEnum ();145 sub_analysis_type=NoneAnalysisEnum; 146 146 147 147 if (sub_analysis_type!=0){ -
issm/trunk/src/c/objects/Node.cpp
r3557 r3567 95 95 96 96 /*Diagnostic Horiz*/ 97 if (iomodel->analysis_type==DiagnosticAnalysisEnum () && iomodel->sub_analysis_type==HorizAnalysisEnum()){97 if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==HorizAnalysisEnum){ 98 98 if (strcmp(iomodel->meshtype,"3d")==0){ 99 99 /*We have a 3d mesh, we may have collapsed elements, hence dead grids. Freeze them out: */ … … 115 115 116 116 /*Diagnostic Stokes*/ 117 if (iomodel->analysis_type==DiagnosticAnalysisEnum () && iomodel->sub_analysis_type==StokesAnalysisEnum()){117 if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==StokesAnalysisEnum){ 118 118 /*On a 3d mesh, in stokes formualtions, only stokes grids are free, the others are frozen: */ 119 119 if (!iomodel->borderstokes) ISSMERROR("iomodel->borderstokes is NULL"); … … 132 132 133 133 /*Diagnostic Hutter*/ 134 if (iomodel->analysis_type==DiagnosticAnalysisEnum () && iomodel->sub_analysis_type==HutterAnalysisEnum()){134 if (iomodel->analysis_type==DiagnosticAnalysisEnum && iomodel->sub_analysis_type==HutterAnalysisEnum){ 135 135 /*Spc all nodes that are not Hutter*/ 136 136 if (!iomodel->gridonhutter) ISSMERROR("iomodel->gridonhutter is NULL"); … … 144 144 /*Prognostic/ Melting/ Slopecompute/ Balancedvelocities/ Balancedthickness*/ 145 145 if ( 146 iomodel->analysis_type==PrognosticAnalysisEnum ()||147 iomodel->analysis_type==MeltingAnalysisEnum ()||148 iomodel->analysis_type==SlopecomputeAnalysisEnum ()||149 iomodel->analysis_type==BalancedvelocitiesAnalysisEnum ()||150 iomodel->analysis_type==BalancedthicknessAnalysisEnum ()146 iomodel->analysis_type==PrognosticAnalysisEnum || 147 iomodel->analysis_type==MeltingAnalysisEnum || 148 iomodel->analysis_type==SlopecomputeAnalysisEnum || 149 iomodel->analysis_type==BalancedvelocitiesAnalysisEnum || 150 iomodel->analysis_type==BalancedthicknessAnalysisEnum 151 151 ){ 152 152 if (strcmp(iomodel->meshtype,"3d")==0){ … … 295 295 int Node::Enum(void){ 296 296 297 return NodeEnum ();297 return NodeEnum; 298 298 299 299 } … … 354 354 355 355 /*get enum type of Node: */ 356 enum_type=NodeEnum ();356 enum_type=NodeEnum; 357 357 358 358 /*marshall enum: */ -
issm/trunk/src/c/objects/NodeProperties.cpp
r3483 r3567 80 80 81 81 /*get enum type of NodeProperties: */ 82 enum_type=NodePropertiesEnum ();82 enum_type=NodePropertiesEnum; 83 83 84 84 /*marshall enum: */ -
issm/trunk/src/c/objects/Numericalflux.cpp
r3565 r3567 96 96 97 97 /*get enum type of Numericalflux: */ 98 enum_type=NumericalfluxEnum ();98 enum_type=NumericalfluxEnum; 99 99 100 100 /*marshall enum: */ … … 571 571 int Numericalflux::Enum(void){ 572 572 573 return NumericalfluxEnum ();573 return NumericalfluxEnum; 574 574 575 575 } -
issm/trunk/src/c/objects/Numpar.cpp
r3463 r3567 60 60 61 61 /*get enum type of Numpar: */ 62 enum_type=NumparEnum ();62 enum_type=NumparEnum; 63 63 64 64 /*marshall enum: */ … … 204 204 int Numpar::Enum(void){ 205 205 206 return NumparEnum ();206 return NumparEnum; 207 207 208 208 } -
issm/trunk/src/c/objects/Param.cpp
r3463 r3567 233 233 234 234 /*get enum type of Param: */ 235 enum_type=ParamEnum ();235 enum_type=ParamEnum; 236 236 237 237 /*marshall enum: */ … … 437 437 int Param::Enum(void){ 438 438 439 return ParamEnum ();439 return ParamEnum; 440 440 } 441 441 /*}}}*/ -
issm/trunk/src/c/objects/Pengrid.cpp
r3463 r3567 64 64 65 65 /*get enum type of Pengrid: */ 66 enum_type=PengridEnum ();66 enum_type=PengridEnum; 67 67 68 68 /*marshall enum: */ … … 219 219 int Pengrid::Enum(void){ 220 220 221 return PengridEnum ();221 return PengridEnum; 222 222 } 223 223 /*}}}1*/ … … 255 255 void Pengrid::PenaltyConstrain(int* punstable,void* vinputs,int analysis_type,int sub_analysis_type){ 256 256 257 if ((analysis_type==DiagnosticAnalysisEnum ()) && ((sub_analysis_type==StokesAnalysisEnum()))){257 if ((analysis_type==DiagnosticAnalysisEnum) && ((sub_analysis_type==StokesAnalysisEnum))){ 258 258 259 259 /*No penalty to check*/ … … 261 261 262 262 } 263 else if (analysis_type==ThermalAnalysisEnum ()){263 else if (analysis_type==ThermalAnalysisEnum){ 264 264 265 265 PenaltyConstrainThermal(punstable,vinputs,analysis_type,sub_analysis_type); 266 266 267 267 } 268 else if (analysis_type==MeltingAnalysisEnum ()){268 else if (analysis_type==MeltingAnalysisEnum){ 269 269 270 270 /*No penalty to check*/ … … 362 362 void Pengrid::PenaltyCreateKMatrix(Mat Kgg,void* inputs,double kmax,int analysis_type,int sub_analysis_type){ 363 363 364 if ((analysis_type==DiagnosticAnalysisEnum ()) && ((sub_analysis_type==StokesAnalysisEnum()))){364 if ((analysis_type==DiagnosticAnalysisEnum) && ((sub_analysis_type==StokesAnalysisEnum))){ 365 365 366 366 PenaltyCreateKMatrixDiagnosticStokes( Kgg,inputs,kmax,analysis_type,sub_analysis_type); 367 367 } 368 else if (analysis_type==ThermalAnalysisEnum ()){368 else if (analysis_type==ThermalAnalysisEnum){ 369 369 370 370 PenaltyCreateKMatrixThermal( Kgg,inputs,kmax,analysis_type,sub_analysis_type); 371 371 372 372 } 373 else if (analysis_type==MeltingAnalysisEnum ()){373 else if (analysis_type==MeltingAnalysisEnum){ 374 374 375 375 PenaltyCreateKMatrixMelting( Kgg,inputs,kmax,analysis_type,sub_analysis_type); … … 500 500 void Pengrid::PenaltyCreatePVector(Vec pg,void* inputs,double kmax,int analysis_type,int sub_analysis_type){ 501 501 502 if (analysis_type==ThermalAnalysisEnum ()){502 if (analysis_type==ThermalAnalysisEnum){ 503 503 504 504 PenaltyCreatePVectorThermal( pg,inputs,kmax,analysis_type,sub_analysis_type); 505 505 506 506 } 507 else if (analysis_type==MeltingAnalysisEnum ()){507 else if (analysis_type==MeltingAnalysisEnum){ 508 508 509 509 PenaltyCreatePVectorMelting( pg,inputs,kmax,analysis_type,sub_analysis_type); 510 510 511 511 } 512 else if (analysis_type==DiagnosticAnalysisEnum ()){512 else if (analysis_type==DiagnosticAnalysisEnum){ 513 513 514 514 /*No loads applied, do nothing: */ -
issm/trunk/src/c/objects/Penpair.cpp
r3463 r3567 61 61 62 62 /*get enum type of Penpair: */ 63 enum_type=PenpairEnum ();63 enum_type=PenpairEnum; 64 64 65 65 /*marshall enum: */ … … 198 198 int Penpair::Enum(void){ 199 199 200 return PenpairEnum ();200 return PenpairEnum; 201 201 } 202 202 /*}}}1*/ -
issm/trunk/src/c/objects/Penta.cpp
r3529 r3567 81 81 //elements of type 3 are macayeal type penta. we collapse the formulation on their base. 82 82 if(iomodel->elements_type){ 83 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum ()){83 if (*(iomodel->elements_type+2*i+0)==MacAyealFormulationEnum){ 84 84 this->properties.collapse=1; 85 85 } … … 196 196 197 197 /*get enum type of Penta: */ 198 enum_type=PentaEnum ();198 enum_type=PentaEnum; 199 199 200 200 /*marshall enum: */ … … 387 387 388 388 /*Check analysis_types*/ 389 if (analysis_type!=DiagnosticAnalysisEnum () || sub_analysis_type!=StokesAnalysisEnum()) ISSMERROR("Not supported yet!");389 if (analysis_type!=DiagnosticAnalysisEnum || sub_analysis_type!=StokesAnalysisEnum) ISSMERROR("Not supported yet!"); 390 390 391 391 /*recover pointers: */ … … 555 555 556 556 /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */ 557 if (analysis_type==ControlAnalysisEnum ()){557 if (analysis_type==ControlAnalysisEnum){ 558 558 559 559 CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type); 560 560 561 561 } 562 else if (analysis_type==DiagnosticAnalysisEnum ()){563 564 if (sub_analysis_type==HorizAnalysisEnum ()){562 else if (analysis_type==DiagnosticAnalysisEnum){ 563 564 if (sub_analysis_type==HorizAnalysisEnum){ 565 565 566 566 CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type); 567 567 } 568 else if (sub_analysis_type==VertAnalysisEnum ()){568 else if (sub_analysis_type==VertAnalysisEnum){ 569 569 570 570 CreateKMatrixDiagnosticVert( Kgg,inputs,analysis_type,sub_analysis_type); 571 571 } 572 else if (sub_analysis_type==StokesAnalysisEnum ()){572 else if (sub_analysis_type==StokesAnalysisEnum){ 573 573 574 574 CreateKMatrixDiagnosticStokes( Kgg,inputs,analysis_type,sub_analysis_type); … … 577 577 else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet")); 578 578 } 579 else if (analysis_type==SlopecomputeAnalysisEnum ()){579 else if (analysis_type==SlopecomputeAnalysisEnum){ 580 580 581 581 CreateKMatrixSlopeCompute( Kgg,inputs,analysis_type,sub_analysis_type); 582 582 } 583 else if (analysis_type==PrognosticAnalysisEnum ()){583 else if (analysis_type==PrognosticAnalysisEnum){ 584 584 585 585 CreateKMatrixPrognostic( Kgg,inputs,analysis_type,sub_analysis_type); 586 586 } 587 else if (analysis_type==BalancedthicknessAnalysisEnum ()){587 else if (analysis_type==BalancedthicknessAnalysisEnum){ 588 588 589 589 CreateKMatrixBalancedthickness( Kgg,inputs,analysis_type,sub_analysis_type); 590 590 } 591 else if (analysis_type==BalancedvelocitiesAnalysisEnum ()){591 else if (analysis_type==BalancedvelocitiesAnalysisEnum){ 592 592 593 593 CreateKMatrixBalancedvelocities( Kgg,inputs,analysis_type,sub_analysis_type); 594 594 } 595 else if (analysis_type==ThermalAnalysisEnum ()){595 else if (analysis_type==ThermalAnalysisEnum){ 596 596 597 597 CreateKMatrixThermal( Kgg,inputs,analysis_type,sub_analysis_type); 598 598 } 599 else if (analysis_type==MeltingAnalysisEnum ()){599 else if (analysis_type==MeltingAnalysisEnum){ 600 600 601 601 CreateKMatrixMelting( Kgg,inputs,analysis_type,sub_analysis_type); … … 1683 1683 1684 1684 /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */ 1685 if (analysis_type==ControlAnalysisEnum ()){1685 if (analysis_type==ControlAnalysisEnum){ 1686 1686 1687 1687 CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type); 1688 1688 } 1689 else if (analysis_type==DiagnosticAnalysisEnum ()){1690 1691 if (sub_analysis_type==HorizAnalysisEnum ()){1689 else if (analysis_type==DiagnosticAnalysisEnum){ 1690 1691 if (sub_analysis_type==HorizAnalysisEnum){ 1692 1692 1693 1693 CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type); 1694 1694 } 1695 else if (sub_analysis_type==VertAnalysisEnum ()){1695 else if (sub_analysis_type==VertAnalysisEnum){ 1696 1696 1697 1697 CreatePVectorDiagnosticVert( pg,inputs,analysis_type,sub_analysis_type); 1698 1698 } 1699 else if (sub_analysis_type==StokesAnalysisEnum ()){1699 else if (sub_analysis_type==StokesAnalysisEnum){ 1700 1700 1701 1701 CreatePVectorDiagnosticStokes( pg,inputs,analysis_type,sub_analysis_type); … … 1703 1703 else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet")); 1704 1704 } 1705 else if (analysis_type==SlopecomputeAnalysisEnum ()){1705 else if (analysis_type==SlopecomputeAnalysisEnum){ 1706 1706 1707 1707 CreatePVectorSlopeCompute( pg,inputs,analysis_type,sub_analysis_type); 1708 1708 } 1709 else if (analysis_type==PrognosticAnalysisEnum ()){1709 else if (analysis_type==PrognosticAnalysisEnum){ 1710 1710 1711 1711 CreatePVectorPrognostic( pg,inputs,analysis_type,sub_analysis_type); 1712 1712 } 1713 else if (analysis_type==BalancedthicknessAnalysisEnum ()){1713 else if (analysis_type==BalancedthicknessAnalysisEnum){ 1714 1714 1715 1715 CreatePVectorBalancedthickness( pg,inputs,analysis_type,sub_analysis_type); 1716 1716 } 1717 else if (analysis_type==BalancedvelocitiesAnalysisEnum ()){1717 else if (analysis_type==BalancedvelocitiesAnalysisEnum){ 1718 1718 1719 1719 CreatePVectorBalancedvelocities( pg,inputs,analysis_type,sub_analysis_type); 1720 1720 } 1721 else if (analysis_type==ThermalAnalysisEnum ()){1721 else if (analysis_type==ThermalAnalysisEnum){ 1722 1722 1723 1723 CreatePVectorThermal( pg,inputs,analysis_type,sub_analysis_type); 1724 1724 } 1725 else if (analysis_type==MeltingAnalysisEnum ()){1725 else if (analysis_type==MeltingAnalysisEnum){ 1726 1726 1727 1727 CreatePVectorMelting( pg,inputs,analysis_type,sub_analysis_type); … … 2620 2620 int Penta::Enum(void){ 2621 2621 2622 return PentaEnum ();2622 return PentaEnum; 2623 2623 2624 2624 } … … 4059 4059 if (!this->properties.onbed) return; 4060 4060 4061 if (sub_analysis_type==HorizAnalysisEnum ()){4061 if (sub_analysis_type==HorizAnalysisEnum){ 4062 4062 4063 4063 /*MacAyeal or Pattyn*/ … … 4067 4067 return; 4068 4068 } 4069 else if (sub_analysis_type==StokesAnalysisEnum ()){4069 else if (sub_analysis_type==StokesAnalysisEnum){ 4070 4070 4071 4071 /*Stokes*/ -
issm/trunk/src/c/objects/Result.cpp
r3463 r3567 202 202 int Result::Enum(void){ 203 203 204 return ResultEnum ();204 return ResultEnum; 205 205 206 206 } -
issm/trunk/src/c/objects/Rgb.cpp
r3463 r3567 52 52 53 53 /*get enum type of Rgb: */ 54 enum_type=RgbEnum ();54 enum_type=RgbEnum; 55 55 56 56 /*marshall enum: */ … … 130 130 int Rgb::Enum(void){ 131 131 132 return RgbEnum ();132 return RgbEnum; 133 133 134 134 } -
issm/trunk/src/c/objects/Riftfront.cpp
r3463 r3567 182 182 183 183 /*get enum type of Riftfront: */ 184 enum_type=RiftfrontEnum ();184 enum_type=RiftfrontEnum; 185 185 186 186 /*marshall enum: */ … … 453 453 int Riftfront::Enum(void){ 454 454 455 return RiftfrontEnum ();455 return RiftfrontEnum; 456 456 457 457 } … … 729 729 730 730 /*Ok, now compute the pressure (in norm) that is being applied to the flanks, depending on the type of fill: */ 731 if(fill==WaterEnum ()){731 if(fill==WaterEnum){ 732 732 if(shelf){ 733 733 /*We are on an ice shelf, hydrostatic equilibrium is used to determine the pressure for water fill: */ … … 739 739 } 740 740 } 741 else if(fill==AirEnum ()){741 else if(fill==AirEnum){ 742 742 pressure=rho_ice*gravity*pow(thickness,(double)2)/(double)2; //icefront on an ice sheet, pressure imbalance ice vs air. 743 743 } 744 else if(fill==IceEnum ()){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice)744 else if(fill==IceEnum){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice) 745 745 pressure=0; 746 746 } 747 else if(fill==MelangeEnum ()){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice)747 else if(fill==MelangeEnum){ //icefront finding itself against another icefront (pressure imbalance is fully compensated, ice vs ice) 748 748 749 749 if(!shelf) ISSMERROR(exprintf("%s%i%s","fill type ",fill," not supported on ice sheets yet.")); -
issm/trunk/src/c/objects/Sing.cpp
r3529 r3567 193 193 194 194 /*get enum type of Sing: */ 195 enum_type=SingEnum ();195 enum_type=SingEnum; 196 196 197 197 /*marshall enum: */ … … 348 348 349 349 /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */ 350 if ((analysis_type==DiagnosticAnalysisEnum ()) && (sub_analysis_type==HutterAnalysisEnum())){350 if ((analysis_type==DiagnosticAnalysisEnum) && (sub_analysis_type==HutterAnalysisEnum)){ 351 351 352 352 CreateKMatrixDiagnosticHutter( Kgg,inputs,analysis_type,sub_analysis_type); … … 385 385 386 386 /*Just branch to the correct load generator, according to the type of analysis we are carrying out: */ 387 if ((analysis_type==DiagnosticAnalysisEnum ()) && (sub_analysis_type==HutterAnalysisEnum())){387 if ((analysis_type==DiagnosticAnalysisEnum) && (sub_analysis_type==HutterAnalysisEnum)){ 388 388 389 389 CreatePVectorDiagnosticHutter( pg,inputs,analysis_type,sub_analysis_type); … … 472 472 int Sing::Enum(void){ 473 473 474 return SingEnum ();474 return SingEnum; 475 475 476 476 } -
issm/trunk/src/c/objects/Spc.cpp
r3463 r3567 53 53 54 54 /*get enum type of Spc: */ 55 enum_type=SpcEnum ();55 enum_type=SpcEnum; 56 56 57 57 /*marshall enum: */ … … 130 130 int Spc::Enum(void){ 131 131 132 return SpcEnum ();132 return SpcEnum; 133 133 134 134 } -
issm/trunk/src/c/objects/Tria.cpp
r3556 r3567 77 77 /*hooks: */ 78 78 //go recover node ids, needed to initialize the node hook. 79 if (iomodel->analysis_type==Prognostic2AnalysisEnum ()){79 if (iomodel->analysis_type==Prognostic2AnalysisEnum){ 80 80 /*Discontinuous Galerkin*/ 81 81 tria_node_ids[0]=3*i+1; … … 209 209 210 210 /*get enum type of Tria: */ 211 enum_type=TriaEnum ();211 enum_type=TriaEnum; 212 212 213 213 /*marshall enum: */ … … 569 569 570 570 /*Just branch to the correct element stiffness matrix generator, according to the type of analysis we are carrying out: */ 571 if (analysis_type==ControlAnalysisEnum ()){571 if (analysis_type==ControlAnalysisEnum){ 572 572 573 573 CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type); 574 574 } 575 else if (analysis_type==DiagnosticAnalysisEnum ()){575 else if (analysis_type==DiagnosticAnalysisEnum){ 576 576 577 if (sub_analysis_type==HorizAnalysisEnum ()){577 if (sub_analysis_type==HorizAnalysisEnum){ 578 578 579 579 CreateKMatrixDiagnosticHoriz( Kgg,inputs,analysis_type,sub_analysis_type); … … 582 582 583 583 } 584 else if (analysis_type==SlopecomputeAnalysisEnum ()){584 else if (analysis_type==SlopecomputeAnalysisEnum){ 585 585 586 586 CreateKMatrixSlopeCompute( Kgg,inputs,analysis_type,sub_analysis_type); 587 587 } 588 else if (analysis_type==PrognosticAnalysisEnum ()){588 else if (analysis_type==PrognosticAnalysisEnum){ 589 589 590 590 CreateKMatrixPrognostic( Kgg,inputs,analysis_type,sub_analysis_type); 591 591 } 592 else if (analysis_type==Prognostic2AnalysisEnum ()){592 else if (analysis_type==Prognostic2AnalysisEnum){ 593 593 594 594 CreateKMatrixPrognostic2(Kgg,inputs,analysis_type,sub_analysis_type); 595 595 } 596 else if (analysis_type==BalancedthicknessAnalysisEnum ()){596 else if (analysis_type==BalancedthicknessAnalysisEnum){ 597 597 598 598 CreateKMatrixBalancedthickness( Kgg,inputs,analysis_type,sub_analysis_type); 599 599 } 600 else if (analysis_type==BalancedvelocitiesAnalysisEnum ()){600 else if (analysis_type==BalancedvelocitiesAnalysisEnum){ 601 601 602 602 CreateKMatrixBalancedvelocities( Kgg,inputs,analysis_type,sub_analysis_type); … … 2177 2177 2178 2178 /*Just branch to the correct load generator, according to the type of analysis we are carrying out: */ 2179 if (analysis_type==ControlAnalysisEnum ()){2179 if (analysis_type==ControlAnalysisEnum){ 2180 2180 2181 2181 CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type); 2182 2182 2183 2183 } 2184 else if (analysis_type==DiagnosticAnalysisEnum ()){2185 if (sub_analysis_type==HorizAnalysisEnum ()){2184 else if (analysis_type==DiagnosticAnalysisEnum){ 2185 if (sub_analysis_type==HorizAnalysisEnum){ 2186 2186 2187 2187 CreatePVectorDiagnosticHoriz( pg,inputs,analysis_type,sub_analysis_type); … … 2189 2189 else ISSMERROR(exprintf("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet")); 2190 2190 } 2191 else if (analysis_type==SlopecomputeAnalysisEnum ()){2191 else if (analysis_type==SlopecomputeAnalysisEnum){ 2192 2192 2193 2193 CreatePVectorSlopeCompute( pg,inputs,analysis_type,sub_analysis_type); 2194 2194 } 2195 else if (analysis_type==PrognosticAnalysisEnum ()){2195 else if (analysis_type==PrognosticAnalysisEnum){ 2196 2196 2197 2197 CreatePVectorPrognostic( pg,inputs,analysis_type,sub_analysis_type); 2198 2198 } 2199 else if (analysis_type==Prognostic2AnalysisEnum ()){2199 else if (analysis_type==Prognostic2AnalysisEnum){ 2200 2200 2201 2201 CreatePVectorPrognostic2( pg,inputs,analysis_type,sub_analysis_type); 2202 2202 } 2203 else if (analysis_type==BalancedthicknessAnalysisEnum ()){2203 else if (analysis_type==BalancedthicknessAnalysisEnum){ 2204 2204 2205 2205 CreatePVectorBalancedthickness( pg,inputs,analysis_type,sub_analysis_type); 2206 2206 } 2207 else if (analysis_type==BalancedvelocitiesAnalysisEnum ()){2207 else if (analysis_type==BalancedvelocitiesAnalysisEnum){ 2208 2208 2209 2209 CreatePVectorBalancedvelocities( pg,inputs,analysis_type,sub_analysis_type); … … 3004 3004 for(i=0;i<numdof;i++) pe_g[i]=0.0; 3005 3005 3006 if ( (sub_analysis_type==SurfaceXAnalysisEnum ()) || (sub_analysis_type==SurfaceYAnalysisEnum())){3006 if ( (sub_analysis_type==SurfaceXAnalysisEnum) || (sub_analysis_type==SurfaceYAnalysisEnum)){ 3007 3007 for(i=0;i<numdof;i++) param[i]=this->properties.s[i]; 3008 3008 } 3009 if ( (sub_analysis_type==BedXAnalysisEnum ()) || (sub_analysis_type==BedYAnalysisEnum())){3009 if ( (sub_analysis_type==BedXAnalysisEnum) || (sub_analysis_type==BedYAnalysisEnum)){ 3010 3010 for(i=0;i<numdof;i++) param[i]=this->properties.b[i]; 3011 3011 } … … 3032 3032 3033 3033 /*Build pe_g_gaussian vector: */ 3034 if ( (sub_analysis_type==SurfaceXAnalysisEnum ()) || (sub_analysis_type==BedXAnalysisEnum())){3034 if ( (sub_analysis_type==SurfaceXAnalysisEnum) || (sub_analysis_type==BedXAnalysisEnum)){ 3035 3035 for(i=0;i<numdof;i++) pe_g_gaussian[i]=Jdet*gauss_weight*slope[0]*l1l2l3[i]; 3036 3036 } 3037 if ( (sub_analysis_type==SurfaceYAnalysisEnum ()) || (sub_analysis_type==BedYAnalysisEnum())){3037 if ( (sub_analysis_type==SurfaceYAnalysisEnum) || (sub_analysis_type==BedYAnalysisEnum)){ 3038 3038 for(i=0;i<numdof;i++) pe_g_gaussian[i]=Jdet*gauss_weight*slope[1]*l1l2l3[i]; 3039 3039 } … … 3613 3613 int Tria::Enum(void){ 3614 3614 3615 return TriaEnum ();3615 return TriaEnum; 3616 3616 3617 3617 } -
issm/trunk/src/c/objects/Vertex.cpp
r3479 r3567 116 116 int Vertex::Enum(void){ 117 117 118 return VertexEnum ();118 return VertexEnum; 119 119 120 120 } … … 138 138 139 139 /*get enum type of Vertex: */ 140 enum_type=VertexEnum ();140 enum_type=VertexEnum; 141 141 142 142 /*marshall enum: */ -
issm/trunk/src/c/parallel/ControlInitialization.cpp
r3529 r3567 60 60 61 61 /*recover fem models: */ 62 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());63 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum (),VertAnalysisEnum());64 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());65 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());66 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum ());62 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 63 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum); 64 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 65 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 66 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum); 67 67 68 68 //specific parameters for specific models … … 81 81 //compute slopes 82 82 if(verbose)_printf_("%s\n","computing bed slope (x and y derivatives)..."); 83 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum (),BedXAnalysisEnum());84 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum (),BedYAnalysisEnum());83 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum,BedXAnalysisEnum); 84 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum,BedYAnalysisEnum); 85 85 FieldExtrudex( slopex, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopex",0); 86 86 FieldExtrudex( slopey, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopey",0); … … 95 95 //horizontal velocity 96 96 if(verbose)_printf_("%s\n"," computing horizontal velocities..."); 97 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_dh,inputs,DiagnosticAnalysisEnum (),HorizAnalysisEnum());97 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_dh,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum); 98 98 if(verbose)_printf_("%s\n"," extruding horizontal velocities..."); 99 99 VecDuplicatePatch(&ug_horiz,ug); FieldExtrudex( ug_horiz,fem_dh->elements,fem_dh->nodes, fem_dh->vertices,fem_dh->loads,fem_dh-> materials,fem_dh->parameters,"velocity",1); … … 102 102 if(verbose)_printf_("%s\n"," computing vertical velocities..."); 103 103 inputs->Add("velocity",ug_horiz,numberofdofspernode_dh,numberofnodes); 104 diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum (),VertAnalysisEnum());104 diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum,VertAnalysisEnum); 105 105 106 106 //Create 3d u_g … … 113 113 //Create 4d u_g 114 114 if(verbose)_printf_("%s\n"," computing pressure according to Pattyn..."); 115 ComputePressurex( &pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads, fem_dh->materials, fem_dh->parameters,inputs,DiagnosticAnalysisEnum (),HorizAnalysisEnum());115 ComputePressurex( &pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads, fem_dh->materials, fem_dh->parameters,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum); 116 116 VecScale(pg,1.0/stokesreconditioning); 117 117 ug_stokes=NewVec(fem_ds->nodesets->GetGSize()); … … 133 133 if(verbose)_printf_("%s\n"," computing stokes velocities and pressure ..."); 134 134 VecFree(&ug); 135 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum (),StokesAnalysisEnum());135 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum,StokesAnalysisEnum); 136 136 137 137 //Add in inputs -
issm/trunk/src/c/parallel/ControlRestart.cpp
r3332 r3567 34 34 for(i=0;i<numberofnodes;i++) param_g_copy[i]=param_g[i]; 35 35 36 temporary_results=new DataSet(ResultsEnum ());36 temporary_results=new DataSet(ResultsEnum); 37 37 result=new Result(temporary_results->Size()+1,0,1,"param_g",param_g_copy,numberofnodes); 38 38 temporary_results->AddObject(result); … … 42 42 43 43 //process results 44 ProcessResults(&results,temporary_results,model,ControlAnalysisEnum ());44 ProcessResults(&results,temporary_results,model,ControlAnalysisEnum); 45 45 46 46 //Write results on disk -
issm/trunk/src/c/parallel/ProcessResults.cpp
r3556 r3567 109 109 110 110 /*Initialize new results: */ 111 newresults=new DataSet(ResultsEnum ());111 newresults=new DataSet(ResultsEnum); 112 112 113 113 /*some flags needed: */ … … 118 118 119 119 /*Recover femmodels first: */ 120 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());120 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 121 121 fem_c=fem_dh; 122 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum (),VertAnalysisEnum());123 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());124 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());125 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum ());126 if(analysis_type==PrognosticAnalysisEnum ()){127 fem_p=model->GetFormulation(PrognosticAnalysisEnum ());122 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum); 123 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 124 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 125 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum); 126 if(analysis_type==PrognosticAnalysisEnum){ 127 fem_p=model->GetFormulation(PrognosticAnalysisEnum); 128 128 } 129 if(analysis_type==Prognostic2AnalysisEnum ()){130 fem_p=model->GetFormulation(Prognostic2AnalysisEnum ());129 if(analysis_type==Prognostic2AnalysisEnum){ 130 fem_p=model->GetFormulation(Prognostic2AnalysisEnum); 131 131 } 132 if(analysis_type==TransientAnalysisEnum ()){133 fem_p=model->GetFormulation(PrognosticAnalysisEnum ());132 if(analysis_type==TransientAnalysisEnum){ 133 fem_p=model->GetFormulation(PrognosticAnalysisEnum); 134 134 } 135 if(analysis_type==BalancedthicknessAnalysisEnum ()){136 fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum ());135 if(analysis_type==BalancedthicknessAnalysisEnum){ 136 fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum); 137 137 } 138 if(analysis_type==BalancedvelocitiesAnalysisEnum ()){139 fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum ());138 if(analysis_type==BalancedvelocitiesAnalysisEnum){ 139 fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum); 140 140 } 141 fem_t=model->GetFormulation(ThermalAnalysisEnum ());141 fem_t=model->GetFormulation(ThermalAnalysisEnum); 142 142 143 143 for(n=0;n<results->Size();n++){ -
issm/trunk/src/c/parallel/balancedthickness.cpp
r3332 r3567 76 76 77 77 _printf_("read and create finite element model:\n"); 78 model->AddFormulation(fid,BalancedthicknessAnalysisEnum ());78 model->AddFormulation(fid,BalancedthicknessAnalysisEnum); 79 79 80 80 /*recover parameters: */ … … 84 84 _printf_("initialize inputs:\n"); 85 85 86 model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedthicknessAnalysisEnum ());87 model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedthicknessAnalysisEnum ());88 model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedthicknessAnalysisEnum ());86 model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedthicknessAnalysisEnum); 87 model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedthicknessAnalysisEnum); 88 model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedthicknessAnalysisEnum); 89 89 model->FindParam(&numberofnodes,"numberofnodes"); 90 90 … … 95 95 96 96 _printf_("initialize results:\n"); 97 results=new DataSet(ResultsEnum ());97 results=new DataSet(ResultsEnum); 98 98 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 99 99 … … 115 115 #ifdef _HAVE_DAKOTA_ 116 116 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 117 Qmux(model,inputs,BalancedthicknessAnalysisEnum (),NoneAnalysisEnum());117 Qmux(model,inputs,BalancedthicknessAnalysisEnum,NoneAnalysisEnum); 118 118 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 119 119 #else … … 127 127 128 128 _printf_("process results:\n"); 129 ProcessResults(&processedresults,results,model,BalancedthicknessAnalysisEnum ());129 ProcessResults(&processedresults,results,model,BalancedthicknessAnalysisEnum); 130 130 131 131 _printf_("write results to disk:\n"); -
issm/trunk/src/c/parallel/balancedthickness_core.cpp
r3446 r3567 35 35 36 36 /*recover fem model: */ 37 fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum ());37 fem_p=model->GetFormulation(BalancedthicknessAnalysisEnum); 38 38 39 39 //first recover parameters common to all solutions … … 48 48 49 49 _printf_("call computational core:\n"); 50 diagnostic_core_linear(&h_g,fem_p,inputs,BalancedthicknessAnalysisEnum (),NoneAnalysisEnum());50 diagnostic_core_linear(&h_g,fem_p,inputs,BalancedthicknessAnalysisEnum,NoneAnalysisEnum); 51 51 52 52 _printf_("extrude computed thickness on all layers:\n"); -
issm/trunk/src/c/parallel/balancedvelocities.cpp
r3332 r3567 73 73 74 74 _printf_("read and create finite element model:\n"); 75 model->AddFormulation(fid,BalancedvelocitiesAnalysisEnum ());75 model->AddFormulation(fid,BalancedvelocitiesAnalysisEnum); 76 76 77 77 /*recover parameters: */ … … 81 81 _printf_("initialize inputs:\n"); 82 82 83 model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedvelocitiesAnalysisEnum ());84 model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedvelocitiesAnalysisEnum ());85 model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedvelocitiesAnalysisEnum ());83 model->FindParam(&u_g_serial,NULL,NULL,"u_g",BalancedvelocitiesAnalysisEnum); 84 model->FindParam(&melting_g,NULL,NULL,"m_g",BalancedvelocitiesAnalysisEnum); 85 model->FindParam(&accumulation_g,NULL,NULL,"a_g",BalancedvelocitiesAnalysisEnum); 86 86 model->FindParam(&numberofnodes,"numberofnodes"); 87 87 … … 92 92 93 93 _printf_("initialize results:\n"); 94 results=new DataSet(ResultsEnum ());94 results=new DataSet(ResultsEnum); 95 95 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 96 96 … … 112 112 #ifdef _HAVE_DAKOTA_ 113 113 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 114 Qmux(model,inputs,BalancedvelocitiesAnalysisEnum (),NoneAnalysisEnum());114 Qmux(model,inputs,BalancedvelocitiesAnalysisEnum,NoneAnalysisEnum); 115 115 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 116 116 #else … … 124 124 125 125 _printf_("process results:\n"); 126 ProcessResults(&processedresults,results,model,BalancedvelocitiesAnalysisEnum ());126 ProcessResults(&processedresults,results,model,BalancedvelocitiesAnalysisEnum); 127 127 128 128 _printf_("write results to disk:\n"); -
issm/trunk/src/c/parallel/balancedvelocities_core.cpp
r3446 r3567 34 34 35 35 /*recover fem model: */ 36 fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum ());36 fem_p=model->GetFormulation(BalancedvelocitiesAnalysisEnum); 37 37 38 38 //first recover parameters common to all solutions … … 47 47 48 48 _printf_("call computational core:\n"); 49 diagnostic_core_linear(&v_g,fem_p,inputs,BalancedvelocitiesAnalysisEnum (),NoneAnalysisEnum());49 diagnostic_core_linear(&v_g,fem_p,inputs,BalancedvelocitiesAnalysisEnum,NoneAnalysisEnum); 50 50 51 51 _printf_("extrude computed thickness on all layers:\n"); -
issm/trunk/src/c/parallel/control_core.cpp
r3446 r3567 97 97 * distribution: */ 98 98 if (control_steady){ 99 steadystate_results=new DataSet(ResultsEnum ());99 steadystate_results=new DataSet(ResultsEnum); 100 100 steadystate_core(steadystate_results,model,inputs); 101 101 VecFree(&t_g); steadystate_results->FindResult(&t_g,"t_g"); … … 108 108 109 109 _printf_("%s\n"," computing gradJ..."); 110 gradjcompute_results=new DataSet(ResultsEnum ());110 gradjcompute_results=new DataSet(ResultsEnum); 111 111 gradjcompute_core(gradjcompute_results,model, inputs); 112 112 gradjcompute_results->FindResult(&grad_g,"grad_g"); … … 205 205 if (control_steady){ 206 206 inputs->Add(control_type,param_g,1,numberofnodes); 207 steadystate_results=new DataSet(ResultsEnum ());207 steadystate_results=new DataSet(ResultsEnum); 208 208 steadystate_core(steadystate_results,model,inputs); 209 209 … … 216 216 else{ 217 217 inputs->Add(control_type,param_g,1,numberofnodes); 218 diagnostic_results=new DataSet(ResultsEnum ());218 diagnostic_results=new DataSet(ResultsEnum); 219 219 diagnostic_core(diagnostic_results,model, inputs); 220 220 -
issm/trunk/src/c/parallel/diagnostic.cpp
r3332 r3567 76 76 _printf_("read and create finite element model:\n"); 77 77 _printf_("\n reading diagnostic horiz model data:\n"); 78 model->AddFormulation(fid,DiagnosticAnalysisEnum (),HorizAnalysisEnum());78 model->AddFormulation(fid,DiagnosticAnalysisEnum,HorizAnalysisEnum); 79 79 80 80 _printf_("\n reading diagnostic vert model data:\n"); 81 model->AddFormulation(fid,DiagnosticAnalysisEnum (),VertAnalysisEnum());81 model->AddFormulation(fid,DiagnosticAnalysisEnum,VertAnalysisEnum); 82 82 83 83 _printf_("\n reading diagnostic stokes model data:\n"); 84 model->AddFormulation(fid,DiagnosticAnalysisEnum (),StokesAnalysisEnum());84 model->AddFormulation(fid,DiagnosticAnalysisEnum,StokesAnalysisEnum); 85 85 86 86 _printf_("\n reading diagnostic hutter model data:\n"); 87 model->AddFormulation(fid,DiagnosticAnalysisEnum (),HutterAnalysisEnum());87 model->AddFormulation(fid,DiagnosticAnalysisEnum,HutterAnalysisEnum); 88 88 89 89 _printf_("\n reading surface and bed slope computation model data:\n"); 90 model->AddFormulation(fid,SlopecomputeAnalysisEnum ());90 model->AddFormulation(fid,SlopecomputeAnalysisEnum); 91 91 92 92 /*get parameters: */ … … 96 96 97 97 _printf_("initialize inputs:\n"); 98 model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum (),HorizAnalysisEnum());98 model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum,HorizAnalysisEnum); 99 99 model->FindParam(&numberofnodes,"numberofnodes"); 100 100 … … 103 103 104 104 if(control_analysis){ 105 model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum (),HorizAnalysisEnum());106 model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum (),HorizAnalysisEnum());105 model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum,HorizAnalysisEnum); 106 model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum,HorizAnalysisEnum); 107 107 inputs->Add("velocity_obs",u_g_obs,2,numberofnodes); 108 108 inputs->Add("weights",weights,1,numberofnodes); … … 110 110 111 111 _printf_("initialize results:\n"); 112 results=new DataSet(ResultsEnum ());112 results=new DataSet(ResultsEnum); 113 113 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 114 114 … … 127 127 128 128 _printf_("process results:\n"); 129 ProcessResults(&processed_results,results,model,DiagnosticAnalysisEnum ());129 ProcessResults(&processed_results,results,model,DiagnosticAnalysisEnum); 130 130 } 131 131 else{ … … 145 145 146 146 _printf_("process results:\n"); 147 ProcessResults(&processed_results,results,model,ControlAnalysisEnum ());147 ProcessResults(&processed_results,results,model,ControlAnalysisEnum); 148 148 } 149 149 … … 157 157 #ifdef _HAVE_DAKOTA_ 158 158 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 159 Qmux(model,inputs,DiagnosticAnalysisEnum (),NoneAnalysisEnum());159 Qmux(model,inputs,DiagnosticAnalysisEnum,NoneAnalysisEnum); 160 160 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 161 161 #else -
issm/trunk/src/c/parallel/diagnostic_core.cpp
r3529 r3567 70 70 71 71 /*recover fem models: */ 72 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());73 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum (),VertAnalysisEnum());74 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());75 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());76 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum ());72 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 73 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum); 74 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 75 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 76 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum); 77 77 78 78 //specific parameters for specific models … … 83 83 //for qmu analysis, be sure the velocity input we are starting from is the one in the parameters: */ 84 84 if(qmu_analysis){ 85 model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum (),HorizAnalysisEnum());85 model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum,HorizAnalysisEnum); 86 86 inputs->Add("velocity",u_g_initial,3,numberofnodes); 87 87 } … … 90 90 91 91 if(verbose)_printf_("%s\n","computing surface slope (x and y derivatives)..."); 92 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum (),SurfaceXAnalysisEnum());93 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum (),SurfaceYAnalysisEnum());92 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceXAnalysisEnum); 93 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceYAnalysisEnum); 94 94 95 95 if (dim==3){ … … 106 106 107 107 if(verbose)_printf_("%s\n"," computing hutter velocities..."); 108 diagnostic_core_linear(&ug,fem_dhu,inputs,DiagnosticAnalysisEnum (),HutterAnalysisEnum());108 diagnostic_core_linear(&ug,fem_dhu,inputs,DiagnosticAnalysisEnum,HutterAnalysisEnum); 109 109 110 110 if(verbose)_printf_("%s\n"," computing pressure according to MacAyeal..."); 111 ComputePressurex(&pg,fem_dhu->elements,fem_dhu->nodes, fem_dhu->vertices,fem_dhu->loads,fem_dhu->materials,fem_dhu->parameters,inputs,DiagnosticAnalysisEnum (),HutterAnalysisEnum());111 ComputePressurex(&pg,fem_dhu->elements,fem_dhu->nodes, fem_dhu->vertices,fem_dhu->loads,fem_dhu->materials,fem_dhu->parameters,inputs,DiagnosticAnalysisEnum,HutterAnalysisEnum); 112 112 113 113 if(verbose)_printf_("%s\n"," update boundary conditions for macyeal pattyn using hutter results..."); … … 123 123 124 124 if(verbose)_printf_("%s\n"," computing horizontal velocities..."); 125 diagnostic_core_nonlinear(&ug,NULL,NULL,fem_dh->loads,fem_dh,inputs,DiagnosticAnalysisEnum (),HorizAnalysisEnum());125 diagnostic_core_nonlinear(&ug,NULL,NULL,fem_dh->loads,fem_dh,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum); 126 126 127 127 if(dim==2){ 128 128 if(verbose)_printf_("%s\n"," computing pressure according to MacAyeal..."); 129 ComputePressurex(&pg,fem_dh->elements,fem_dh->nodes, fem_dh->vertices,fem_dh->loads,fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum (),HorizAnalysisEnum());129 ComputePressurex(&pg,fem_dh->elements,fem_dh->nodes, fem_dh->vertices,fem_dh->loads,fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum); 130 130 } 131 131 … … 140 140 if(verbose)_printf_("%s\n"," computing vertical velocities..."); 141 141 inputs->Add("velocity",ug_horiz,numberofdofspernode_dh,numberofnodes); 142 diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum (),VertAnalysisEnum());142 diagnostic_core_linear(&ug_vert,fem_dv,inputs,DiagnosticAnalysisEnum,VertAnalysisEnum); 143 143 144 144 if(verbose)_printf_("%s\n"," combining horizontal and vertical velocities..."); … … 150 150 151 151 if(verbose)_printf_("%s\n"," computing pressure according to Pattyn..."); 152 ComputePressurex(&pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads, fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum (),HorizAnalysisEnum());152 ComputePressurex(&pg,fem_dh->elements, fem_dh->nodes, fem_dh->vertices,fem_dh->loads, fem_dh->materials,fem_dh->parameters,inputs,DiagnosticAnalysisEnum,HorizAnalysisEnum); 153 153 154 154 if (isstokes){ … … 158 158 159 159 if(verbose)_printf_("%s\n","computing bed slope (x and y derivatives)..."); 160 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum (),BedXAnalysisEnum());161 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum (),BedYAnalysisEnum());160 diagnostic_core_linear(&slopex,fem_sl,inputs,SlopecomputeAnalysisEnum,BedXAnalysisEnum); 161 diagnostic_core_linear(&slopey,fem_sl,inputs,SlopecomputeAnalysisEnum,BedYAnalysisEnum); 162 162 FieldExtrudex( slopex, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopex",0); 163 163 FieldExtrudex( slopey, fem_sl->elements,fem_sl->nodes,fem_sl->vertices,fem_sl->loads,fem_sl->materials,fem_sl->parameters,"slopey",0); … … 182 182 if(verbose)_printf_("%s\n"," computing stokes velocities and pressure ..."); 183 183 VecFree(&ug); 184 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum (),StokesAnalysisEnum());184 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum,StokesAnalysisEnum); 185 185 186 186 //decondition" pressure -
issm/trunk/src/c/parallel/gradjcompute_core.cpp
r3446 r3567 62 62 63 63 _printf_("%s\n"," recover solution for this stiffness and right hand side:"); 64 diagnostic_core_nonlinear(&u_g,&K_ff0,&K_fs0,NULL, femmodel,inputs,DiagnosticAnalysisEnum (),sub_analysis_type);64 diagnostic_core_nonlinear(&u_g,&K_ff0,&K_fs0,NULL, femmodel,inputs,DiagnosticAnalysisEnum,sub_analysis_type); 65 65 VecToMPISerial(&u_g_double,u_g); VecFree(&u_g); 66 66 inputs->Add("velocity",u_g_double,numberofdofspernode,numberofnodes); … … 95 95 96 96 if(control_steady){ 97 diagnostic_results=new DataSet(ResultsEnum ());97 diagnostic_results=new DataSet(ResultsEnum); 98 98 diagnostic_core(diagnostic_results,model, inputs); 99 99 -
issm/trunk/src/c/parallel/objectivefunctionC.cpp
r3563 r3567 90 90 /*Run diagnostic with updated parameters.*/ 91 91 if(!control_steady){ 92 diagnostic_core_nonlinear(&u_g,NULL,NULL,NULL,femmodel,inputs,DiagnosticAnalysisEnum (),sub_analysis_type);92 diagnostic_core_nonlinear(&u_g,NULL,NULL,NULL,femmodel,inputs,DiagnosticAnalysisEnum,sub_analysis_type); 93 93 VecToMPISerial(&u_g_double,u_g); VecFree(&u_g); 94 94 inputs->Add("velocity",u_g_double,numberofdofspernode,numberofnodes); … … 96 96 else{ 97 97 //We need a 3D velocity!! (vz is required for the next thermal run) 98 diagnostic_results=new DataSet(ResultsEnum ());98 diagnostic_results=new DataSet(ResultsEnum); 99 99 diagnostic_core(diagnostic_results,model, inputs); 100 100 -
issm/trunk/src/c/parallel/prognostic.cpp
r3381 r3567 74 74 75 75 _printf_("read and create finite element model:\n"); 76 model->AddFormulation(fid,PrognosticAnalysisEnum ());76 model->AddFormulation(fid,PrognosticAnalysisEnum); 77 77 78 78 /*recover parameters: */ … … 82 82 _printf_("initialize inputs:\n"); 83 83 84 model->FindParam(&u_g_serial,NULL,NULL,"u_g",PrognosticAnalysisEnum ());85 model->FindParam(&h_g_initial,NULL,NULL,"h_g",PrognosticAnalysisEnum ());86 model->FindParam(&melting_g,NULL,NULL,"m_g",PrognosticAnalysisEnum ());87 model->FindParam(&accumulation_g,NULL,NULL,"a_g",PrognosticAnalysisEnum ());84 model->FindParam(&u_g_serial,NULL,NULL,"u_g",PrognosticAnalysisEnum); 85 model->FindParam(&h_g_initial,NULL,NULL,"h_g",PrognosticAnalysisEnum); 86 model->FindParam(&melting_g,NULL,NULL,"m_g",PrognosticAnalysisEnum); 87 model->FindParam(&accumulation_g,NULL,NULL,"a_g",PrognosticAnalysisEnum); 88 88 model->FindParam(&dt,"dt"); 89 89 model->FindParam(&yts,"yts"); … … 98 98 99 99 _printf_("initialize results:\n"); 100 results=new DataSet(ResultsEnum ());100 results=new DataSet(ResultsEnum); 101 101 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 102 102 … … 118 118 #ifdef _HAVE_DAKOTA_ 119 119 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 120 Qmux(model,inputs,PrognosticAnalysisEnum (),NoneAnalysisEnum());120 Qmux(model,inputs,PrognosticAnalysisEnum,NoneAnalysisEnum); 121 121 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 122 122 #else … … 130 130 131 131 _printf_("process results:\n"); 132 ProcessResults(&processedresults,results,model,PrognosticAnalysisEnum ());132 ProcessResults(&processedresults,results,model,PrognosticAnalysisEnum); 133 133 134 134 _printf_("write results to disk:\n"); -
issm/trunk/src/c/parallel/prognostic2.cpp
r3381 r3567 74 74 75 75 _printf_("read and create finite element model:\n"); 76 model->AddFormulation(fid,Prognostic2AnalysisEnum ());76 model->AddFormulation(fid,Prognostic2AnalysisEnum); 77 77 78 78 /*recover parameters: */ … … 82 82 _printf_("initialize inputs:\n"); 83 83 84 model->FindParam(&vx_g,NULL,NULL,"vx_g",Prognostic2AnalysisEnum ());85 model->FindParam(&vy_g,NULL,NULL,"vy_g",Prognostic2AnalysisEnum ());86 model->FindParam(&h_g_initial,NULL,NULL,"h_g",Prognostic2AnalysisEnum ());87 model->FindParam(&m_g,NULL,NULL,"m_g",Prognostic2AnalysisEnum ());88 model->FindParam(&a_g,NULL,NULL,"a_g",Prognostic2AnalysisEnum ());84 model->FindParam(&vx_g,NULL,NULL,"vx_g",Prognostic2AnalysisEnum); 85 model->FindParam(&vy_g,NULL,NULL,"vy_g",Prognostic2AnalysisEnum); 86 model->FindParam(&h_g_initial,NULL,NULL,"h_g",Prognostic2AnalysisEnum); 87 model->FindParam(&m_g,NULL,NULL,"m_g",Prognostic2AnalysisEnum); 88 model->FindParam(&a_g,NULL,NULL,"a_g",Prognostic2AnalysisEnum); 89 89 model->FindParam(&dt,"dt"); 90 90 model->FindParam(&yts,"yts"); … … 100 100 101 101 _printf_("initialize results:\n"); 102 results=new DataSet(ResultsEnum ());102 results=new DataSet(ResultsEnum); 103 103 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 104 104 … … 120 120 #ifdef _HAVE_DAKOTA_ 121 121 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 122 Qmux(model,inputs,Prognostic2AnalysisEnum (),NoneAnalysisEnum());122 Qmux(model,inputs,Prognostic2AnalysisEnum,NoneAnalysisEnum); 123 123 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 124 124 #else … … 132 132 133 133 _printf_("process results:\n"); 134 ProcessResults(&processedresults,results,model,Prognostic2AnalysisEnum ());134 ProcessResults(&processedresults,results,model,Prognostic2AnalysisEnum); 135 135 136 136 _printf_("write results to disk:\n"); -
issm/trunk/src/c/parallel/prognostic2_core.cpp
r3556 r3567 35 35 36 36 /*recover fem model: */ 37 fem_p=model->GetFormulation(Prognostic2AnalysisEnum ());37 fem_p=model->GetFormulation(Prognostic2AnalysisEnum); 38 38 39 39 //first recover parameters common to all solutions … … 54 54 55 55 _printf_("call computational core:\n"); 56 diagnostic_core_linear(&h_g,fem_p,inputs,Prognostic2AnalysisEnum (),NoneAnalysisEnum());56 diagnostic_core_linear(&h_g,fem_p,inputs,Prognostic2AnalysisEnum,NoneAnalysisEnum); 57 57 58 58 _printf_("Averaging over vertices:\n"); -
issm/trunk/src/c/parallel/prognostic_core.cpp
r3446 r3567 34 34 35 35 /*recover fem model: */ 36 fem_p=model->GetFormulation(PrognosticAnalysisEnum ());36 fem_p=model->GetFormulation(PrognosticAnalysisEnum); 37 37 38 38 //first recover parameters common to all solutions … … 47 47 48 48 _printf_("call computational core:\n"); 49 diagnostic_core_linear(&h_g,fem_p,inputs,PrognosticAnalysisEnum (),NoneAnalysisEnum());49 diagnostic_core_linear(&h_g,fem_p,inputs,PrognosticAnalysisEnum,NoneAnalysisEnum); 50 50 51 51 _printf_("extrude computed thickness on all layers:\n"); -
issm/trunk/src/c/parallel/slopecompute.cpp
r3332 r3567 68 68 69 69 _printf_("read and create finite element model:\n"); 70 model->AddFormulation(fid,SlopecomputeAnalysisEnum ());70 model->AddFormulation(fid,SlopecomputeAnalysisEnum); 71 71 72 72 /*recover parameters: */ … … 80 80 81 81 _printf_("initialize results:\n"); 82 results=new DataSet(ResultsEnum ());82 results=new DataSet(ResultsEnum); 83 83 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 84 84 … … 94 94 95 95 _printf_("process results:\n"); 96 ProcessResults(&processedresults,results,model,SlopecomputeAnalysisEnum ());96 ProcessResults(&processedresults,results,model,SlopecomputeAnalysisEnum); 97 97 98 98 _printf_("write results to disk:\n"); -
issm/trunk/src/c/parallel/slopecompute_core.cpp
r3446 r3567 31 31 32 32 /*recover fem model: */ 33 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum ());33 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum); 34 34 35 35 //first recover parameters common to all solutions … … 39 39 40 40 _printf_("call computational core:\n"); 41 diagnostic_core_linear(&sx_g,fem_sl,inputs,SlopecomputeAnalysisEnum (),SurfaceXAnalysisEnum());42 diagnostic_core_linear(&sy_g,fem_sl,inputs,SlopecomputeAnalysisEnum (),SurfaceYAnalysisEnum());41 diagnostic_core_linear(&sx_g,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceXAnalysisEnum); 42 diagnostic_core_linear(&sy_g,fem_sl,inputs,SlopecomputeAnalysisEnum,SurfaceYAnalysisEnum); 43 43 44 44 _printf_("extrude computed slope on all layers:\n"); -
issm/trunk/src/c/parallel/steadystate.cpp
r3332 r3567 79 79 80 80 _printf_("\n reading diagnostic horiz model data:\n"); 81 model->AddFormulation(fid,DiagnosticAnalysisEnum (),HorizAnalysisEnum());81 model->AddFormulation(fid,DiagnosticAnalysisEnum,HorizAnalysisEnum); 82 82 83 83 _printf_("\n reading diagnostic vert model data:\n"); 84 model->AddFormulation(fid,DiagnosticAnalysisEnum (),VertAnalysisEnum());84 model->AddFormulation(fid,DiagnosticAnalysisEnum,VertAnalysisEnum); 85 85 86 86 _printf_("\n reading diagnostic stokes model data:\n"); 87 model->AddFormulation(fid,DiagnosticAnalysisEnum (),StokesAnalysisEnum());87 model->AddFormulation(fid,DiagnosticAnalysisEnum,StokesAnalysisEnum); 88 88 89 89 _printf_("\n reading diagnostic hutter model data:\n"); 90 model->AddFormulation(fid,DiagnosticAnalysisEnum (),HutterAnalysisEnum());90 model->AddFormulation(fid,DiagnosticAnalysisEnum,HutterAnalysisEnum); 91 91 92 92 _printf_("\n reading surface and bed slope computation model data:\n"); 93 model->AddFormulation(fid,SlopecomputeAnalysisEnum ());93 model->AddFormulation(fid,SlopecomputeAnalysisEnum); 94 94 95 95 _printf_("\n read and create thermal finite element model:\n"); 96 model->AddFormulation(fid,ThermalAnalysisEnum ());96 model->AddFormulation(fid,ThermalAnalysisEnum); 97 97 _printf_("\n read and create melting finite element model:\n"); 98 model->AddFormulation(fid,MeltingAnalysisEnum ());98 model->AddFormulation(fid,MeltingAnalysisEnum); 99 99 100 100 /*recover parameters: */ … … 104 104 105 105 _printf_("initialize inputs:\n"); 106 model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum (),HorizAnalysisEnum());107 model->FindParam(&p_g_initial,NULL,NULL,"p_g",ThermalAnalysisEnum ());106 model->FindParam(&u_g_initial,NULL,NULL,"u_g",DiagnosticAnalysisEnum,HorizAnalysisEnum); 107 model->FindParam(&p_g_initial,NULL,NULL,"p_g",ThermalAnalysisEnum); 108 108 model->FindParam(&dt,"dt"); 109 109 model->FindParam(&numberofnodes,"numberofnodes"); … … 115 115 116 116 if(control_analysis){ 117 model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum (),HorizAnalysisEnum());118 model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum (),HorizAnalysisEnum());117 model->FindParam(&u_g_obs,NULL,NULL,"u_g_obs",DiagnosticAnalysisEnum,HorizAnalysisEnum); 118 model->FindParam(&weights,NULL,NULL,"weights",DiagnosticAnalysisEnum,HorizAnalysisEnum); 119 119 inputs->Add("velocity_obs",u_g_obs,2,numberofnodes); 120 120 inputs->Add("weights",weights,1,numberofnodes); … … 122 122 123 123 _printf_("initialize results:\n"); 124 results=new DataSet(ResultsEnum ());124 results=new DataSet(ResultsEnum); 125 125 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 126 126 … … 140 140 141 141 _printf_("process results:\n"); 142 ProcessResults(&processed_results,results,model,SteadystateAnalysisEnum ());142 ProcessResults(&processed_results,results,model,SteadystateAnalysisEnum); 143 143 } 144 144 else{ 145 145 /*change control_steady to 1 to know we are doing steadystate*/ 146 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());147 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());146 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 147 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 148 148 param=(Param*)fem_dh->parameters->FindParamObject("control_steady"); 149 149 param->SetDouble(1); … … 165 165 166 166 _printf_("process results:\n"); 167 ProcessResults(&processed_results,results,model,ControlAnalysisEnum ());167 ProcessResults(&processed_results,results,model,ControlAnalysisEnum); 168 168 } 169 169 … … 177 177 #ifdef _HAVE_DAKOTA_ 178 178 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 179 Qmux(model,inputs,SteadystateAnalysisEnum (),NoneAnalysisEnum());179 Qmux(model,inputs,SteadystateAnalysisEnum,NoneAnalysisEnum); 180 180 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 181 181 #else -
issm/trunk/src/c/parallel/steadystate_core.cpp
r3446 r3567 51 51 52 52 /*recover fem models: */ 53 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());54 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum (),VertAnalysisEnum());55 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());56 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());57 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum ());58 fem_t=model->GetFormulation(ThermalAnalysisEnum ());59 fem_m=model->GetFormulation(MeltingAnalysisEnum ());53 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 54 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum); 55 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 56 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 57 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum); 58 fem_t=model->GetFormulation(ThermalAnalysisEnum); 59 fem_m=model->GetFormulation(MeltingAnalysisEnum); 60 60 61 61 … … 81 81 inputs->Add("velocity",u_g,ndof,numberofnodes); 82 82 } 83 results_thermal=new DataSet(ResultsEnum ());83 results_thermal=new DataSet(ResultsEnum); 84 84 thermal_core(results_thermal,model,inputs); 85 85 … … 98 98 99 99 //now compute diagnostic velocity using the steady state temperature. 100 results_diagnostic=new DataSet(ResultsEnum ());100 results_diagnostic=new DataSet(ResultsEnum); 101 101 diagnostic_core(results_diagnostic,model, inputs); 102 102 -
issm/trunk/src/c/parallel/thermal.cpp
r3332 r3567 73 73 74 74 _printf_("read and create thermal finite element model:\n"); 75 model->AddFormulation(fid,ThermalAnalysisEnum ());75 model->AddFormulation(fid,ThermalAnalysisEnum); 76 76 _printf_("read and create melting finite element model:\n"); 77 model->AddFormulation(fid,MeltingAnalysisEnum ());77 model->AddFormulation(fid,MeltingAnalysisEnum); 78 78 79 79 /*recover parameters: */ … … 82 82 83 83 _printf_("initialize inputs:\n"); 84 model->FindParam(&u_g,NULL,NULL,"u_g",ThermalAnalysisEnum ());85 model->FindParam(&p_g,NULL,NULL,"p_g",ThermalAnalysisEnum ());84 model->FindParam(&u_g,NULL,NULL,"u_g",ThermalAnalysisEnum); 85 model->FindParam(&p_g,NULL,NULL,"p_g",ThermalAnalysisEnum); 86 86 model->FindParam(&numberofnodes,"numberofnodes"); 87 87 model->FindParam(&dt,"dt"); … … 95 95 96 96 _printf_("initialize results:\n"); 97 results=new DataSet(ResultsEnum ());97 results=new DataSet(ResultsEnum); 98 98 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 99 99 … … 112 112 113 113 _printf_("process results:\n"); 114 ProcessResults(&processed_results,results,model,ThermalAnalysisEnum ());114 ProcessResults(&processed_results,results,model,ThermalAnalysisEnum); 115 115 116 116 _printf_("write results to disk:\n"); … … 123 123 #ifdef _HAVE_DAKOTA_ 124 124 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 125 Qmux(model,inputs,ThermalAnalysisEnum (),NoneAnalysisEnum());125 Qmux(model,inputs,ThermalAnalysisEnum,NoneAnalysisEnum); 126 126 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 127 127 #else -
issm/trunk/src/c/parallel/thermal_core.cpp
r3332 r3567 45 45 46 46 /*recover fem models: */ 47 fem_t=model->GetFormulation(ThermalAnalysisEnum ());48 fem_m=model->GetFormulation(MeltingAnalysisEnum ());47 fem_t=model->GetFormulation(ThermalAnalysisEnum); 48 fem_m=model->GetFormulation(MeltingAnalysisEnum); 49 49 50 50 //first recover parameters common to all solutions … … 65 65 66 66 if(verbose)_printf_("computing temperatures:\n"); 67 thermal_core_nonlinear(&t_g[0],&melting_offset,fem_t,inputs,ThermalAnalysisEnum (),NoneAnalysisEnum());67 thermal_core_nonlinear(&t_g[0],&melting_offset,fem_t,inputs,ThermalAnalysisEnum,NoneAnalysisEnum); 68 68 inputs->Add("temperature",t_g[0],1,numberofnodes); 69 69 inputs->Add("melting_offset",melting_offset); 70 70 71 71 if(verbose)_printf_("computing melting:\n"); 72 diagnostic_core_linear(&m_g[0],fem_m,inputs,MeltingAnalysisEnum (),NoneAnalysisEnum());72 diagnostic_core_linear(&m_g[0],fem_m,inputs,MeltingAnalysisEnum,NoneAnalysisEnum); 73 73 } 74 74 else{ … … 96 96 if(verbose)_printf_("computing temperatures:\n"); 97 97 inputs->Add("temperature",t_g[i],1,numberofnodes); 98 thermal_core_nonlinear(&t_g[i+1],&melting_offset,fem_t,inputs,ThermalAnalysisEnum (),NoneAnalysisEnum());98 thermal_core_nonlinear(&t_g[i+1],&melting_offset,fem_t,inputs,ThermalAnalysisEnum,NoneAnalysisEnum); 99 99 100 100 if(verbose)_printf_("computing melting:\n"); 101 101 inputs->Add("temperature",t_g[i+1],1,numberofnodes); 102 102 inputs->Add("melting_offset",melting_offset); 103 diagnostic_core_linear(&m_g[i+1],fem_m,inputs,MeltingAnalysisEnum (),NoneAnalysisEnum());103 diagnostic_core_linear(&m_g[i+1],fem_m,inputs,MeltingAnalysisEnum,NoneAnalysisEnum); 104 104 } 105 105 } -
issm/trunk/src/c/parallel/transient.cpp
r3332 r3567 76 76 _printf_("read and create finite element model:\n"); 77 77 _printf_("\n reading diagnostic horiz model data:\n"); 78 model->AddFormulation(fid,DiagnosticAnalysisEnum (),HorizAnalysisEnum());78 model->AddFormulation(fid,DiagnosticAnalysisEnum,HorizAnalysisEnum); 79 79 80 80 _printf_("\n reading diagnostic vert model data:\n"); 81 model->AddFormulation(fid,DiagnosticAnalysisEnum (),VertAnalysisEnum());81 model->AddFormulation(fid,DiagnosticAnalysisEnum,VertAnalysisEnum); 82 82 83 83 _printf_("\n reading diagnostic stokes model data:\n"); 84 model->AddFormulation(fid,DiagnosticAnalysisEnum (),StokesAnalysisEnum());84 model->AddFormulation(fid,DiagnosticAnalysisEnum,StokesAnalysisEnum); 85 85 86 86 _printf_("\n reading diagnostic hutter model data:\n"); 87 model->AddFormulation(fid,DiagnosticAnalysisEnum (),HutterAnalysisEnum());87 model->AddFormulation(fid,DiagnosticAnalysisEnum,HutterAnalysisEnum); 88 88 89 89 _printf_("\n reading surface and bed slope computation model data:\n"); 90 model->AddFormulation(fid,SlopecomputeAnalysisEnum ());90 model->AddFormulation(fid,SlopecomputeAnalysisEnum); 91 91 92 92 _printf_("\n reading prognositc model data:\n"); 93 model->AddFormulation(fid,PrognosticAnalysisEnum ());93 model->AddFormulation(fid,PrognosticAnalysisEnum); 94 94 95 95 /*Do we run in 3d?, in which case we need thermal and melting also:*/ … … 97 97 if(dim==3){ 98 98 _printf_("read and create thermal finite element model:\n"); 99 model->AddFormulation(fid,ThermalAnalysisEnum (),TransientAnalysisEnum());99 model->AddFormulation(fid,ThermalAnalysisEnum,TransientAnalysisEnum); 100 100 _printf_("read and create melting finite element model:\n"); 101 model->AddFormulation(fid,MeltingAnalysisEnum (),TransientAnalysisEnum());101 model->AddFormulation(fid,MeltingAnalysisEnum,TransientAnalysisEnum); 102 102 } 103 103 … … 107 107 108 108 _printf_("initialize inputs:\n"); 109 model->FindParam(&u_g,NULL,NULL,"u_g",PrognosticAnalysisEnum ());110 model->FindParam(&m_g,NULL,NULL,"m_g",PrognosticAnalysisEnum ());111 model->FindParam(&a_g,NULL,NULL,"a_g",PrognosticAnalysisEnum ());109 model->FindParam(&u_g,NULL,NULL,"u_g",PrognosticAnalysisEnum); 110 model->FindParam(&m_g,NULL,NULL,"m_g",PrognosticAnalysisEnum); 111 model->FindParam(&a_g,NULL,NULL,"a_g",PrognosticAnalysisEnum); 112 112 model->FindParam(&numberofnodes,"numberofnodes"); 113 113 model->FindParam(&dt,"dt"); … … 122 122 123 123 _printf_("initialize results:\n"); 124 results=new DataSet(ResultsEnum ());124 results=new DataSet(ResultsEnum); 125 125 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); 126 126 … … 139 139 140 140 _printf_("process results:\n"); 141 ProcessResults(&processed_results,results,model,TransientAnalysisEnum ());141 ProcessResults(&processed_results,results,model,TransientAnalysisEnum); 142 142 143 143 _printf_("write results to disk:\n"); … … 150 150 #ifdef _HAVE_DAKOTA_ 151 151 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); 152 Qmux(model,inputs,TransientAnalysisEnum (),NoneAnalysisEnum());152 Qmux(model,inputs,TransientAnalysisEnum,NoneAnalysisEnum); 153 153 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); 154 154 #else -
issm/trunk/src/c/parallel/transient_core_2d.cpp
r3446 r3567 63 63 64 64 /*recover fem models: */ 65 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());66 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum (),VertAnalysisEnum());67 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());68 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());69 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum ());70 fem_p=model->GetFormulation(PrognosticAnalysisEnum ());65 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 66 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum); 67 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 68 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 69 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum); 70 fem_p=model->GetFormulation(PrognosticAnalysisEnum); 71 71 72 72 … … 119 119 //Deal with velocities. 120 120 _printf_("%s\n","computing new velocity"); 121 diagnostic_results=new DataSet(ResultsEnum ());121 diagnostic_results=new DataSet(ResultsEnum); 122 122 diagnostic_core(diagnostic_results,model, inputs); 123 123 … … 131 131 132 132 inputs->Add("velocity",u_g,2,numberofnodes); 133 prognostic_results=new DataSet(ResultsEnum ());133 prognostic_results=new DataSet(ResultsEnum); 134 134 prognostic_core(prognostic_results,model, inputs); 135 135 -
issm/trunk/src/c/parallel/transient_core_3d.cpp
r3544 r3567 67 67 68 68 /*recover fem models: */ 69 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum (),HorizAnalysisEnum());70 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum (),VertAnalysisEnum());71 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum (),StokesAnalysisEnum());72 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum (),HutterAnalysisEnum());73 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum ());74 fem_p=model->GetFormulation(PrognosticAnalysisEnum ());75 fem_t=model->GetFormulation(ThermalAnalysisEnum ());76 fem_m=model->GetFormulation(MeltingAnalysisEnum ());69 fem_dh=model->GetFormulation(DiagnosticAnalysisEnum,HorizAnalysisEnum); 70 fem_dv=model->GetFormulation(DiagnosticAnalysisEnum,VertAnalysisEnum); 71 fem_ds=model->GetFormulation(DiagnosticAnalysisEnum,StokesAnalysisEnum); 72 fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum,HutterAnalysisEnum); 73 fem_sl=model->GetFormulation(SlopecomputeAnalysisEnum); 74 fem_p=model->GetFormulation(PrognosticAnalysisEnum); 75 fem_t=model->GetFormulation(ThermalAnalysisEnum); 76 fem_m=model->GetFormulation(MeltingAnalysisEnum); 77 77 78 78 … … 134 134 //Deal with temperature first 135 135 if(verbose)_printf_("%s\n","computing temperature"); 136 thermal_core_nonlinear(&t_g,&melting_offset,fem_t,inputs,ThermalAnalysisEnum (),TransientAnalysisEnum());136 thermal_core_nonlinear(&t_g,&melting_offset,fem_t,inputs,ThermalAnalysisEnum,TransientAnalysisEnum); 137 137 if(verbose)_printf_("%s\n","computing melting"); 138 138 inputs->Add("temperature",t_g,1,numberofnodes); 139 139 inputs->Add("melting_offset",melting_offset); 140 diagnostic_core_linear(&m_g,fem_m,inputs,MeltingAnalysisEnum (),TransientAnalysisEnum());140 diagnostic_core_linear(&m_g,fem_m,inputs,MeltingAnalysisEnum,TransientAnalysisEnum); 141 141 142 142 //Compute depth averaged temperature and add to inputs … … 148 148 149 149 //Deal with velocities. 150 diagnostic_results=new DataSet(ResultsEnum ());150 diagnostic_results=new DataSet(ResultsEnum); 151 151 diagnostic_core(diagnostic_results,model, inputs); 152 152 … … 160 160 161 161 inputs->Add("velocity",u_g,3,numberofnodes); 162 prognostic_results=new DataSet(ResultsEnum ()); prognostic_core(prognostic_results,model, inputs);162 prognostic_results=new DataSet(ResultsEnum); prognostic_core(prognostic_results,model, inputs); 163 163 164 164 //extract h_g prognostic_results, and erase prognostic_results; … … 198 198 199 199 model->FindParam(&outputfilename,"outputfilename"); 200 ProcessResults(&processed_results,results,model,TransientAnalysisEnum ());200 ProcessResults(&processed_results,results,model,TransientAnalysisEnum); 201 201 OutputResults(processed_results,outputfilename); 202 202 -
issm/trunk/src/c/shared/Dofs/DistributeNumDofs.cpp
r3354 r3567 13 13 14 14 /*ok, according to analysis type: */ 15 if (analysis_type==ControlAnalysisEnum ()){16 if (sub_analysis_type==HorizAnalysisEnum ()){15 if (analysis_type==ControlAnalysisEnum){ 16 if (sub_analysis_type==HorizAnalysisEnum){ 17 17 numdofs=2; 18 18 } 19 else if (sub_analysis_type==VertAnalysisEnum ()){19 else if (sub_analysis_type==VertAnalysisEnum){ 20 20 numdofs=1; 21 21 } 22 else if (sub_analysis_type==StokesAnalysisEnum ()){22 else if (sub_analysis_type==StokesAnalysisEnum){ 23 23 numdofs=4; 24 24 } 25 else if (sub_analysis_type==HutterAnalysisEnum ()){25 else if (sub_analysis_type==HutterAnalysisEnum){ 26 26 numdofs=2; 27 27 } 28 28 } 29 else if (analysis_type==DiagnosticAnalysisEnum ()){30 if (sub_analysis_type==HorizAnalysisEnum ()){29 else if (analysis_type==DiagnosticAnalysisEnum){ 30 if (sub_analysis_type==HorizAnalysisEnum){ 31 31 numdofs=2; 32 32 } 33 else if (sub_analysis_type==VertAnalysisEnum ()){33 else if (sub_analysis_type==VertAnalysisEnum){ 34 34 numdofs=1; 35 35 } 36 else if (sub_analysis_type==StokesAnalysisEnum ()){36 else if (sub_analysis_type==StokesAnalysisEnum){ 37 37 numdofs=4; 38 38 } 39 else if (sub_analysis_type==HutterAnalysisEnum ()){39 else if (sub_analysis_type==HutterAnalysisEnum){ 40 40 numdofs=2; 41 41 } 42 42 } 43 else if (analysis_type==SlopecomputeAnalysisEnum ()){43 else if (analysis_type==SlopecomputeAnalysisEnum){ 44 44 numdofs=1; 45 45 } 46 else if (analysis_type==ThermalAnalysisEnum ()){46 else if (analysis_type==ThermalAnalysisEnum){ 47 47 numdofs=1; 48 48 } 49 else if (analysis_type==MeltingAnalysisEnum ()){49 else if (analysis_type==MeltingAnalysisEnum){ 50 50 numdofs=1; 51 51 } 52 else if (analysis_type==PrognosticAnalysisEnum ()){52 else if (analysis_type==PrognosticAnalysisEnum){ 53 53 numdofs=1; 54 54 } 55 else if (analysis_type==Prognostic2AnalysisEnum ()){55 else if (analysis_type==Prognostic2AnalysisEnum){ 56 56 numdofs=1; 57 57 } 58 else if (analysis_type==BalancedthicknessAnalysisEnum ()){58 else if (analysis_type==BalancedthicknessAnalysisEnum){ 59 59 numdofs=1; 60 60 } 61 else if (analysis_type==BalancedvelocitiesAnalysisEnum ()){61 else if (analysis_type==BalancedvelocitiesAnalysisEnum){ 62 62 numdofs=1; 63 63 } -
issm/trunk/src/m/enum/AdjointAnalysisEnum.m
r1714 r3567 7 7 % macro=AdjointAnalysisEnum() 8 8 9 macro= 221;9 macro=16; -
issm/trunk/src/m/enum/AirEnum.m
r1749 r3567 7 7 % macro=AirEnum() 8 8 9 macro= 603;9 macro=75; -
issm/trunk/src/m/enum/AnalysisEnum.m
r1714 r3567 7 7 % macro=AnalysisEnum() 8 8 9 macro= 200;9 macro=9; -
issm/trunk/src/m/enum/AnalysisTypeFromEnum.m
r3362 r3567 11 11 %initialize output 12 12 string='not found'; 13 14 13 15 14 if enum==AnalysisEnum(), -
issm/trunk/src/m/enum/BalancedthicknessAnalysisEnum.m
r3354 r3567 7 7 % macro=BalancedthicknessAnalysisEnum() 8 8 9 macro=2 52;9 macro=29; -
issm/trunk/src/m/enum/BalancedvelocitiesAnalysisEnum.m
r3354 r3567 7 7 % macro=BalancedvelocitiesAnalysisEnum() 8 8 9 macro= 253;9 macro=30; -
issm/trunk/src/m/enum/BeamEnum.m
r3417 r3567 7 7 % macro=BeamEnum() 8 8 9 macro=4 16;9 macro=49; -
issm/trunk/src/m/enum/BedXAnalysisEnum.m
r1714 r3567 7 7 % macro=BedXAnalysisEnum() 8 8 9 macro=2 43;9 macro=25; -
issm/trunk/src/m/enum/BedYAnalysisEnum.m
r1714 r3567 7 7 % macro=BedYAnalysisEnum() 8 8 9 macro=2 44;9 macro=26; -
issm/trunk/src/m/enum/ConstraintsEnum.m
r1714 r3567 7 7 % macro=ConstraintsEnum() 8 8 9 macro= 103;9 macro=3; -
issm/trunk/src/m/enum/ControlAnalysisEnum.m
r1714 r3567 7 7 % macro=ControlAnalysisEnum() 8 8 9 macro= 220;9 macro=15; -
issm/trunk/src/m/enum/DatasetsEnum.m
r1714 r3567 7 7 % macro=DatasetsEnum() 8 8 9 macro= 100;9 macro=0; -
issm/trunk/src/m/enum/DiagnosticAnalysisEnum.m
r1714 r3567 7 7 % macro=DiagnosticAnalysisEnum() 8 8 9 macro= 210;9 macro=10; -
issm/trunk/src/m/enum/DofIndexingEnum.m
r3417 r3567 7 7 % macro=DofIndexingEnum() 8 8 9 macro= 417;9 macro=51; -
issm/trunk/src/m/enum/DofVecEnum.m
r2311 r3567 7 7 % macro=DofVecEnum() 8 8 9 macro= 495;9 macro=69; -
issm/trunk/src/m/enum/ElementEnum.m
r1714 r3567 7 7 % macro=ElementEnum() 8 8 9 macro=4 10;9 macro=44; -
issm/trunk/src/m/enum/ElementPropertiesEnum.m
r3383 r3567 7 7 % macro=ElementPropertiesEnum() 8 8 9 macro=4 12;9 macro=46; -
issm/trunk/src/m/enum/ElementsEnum.m
r1714 r3567 7 7 % macro=ElementsEnum() 8 8 9 macro=1 01;9 macro=1; -
issm/trunk/src/m/enum/FormulationEnum.m
r1714 r3567 7 7 % macro=FormulationEnum() 8 8 9 macro=3 00;9 macro=36; -
issm/trunk/src/m/enum/GeographyEnum.m
r1729 r3567 7 7 % macro=GeographyEnum() 8 8 9 macro= 500;9 macro=70; -
issm/trunk/src/m/enum/GradientAnalysisEnum.m
r1714 r3567 7 7 % macro=GradientAnalysisEnum() 8 8 9 macro= 223;9 macro=18; -
issm/trunk/src/m/enum/HookEnum.m
r3383 r3567 7 7 % macro=HookEnum() 8 8 9 macro=4 09;9 macro=43; -
issm/trunk/src/m/enum/HorizAnalysisEnum.m
r1714 r3567 7 7 % macro=HorizAnalysisEnum() 8 8 9 macro= 211;9 macro=11; -
issm/trunk/src/m/enum/HutterAnalysisEnum.m
r1714 r3567 7 7 % macro=HutterAnalysisEnum() 8 8 9 macro= 213;9 macro=13; -
issm/trunk/src/m/enum/HutterFormulationEnum.m
r1714 r3567 7 7 % macro=HutterFormulationEnum() 8 8 9 macro=3 02;9 macro=38; -
issm/trunk/src/m/enum/IceEnum.m
r1749 r3567 7 7 % macro=IceEnum() 8 8 9 macro= 602;9 macro=74; -
issm/trunk/src/m/enum/IceSheetEnum.m
r1729 r3567 7 7 % macro=IceSheetEnum() 8 8 9 macro= 502;9 macro=71; -
issm/trunk/src/m/enum/IceShelfEnum.m
r1729 r3567 7 7 % macro=IceShelfEnum() 8 8 9 macro= 502;9 macro=72; -
issm/trunk/src/m/enum/IcefrontEnum.m
r1714 r3567 7 7 % macro=IcefrontEnum() 8 8 9 macro= 431;9 macro=55; -
issm/trunk/src/m/enum/InputEnum.m
r1714 r3567 7 7 % macro=InputEnum() 8 8 9 macro= 450;9 macro=64; -
issm/trunk/src/m/enum/InverseAnalysisEnum.m
r1714 r3567 7 7 % macro=InverseAnalysisEnum() 8 8 9 macro= 222;9 macro=17; -
issm/trunk/src/m/enum/LoadEnum.m
r1714 r3567 7 7 % macro=LoadEnum() 8 8 9 macro= 430;9 macro=54; -
issm/trunk/src/m/enum/LoadsEnum.m
r1714 r3567 7 7 % macro=LoadsEnum() 8 8 9 macro= 104;9 macro=4; -
issm/trunk/src/m/enum/MacAyealFormulationEnum.m
r1714 r3567 7 7 % macro=MacAyealFormulationEnum() 8 8 9 macro=3 03;9 macro=39; -
issm/trunk/src/m/enum/MaterialEnum.m
r1714 r3567 7 7 % macro=MaterialEnum() 8 8 9 macro= 440;9 macro=60; -
issm/trunk/src/m/enum/MaterialsEnum.m
r1714 r3567 7 7 % macro=MaterialsEnum() 8 8 9 macro= 105;9 macro=5; -
issm/trunk/src/m/enum/MaticeEnum.m
r1714 r3567 7 7 % macro=MaticeEnum() 8 8 9 macro= 441;9 macro=61; -
issm/trunk/src/m/enum/MatparEnum.m
r1714 r3567 7 7 % macro=MatparEnum() 8 8 9 macro= 442;9 macro=62; -
issm/trunk/src/m/enum/MelangeEnum.m
r1749 r3567 7 7 % macro=MelangeEnum() 8 8 9 macro= 604;9 macro=76; -
issm/trunk/src/m/enum/MeltingAnalysisEnum.m
r1714 r3567 7 7 % macro=MeltingAnalysisEnum() 8 8 9 macro= 260;9 macro=31; -
issm/trunk/src/m/enum/Mesh2gridAnalysisEnum.m
r1714 r3567 7 7 % macro=Mesh2gridAnalysisEnum() 8 8 9 macro= 270;9 macro=32; -
issm/trunk/src/m/enum/NodeEnum.m
r3417 r3567 7 7 % macro=NodeEnum() 8 8 9 macro= 421;9 macro=52; -
issm/trunk/src/m/enum/NodePropertiesEnum.m
r3417 r3567 7 7 % macro=NodePropertiesEnum() 8 8 9 macro=4 13;9 macro=47; -
issm/trunk/src/m/enum/NodesEnum.m
r1714 r3567 7 7 % macro=NodesEnum() 8 8 9 macro= 102;9 macro=2; -
issm/trunk/src/m/enum/NoneAnalysisEnum.m
r1714 r3567 7 7 % macro=NoneAnalysisEnum() 8 8 9 macro= 290;9 macro=35; -
issm/trunk/src/m/enum/NoneFormulationEnum.m
r1714 r3567 7 7 % macro=NoneFormulationEnum() 8 8 9 macro=3 01;9 macro=37; -
issm/trunk/src/m/enum/NumericalfluxEnum.m
r3359 r3567 7 7 % macro=NumericalfluxEnum() 8 8 9 macro= 435;9 macro=59; -
issm/trunk/src/m/enum/NumparEnum.m
r2333 r3567 7 7 % macro=NumparEnum() 8 8 9 macro= 443;9 macro=63; -
issm/trunk/src/m/enum/ObjectEnum.m
r1714 r3567 7 7 % macro=ObjectEnum() 8 8 9 macro=4 00;9 macro=42; -
issm/trunk/src/m/enum/ParamEnum.m
r1714 r3567 7 7 % macro=ParamEnum() 8 8 9 macro= 460;9 macro=65; -
issm/trunk/src/m/enum/ParametersAnalysisEnum.m
r1714 r3567 7 7 % macro=ParametersAnalysisEnum() 8 8 9 macro= 280;9 macro=33; -
issm/trunk/src/m/enum/ParametersEnum.m
r1714 r3567 7 7 % macro=ParametersEnum() 8 8 9 macro= 106;9 macro=6; -
issm/trunk/src/m/enum/PattynFormulationEnum.m
r1714 r3567 7 7 % macro=PattynFormulationEnum() 8 8 9 macro= 304;9 macro=40; -
issm/trunk/src/m/enum/PengridEnum.m
r1714 r3567 7 7 % macro=PengridEnum() 8 8 9 macro= 434;9 macro=58; -
issm/trunk/src/m/enum/PenpairEnum.m
r1714 r3567 7 7 % macro=PenpairEnum() 8 8 9 macro= 433;9 macro=57; -
issm/trunk/src/m/enum/PentaEnum.m
r3417 r3567 7 7 % macro=PentaEnum() 8 8 9 macro=4 14;9 macro=48; -
issm/trunk/src/m/enum/Prognostic2AnalysisEnum.m
r3354 r3567 7 7 % macro=Prognostic2AnalysisEnum() 8 8 9 macro=2 51;9 macro=28; -
issm/trunk/src/m/enum/PrognosticAnalysisEnum.m
r1714 r3567 7 7 % macro=PrognosticAnalysisEnum() 8 8 9 macro=2 50;9 macro=27; -
issm/trunk/src/m/enum/ResultEnum.m
r1714 r3567 7 7 % macro=ResultEnum() 8 8 9 macro= 470;9 macro=66; -
issm/trunk/src/m/enum/ResultsEnum.m
r1714 r3567 7 7 % macro=ResultsEnum() 8 8 9 macro= 107;9 macro=7; -
issm/trunk/src/m/enum/RgbEnum.m
r1714 r3567 7 7 % macro=RgbEnum() 8 8 9 macro= 480;9 macro=67; -
issm/trunk/src/m/enum/RiftfrontEnum.m
r1714 r3567 7 7 % macro=RiftfrontEnum() 8 8 9 macro= 432;9 macro=56; -
issm/trunk/src/m/enum/SingEnum.m
r3417 r3567 7 7 % macro=SingEnum() 8 8 9 macro= 415;9 macro=50; -
issm/trunk/src/m/enum/SlopecomputeAnalysisEnum.m
r3087 r3567 7 7 % macro=SlopecomputeAnalysisEnum() 8 8 9 macro=2 40;9 macro=22; -
issm/trunk/src/m/enum/SpcEnum.m
r1714 r3567 7 7 % macro=SpcEnum() 8 8 9 macro= 490;9 macro=68; -
issm/trunk/src/m/enum/SteadyAnalysisEnum.m
r3362 r3567 7 7 % macro=SteadyAnalysisEnum() 8 8 9 macro=2 32;9 macro=21; -
issm/trunk/src/m/enum/SteadystateAnalysisEnum.m
r1865 r3567 7 7 % macro=SteadystateAnalysisEnum() 8 8 9 macro= 281;9 macro=34; -
issm/trunk/src/m/enum/StokesAnalysisEnum.m
r1714 r3567 7 7 % macro=StokesAnalysisEnum() 8 8 9 macro= 212;9 macro=12; -
issm/trunk/src/m/enum/StokesFormulationEnum.m
r1714 r3567 7 7 % macro=StokesFormulationEnum() 8 8 9 macro= 305;9 macro=41; -
issm/trunk/src/m/enum/SurfaceXAnalysisEnum.m
r1714 r3567 7 7 % macro=SurfaceXAnalysisEnum() 8 8 9 macro=2 41;9 macro=23; -
issm/trunk/src/m/enum/SurfaceYAnalysisEnum.m
r1714 r3567 7 7 % macro=SurfaceYAnalysisEnum() 8 8 9 macro=24 2;9 macro=24; -
issm/trunk/src/m/enum/ThermalAnalysisEnum.m
r1714 r3567 7 7 % macro=ThermalAnalysisEnum() 8 8 9 macro= 230;9 macro=19; -
issm/trunk/src/m/enum/TransientAnalysisEnum.m
r1901 r3567 7 7 % macro=TransientAnalysisEnum() 8 8 9 macro=2 31;9 macro=20; -
issm/trunk/src/m/enum/TriaEnum.m
r1714 r3567 7 7 % macro=TriaEnum() 8 8 9 macro=4 11;9 macro=45; -
issm/trunk/src/m/enum/VertAnalysisEnum.m
r1714 r3567 7 7 % macro=VertAnalysisEnum() 8 8 9 macro= 214;9 macro=14; -
issm/trunk/src/m/enum/VertexEnum.m
r3417 r3567 7 7 % macro=VertexEnum() 8 8 9 macro= 420;9 macro=53; -
issm/trunk/src/m/enum/VerticesEnum.m
r3417 r3567 7 7 % macro=VerticesEnum() 8 8 9 macro= 108;9 macro=8; -
issm/trunk/src/m/enum/WaterEnum.m
r1750 r3567 7 7 % macro=WaterEnum() 8 8 9 macro= 601;9 macro=73; -
issm/trunk/src/mex/TriMeshProcessRifts/TriMeshProcessRifts.cpp
r1904 r3567 287 287 /*Friction fraction, fractionincrement and fill: */ 288 288 mxSetField(pmxa_array,i,"friction",mxCreateDoubleScalar(0)); 289 mxSetField(pmxa_array,i,"fill",mxCreateDoubleScalar(IceEnum ())); //default is ice289 mxSetField(pmxa_array,i,"fill",mxCreateDoubleScalar(IceEnum)); //default is ice 290 290 mxSetField(pmxa_array,i,"fraction",mxCreateDoubleScalar(0)); //default is ice 291 291 mxSetField(pmxa_array,i,"fractionincrement",mxCreateDoubleScalar(0.1)); … … 295 295 /*output NaN :*/ 296 296 pNaN=(double*)xmalloc(sizeof(double)); 297 *pNaN= 0.0/0.0;297 *pNaN=NAN; 298 298 pmxa_array= mxCreateDoubleMatrix(0,0,mxREAL); 299 299 mxSetM(pmxa_array,1);
Note:
See TracChangeset
for help on using the changeset viewer.