382 void cmatrixrighttrsm(
const ae_int_t m,
const ae_int_t n,
const complex_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const complex_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
383 void smp_cmatrixrighttrsm(
const ae_int_t m,
const ae_int_t n,
const complex_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const complex_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
389 void cmatrixlefttrsm(
const ae_int_t m,
const ae_int_t n,
const complex_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const complex_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
390 void smp_cmatrixlefttrsm(
const ae_int_t m,
const ae_int_t n,
const complex_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const complex_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
396 void rmatrixrighttrsm(
const ae_int_t m,
const ae_int_t n,
const real_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const real_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
397 void smp_rmatrixrighttrsm(
const ae_int_t m,
const ae_int_t n,
const real_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const real_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
403 void rmatrixlefttrsm(
const ae_int_t m,
const ae_int_t n,
const real_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const real_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
404 void smp_rmatrixlefttrsm(
const ae_int_t m,
const ae_int_t n,
const real_2d_array &a,
const ae_int_t i1,
const ae_int_t j1,
const bool isupper,
const bool isunit,
const ae_int_t optype,
const real_2d_array &x,
const ae_int_t i2,
const ae_int_t j2);
410 void cmatrixsyrk(
const ae_int_t n,
const ae_int_t k,
const double alpha,
const complex_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const double beta,
const complex_2d_array &
c,
const ae_int_t ic,
const ae_int_t jc,
const bool isupper);
411 void smp_cmatrixsyrk(
const ae_int_t n,
const ae_int_t k,
const double alpha,
const complex_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const double beta,
const complex_2d_array &c,
const ae_int_t ic,
const ae_int_t jc,
const bool isupper);
417 void rmatrixsyrk(
const ae_int_t n,
const ae_int_t k,
const double alpha,
const real_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const double beta,
const real_2d_array &c,
const ae_int_t ic,
const ae_int_t jc,
const bool isupper);
418 void smp_rmatrixsyrk(
const ae_int_t n,
const ae_int_t k,
const double alpha,
const real_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const double beta,
const real_2d_array &c,
const ae_int_t ic,
const ae_int_t jc,
const bool isupper);
424 void cmatrixgemm(
const ae_int_t m,
const ae_int_t n,
const ae_int_t k,
const alglib::complex alpha,
const complex_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const complex_2d_array &b,
const ae_int_t ib,
const ae_int_t jb,
const ae_int_t optypeb,
const alglib::complex beta,
const complex_2d_array &c,
const ae_int_t ic,
const ae_int_t jc);
425 void smp_cmatrixgemm(
const ae_int_t m,
const ae_int_t n,
const ae_int_t k,
const alglib::complex alpha,
const complex_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const complex_2d_array &b,
const ae_int_t ib,
const ae_int_t jb,
const ae_int_t optypeb,
const alglib::complex beta,
const complex_2d_array &c,
const ae_int_t ic,
const ae_int_t jc);
431 void rmatrixgemm(
const ae_int_t m,
const ae_int_t n,
const ae_int_t k,
const double alpha,
const real_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const real_2d_array &b,
const ae_int_t ib,
const ae_int_t jb,
const ae_int_t optypeb,
const double beta,
const real_2d_array &c,
const ae_int_t ic,
const ae_int_t jc);
432 void smp_rmatrixgemm(
const ae_int_t m,
const ae_int_t n,
const ae_int_t k,
const double alpha,
const real_2d_array &a,
const ae_int_t ia,
const ae_int_t ja,
const ae_int_t optypea,
const real_2d_array &b,
const ae_int_t ib,
const ae_int_t jb,
const ae_int_t optypeb,
const double beta,
const real_2d_array &c,
const ae_int_t ic,
const ae_int_t jc);
1268 bool rmatrixbdsvd(
real_1d_array &d,
const real_1d_array &e,
const ae_int_t n,
const bool isupper,
const bool isfractionalaccuracyrequired,
real_2d_array &u,
const ae_int_t nru,
real_2d_array &c,
const ae_int_t ncc,
real_2d_array &vt,
const ae_int_t ncvt);
4630 ae_bool isfractionalaccuracyrequired,
4642 ae_bool isfractionalaccuracyrequired,
5082 fblslincgstate* state,
void smp_rmatrixlefttrsm(const ae_int_t m, const ae_int_t n, const real_2d_array &a, const ae_int_t i1, const ae_int_t j1, const bool isupper, const bool isunit, const ae_int_t optype, const real_2d_array &x, const ae_int_t i2, const ae_int_t j2)
void hpdmatrixrndcond(ae_int_t n, double c, ae_matrix *a, ae_state *_state)
void cmatrixlu(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *pivots, ae_state *_state)
ae_bool spdmatrixcholesky(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
void smp_cmatrixgemm(const ae_int_t m, const ae_int_t n, const ae_int_t k, const alglib::complex alpha, const complex_2d_array &a, const ae_int_t ia, const ae_int_t ja, const ae_int_t optypea, const complex_2d_array &b, const ae_int_t ib, const ae_int_t jb, const ae_int_t optypeb, const alglib::complex beta, const complex_2d_array &c, const ae_int_t ic, const ae_int_t jc)
struct alglib_impl::ae_state ae_state
double beta(double a, double b, ae_state *_state)
double rmatrixrcond1(ae_matrix *a, ae_int_t n, ae_state *_state)
double rmatrixrcondinf(ae_matrix *a, ae_int_t n, ae_state *_state)
ae_int_t ablasblocksize(ae_matrix *a, ae_state *_state)
double sparseget(sparsematrix *s, ae_int_t i, ae_int_t j, ae_state *_state)
double cmatrixtrrcondinf(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_bool isunit, ae_state *_state)
void smp_rmatrixgemm(const ae_int_t m, const ae_int_t n, const ae_int_t k, const double alpha, const real_2d_array &a, const ae_int_t ia, const ae_int_t ja, const ae_int_t optypea, const real_2d_array &b, const ae_int_t ib, const ae_int_t jb, const ae_int_t optypeb, const double beta, const real_2d_array &c, const ae_int_t ic, const ae_int_t jc)
alglib_impl::matinvreport * p_struct
void rmatrixlqunpackq(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_int_t qrows, ae_matrix *q, ae_state *_state)
ae_bool _fblslincgstate_init(void *_p, ae_state *_state, ae_bool make_automatic)
ae_bool bidiagonalsvddecomposition(ae_vector *d, ae_vector *e, ae_int_t n, ae_bool isupper, ae_bool isfractionalaccuracyrequired, ae_matrix *u, ae_int_t nru, ae_matrix *c, ae_int_t ncc, ae_matrix *vt, ae_int_t ncvt, ae_state *_state)
double spdmatrixcholeskydet(ae_matrix *a, ae_int_t n, ae_state *_state)
void sparsecopy(sparsematrix *s0, sparsematrix *s1, ae_state *_state)
void ablascomplexsplitlength(ae_matrix *a, ae_int_t n, ae_int_t *n1, ae_int_t *n2, ae_state *_state)
ae_bool normestimatoriteration(normestimatorstate *state, ae_state *_state)
ae_int_t sparsegetnrows(sparsematrix *s, ae_state *_state)
void rmatrixqr(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_state *_state)
void rmatrixgemm(ae_int_t m, ae_int_t n, ae_int_t k, double alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_int_t optypeb, double beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_state *_state)
void cmatrixlq(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_state *_state)
void sparsemv(sparsematrix *s, ae_vector *x, ae_vector *y, ae_state *_state)
void rmatrixbdmultiplybyp(ae_matrix *qp, ae_int_t m, ae_int_t n, ae_vector *taup, ae_matrix *z, ae_int_t zrows, ae_int_t zcolumns, ae_bool fromtheright, ae_bool dotranspose, ae_state *_state)
void rmatrixbdunpackdiagonals(ae_matrix *b, ae_int_t m, ae_int_t n, ae_bool *isupper, ae_vector *d, ae_vector *e, ae_state *_state)
void hmatrixtdunpackq(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_vector *tau, ae_matrix *q, ae_state *_state)
void smp_cmatrixsyrk(const ae_int_t n, const ae_int_t k, const double alpha, const complex_2d_array &a, const ae_int_t ia, const ae_int_t ja, const ae_int_t optypea, const double beta, const complex_2d_array &c, const ae_int_t ic, const ae_int_t jc, const bool isupper)
void sparseresizematrix(sparsematrix *s, ae_state *_state)
void normestimatorsetseed(normestimatorstate *state, ae_int_t seedval, ae_state *_state)
ae_bool smatrixevdr(ae_matrix *a, ae_int_t n, ae_int_t zneeded, ae_bool isupper, double b1, double b2, ae_int_t *m, ae_vector *w, ae_matrix *z, ae_state *_state)
void _normestimatorstate_clear(void *_p)
ae_bool hmatrixevdi(ae_matrix *a, ae_int_t n, ae_int_t zneeded, ae_bool isupper, ae_int_t i1, ae_int_t i2, ae_vector *w, ae_matrix *z, ae_state *_state)
ae_bool _sparsematrix_init(void *_p, ae_state *_state, ae_bool make_automatic)
void _sparsematrix_clear(void *_p)
void cmatrixrighttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
void _pexec_rmatrixlefttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
void normestimatorresults(normestimatorstate *state, double *nrm, ae_state *_state)
void _pexec_rmatrixsyrk(ae_int_t n, ae_int_t k, double alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, double beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_bool isupper, ae_state *_state)
ae_bool hpdmatrixcholesky(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
ae_int_t ablascomplexblocksize(ae_matrix *a, ae_state *_state)
ae_bool smatrixgevdreduce(ae_matrix *a, ae_int_t n, ae_bool isuppera, ae_matrix *b, ae_bool isupperb, ae_int_t problemtype, ae_matrix *r, ae_bool *isupperr, ae_state *_state)
alglib_impl::sparsematrix * p_struct
void smatrixrndmultiply(ae_matrix *a, ae_int_t n, ae_state *_state)
ae_bool rmatrixschur(ae_matrix *a, ae_int_t n, ae_matrix *s, ae_state *_state)
double rmatrixdet(ae_matrix *a, ae_int_t n, ae_state *_state)
void rmatrixlup(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *pivots, ae_state *_state)
void cmatrixlqunpackq(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_int_t qrows, ae_matrix *q, ae_state *_state)
void rmatrixlefttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
void _pexec_rmatrixgemm(ae_int_t m, ae_int_t n, ae_int_t k, double alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_int_t optypeb, double beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_state *_state)
void sparsemtv(sparsematrix *s, ae_vector *x, ae_vector *y, ae_state *_state)
double cmatrixlurcond1(ae_matrix *lua, ae_int_t n, ae_state *_state)
void cmatrixqr(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_state *_state)
ae_bool _matinvreport_init_copy(void *_dst, void *_src, ae_state *_state, ae_bool make_automatic)
void cmatrixmv(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t opa, ae_vector *x, ae_int_t ix, ae_vector *y, ae_int_t iy, ae_state *_state)
void rmatrixinvupdaterow(ae_matrix *inva, ae_int_t n, ae_int_t updrow, ae_vector *v, ae_state *_state)
void sparsesmm(sparsematrix *s, ae_bool isupper, ae_matrix *a, ae_int_t k, ae_matrix *b, ae_state *_state)
double rmatrixludet(ae_matrix *a, ae_vector *pivots, ae_int_t n, ae_state *_state)
void sparsemtm(sparsematrix *s, ae_matrix *a, ae_int_t k, ae_matrix *b, ae_state *_state)
void rmatrixhessenberg(ae_matrix *a, ae_int_t n, ae_vector *tau, ae_state *_state)
ae_bool sparserewriteexisting(sparsematrix *s, ae_int_t i, ae_int_t j, double v, ae_state *_state)
void sparseconverttocrs(sparsematrix *s, ae_state *_state)
double spdmatrixdet(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
void cmatrixrndcond(ae_int_t n, double c, ae_matrix *a, ae_state *_state)
void _fblslincgstate_destroy(void *_p)
ae_int_t sparsegetmatrixtype(sparsematrix *s, ae_state *_state)
void rmatrixenforcesymmetricity(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
void cmatrixplu(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *pivots, ae_state *_state)
void fblssolvels(ae_matrix *a, ae_vector *b, ae_int_t m, ae_int_t n, ae_vector *tmp0, ae_vector *tmp1, ae_vector *tmp2, ae_state *_state)
ql0001_ & k(htemp+1),(cvec+1),(atemp+1),(bj+1),(bl+1),(bu+1),(x+1),(clamda+1), &iout, infoqp, &zero,(w+1), &lenw,(iw+1), &leniw, &glob_grd.epsmac
ae_complex cmatrixdet(ae_matrix *a, ae_int_t n, ae_state *_state)
void cmatrixrndorthogonalfromtheright(ae_matrix *a, ae_int_t m, ae_int_t n, ae_state *_state)
void sparseconverttohash(sparsematrix *s, ae_state *_state)
double hpdmatrixrcond(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
void cmatrixlqunpackl(ae_matrix *a, ae_int_t m, ae_int_t n, ae_matrix *l, ae_state *_state)
void rmatrixsyrk(ae_int_t n, ae_int_t k, double alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, double beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_bool isupper, ae_state *_state)
void spdmatrixrndcond(ae_int_t n, double c, ae_matrix *a, ae_state *_state)
void rmatrixrndorthogonal(ae_int_t n, ae_matrix *a, ae_state *_state)
ae_bool sparseishash(sparsematrix *s, ae_state *_state)
void hmatrixrndcond(ae_int_t n, double c, ae_matrix *a, ae_state *_state)
void normestimatorrestart(normestimatorstate *state, ae_state *_state)
void rmatrixrndcond(ae_int_t n, double c, ae_matrix *a, ae_state *_state)
void sparseset(sparsematrix *s, ae_int_t i, ae_int_t j, double v, ae_state *_state)
void rmatrixrndorthogonalfromtheleft(ae_matrix *a, ae_int_t m, ae_int_t n, ae_state *_state)
void rmatrixlqunpackl(ae_matrix *a, ae_int_t m, ae_int_t n, ae_matrix *l, ae_state *_state)
void rmatrixtrinverse(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_bool isunit, ae_int_t *info, matinvreport *rep, ae_state *_state)
void cmatrixqrunpackq(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_int_t qcolumns, ae_matrix *q, ae_state *_state)
ae_bool _normestimatorstate_init_copy(void *_dst, void *_src, ae_state *_state, ae_bool make_automatic)
void hpdmatrixinverse(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_int_t *info, matinvreport *rep, ae_state *_state)
void smp_cmatrixlefttrsm(const ae_int_t m, const ae_int_t n, const complex_2d_array &a, const ae_int_t i1, const ae_int_t j1, const bool isupper, const bool isunit, const ae_int_t optype, const complex_2d_array &x, const ae_int_t i2, const ae_int_t j2)
void rmatrixcopy(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_state *_state)
void sparsemm(sparsematrix *s, ae_matrix *a, ae_int_t k, ae_matrix *b, ae_state *_state)
double cmatrixtrrcond1(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_bool isunit, ae_state *_state)
void rmatrixbdunpackq(ae_matrix *qp, ae_int_t m, ae_int_t n, ae_vector *tauq, ae_int_t qcolumns, ae_matrix *q, ae_state *_state)
void smatrixrndcond(ae_int_t n, double c, ae_matrix *a, ae_state *_state)
void spdmatrixinverse(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_int_t *info, matinvreport *rep, ae_state *_state)
void cmatrixtrinverse(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_bool isunit, ae_int_t *info, matinvreport *rep, ae_state *_state)
void sparsemv2(sparsematrix *s, ae_vector *x, ae_vector *y0, ae_vector *y1, ae_state *_state)
void rmatrixplu(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *pivots, ae_state *_state)
double cmatrixrcondinf(ae_matrix *a, ae_int_t n, ae_state *_state)
void rmatrixinvupdatesimple(ae_matrix *inva, ae_int_t n, ae_int_t updrow, ae_int_t updcolumn, double updval, ae_state *_state)
void hmatrixtd(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_vector *tau, ae_vector *d, ae_vector *e, ae_state *_state)
void sparsegetrow(sparsematrix *s, ae_int_t i, ae_vector *irow, ae_state *_state)
void _pexec_cmatrixrighttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
void _pexec_cmatrixgemm(ae_int_t m, ae_int_t n, ae_int_t k, ae_complex alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_int_t optypeb, ae_complex beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_state *_state)
void sparsecreatecrs(ae_int_t m, ae_int_t n, ae_vector *ner, sparsematrix *s, ae_state *_state)
void rmatrixhessenbergunpackh(ae_matrix *a, ae_int_t n, ae_matrix *h, ae_state *_state)
void rmatrixinvupdateuv(ae_matrix *inva, ae_int_t n, ae_vector *u, ae_vector *v, ae_state *_state)
void rmatrixlu(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *pivots, ae_state *_state)
void _pexec_cmatrixlefttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
void rmatrixtranspose(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_state *_state)
double rmatrixtrrcond1(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_bool isunit, ae_state *_state)
void cmatrixrank1(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_vector *u, ae_int_t iu, ae_vector *v, ae_int_t iv, ae_state *_state)
void cmatrixinverse(ae_matrix *a, ae_int_t n, ae_int_t *info, matinvreport *rep, ae_state *_state)
void hpdmatrixcholeskyinverse(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_int_t *info, matinvreport *rep, ae_state *_state)
void sparsecopytohash(sparsematrix *s0, sparsematrix *s1, ae_state *_state)
void rmatrixinverse(ae_matrix *a, ae_int_t n, ae_int_t *info, matinvreport *rep, ae_state *_state)
void sparseadd(sparsematrix *s, ae_int_t i, ae_int_t j, double v, ae_state *_state)
void rmatrixinvupdatecolumn(ae_matrix *inva, ae_int_t n, ae_int_t updcolumn, ae_vector *u, ae_state *_state)
void rmatrixrank1(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_vector *u, ae_int_t iu, ae_vector *v, ae_int_t iv, ae_state *_state)
ae_bool fblscgiteration(fblslincgstate *state, ae_state *_state)
void rmatrixbdmultiplybyq(ae_matrix *qp, ae_int_t m, ae_int_t n, ae_vector *tauq, ae_matrix *z, ae_int_t zrows, ae_int_t zcolumns, ae_bool fromtheright, ae_bool dotranspose, ae_state *_state)
void fblscgcreate(ae_vector *x, ae_vector *b, ae_int_t n, fblslincgstate *state, ae_state *_state)
void rmatrixlq(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_state *_state)
ae_bool _fblslincgstate_init_copy(void *_dst, void *_src, ae_state *_state, ae_bool make_automatic)
ae_bool sparseiscrs(sparsematrix *s, ae_state *_state)
ae_bool sparseenumerate(sparsematrix *s, ae_int_t *t0, ae_int_t *t1, ae_int_t *i, ae_int_t *j, double *v, ae_state *_state)
double rcondthreshold(ae_state *_state)
void _fblslincgstate_clear(void *_p)
double cmatrixlurcondinf(ae_matrix *lua, ae_int_t n, ae_state *_state)
ae_int_t ablasmicroblocksize(ae_state *_state)
ae_bool hmatrixevd(ae_matrix *a, ae_int_t n, ae_int_t zneeded, ae_bool isupper, ae_vector *d, ae_matrix *z, ae_state *_state)
ae_bool smatrixgevd(ae_matrix *a, ae_int_t n, ae_bool isuppera, ae_matrix *b, ae_bool isupperb, ae_int_t zneeded, ae_int_t problemtype, ae_vector *d, ae_matrix *z, ae_state *_state)
void smp_rmatrixrighttrsm(const ae_int_t m, const ae_int_t n, const real_2d_array &a, const ae_int_t i1, const ae_int_t j1, const bool isupper, const bool isunit, const ae_int_t optype, const real_2d_array &x, const ae_int_t i2, const ae_int_t j2)
ae_int_t sparsegetncols(sparsematrix *s, ae_state *_state)
void _matinvreport_clear(void *_p)
void normestimatorcreate(ae_int_t m, ae_int_t n, ae_int_t nstart, ae_int_t nits, normestimatorstate *state, ae_state *_state)
ae_bool smatrixtdevdi(ae_vector *d, ae_vector *e, ae_int_t n, ae_int_t zneeded, ae_int_t i1, ae_int_t i2, ae_matrix *z, ae_state *_state)
double rmatrixlurcond1(ae_matrix *lua, ae_int_t n, ae_state *_state)
ae_complex cmatrixludet(ae_matrix *a, ae_vector *pivots, ae_int_t n, ae_state *_state)
double rmatrixtrrcondinf(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_bool isunit, ae_state *_state)
struct alglib_impl::ae_vector ae_vector
void ablassplitlength(ae_matrix *a, ae_int_t n, ae_int_t *n1, ae_int_t *n2, ae_state *_state)
double sparsegetaveragelengthofchain(sparsematrix *s, ae_state *_state)
void rmatrixmv(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t opa, ae_vector *x, ae_int_t ix, ae_vector *y, ae_int_t iy, ae_state *_state)
void rmatrixqrunpackr(ae_matrix *a, ae_int_t m, ae_int_t n, ae_matrix *r, ae_state *_state)
void smp_rmatrixsyrk(const ae_int_t n, const ae_int_t k, const double alpha, const real_2d_array &a, const ae_int_t ia, const ae_int_t ja, const ae_int_t optypea, const double beta, const real_2d_array &c, const ae_int_t ic, const ae_int_t jc, const bool isupper)
double spdmatrixrcond(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
void _matinvreport_destroy(void *_p)
ae_bool rmatrixevd(ae_matrix *a, ae_int_t n, ae_int_t vneeded, ae_vector *wr, ae_vector *wi, ae_matrix *vl, ae_matrix *vr, ae_state *_state)
void cmatrixgemm(ae_int_t m, ae_int_t n, ae_int_t k, ae_complex alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_int_t optypeb, ae_complex beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_state *_state)
void rmatrixqrunpackq(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tau, ae_int_t qcolumns, ae_matrix *q, ae_state *_state)
double cmatrixrcond1(ae_matrix *a, ae_int_t n, ae_state *_state)
struct alglib_impl::ae_matrix ae_matrix
ae_bool rmatrixbdsvd(ae_vector *d, ae_vector *e, ae_int_t n, ae_bool isupper, ae_bool isfractionalaccuracyrequired, ae_matrix *u, ae_int_t nru, ae_matrix *c, ae_int_t ncc, ae_matrix *vt, ae_int_t ncvt, ae_state *_state)
ae_bool smatrixtdevd(ae_vector *d, ae_vector *e, ae_int_t n, ae_int_t zneeded, ae_matrix *z, ae_state *_state)
void cmatrixlefttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
ae_bool rmatrixsvd(ae_matrix *a, ae_int_t m, ae_int_t n, ae_int_t uneeded, ae_int_t vtneeded, ae_int_t additionalmemory, ae_vector *w, ae_matrix *u, ae_matrix *vt, ae_state *_state)
void rmatrixhessenbergunpackq(ae_matrix *a, ae_int_t n, ae_vector *tau, ae_matrix *q, ae_state *_state)
void cmatrixsyrk(ae_int_t n, ae_int_t k, double alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, double beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_bool isupper, ae_state *_state)
void rmatrixluinverse(ae_matrix *a, ae_vector *pivots, ae_int_t n, ae_int_t *info, matinvreport *rep, ae_state *_state)
double spdmatrixcholeskyrcond(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
void cmatrixtranspose(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_state *_state)
ae_bool smatrixevdi(ae_matrix *a, ae_int_t n, ae_int_t zneeded, ae_bool isupper, ae_int_t i1, ae_int_t i2, ae_vector *w, ae_matrix *z, ae_state *_state)
ae_bool spdmatrixcholeskyrec(ae_matrix *a, ae_int_t offs, ae_int_t n, ae_bool isupper, ae_vector *tmp, ae_state *_state)
void cmatrixrndorthogonalfromtheleft(ae_matrix *a, ae_int_t m, ae_int_t n, ae_state *_state)
ae_bool hmatrixevdr(ae_matrix *a, ae_int_t n, ae_int_t zneeded, ae_bool isupper, double b1, double b2, ae_int_t *m, ae_vector *w, ae_matrix *z, ae_state *_state)
void rmatrixbdunpackpt(ae_matrix *qp, ae_int_t m, ae_int_t n, ae_vector *taup, ae_int_t ptrows, ae_matrix *pt, ae_state *_state)
void cmatrixqrunpackr(ae_matrix *a, ae_int_t m, ae_int_t n, ae_matrix *r, ae_state *_state)
void _pexec_rmatrixrighttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
void _normestimatorstate_destroy(void *_p)
void sparsefree(sparsematrix *s, ae_state *_state)
void rmatrixlqbasecase(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *work, ae_vector *t, ae_vector *tau, ae_state *_state)
void cmatrixrndorthogonal(ae_int_t n, ae_matrix *a, ae_state *_state)
void sparsecreate(ae_int_t m, ae_int_t n, ae_int_t k, sparsematrix *s, ae_state *_state)
ae_bool _normestimatorstate_init(void *_p, ae_state *_state, ae_bool make_automatic)
void rmatrixbd(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *tauq, ae_vector *taup, ae_state *_state)
void _sparsematrix_destroy(void *_p)
ae_bool smatrixtdevdr(ae_vector *d, ae_vector *e, ae_int_t n, ae_int_t zneeded, double a, double b, ae_int_t *m, ae_matrix *z, ae_state *_state)
void smp_cmatrixrighttrsm(const ae_int_t m, const ae_int_t n, const complex_2d_array &a, const ae_int_t i1, const ae_int_t j1, const bool isupper, const bool isunit, const ae_int_t optype, const complex_2d_array &x, const ae_int_t i2, const ae_int_t j2)
void cmatrixlup(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *pivots, ae_state *_state)
double rmatrixlurcondinf(ae_matrix *lua, ae_int_t n, ae_state *_state)
void cmatrixluinverse(ae_matrix *a, ae_vector *pivots, ae_int_t n, ae_int_t *info, matinvreport *rep, ae_state *_state)
void rmatrixqrbasecase(ae_matrix *a, ae_int_t m, ae_int_t n, ae_vector *work, ae_vector *t, ae_vector *tau, ae_state *_state)
void smatrixtd(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_vector *tau, ae_vector *d, ae_vector *e, ae_state *_state)
alglib_impl::ae_int_t ae_int_t
void sparsesmv(sparsematrix *s, ae_bool isupper, ae_vector *x, ae_vector *y, ae_state *_state)
void _pexec_cmatrixsyrk(ae_int_t n, ae_int_t k, double alpha, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_int_t optypea, double beta, ae_matrix *c, ae_int_t ic, ae_int_t jc, ae_bool isupper, ae_state *_state)
ae_bool _sparsematrix_init_copy(void *_dst, void *_src, ae_state *_state, ae_bool make_automatic)
void hmatrixrndmultiply(ae_matrix *a, ae_int_t n, ae_state *_state)
double hpdmatrixcholeskyrcond(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_state *_state)
void fblscholeskysolve(ae_matrix *cha, double sqrtscalea, ae_int_t n, ae_bool isupper, ae_vector *xb, ae_vector *tmp, ae_state *_state)
double sparsegetdiagonal(sparsematrix *s, ae_int_t i, ae_state *_state)
void spdmatrixcholeskyinverse(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_int_t *info, matinvreport *rep, ae_state *_state)
void sparsecopytocrs(sparsematrix *s0, sparsematrix *s1, ae_state *_state)
void normestimatorestimatesparse(normestimatorstate *state, sparsematrix *a, ae_state *_state)
void smatrixtdunpackq(ae_matrix *a, ae_int_t n, ae_bool isupper, ae_vector *tau, ae_matrix *q, ae_state *_state)
void rmatrixrighttrsm(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t i1, ae_int_t j1, ae_bool isupper, ae_bool isunit, ae_int_t optype, ae_matrix *x, ae_int_t i2, ae_int_t j2, ae_state *_state)
ae_bool smatrixevd(ae_matrix *a, ae_int_t n, ae_int_t zneeded, ae_bool isupper, ae_vector *d, ae_matrix *z, ae_state *_state)
void fblssolvecgx(ae_matrix *a, ae_int_t m, ae_int_t n, double alpha, ae_vector *b, ae_vector *x, ae_vector *buf, ae_state *_state)
ae_bool _matinvreport_init(void *_p, ae_state *_state, ae_bool make_automatic)
void cmatrixcopy(ae_int_t m, ae_int_t n, ae_matrix *a, ae_int_t ia, ae_int_t ja, ae_matrix *b, ae_int_t ib, ae_int_t jb, ae_state *_state)
alglib_impl::normestimatorstate * p_struct
void rmatrixrndorthogonalfromtheright(ae_matrix *a, ae_int_t m, ae_int_t n, ae_state *_state)
void sparsemm2(sparsematrix *s, ae_matrix *a, ae_int_t k, ae_matrix *b0, ae_matrix *b1, ae_state *_state)