From 0adef2be78b0c8e7fd69c7d0a44ec822c0efc2b4 Mon Sep 17 00:00:00 2001
From: deseilligny <marc.pierrot-deseilligny@ign.fr>
Date: Mon, 13 May 2024 21:56:29 +0200
Subject: [PATCH] Pb compile

---
 MMVII/src/UtiMaths/Interpolators.cpp | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/MMVII/src/UtiMaths/Interpolators.cpp b/MMVII/src/UtiMaths/Interpolators.cpp
index cce2539056..e8217dcbe2 100755
--- a/MMVII/src/UtiMaths/Interpolators.cpp
+++ b/MMVII/src/UtiMaths/Interpolators.cpp
@@ -468,6 +468,7 @@ void cTabulatedInterpolator::SetDiff(const cTabulatedInterpolator & anInt)
 void cTabulatedInterpolator::DoNormalize(bool ForDerive)
 {
       tREAL8 aSomWDif = 0; // observation of inital deviation, for eventual show/debug
+      tREAL8 aCheckS=0.0;   // Check sum, useful 4 derive
 			   
       int aAmpl = (mSzTot/mNbTabul+3) * mNbTabul;
       // Avoid  parse twice the phase
@@ -490,7 +491,6 @@ void cTabulatedInterpolator::DoNormalize(bool ForDerive)
 
 
 	  // divide/substratc to all value same phase
-	  tREAL8 aCheckS=0.0;
           for (int aKSigned=aKPhase-aAmpl ; aKSigned<aAmpl ; aKSigned+=mNbTabul)
 	  {
                int aKAbs =  std::abs(aKSigned);
@@ -504,7 +504,7 @@ void cTabulatedInterpolator::DoNormalize(bool ForDerive)
 		   }
 	           else
 		   {
-                      // if the phase = mNbTabul we dont want to divide twice the result
+                      // if the 2*phase = mNbTabul we dont want to divide twice the result
                       if (  ((2*aKPhase)!=mNbTabul) || (aKSigned<=0))
 		      {
                            mDIm->SetV(aKAbs,mDIm->GetV(aKAbs)/aSumV);
@@ -520,7 +520,7 @@ void cTabulatedInterpolator::DoNormalize(bool ForDerive)
       if (0)
       {
           aSomWDif /= mNbTabul;
-          StdOut() << "SSS= " << aSomWDif  << " NbT=" << mNbTabul << "\n";
+          StdOut() << "SSS= " << aSomWDif  << " NbT=" << mNbTabul   << " CheckS=" << aCheckS << "\n";
           getchar();
       }
 }