ref: 590521f8c2da89fa1bd1a4a67bf80e91a008d36d
parent: a3cce095457eae34469dcdb2d099df33082e0c50
author: menno <menno>
date: Mon Feb 7 13:57:47 EST 2000
Same change to short blocks as oxy2000 made to long blocks
--- a/psych.c
+++ b/psych.c
@@ -52,9 +52,9 @@
Source file:
-$Id: psych.c,v 1.40 2000/02/07 18:17:30 menno Exp $
-$Id: psych.c,v 1.40 2000/02/07 18:17:30 menno Exp $
-$Id: psych.c,v 1.40 2000/02/07 18:17:30 menno Exp $
+$Id: psych.c,v 1.41 2000/02/07 18:57:47 menno Exp $
+$Id: psych.c,v 1.41 2000/02/07 18:57:47 menno Exp $
+$Id: psych.c,v 1.41 2000/02/07 18:57:47 menno Exp $
**********************************************************************/
@@ -761,8 +761,11 @@
fp = psy_var_short->f_pred[i][w];
if( fabs(r) + fabs(rp) != 0.0 )
+/* Replacing with something faster
psy_var_short->c[i][w] = sqrt( psy_sqr(r*cos(f) - rp*cos(fp))
+psy_sqr(r*sin(f) - rp*sin(fp)) )/ ( r + fabs(rp) ) ;
+*/
+ psy_var_short->c[i][w] = sqrt( r*r+rp*rp-2*r*rp*cos(f+fp) )/ ( r + fabs(rp) );
else
psy_var_short->c[i][w] = 0.0; /* tmp */
}