Changeset 22015

Timestamp:

08/25/17 21:33:42 (8 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing AD

File:

• issm/trunk-jpl/src/c/cores/love_core.cpp (modified) (5 diffs)

issm/trunk-jpl/src/c/cores/love_core.cpp

@@ -17 +17 @@
 
         /*love parameters: */
-        IssmDouble* frequencies=NULL;
-        int nfreq,dummy; 
-        int sh_nmin,sh_nmax;
-        IssmDouble g0,r0,mu0;
-        bool allow_layer_deletion;
-        int forcing_type;
-        bool verbosemod =(int) VerboseModule();
+        IssmDouble *frequencies = NULL;
+        int         nfreq,dummy;
+        int         sh_nmin,sh_nmax;
+        IssmDouble  g0,r0,mu0;
+        bool        allow_layer_deletion;
+        int         forcing_type;
+        bool        verbosemod = (int)VerboseModule();
 
         /*parameters: */
@@ -29 +29 @@
         
         if(VerboseSolution()) _printf0_("   computing LOVE numbers\n");
-
 
         /*Recover some parameters: */
@@ -60 +59 @@
         IssmDouble*  LoveHi = xNewZeroInit<IssmDouble>(nfreq*(sh_nmax+1));
         IssmDouble*  LoveLi = xNewZeroInit<IssmDouble>(nfreq*(sh_nmax+1));
-
-
         
         /*call the main module: */
@@ -71 +68 @@
 
         /*Add love matrices to results:*/
-        femmodel->results->AddObject(new GenericExternalResult<IssmPDouble*>(femmodel->results->Size()+1,LoveKrEnum,LoveKr,sh_nmax+1,nfreq,0,0));
-        femmodel->results->AddObject(new GenericExternalResult<IssmPDouble*>(femmodel->results->Size()+1,LoveHrEnum,LoveHr,sh_nmax+1,nfreq,0,0));
-        femmodel->results->AddObject(new GenericExternalResult<IssmPDouble*>(femmodel->results->Size()+1,LoveLrEnum,LoveLr,sh_nmax+1,nfreq,0,0));
-        femmodel->results->AddObject(new GenericExternalResult<IssmPDouble*>(femmodel->results->Size()+1,LoveKiEnum,LoveKi,sh_nmax+1,nfreq,0,0));
-        femmodel->results->AddObject(new GenericExternalResult<IssmPDouble*>(femmodel->results->Size()+1,LoveHiEnum,LoveHi,sh_nmax+1,nfreq,0,0));
-        femmodel->results->AddObject(new GenericExternalResult<IssmPDouble*>(femmodel->results->Size()+1,LoveLiEnum,LoveLi,sh_nmax+1,nfreq,0,0));
+        femmodel->results->AddObject(new GenericExternalResult<IssmDouble*>(femmodel->results->Size()+1,LoveKrEnum,LoveKr,sh_nmax+1,nfreq,0,0));
+        femmodel->results->AddObject(new GenericExternalResult<IssmDouble*>(femmodel->results->Size()+1,LoveHrEnum,LoveHr,sh_nmax+1,nfreq,0,0));
+        femmodel->results->AddObject(new GenericExternalResult<IssmDouble*>(femmodel->results->Size()+1,LoveLrEnum,LoveLr,sh_nmax+1,nfreq,0,0));
+        femmodel->results->AddObject(new GenericExternalResult<IssmDouble*>(femmodel->results->Size()+1,LoveKiEnum,LoveKi,sh_nmax+1,nfreq,0,0));
+        femmodel->results->AddObject(new GenericExternalResult<IssmDouble*>(femmodel->results->Size()+1,LoveHiEnum,LoveHi,sh_nmax+1,nfreq,0,0));
+        femmodel->results->AddObject(new GenericExternalResult<IssmDouble*>(femmodel->results->Size()+1,LoveLiEnum,LoveLi,sh_nmax+1,nfreq,0,0));
 
         /*Free ressources:*/
@@ -86 +83 @@
         xDelete<IssmDouble>(LoveHi);
         xDelete<IssmDouble>(LoveLi);
-        
-
-
 }