31 static double TeEvalCS = -1.;
32 static long int nzoneEval=0;
120 ASSERT( fabs(OTSLevel2) < 1e37 );
131 if( OTSLevel2/
SDIV(sumots) > 1e-4 )
157 enum {DEBUG_LOC=
false};
161 fprintf(
ioQQQ,
"%li\t%.2e\t%.2e\n",
210 return( ColStrGBar_v );
222 return( ColStrGBar_v );
236 yy = (1.398813573838321 + xx*(0.02943050869177121 + xx*
237 (-0.4439783893114510 + xx*(0.2316073358577902 + xx*(0.001870493481643103 +
238 xx*(-0.008227246351067403))))))/(1.0 + xx*(-0.6064792600526370 +
239 xx*(0.1958559534507252 + xx*(-0.02110452007196644 +
240 xx*(0.01348743933722316 + xx*(-0.0001944731334371711))))));
245 yy = (1.359675968512206 + xx*(0.04636500015069853 + xx*
246 (-0.4491620298246676 + xx*(0.2498199231048967 + xx*(0.005053803073345794 +
247 xx*(-0.01015647880244268))))))/(1.0 + xx*(-0.5904799485819767 +
248 xx*(0.1877833737815317 + xx*(-0.01536634911179847 +
249 xx*(0.01530712091180953 + xx*(-0.0001909176790831023))))));
260 e1 = log(1.0+1.0/y) - 0.4/
POW2(y + 1.0);
265 x = (double)t->
Hi->
nelem - 3.0;
287 gb = a + (b*y - c*y*y + d)*e1 + c*y;
325 ColStrGBar_v =
MAX2(ColStrGBar_v,1e-10);
326 return( ColStrGBar_v );