pyexpint/expint.cpp
/**
 * Algorithm for the generalized exponential integral with double nu
 */
Algorithm for the generalized exponential integral with double nu
Algorithm for the generalized exponential integral with double nu
double expint_v(const double v, const double x)
{
  int iflag = 0, vint;
  double result = 0.0;
  vint = int(v);
  if (v == vint)
    return expint_n(vint, x);
  expint_v_special_cases(v, x, &iflag, &result);
  if (iflag) return result;
  if (x <= 1.0)
    return expint_small_x(v, x);
  else if (v >= x)
    return expint_large_v(v, x);
  else
    return expint_large_x(v, x);
}
double expint_v ( const double v , const double x ) { int iflag = 0 , vint ; double result = 0.0 ; vint = int ( v ) ; if ( v == vint ) return expint_n ( vint , x ) ; expint_v_special_cases ( v , x , & iflag , & result ) ; if ( iflag ) return result ; if ( x <= 1.0 ) return expint_small_x ( v , x ) ; else if ( v >= x ) return expint_large_v ( v , x ) ; else return expint_large_x ( v , x ) ; }
Algorithm for the generalized exponential integral with double nu
Algorithm for the generalized exponential integral with double nu
// v is integer // special cases // general cases