Changeset 87 for altifloat

Timestamp:

04/07/13 10:33:17 (11 years ago)

Author:

leila_ocean

Message:

correction of scalar product

File:

• altifloat/src/r_float_d.h (modified) (2 diffs)

altifloat/src/r_float_d.h

@@ -1 +1 @@
 forward(YREAL xd, YREAL yd, YREAL ud, YREAL vd)
 {
-        YS1=xd+rdtflo*ud+xd*rdtflo*(ur_dx+ur_dy);//mot cle YS1: yao sortie 1
+        
+        //what we had up to 01/04
+        /*YS1=xd+rdtflo*ud+xd*rdtflo*(ur_dx+ur_dy);//mot cle YS1: yao sortie 1
         printf("xd=%f 2term=%f 3term=%f\n",xd,rdtflo*ud,xd*rdtflo*(ur_dx+ur_dy));
-        YS2=yd+rdtflo*vd+yd*rdtflo*(vr_dx+vr_dy);
+        YS2=yd+rdtflo*vd+yd*rdtflo*(vr_dx+vr_dy);*/
+        
+        //what we should have corrected on 04/04
+        YS1=xd+rdtflo*ud+rdtflo*(xd*ur_dx+yd*ur_dy);
+        YS2=yd+rdtflo*vd+rdtflo*(xd*vr_dx+yd*vr_dy);
 }
 
 backward(YREAL xd, YREAL yd, YREAL ud, YREAL vd)
 {
-        YJ1I1=1+rdtflo*(ur_dx+ur_dy);   //mot cle, les entrees du gradient du out par rapport au in, dYs1/dx,y,u,v
+        /*YJ1I1=1+rdtflo*(ur_dx+ur_dy);   //mot cle, les entrees du gradient du out par rapport au in, dYs1/dx,y,u,v
         YJ1I2=0;
         YJ1I3=rdtflo;
@@ -16 +22 @@
         YJ2I2=1+rdtflo*(vr_dx+vr_dy);
         YJ2I3=0;
+        YJ2I4=rdtflo;*/
+        
+        
+        //what we sould have corrected on 04/04
+        YJ1I1=1+rdtflo*ur_dx;   
+    YJ1I2=rdtflo*ur_dy;
+        YJ1I3=rdtflo;
+        YJ1I4=0;
+        
+        YJ2I1=rdtflo*vr_dx; 
+        YJ2I2=1+rdtflo*vr_dy;
+        YJ2I3=0;
         YJ2I4=rdtflo;