Changeset 20645 for issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
- Timestamp:
- 05/25/16 10:52:44 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp ¶
r20644 r20645 41 41 42 42 /*Fetch parameters: */ 43 iomodel-> Constant(&g,ConstantsGEnum);44 iomodel-> Constant(&rho_ice,MaterialsRhoIceEnum);45 iomodel-> Constant(&FSreconditioning,StressbalanceFSreconditioningEnum);46 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);47 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);48 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);49 iomodel-> Constant(&isHO,FlowequationIsHOEnum);50 iomodel-> Constant(&isFS,FlowequationIsFSEnum);43 iomodel->FindConstant(&g,ConstantsGEnum); 44 iomodel->FindConstant(&rho_ice,MaterialsRhoIceEnum); 45 iomodel->FindConstant(&FSreconditioning,StressbalanceFSreconditioningEnum); 46 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 47 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 48 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 49 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 50 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 51 51 52 52 /*Now, is the flag macayaealHO on? otherwise, do nothing: */ … … 63 63 64 64 /*Get finite element type*/ 65 if(isSSA) iomodel-> Constant(&finiteelement,FlowequationFeSSAEnum);65 if(isSSA) iomodel->FindConstant(&finiteelement,FlowequationFeSSAEnum); 66 66 else if(isL1L2) finiteelement = P1Enum; 67 else if(isHO) iomodel-> Constant(&finiteelement,FlowequationFeHOEnum);68 else if(isFS){ iomodel-> Constant(&finiteelement,FlowequationFeFSEnum);67 else if(isHO) iomodel->FindConstant(&finiteelement,FlowequationFeHOEnum); 68 else if(isFS){ iomodel->FindConstant(&finiteelement,FlowequationFeFSEnum); 69 69 /*Deduce velocity interpolation from finite element*/ 70 70 switch(finiteelement){ … … 451 451 452 452 /*Fetch parameters: */ 453 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);454 iomodel-> Constant(&isFS,FlowequationIsFSEnum);455 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);456 iomodel-> Constant(&isHO,FlowequationIsHOEnum);457 iomodel-> Constant(&numrifts,RiftsNumriftsEnum);453 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 454 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 455 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 456 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 457 iomodel->FindConstant(&numrifts,RiftsNumriftsEnum); 458 458 459 459 /*Now, is the flag macayaealHO on? otherwise, do nothing: */ … … 507 507 508 508 /*Fetch parameters: */ 509 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);510 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);511 iomodel-> Constant(&isHO,FlowequationIsHOEnum);512 iomodel-> Constant(&isFS,FlowequationIsFSEnum);509 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 510 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 511 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 512 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 513 513 514 514 /*Now, check that we have non SIA elements */ … … 527 527 if(isSSA){ 528 528 approximation=SSAApproximationEnum; 529 iomodel-> Constant(&finiteelement,FlowequationFeSSAEnum);529 iomodel->FindConstant(&finiteelement,FlowequationFeSSAEnum); 530 530 } 531 531 else if(isL1L2){ … … 535 535 else if(isHO){ 536 536 approximation = HOApproximationEnum; 537 iomodel-> Constant(&finiteelement,FlowequationFeHOEnum);537 iomodel->FindConstant(&finiteelement,FlowequationFeHOEnum); 538 538 } 539 539 else if(isFS){ 540 540 approximation = FSApproximationEnum; 541 iomodel-> Constant(&finiteelement,FlowequationFeFSEnum);541 iomodel->FindConstant(&finiteelement,FlowequationFeFSEnum); 542 542 } 543 543 iomodel->FetchData(3,FlowequationBorderSSAEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum); … … 681 681 682 682 /*Fetch constants needed: */ 683 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);684 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);685 iomodel-> Constant(&isHO,FlowequationIsHOEnum);686 iomodel-> Constant(&isFS,FlowequationIsFSEnum);687 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);688 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);689 iomodel-> Constant(&materials_type,MaterialsEnum);690 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);691 iomodel-> Constant(&frictionlaw,FrictionLawEnum);683 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 684 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 685 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 686 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 687 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 688 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 689 iomodel->FindConstant(&materials_type,MaterialsEnum); 690 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 691 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 692 692 693 693 /*return if no processing required*/ … … 706 706 /*Get finite element type*/ 707 707 if(!iscoupling){ 708 if(isSSA) iomodel-> Constant(&finiteelement,FlowequationFeSSAEnum);708 if(isSSA) iomodel->FindConstant(&finiteelement,FlowequationFeSSAEnum); 709 709 else if(isL1L2) finiteelement = P1Enum; 710 else if(isHO) iomodel-> Constant(&finiteelement,FlowequationFeHOEnum);711 else if(isFS) iomodel-> Constant(&finiteelement,FlowequationFeFSEnum);710 else if(isHO) iomodel->FindConstant(&finiteelement,FlowequationFeHOEnum); 711 else if(isFS) iomodel->FindConstant(&finiteelement,FlowequationFeFSEnum); 712 712 for(int i=0;i<iomodel->numberofelements;i++){ 713 713 finiteelement_list[i]=finiteelement; … … 773 773 } 774 774 /*LATH parameters*/ 775 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);775 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 776 776 if(fe_FS==LATaylorHoodEnum || fe_FS==LACrouzeixRaviartEnum){ 777 777 InputUpdateFromConstantx(elements,0.,SigmaNNEnum); … … 790 790 break; 791 791 case 3: 792 iomodel-> Constant(&FrictionCoupling,FrictionCouplingEnum);792 iomodel->FindConstant(&FrictionCoupling,FrictionCouplingEnum); 793 793 iomodel->FetchDataToInput(elements,FrictionCEnum); 794 794 iomodel->FetchDataToInput(elements,FrictionAsEnum); … … 862 862 863 863 /*XTH LATH parameters*/ 864 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);864 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 865 865 if(fe_FS==XTaylorHoodEnum || fe_FS==LATaylorHoodEnum || fe_FS==LACrouzeixRaviartEnum){ 866 866 parameters->AddObject(iomodel->CopyConstantObject(AugmentedLagrangianREnum)); … … 869 869 } 870 870 871 iomodel-> Constant(&materials_type,MaterialsEnum);871 iomodel->FindConstant(&materials_type,MaterialsEnum); 872 872 if(materials_type==MatdamageiceEnum){ 873 873 parameters->AddObject(iomodel->CopyConstantObject(DamageLawEnum)); … … 884 884 /*Deal with friction parameters*/ 885 885 int frictionlaw; 886 iomodel-> Constant(&frictionlaw,FrictionLawEnum);886 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 887 887 if(frictionlaw==4 || frictionlaw==6) parameters->AddObject(iomodel->CopyConstantObject(FrictionGammaEnum)); 888 888 if(frictionlaw==3) parameters->AddObject(iomodel->CopyConstantObject(FrictionCouplingEnum));
Note:
See TracChangeset
for help on using the changeset viewer.