Changeset 20022

Timestamp:

01/29/16 17:15:42 (9 years ago)

Author:

Eric.Larour

Message:

CHG: precompute legendre polynomial values.

File:

• issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp (modified) (1 diff)

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

@@ -66 +66 @@
         iomodel->Constant(&legendre_precompute,SealevelriseLegendrePrecomputeEnum);
         if(legendre_precompute){
-                iomodel->FetchData(&legendre_coefficients,&M,&nl,SealevelriseLoveKEnum);
-                parameters->AddObject(new DoubleMatParam(SealevelriseLegendreCoefficientsEnum,legendre_coefficients,M,nl));
+
+                /*compute values at which to evaluate the legendre polynomical values:*/
+                const IssmDouble degacc=.01; M=reCast<int>(180/degacc+1);
+                IssmDouble* x=xNew<IssmDouble>(M);
+                for(int i=0;i<=M;i++){ //watch out the <=
+                        x[i]=-cos( (reCast<IssmDouble>(i)*degacc) * PI / 180.0);
+                }
+                /*compute legendre coefficient matrix:*/
+                legendre_coefficients=p_polynomial_value(M,nl-1,x);
+                parameters->AddObject(new DoubleMatParam(SealevelriseLegendreCoefficientsEnum,legendre_coefficients,nl,M));
+
+                /*free ressources:*/
+                xDelete<IssmDouble>(x);
         }