VVIObjFDetails Namespace Reference

Functions
template<typename F >
int	dzero (float a, float b, float &x0, float &rv, float eps, int mxf, F func)
	Private version of the exponential integral. More...
float	expint (float x)
	Private version of the cosine and sine integral. More...
void	sincosint (float x, float &sint, float &cint)

Function Documentation

template<typename F >

int VVIObjFDetails::dzero	(	float	a,
		float	b,
		float &	x0,
		float &	rv,
		float	eps,
		int	mxf,
		F	func
	)

Private version of the exponential integral.

Definition at line 321 of file VVIObjF.cc.

References funct::abs(), f, python.connectstrParser::f1, python.connectstrParser::f2, python.connectstrParser::f3, benchmark_cfg::fa, benchmark_cfg::fb, alignCSCRings::ff, bookConverter::max, and min().

Referenced by VVIObjF::VVIObjF().

                                               {
    /* System generated locals */
    float d__1, d__2, d__3, d__4;
    
    // Local variables
    float f1, f2, f3, u1, u2, x1, x2, u3, u4, x3, ca, cb, cc, fa, fb, ee, ff;
    int mc;
    float xa, xb, fx, xx, su4;
    
    xa = std::min(a,b);
    xb = std::max(a,b);
    fa = func(xa);
    fb = func(xb);
    if (fa * fb > 0.f) {
      rv = (xb - xa) * -2.f;
      x0 = 0.f;
      return 1;
    }
    mc = 0;
  L1:
    x0 = (xa + xb) * 0.5f;
    rv = x0 - xa;
    ee = eps * (std::abs(x0) + 1.f);
    if (rv <= ee) {
      rv = ee;
      ff = func(x0);
      return 0;
    }
    f1 = fa;
    x1 = xa;
    f2 = fb;
    x2 = xb;
  L2:
    fx = func(x0);
    ++mc;
    if (mc > mxf) {
      rv = (d__1 = xb - xa, fabs(d__1)) * -0.5f;
      x0 = 0.;
      return 0;
    }
    if (fx * fa > 0.f) {
      xa = x0;
      fa = fx;
    } else {
      xb = x0;
      fb = fx;
    }
  L3:
    u1 = f1 - f2;
    u2 = x1 - x2;
    u3 = f2 - fx;
    u4 = x2 - x0;
    if (u2 == 0.f || u4 == 0.f) {goto L1;}
    f3 = fx;
    x3 = x0;
    u1 /= u2;
    u2 = u3 / u4;
    ca = u1 - u2;
    cb = (x1 + x2) * u2 - (x2 + x0) * u1;
    cc = (x1 - x0) * f1 - x1 * (ca * x1 + cb);
    if (ca == 0.f) {
      if (cb == 0.f) {goto L1;}
      x0 = -cc / cb;
    } else {
      u3 = cb / (ca * 2.f);
      u4 = u3 * u3 - cc / ca;
      if (u4 < 0.f) {goto L1;}
      su4 = std::abs(u4);
      if (x0 + u3 < 0.f) {su4 = -su4;}
      x0 = -u3 + su4;
    }
    if (x0 < xa || x0 > xb) {goto L1;}
    // Computing MIN
    d__3 = (d__1 = x0 - x3, std::abs(d__1));
    d__4 = (d__2 = x0 - x2, std::abs(d__2));
    rv = std::min(d__3,d__4);
    ee = eps * (std::abs(x0) + 1);
    if (rv > ee) {
      f1 = f2;
      x1 = x2;
      f2 = f3;
      x2 = x3;
      goto L2;
    }
    fx = func(x0);
    if (fx == 0.f) {
      rv = ee;
      ff = func(x0);
      return 0;
    }
    if (fx * fa < 0.f) {
      xx = x0 - ee;
      if (xx <= xa) {
    rv = ee;
    ff = func(x0);
    return 0;
      }
      ff = func(xx);
      fb = ff;
      xb = xx;
    } else {
      xx = x0 + ee;
      if (xx >= xb) {
    rv = ee;
    ff = func(x0);
    return 0;
      }
      ff = func(xx);
      fa = ff;
      xa = xx;
    }
    if (fx * ff > 0.f) {
      mc += 2;
      if (mc > mxf) {
    rv = (d__1 = xb - xa, std::abs(d__1)) * -0.5f;
    x0 = 0.f;
    return 0;
      }
      f1 = f3;
      x1 = x3;
      f2 = fx;
      x2 = x0;
      x0 = xx;
      fx = ff;
      goto L3;
    }
    /* L4: */
    rv = ee;
    ff = func(x0);
    return 0;
  } // dzero

float VVIObjFDetails::expint

(

float

x

)

Private version of the cosine and sine integral.

Definition at line 198 of file VVIObjF.cc.

References myMath::fast_expf(), myMath::fast_logf(), p1, p2, p3, p4, q1, q2, findQualityFiles::v, x, and detailsBasic3DVector::y.

