Changeset 20018

Timestamp:

01/29/16 14:47:47 (9 years ago)

Author:

Eric.Larour

Message:

CHG: new implementation of legendre function.

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -3648 +3648 @@
                                 else{
                                         for(int n=0;n<nl;n++){
-                                                Pn=legendre(Pn1,Pn2,reCast<int,IssmDouble>(cosalpha),n); Pn1=Pn2; Pn2=Pn;
+                                                if(n==0)Pn=1;
+                                                else if(n==1)Pn=cosalpha;
+                                                else Pn= ( (2*n-1)*cosalpha*Pn1 - (n-1)*Pn2 ) /n;
+                                                Pn2=Pn1; Pn1=Pn;
                                                 G_elastic += deltalove[n]*Pn;
                                         }
@@ -3809 +3812 @@
 
                 /*initialize G_elastic:*/
-                G_elastic=xNew<IssmDouble>(gsize);
+                G_elastic=xNewZeroInit<IssmDouble>(gsize);
         }
         if(compute_G_rigid){
                 /*Initialize G_rigid and G_elastic:*/
-                G_rigid=xNew<IssmDouble>(gsize);
+                G_rigid=xNewZeroInit<IssmDouble>(gsize);
         }
 
@@ -3846 +3849 @@
                                 else{
                                         for(int n=0;n<nl;n++){
-                                                Pn=legendre(Pn1,Pn2,reCast<int,IssmDouble>(cosalpha),n); Pn1=Pn2; Pn2=Pn;
+                                                if(n==0)Pn=1;
+                                                else if(n==1)Pn=cosalpha;
+                                                else Pn= ( (2*n-1)*cosalpha*Pn1 - (n-1)*Pn2 ) /n;
+                                                Pn2=Pn1; Pn1=Pn;
                                                 G_elastic[i] += deltalove[n]*Pn;
                                         }