80 #if defined (HAVE_CONFIG_H)
89 #define INFINITE lo_ieee_isinf
90 #define RUNI oct_randu()
91 #define RNOR oct_randn()
92 #define REXP oct_rande()
99 const double d = (a < 1. ? 1.+a :
a) - 1./3.;
100 const double c = 1./sqrt (9.*d);
105 for (i=0; i < n; i++)
110 for (i=0; i < n; i++)
121 if (u >= 1.-0.0331*xsq*xsq &&
std::log (u) >= 0.5*xsq + d*(1-v+
std::log (v)))
129 for (i = 0; i < n; i++)
130 r[i] *= exp (-
REXP/a);
145 #define RUNI oct_float_randu()
146 #define RNOR oct_float_randn()
147 #define REXP oct_float_rande()
154 const float d = (a < 1. ? 1.+a :
a) - 1./3.;
155 const float c = 1./sqrt (9.*d);
160 for (i=0; i < n; i++)
165 for (i=0; i < n; i++)
176 if (u >= 1.-0.0331*xsq*xsq &&
std::log (u) >= 0.5*xsq + d*(1-v+
std::log (v)))
184 for (i = 0; i < n; i++)
185 r[i] *= exp (-
REXP/a);
OCTAVE_EXPORT octave_value_list or N dimensional array whose elements are all equal to the base of natural logarithms The constant ex $e satisfies the equation log(e)
void oct_fill_randg(double a, octave_idx_type n, double *r)
F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &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 F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T F77_REAL F77_REAL &F77_RET_T F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
calling an anonymous function involves an overhead quite comparable to the overhead of an m file function Passing a handle to a built in function is because the interpreter is not involved in the internal loop For a
float oct_float_randg(float a)
the sparsity preserving column transformation such that that defines the pivoting threshold can be given in which case it defines the c
double oct_randg(double a)
void oct_fill_float_randg(float a, octave_idx_type n, float *r)
=val(i)}if ode{val(i)}occurs in table i
OCTAVE_EXPORT octave_value_list or N dimensional array whose elements are all equal to the IEEE symbol NaN(Not a Number).NaN is the result of operations which do not produce a well defined 0 result.Common operations which produce a NaN are arithmetic with infinity ex($\infty-\infty $)
F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &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 F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T F77_REAL F77_REAL &F77_RET_T F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE * x