24 #if ! defined (octave_chol_h) 25 #define octave_chol_h 1 27 #include "octave-config.h" 39 typedef typename T::column_vector_type
VT;
40 typedef typename T::real_elt_type
COND_T;
42 chol (
void) : chol_mat (), xrcond (0) { }
44 chol (
const T&
a,
bool upper =
true,
bool calc_cond =
false)
45 : chol_mat (), xrcond (0)
47 init (
a, upper, calc_cond);
51 bool calc_cond =
false)
52 : chol_mat (), xrcond (0)
54 info = init (
a, upper, calc_cond);
58 : chol_mat (
a.chol_mat), xrcond (
a.xrcond) { }
64 chol_mat =
a.chol_mat;
76 T inverse (
void)
const;
78 void set (
const T& R);
80 void update (
const VT&
u);
101 template <
typename T>
chol(const T &a, bool upper=true, bool calc_cond=false)
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
static void downdate(double *c, int n, int d, int *nans, int nnans)
chol(const T &a, octave_idx_type &info, bool upper=true, bool calc_cond=false)
T chol_matrix(void) const