Index: /issm/trunk-jpl/src/c/classes/Cfsurfacelogvel.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Cfsurfacelogvel.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/Cfsurfacelogvel.cpp (revision 25446) @@ -168,7 +168,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=topelement->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/classes/Cfsurfacesquare.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Cfsurfacesquare.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/Cfsurfacesquare.cpp (revision 25446) @@ -172,7 +172,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=topelement->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp (revision 25446) @@ -736,6 +736,5 @@ Gauss* gauss=this->NewGauss(5); - for(int ig=gauss->begin();igend();ig++){ - gauss->GaussPoint(ig); + while(gauss->next()){ this->JacobianDeterminant(&Jdet,xyz_list,gauss); Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/FemModel.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/FemModel.cpp (revision 25446) @@ -990,7 +990,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ @@ -2080,7 +2078,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ @@ -2138,7 +2134,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ @@ -2697,7 +2691,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ @@ -2752,7 +2744,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ @@ -2871,7 +2861,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ @@ -4148,6 +4136,5 @@ Gauss* gauss=element->NewGauss(2); error=0.; - for(int ig=gauss->begin();igend();ig++){ - gauss->GaussPoint(ig); + while(gauss->next()){ element->JacobianDeterminant(&Jdet,xyz_list,gauss); element->NodalFunctions(basis,gauss); @@ -4212,5 +4199,5 @@ /*Get the gradient of thickness at the center point (in fact, GradH is constante over the element)*/ Gauss* gauss=element->NewGauss(1); - gauss->GaussPoint(gauss->begin()); + gauss->GaussPoint(0); element->ValueP1DerivativesOnGauss(GradH,H,xyz_list,gauss); @@ -4295,11 +4282,10 @@ /*Get the gradient of thickness*/ Gauss* gaussH=element->NewGauss(1); - gaussH->GaussPoint(gaussH->begin()); + gaussH->GaussPoint(0); element->ValueP1DerivativesOnGauss(GradH,H,xyz_list,gaussH); /*Integrate*/ Gauss* gauss=element->NewGauss(2); error=0.; - for(int ig=gauss->begin();igend();ig++){ - gauss->GaussPoint(ig); + while(gauss->next()){ element->JacobianDeterminant(&Jdet,xyz_list,gauss); element->NodalFunctions(basis,gauss); Index: /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp (revision 25446) @@ -403,6 +403,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentJacobianDeterminant(&Jdet,&xyz_list[0][0],gauss); @@ -534,6 +533,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentJacobianDeterminant(&Jdet,&xyz_list[0][0],gauss); Index: /issm/trunk-jpl/src/c/classes/Loads/Neumannflux.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Loads/Neumannflux.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/Loads/Neumannflux.cpp (revision 25446) @@ -368,7 +368,5 @@ /* Start looping on the number of gaussian points: */ GaussTria* gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentJacobianDeterminant(&Jdet,&xyz_list[0][0],gauss); @@ -412,7 +410,5 @@ /* Start looping on the number of gaussian points: */ GaussTria* gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentJacobianDeterminant(&Jdet,&xyz_list[0][0],gauss); Index: /issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/Loads/Numericalflux.cpp (revision 25446) @@ -541,7 +541,5 @@ /* Start looping on the number of gaussian points: */ gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentNodalFunctions(&basis[0],gauss,index1,index2,tria->FiniteElement()); @@ -598,7 +596,5 @@ int index2=tria->GetVertexIndex(vertices[1]); GaussTria* gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentNodalFunctions(&basis_plus[0] ,gauss,index1,index2,tria->FiniteElement()); @@ -713,7 +709,5 @@ /* Start looping on the number of gaussian points: */ gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentNodalFunctions(&basis[0],gauss,index1,index2,tria->FiniteElement()); @@ -771,7 +765,5 @@ int index2=tria->GetVertexIndex(vertices[1]); GaussTria* gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentNodalFunctions(&basis_plus[0] ,gauss,index1,index2,tria->FiniteElement()); @@ -890,7 +882,5 @@ /* Start looping on the number of gaussian points: */ gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentNodalFunctions(&basis[0],gauss,index1,index2,tria->FiniteElement()); @@ -971,7 +961,5 @@ /* Start looping on the number of gaussian points: */ gauss=new GaussTria(index1,index2,2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ tria->GetSegmentNodalFunctions(&basis[0],gauss,index1,index2,tria->FiniteElement()); Index: /issm/trunk-jpl/src/c/classes/gauss/Gauss.h =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/Gauss.h (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/Gauss.h (revision 25446) @@ -12,7 +12,5 @@ virtual ~Gauss(){}; - virtual int begin(void)=0; virtual void Echo(void)=0; - virtual int end(void)=0; virtual void Reset(void)=0; virtual bool next(void)=0; Index: /issm/trunk-jpl/src/c/classes/gauss/GaussPenta.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussPenta.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussPenta.cpp (revision 25446) @@ -466,18 +466,4 @@ /*Methods*/ -int GaussPenta::begin(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - _assert_(coords2); - _assert_(coords3); - _assert_(coords4); - - /*return first gauss index*/ - return 0; -} -/*}}}*/ void GaussPenta::Echo(void){/*{{{*/ @@ -522,18 +508,4 @@ _printf_(" coord4 = " << coord4 << "\n"); -} -/*}}}*/ -int GaussPenta::end(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - _assert_(coords2); - _assert_(coords3); - _assert_(coords4); - - /*return last gauss index +1*/ - return numgauss; } /*}}}*/ Index: /issm/trunk-jpl/src/c/classes/gauss/GaussPenta.h =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussPenta.h (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussPenta.h (revision 25446) @@ -43,7 +43,5 @@ /*Methods*/ bool next(void); - int begin(void); void Echo(void); - int end(void); int Enum(void); void GaussFaceTria(int index1, int index2, int index3, int order); Index: /issm/trunk-jpl/src/c/classes/gauss/GaussSeg.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussSeg.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussSeg.cpp (revision 25446) @@ -77,15 +77,4 @@ /*Methods*/ -int GaussSeg::begin(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - - /*return first gauss index*/ - return 0; -} -/*}}}*/ void GaussSeg::Echo(void){/*{{{*/ @@ -107,15 +96,4 @@ _printf_(" coord1 = " << coord1 << "\n"); -} -/*}}}*/ -int GaussSeg::end(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - - /*return last gauss index +1*/ - return numgauss; } /*}}}*/ Index: /issm/trunk-jpl/src/c/classes/gauss/GaussSeg.h =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussSeg.h (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussSeg.h (revision 25446) @@ -31,7 +31,5 @@ /*Methods*/ bool next(void); - int begin(void); void Echo(void); - int end(void); int Enum(void); void GaussPoint(int ig); Index: /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.cpp (revision 25446) @@ -97,18 +97,4 @@ /*Methods*/ -int GaussTetra::begin(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - _assert_(coords2); - _assert_(coords3); - _assert_(coords4); - - /*return first gauss index*/ - return 0; -} -/*}}}*/ void GaussTetra::Echo(void){/*{{{*/ @@ -153,18 +139,4 @@ _printf_(" coord4 = " << coord4 << "\n"); -} -/*}}}*/ -int GaussTetra::end(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - _assert_(coords2); - _assert_(coords3); - _assert_(coords4); - - /*return last gauss index +1*/ - return numgauss; } /*}}}*/ Index: /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.h =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.h (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussTetra.h (revision 25446) @@ -37,7 +37,5 @@ /*Methods*/ bool next(void); - int begin(void); void Echo(void); - int end(void); int Enum(void); void GaussPoint(int ig); Index: /issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussTria.cpp (revision 25446) @@ -426,17 +426,4 @@ /*Methods*/ -int GaussTria::begin(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - _assert_(coords2); - _assert_(coords3); - - /*return first gauss index*/ - return 0; -} -/*}}}*/ void GaussTria::Echo(void){/*{{{*/ @@ -474,17 +461,4 @@ _printf_(" coord3 = " << coord3 << "\n"); -} -/*}}}*/ -int GaussTria::end(void){/*{{{*/ - - /*Check that this has been initialized*/ - _assert_(numgauss>0); - _assert_(weights); - _assert_(coords1); - _assert_(coords2); - _assert_(coords3); - - /*return last gauss index +1*/ - return numgauss; } /*}}}*/ Index: /issm/trunk-jpl/src/c/classes/gauss/GaussTria.h =================================================================== --- /issm/trunk-jpl/src/c/classes/gauss/GaussTria.h (revision 25445) +++ /issm/trunk-jpl/src/c/classes/gauss/GaussTria.h (revision 25446) @@ -39,7 +39,5 @@ /*Methods*/ bool next(void); - int begin(void); void Echo(void); - int end(void); int Enum(void); void GaussEdgeCenter(int index1,int index2); Index: /issm/trunk-jpl/src/c/modules/DragCoefficientAbsGradientx/DragCoefficientAbsGradientx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/DragCoefficientAbsGradientx/DragCoefficientAbsGradientx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/DragCoefficientAbsGradientx/DragCoefficientAbsGradientx.cpp (revision 25446) @@ -77,7 +77,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=basalelement->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/RheologyBAbsGradientx/RheologyBAbsGradientx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/RheologyBAbsGradientx/RheologyBAbsGradientx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/RheologyBAbsGradientx/RheologyBAbsGradientx.cpp (revision 25446) @@ -59,7 +59,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ @@ -131,7 +129,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/RheologyBbarAbsGradientx/RheologyBbarAbsGradientx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/RheologyBbarAbsGradientx/RheologyBbarAbsGradientx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/RheologyBbarAbsGradientx/RheologyBbarAbsGradientx.cpp (revision 25446) @@ -65,7 +65,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=basalelement->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/SurfaceAbsVelMisfitx/SurfaceAbsVelMisfitx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/SurfaceAbsVelMisfitx/SurfaceAbsVelMisfitx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/SurfaceAbsVelMisfitx/SurfaceAbsVelMisfitx.cpp (revision 25446) @@ -72,7 +72,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=topelement->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/SurfaceAverageVelMisfitx/SurfaceAverageVelMisfitx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/SurfaceAverageVelMisfitx/SurfaceAverageVelMisfitx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/SurfaceAverageVelMisfitx/SurfaceAverageVelMisfitx.cpp (revision 25446) @@ -80,7 +80,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=topelement->NewGauss(3); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/SurfaceLogVelMisfitx/SurfaceLogVelMisfitx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/SurfaceLogVelMisfitx/SurfaceLogVelMisfitx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/SurfaceLogVelMisfitx/SurfaceLogVelMisfitx.cpp (revision 25446) @@ -74,7 +74,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=topelement->NewGauss(4); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/SurfaceLogVxVyMisfitx/SurfaceLogVxVyMisfitx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/SurfaceLogVxVyMisfitx/SurfaceLogVxVyMisfitx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/SurfaceLogVxVyMisfitx/SurfaceLogVxVyMisfitx.cpp (revision 25446) @@ -74,7 +74,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=topelement->NewGauss(4); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/SurfaceRelVelMisfitx/SurfaceRelVelMisfitx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/SurfaceRelVelMisfitx/SurfaceRelVelMisfitx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/SurfaceRelVelMisfitx/SurfaceRelVelMisfitx.cpp (revision 25446) @@ -74,7 +74,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=topelement->NewGauss(4); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/ThicknessAbsMisfitx/ThicknessAbsMisfitx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ThicknessAbsMisfitx/ThicknessAbsMisfitx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/ThicknessAbsMisfitx/ThicknessAbsMisfitx.cpp (revision 25446) @@ -50,7 +50,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/ThicknessAcrossGradientx/ThicknessAcrossGradientx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ThicknessAcrossGradientx/ThicknessAcrossGradientx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/ThicknessAcrossGradientx/ThicknessAcrossGradientx.cpp (revision 25446) @@ -53,7 +53,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */ Index: /issm/trunk-jpl/src/c/modules/ThicknessAlongGradientx/ThicknessAlongGradientx.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ThicknessAlongGradientx/ThicknessAlongGradientx.cpp (revision 25445) +++ /issm/trunk-jpl/src/c/modules/ThicknessAlongGradientx/ThicknessAlongGradientx.cpp (revision 25446) @@ -53,7 +53,5 @@ /* Start looping on the number of gaussian points: */ Gauss* gauss=element->NewGauss(2); - for(int ig=gauss->begin();igend();ig++){ - - gauss->GaussPoint(ig); + while(gauss->next()){ /* Get Jacobian determinant: */