Changeset 343 for trunk/projects/polly/lreg.c

Timestamp:

01/17/07 00:38:49 (4 years ago)

Author:

anthony_rowe

Message:

added ilp functions

Files:

• trunk/projects/polly/lreg.c (modified) (2 diffs)

trunk/projects/polly/lreg.c

@@ -114 +114 @@
 
     // standard deviation of the points
-    double stddevPoints = isqrt( (SSyy - b * SSxy)/(size-2) );
+   // printf( "1 sqrt( %f )=",(SSyy - b * SSxy)/(size-2));
+   // double stddevPoints = isqrt( (SSyy - b * SSxy)/(size-2) );
+   // printf( "%f\n",stddevPoints);
 
     // Error of the slope
-    double bError = stddevPoints / isqrt( SSxx );
+   // printf( "2 sqrt( %f )=",SSxx);
+   // double bError = stddevPoints / isqrt( SSxx );
+   // printf( "%f\n",isqrt( SSxx));
 
     double r2Numerator = (size * Sxy) - (Sx * Sy);
@@ -123 +127 @@
     double r2 = (r2Numerator * r2Numerator) / r2Denominator;
 
-    double signB = (b < 0) ? -1.0 : 1.0;
-    printf( "signB = %f r2 = %f\n",signB,r2);
-    double r = signB * isqrt( r2*10000 );
+    //double signB = (b < 0) ? -1.0 : 1.0;
+    //printf( "3 sqrt(%f)=",r2);
+    //double r = signB * isqrt(r2);
+    //printf( "%f\n",r);
 
-    reg_out->r = r; 
+    reg_out->r_sqr = r2; 
     reg_out->a = a; 
     reg_out->b = b; 
-    reg_out->bError = bError; 
-    reg_out->stddevPoints= stddevPoints; 
+    reg_out->bError = 0; 
+    reg_out->stddevPoints= 0; 
+    //reg_out->bError = bError; 
+    //reg_out->stddevPoints= stddevPoints; 
   } // if N > 0