Changeset 19986 for issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp

Timestamp:

01/22/16 10:47:46 (9 years ago)

Author:

Eric.Larour

Message:

CHG: added precomputation of legendre polynomials.

File:

• issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp

@@ -38 +38 @@
 void SealevelriseAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/
 
-        int nl;
+        int         nl;
         IssmDouble* love_h=NULL;
         IssmDouble* love_k=NULL;
+        
+        bool        legendre_precompute=false;
+        IssmDouble* legendre_coefficients=NULL;
+        int         M;
 
+        /*some constant parameters: */
         parameters->AddObject(iomodel->CopyConstantObject(SealevelriseReltolEnum));
         parameters->AddObject(iomodel->CopyConstantObject(SealevelriseAbstolEnum));
@@ -48 +53 @@
         parameters->AddObject(iomodel->CopyConstantObject(SealevelriseElasticEnum));
         parameters->AddObject(iomodel->CopyConstantObject(SealevelriseEustaticEnum));
+        parameters->AddObject(iomodel->CopyConstantObject(SealevelriseLegendrePrecomputeEnum));
+
         
+        /*love numbers: */
         iomodel->FetchData(&love_h,&nl,NULL,SealevelriseLoveHEnum);
         iomodel->FetchData(&love_k,&nl,NULL,SealevelriseLoveKEnum);
-        
         parameters->AddObject(new DoubleVecParam(SealevelriseLoveHEnum,love_h,nl));
         parameters->AddObject(new DoubleVecParam(SealevelriseLoveKEnum,love_k,nl));
+
+        /*legendre coefficients: */
+        iomodel->Constant(&legendre_precompute,SealevelriseLegendrePrecomputeEnum);
+        if(legendre_precompute){
+                iomodel->FetchData(&legendre_coefficients,&M,&nl,SealevelriseLoveKEnum);
+                parameters->AddObject(new DoubleMatParam(SealevelriseLegendreCoefficientsEnum,legendre_coefficients,M,nl));
+        }
 
         /*free ressources: */
         xDelete<IssmDouble>(love_h);
         xDelete<IssmDouble>(love_k);
+        xDelete<IssmDouble>(legendre_coefficients);
 
 }/*}}}*/