23 #if defined (HAVE_CONFIG_H) 31 DEFUN (__expint__, args, ,
37 int nargin = args.length ();
46 int numel_x = args(0).numel ();
61 x = args(0).float_complex_array_value ();
65 static const float eps = std::numeric_limits<float>::epsilon ();
68 const int maxit = 100;
87 while ((
std::abs (Deltaj - cone) >
eps) && (j < maxit))
89 Dj = beta_j + alpha_j * Dj;
92 Cj = beta_j + alpha_j / Cj;
98 alpha_j = (j + 1) / 2;
119 x = args(0).complex_array_value ();
123 static const double eps = std::numeric_limits<double>::epsilon ();
125 const Complex czero (0.0, 0.0);
126 const int maxit = 200;
145 while ((
std::abs (Deltaj - cone) >
eps) && (j < maxit))
147 Dj = beta_j + alpha_j * Dj;
150 Cj = beta_j + alpha_j / Cj;
156 alpha_j = (j + 1) / 2;
OCTINTERP_API void print_usage(void)
F77_RET_T const F77_REAL const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE * f
#define DEFUN(name, args_name, nargout_name, doc)
Macro to define a builtin function.
bool is_single_type(void) const
the element is set to zero In other the statement xample y
std::complex< float > FloatComplex
std::complex< double > Complex
Vector representing the dimensions (size) of an Array.
F77_RET_T const F77_REAL const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T const F77_DBLE * x