Changeset 2148 in Sophya for trunk/SophyaLib/Samba

Timestamp:

Jul 31, 2002, 1:42:54 PM (23 years ago)

Author:

ansari

Message:

Adaptation a fmath.h vide (pas de #define fsqrt() fabsf() , Reza 31/7/2002

File:

• trunk/SophyaLib/Samba/bruit.cc (modified) (9 diffs)

trunk/SophyaLib/Samba/bruit.cc

@@ -37 +37 @@
 //--
 
+// Definition de fonctions pour float - a changer en utilisant cmath
+//                              Reza Juillet 2002
+inline float _myfabsf(float x) { return ((float)(fabs((double)x)) ); }
+inline float _mysqrtf(float x) { return ((float)(sqrt((double)x)) ); }
+ 
 /* --Methode-- */
 //++
@@ -345 +350 @@
 //++
 MemNoise::MemNoise(float sigma, int mem, float tau, int ava)
-    : NoiseGenerator(sqrtf(fabsf(tau))*sigma)
+    : NoiseGenerator(_mysqrtf(_myfabsf(tau))*sigma)
 //
 //--
@@ -468 +473 @@
         if(mStPos[i]==0.)continue;
         //      printf("i,temps,bruit %d,%g,%g \n",i,mTePos[i],mStPos[i]);
-        Somme = mStPos[i]/sqrtf(mTePos[i]);
+        Somme = mStPos[i]/_mysqrtf(mTePos[i]);
             SbaryT = mStPos[i] * mTePos[i];
             Poids = mStPos[i];
@@ -474 +479 @@
                 {
                     //          printf("  j,temps,bruit %d,%g,%g \n",j,mTePos[j],mStPos[j]);
-                Somme += mStPos[j]/sqrtf(mTePos[j]);
+                Somme += mStPos[j]/_mysqrtf(mTePos[j]);
                     Poids += mStPos[j];
                     SbaryT += mStPos[j] * mTePos[j];
                     BaryT = SbaryT / Poids;
-                    Approx = Poids / sqrtf(BaryT);
+                    Approx = Poids / _mysqrtf(BaryT);
                     Erreur = (Approx - Somme) / Somme;
-                    //                              printf("i,j = %d %d Approx = %g \n",i,j,fabsf(Erreur));
+                    //                              printf("i,j = %d %d Approx = %g \n",i,j,_myfabsf(Erreur));
                     //              printf("Temps,Bruit = %g %g \n",BaryT,Poids);
-                    if(fabsf(Erreur)<epsilon)
+                    if(_myfabsf(Erreur)<epsilon)
                         {
                             mStPos[i]=Poids;
@@ -521 +526 @@
         if(mStNeg[i]==0.)continue;
         //      printf("i,temps,bruit %d,%g,%g \n",i,mTeNeg[i],mStNeg[i]);
-        Somme = mStNeg[i]/sqrtf(mTeNeg[i]);
+        Somme = mStNeg[i]/_mysqrtf(mTeNeg[i]);
             SbaryT = mStNeg[i] * mTeNeg[i];
             Poids = mStNeg[i];
@@ -527 +532 @@
                 {
                     //          printf("  j,temps,bruit %d,%g,%g \n",j,mTeNeg[j],mStNeg[j]);
-                Somme += mStNeg[j]/sqrtf(mTeNeg[j]);
+                Somme += mStNeg[j]/_mysqrtf(mTeNeg[j]);
                     Poids += mStNeg[j];
                     SbaryT += mStNeg[j] * mTeNeg[j];
                     BaryT = SbaryT / Poids;
-                    Approx = Poids / sqrtf(BaryT);
+                    Approx = Poids / _mysqrtf(BaryT);
                     Erreur = (Approx - Somme) / Somme;
                     //              printf("i,j = %d %d Approx = %g \n",i,j,Erreur);
                     //              printf("Temps,Bruit = %g %g \n",BaryT,Poids);
-                    if(fabsf(Erreur)<epsilon)
+                    if(_myfabsf(Erreur)<epsilon)
                         {
                             mStNeg[i]=Poids;
@@ -588 +593 @@
 }
 // float ajout=0.;
-for(i=indlim; i < mMemPos-1; i++)  rn += mStPos[i]/sqrtf(mTePos[i]);
+for(i=indlim; i < mMemPos-1; i++)  rn += mStPos[i]/_mysqrtf(mTePos[i]);
 indlim=0;
 for(i=0; i < mMemNeg-1; i++) 
@@ -596 +601 @@
 indlim++;
 }
-for(i=indlim; i < mMemNeg-1; i++) rn += mStNeg[i]/sqrtf(mTeNeg[i]);
+for(i=indlim; i < mMemNeg-1; i++) rn += mStNeg[i]/_mysqrtf(mTeNeg[i]);
 }
 }
@@ -649 +654 @@
 // Je calcule les sigma du MemNoise d'apres la formule de Francois 
 // mn = new MemNoise(sig*fknee);  
-// mn = new MemNoise(1.4142*sig*sqrtf(fknee),1,10.,1000000); // modif appel FC 13 mai 98.  
+// mn = new MemNoise(1.4142*sig*_mysqrtf(fknee),1,10.,1000000); // modif appel FC 13 mai 98.  
 mn = new MemNoise(1.4142*sig*sqrt((double)fknee),mem,tau,12000); // FC 13 mai 98. -> Reza 1/12/99 
 }