Referenced by VVIObjF::VVIObjF().

                      {
  
  // Initialized data
  
  const float zero = 0.;
  const float q2[7] = { .10340013040487,3.319092135933,
             20.449478501379,41.280784189142,32.426421069514,10.041164382905,
             1. };
  const float p3[6] = { -2.3909964453136,-147.98219500504,
             -254.3763397689,-119.55761038372,-19.630408535939,-.9999999999036 
  };
  const float q3[6] = { 177.60070940351,530.68509610812,
             462.23027156148,156.81843364539,21.630408494238,1. };
  const float p4[8] = { -8.6693733995107,-549.14226552109,
             -4210.0161535707,-249301.39345865,-119623.66934925,
             -22174462.775885,3892804.213112,-391546073.8091 };
  const float q4[8] = { 34.171875,-1607.0892658722,35730.029805851,
             -483547.43616216,4285596.2461175,-24903337.574054,89192576.757561,
             -165254299.72521 };
  const float a1[8] = { -2.1808638152072,-21.901023385488,
             9.3081638566217,25.076281129356,-33.184253199722,60.121799083008,
             -43.253113287813,1.0044310922808 };
  const float b1[8] = { 0.,3.9370770185272,300.89264837292,
             -6.2504116167188,1003.6743951673,14.325673812194,2736.2411988933,
             .52746885196291 };
  const float a2[8] = { -3.4833465360285,-18.65454548834,
             -8.2856199414064,-32.34673303054,17.960168876925,1.7565631546961,
             -1.9502232128966,.99999429607471 };
  const float b2[8] = { 0.,69.500065588743,57.283719383732,
             25.777638423844,760.76114800773,28.951672792514,-3.4394226689987,
             1.0008386740264 };
  const float a3[6] = { -27.780928934438,-10.10479081576,
             -9.1483008216736,-5.0223317461851,-3.0000077799358,
             1.0000000000704 };
  const float one = 1.;
  const float b3[6] = { 0.,122.39993926823,2.7276100778779,
             -7.1897518395045,-2.9990118065262,1.999999942826 };
  const float two = 2.;
  const float three = 3.;
  const float x0 = .37250741078137;
  const float xl[6] = { -24.,-12.,-6.,0.,1.,4. };
  const float p1[5] = { 4.293125234321,39.894153870321,
             292.52518866921,425.69682638592,-434.98143832952 };
  const float q1[5] = { 1.,18.899288395003,150.95038744251,
             568.05252718987,753.58564359843 };
  const float p2[7] = { .43096783946939,6.9052252278444,
             23.019255939133,24.378408879132,9.0416155694633,.99997957705159,
             4.656271079751e-7 };
  
  // Local variables 
   float v, y, ap, bp, aq, dp, bq, dq;
  
  if (x <= xl[0]) {
    ap = a3[0] - x;
    for ( int i__ = 2; i__ <= 5; ++i__) {
      ap = a3[i__ - 1] - x + b3[i__ - 1] / ap;
    }
    y = vdt::fast_expf(-x) / x * (one - (a3[5] + b3[5] / ap) / x);
  } else if (x <= xl[1]) {
    ap = a2[0] - x;
    for ( int i__ = 2; i__ <= 7; ++i__) {
      ap = a2[i__ - 1] - x + b2[i__ - 1] / ap;
    }
    y = vdt::fast_expf(-x) / x * (a2[7] + b2[7] / ap);
  } else if (x <= xl[2]) {
    ap = a1[0] - x;
    for ( int i__ = 2; i__ <= 7; ++i__) {
      ap = a1[i__ - 1] - x + b1[i__ - 1] / ap;
    }
    y = vdt::fast_expf(-x) / x * (a1[7] + b1[7] / ap);
  } else if (x < xl[3]) {
    v = -two * (x / three + one);
    bp = zero;
    dp = p4[0];
    for ( int i__ = 2; i__ <= 8; ++i__) {
      ap = bp;
      bp = dp;
      dp = p4[i__ - 1] - ap + v * bp;
    }
    bq = zero;
    dq = q4[0];
    for ( int i__ = 2; i__ <= 8; ++i__) {
      aq = bq;
      bq = dq;
      dq = q4[i__ - 1] - aq + v * bq;
    }
    y = -vdt::fast_logf(-x / x0) + (x + x0) * (dp - ap) / (dq - aq);
  } else if (x == xl[3]) {
    return zero;
  } else if (x < xl[4]) {
    ap = p1[0];
    aq = q1[0];
    for ( int i__ = 2; i__ <= 5; ++i__) {
      ap = p1[i__ - 1] + x * ap;
      aq = q1[i__ - 1] + x * aq;
    }
    y = -vdt::fast_logf(x) + ap / aq;
  } else if (x <= xl[5]) {
    y = one / x;
    ap = p2[0];
    aq = q2[0];
    for ( int i__ = 2; i__ <= 7; ++i__) {
      ap = p2[i__ - 1] + y * ap;
      aq = q2[i__ - 1] + y * aq;
    }
    y = vdt::fast_expf(-x) * ap / aq;
  } else {
    y = one / x;
    ap = p3[0];
    aq = q3[0];
    for ( int i__ = 2; i__ <= 6; ++i__) {
      ap = p3[i__ - 1] + y * ap;
      aq = q3[i__ - 1] + y * aq;
    }
    y = vdt::fast_expf(-x) * y * (one + y * ap / aq);
  }
  return y;
} // expint

void VVIObjFDetails::sincosint	(	float	x,
		float &	sint,
		float &	cint
	)

Definition at line 193 of file VVIObjF.cc.

References sicif().

Referenced by VVIObjF::VVIObjF().

                                                      {
    sicif(x,sint,cint);
  }