47 #define AE_LOCK_ALIGNMENT 16 51 #define AE_OS AE_UNKNOWN 55 #define AE_COMPILER AE_UNKNOWN 58 #define AE_COMPILER AE_GNUC 60 #if defined(__SUNPRO_C)||defined(__SUNPRO_CC) 62 #define AE_COMPILER AE_SUNC 66 #define AE_COMPILER AE_MSVC 87 #if defined(AE_HAVE_STDINT) 110 #if AE_COMPILER==AE_MSVC 111 #include <emmintrin.h> 112 #define AE_HAS_SSE2_INTRINSICS 114 #if AE_COMPILER==AE_GNUC 115 #include <xmmintrin.h> 116 #define AE_HAS_SSE2_INTRINSICS 118 #if AE_COMPILER==AE_SUNC 119 #include <xmmintrin.h> 120 #include <emmintrin.h> 121 #define AE_HAS_SSE2_INTRINSICS 127 #ifdef AE_DEBUG4WINDOWS 146 #define AE_USE_CPP_BOOL 147 #define AE_USE_CPP_ERROR_HANDLING 148 #define AE_USE_CPP_SERIALIZATION 155 #if defined(AE_INT32_T) 158 #if defined(AE_HAVE_STDINT) && !defined(AE_INT32_T) 161 #if !defined(AE_HAVE_STDINT) && !defined(AE_INT32_T) 162 #if AE_COMPILER==AE_MSVC 165 #if (AE_COMPILER==AE_GNUC) || (AE_COMPILER==AE_SUNC) || (AE_COMPILER==AE_UNKNOWN) 170 #if defined(AE_INT64_T) 173 #if defined(AE_HAVE_STDINT) && !defined(AE_INT64_T) 176 #if !defined(AE_HAVE_STDINT) && !defined(AE_INT64_T) 177 #if AE_COMPILER==AE_MSVC 180 #if (AE_COMPILER==AE_GNUC) || (AE_COMPILER==AE_SUNC) || (AE_COMPILER==AE_UNKNOWN) 185 #if !defined(AE_INT_T) 189 #if !defined(AE_USE_CPP_BOOL) 196 #define ae_false false 377 #ifndef AE_USE_CPP_ERROR_HANDLING 378 jmp_buf *
volatile break_jump;
403 void (*thread_exception_handler)(
void*);
418 #ifdef AE_USE_CPP_SERIALIZATION 476 void (*destroy)(
void*);
498 #if AE_OS==AE_WINDOWS 499 volatile ae_int_t *
volatile p_lock;
501 #elif AE_OS==AE_POSIX 502 pthread_mutex_t
mutex;
557 void (*destroy)(
void*
ptr);
570 ae_int_t
ae_sizeof(ae_datatype datatype);
575 #ifndef AE_USE_CPP_ERROR_HANDLING 576 void ae_state_set_break_jump(
ae_state *state, jmp_buf *buf);
625 ae_int_t size_of_object,
628 void (*destroy)(
void*
ptr),
676 #ifdef AE_USE_CPP_SERIALIZATION 804 double ae_v_dotproduct(
const double *
v0, ae_int_t stride0,
const double *
v1, ae_int_t stride1, ae_int_t
n);
805 void ae_v_move(
double *vdst, ae_int_t stride_dst,
const double* vsrc, ae_int_t stride_src, ae_int_t
n);
806 void ae_v_moveneg(
double *vdst, ae_int_t stride_dst,
const double* vsrc, ae_int_t stride_src, ae_int_t
n);
807 void ae_v_moved(
double *vdst, ae_int_t stride_dst,
const double* vsrc, ae_int_t stride_src, ae_int_t
n,
double alpha);
808 void ae_v_add(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n);
809 void ae_v_addd(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n,
double alpha);
810 void ae_v_sub(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n);
811 void ae_v_subd(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n,
double alpha);
812 void ae_v_muld(
double *vdst, ae_int_t stride_dst, ae_int_t
n,
double alpha);
825 #define ae_machineepsilon 5E-16 826 #define ae_maxrealnumber 1E300 827 #define ae_minrealnumber 1E-300 828 #define ae_pi 3.1415926535897932384626433832795 847 #ifdef AE_USE_ALLOC_COUNTER 848 extern ae_int64_t _alloc_counter;
860 #ifdef AE_DEBUG4WINDOWS 861 #define flushconsole(s) fflush(stdout) 862 #define tickcount(s) _tickcount() 865 #ifdef AE_DEBUG4POSIX 866 #define flushconsole(s) fflush(stdout) 867 #define tickcount(s) _tickcount() 871 ae_int_t ae_debugrng();
872 void ae_set_seed(ae_int_t s0, ae_int_t s1);
873 void ae_get_seed(ae_int_t *s0, ae_int_t *s1);
896 ae_int_t
vlen(ae_int_t n1, ae_int_t n2);
908 static void make_assertion(
bool bClause);
909 static void make_assertion(
bool bClause,
const char *msg);
921 complex(
const double &_x,
const double &_y);
924 complex& operator= (
const double& v);
939 std::string tostring(
int dps)
const;
979 double vdotproduct(
const double *
v0, ae_int_t stride0,
const double *
v1, ae_int_t stride1, ae_int_t
n);
980 double vdotproduct(
const double *
v1,
const double *v2, ae_int_t N);
985 void vmove(
double *vdst, ae_int_t stride_dst,
const double* vsrc, ae_int_t stride_src, ae_int_t
n);
986 void vmove(
double *vdst,
const double* vsrc, ae_int_t N);
991 void vmoveneg(
double *vdst, ae_int_t stride_dst,
const double* vsrc, ae_int_t stride_src, ae_int_t
n);
992 void vmoveneg(
double *vdst,
const double *vsrc, ae_int_t N);
997 void vmove(
double *vdst, ae_int_t stride_dst,
const double* vsrc, ae_int_t stride_src, ae_int_t
n,
double alpha);
998 void vmove(
double *vdst,
const double *vsrc, ae_int_t N,
double alpha);
1006 void vadd(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n);
1007 void vadd(
double *vdst,
const double *vsrc, ae_int_t N);
1012 void vadd(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n,
double alpha);
1013 void vadd(
double *vdst,
const double *vsrc, ae_int_t N,
double alpha);
1021 void vsub(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n);
1022 void vsub(
double *vdst,
const double *vsrc, ae_int_t N);
1027 void vsub(
double *vdst, ae_int_t stride_dst,
const double *vsrc, ae_int_t stride_src, ae_int_t
n,
double alpha);
1028 void vsub(
double *vdst,
const double *vsrc, ae_int_t N,
double alpha);
1036 void vmul(
double *vdst, ae_int_t stride_dst, ae_int_t
n,
double alpha);
1037 void vmul(
double *vdst, ae_int_t N,
double alpha);
1060 void setlength(ae_int_t iLen);
1117 const ae_bool& operator()(ae_int_t
i)
const;
1118 ae_bool& operator()(ae_int_t i);
1120 const ae_bool& operator[](ae_int_t i)
const;
1121 ae_bool& operator[](ae_int_t i);
1123 void setcontent(ae_int_t iLen,
const bool *pContent );
1125 const ae_bool* getcontent()
const;
1127 std::string tostring()
const;
1140 const ae_int_t& operator()(ae_int_t
i)
const;
1141 ae_int_t& operator()(ae_int_t i);
1143 const ae_int_t& operator[](ae_int_t i)
const;
1144 ae_int_t& operator[](ae_int_t i);
1146 void setcontent(ae_int_t iLen,
const ae_int_t *pContent );
1148 ae_int_t* getcontent();
1149 const ae_int_t* getcontent()
const;
1151 std::string tostring()
const;
1164 const double& operator()(ae_int_t
i)
const;
1165 double& operator()(ae_int_t i);
1167 const double& operator[](ae_int_t i)
const;
1168 double& operator[](ae_int_t i);
1170 void setcontent(ae_int_t iLen,
const double *pContent );
1171 double* getcontent();
1172 const double* getcontent()
const;
1174 std::string tostring(
int dps)
const;
1197 std::string tostring(
int dps)
const;
1207 void setlength(ae_int_t rows, ae_int_t cols);
1208 ae_int_t rows()
const;
1209 ae_int_t cols()
const;
1210 bool isempty()
const;
1211 ae_int_t getstride()
const;
1265 const ae_bool& operator()(ae_int_t
i, ae_int_t
j)
const;
1266 ae_bool& operator()(ae_int_t i, ae_int_t j);
1268 const ae_bool* operator[](ae_int_t i)
const;
1269 ae_bool* operator[](ae_int_t i);
1271 void setcontent(ae_int_t irows, ae_int_t icols,
const bool *pContent );
1273 std::string tostring()
const ;
1285 const ae_int_t& operator()(ae_int_t
i, ae_int_t
j)
const;
1286 ae_int_t& operator()(ae_int_t i, ae_int_t j);
1288 const ae_int_t* operator[](ae_int_t i)
const;
1289 ae_int_t* operator[](ae_int_t i);
1291 void setcontent(ae_int_t irows, ae_int_t icols,
const ae_int_t *pContent );
1293 std::string tostring()
const;
1305 const double& operator()(ae_int_t
i, ae_int_t
j)
const;
1306 double& operator()(ae_int_t i, ae_int_t j);
1308 const double* operator[](ae_int_t i)
const;
1309 double* operator[](ae_int_t i);
1311 void setcontent(ae_int_t irows, ae_int_t icols,
const double *pContent );
1313 std::string tostring(
int dps)
const;
1331 void setcontent(ae_int_t irows, ae_int_t icols,
const alglib::complex *pContent );
1333 std::string tostring(
int dps)
const;
1382 extern const double fp_nan;
1390 int round(
double x);
1391 int trunc(
double x);
1393 int iceil(
double x);
1395 double sqr(
double x);
1396 int maxint(
int m1,
int m2);
1397 int minint(
int m1,
int m2);
1398 double maxreal(
double m1,
double m2);
1399 double minreal(
double m1,
double m2);
1401 bool fp_eq(
double v1,
double v2);
1402 bool fp_neq(
double v1,
double v2);
1403 bool fp_less(
double v1,
double v2);
1427 #define ALGLIB_INTERCEPTS_ABLAS 1430 void _ialglib_vcopy(ae_int_t
n,
const double *
a, ae_int_t stridea,
double *
b, ae_int_t strideb);
void ae_shared_pool_reset(ae_shared_pool *pool, ae_state *state)
ae_bool _ialglib_i_cmatrixrighttrsmf(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)
void ae_vector_init_from_x(ae_vector *dst, x_vector *src, ae_state *state, ae_bool make_automatic)
void ae_v_cmove(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n)
struct alglib_impl::ae_state ae_state
double beta(double a, double b, ae_state *_state)
void _ialglib_vcopy_dcomplex(ae_int_t n, const double *a, ae_int_t stridea, double *b, ae_int_t strideb, const char *conj)
void *volatile next_entry
ae_bool ae_fp_greater_eq(double v1, double v2)
alglib::complex vdotproduct(const alglib::complex *v1, const alglib::complex *v2, ae_int_t N)
bool fp_isneginf(double x)
void setnworkers(alglib::ae_int_t nworkers)
void ae_serializer_init(ae_serializer *serializer)
const bool operator==(const alglib::complex &lhs, const alglib::complex &rhs)
double ae_c_abs(ae_complex z, ae_state *state)
void ae_v_moved(double *vdst, ae_int_t stride_dst, const double *vsrc, ae_int_t stride_src, ae_int_t n, double alpha)
bool fp_less(double v1, double v2)
void ae_x_set_matrix(x_matrix *dst, ae_matrix *src, ae_state *state)
alglib::complex conj(const alglib::complex &z)
double ae_sin(double x, ae_state *state)
ae_bool ae_shared_pool_init(void *_dst, ae_state *state, ae_bool make_automatic)
double minreal(double m1, double m2)
void ae_serializer_serialize_bool(ae_serializer *serializer, ae_bool v, ae_state *state)
ae_int_t ae_serializer_get_alloc_size(ae_serializer *serializer)
void ae_serializer_sstart_str(ae_serializer *serializer, std::string *buf)
ae_bool ae_is_symmetric(ae_matrix *a)
bool fp_greater(double v1, double v2)
void ae_shared_pool_first_recycled(ae_shared_pool *pool, ae_smart_ptr *pptr, ae_state *state)
void ae_v_muld(double *vdst, ae_int_t stride_dst, ae_int_t n, double alpha)
ae_int_t ae_misalignment(const void *ptr, size_t alignment)
void _ialglib_vzero_complex(ae_int_t n, ae_complex *p, ae_int_t stride)
ae_bool ae_c_neq(ae_complex lhs, ae_complex rhs)
void ae_shared_pool_retrieve(ae_shared_pool *pool, ae_smart_ptr *pptr, ae_state *state)
double ae_fabs(double x, ae_state *state)
void ae_shared_pool_clear(void *_dst)
alglib_impl::ae_matrix mat
double ae_tan(double x, ae_state *state)
bool fp_eq(double v1, double v2)
double ae_pow(double x, double y, ae_state *state)
void * ae_malloc(size_t size, ae_state *state)
void ae_db_free(ae_dyn_block *block)
void _ialglib_mcopyblock(ae_int_t m, ae_int_t n, const double *a, ae_int_t op, ae_int_t stride, double *b)
bool fp_less_eq(double v1, double v2)
void _ialglib_mcopyblock_complex(ae_int_t m, ae_int_t n, const ae_complex *a, ae_int_t op, ae_int_t stride, double *b)
ae_bool ae_isposinf(double x, ae_state *state)
void ae_frame_make(ae_state *state, ae_frame *tmp)
#define AE_LOCK_ALIGNMENT
const bool operator!=(const alglib::complex &lhs, const alglib::complex &rhs)
ae_complex ae_c_d_sub(double lhs, ae_complex rhs)
const ae_int_t endianness
ae_shared_pool_entry *volatile recycled_objects
void x_vector_clear(x_vector *dst)
int maxint(int m1, int m2)
alglib_impl::ae_matrix * p_mat
ae_error_type volatile last_error
ae_bool x_force_symmetric(x_matrix *a)
ae_bool ae_c_eq_d(ae_complex lhs, double rhs)
ae_complex ae_c_conj(ae_complex lhs, ae_state *state)
ae_bool ae_db_realloc(ae_dyn_block *block, ae_int_t size, ae_state *state)
void(* deallocator)(void *)
ae_complex ae_complex_from_d(double v)
ae_int_t ae_sign(double x, ae_state *state)
void ae_serializer_ustart_str(ae_serializer *serializer, const std::string *buf)
void ae_db_swap(ae_dyn_block *block1, ae_dyn_block *block2)
ae_bool ae_force_hermitian(ae_matrix *a)
void ae_swap_vectors(ae_vector *vec1, ae_vector *vec2)
ae_bool ae_isnan_stateless(double x, ae_int_t endianness)
ae_bool ae_isneginf_stateless(double x, ae_int_t endianness)
ae_bool _ialglib_i_rmatrixgemmf(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)
const alglib::complex operator+(const double &lhs, const alglib::complex &rhs)
void ae_v_csubd(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n, double alpha)
void ae_shared_pool_recycle(ae_shared_pool *pool, ae_smart_ptr *pptr, ae_state *state)
ae_complex ae_c_div_d(ae_complex lhs, double rhs)
int minint(int m1, int m2)
void ae_serializer_stop(ae_serializer *serializer)
ae_bool ae_c_neq_d(ae_complex lhs, double rhs)
void ae_state_clear(ae_state *state)
ae_bool ae_isposinf_stateless(double x, ae_int_t endianness)
ae_bool ae_fp_eq(double v1, double v2)
void ae_smart_ptr_assign(ae_smart_ptr *dst, void *new_ptr, ae_bool is_owner, ae_bool is_dynamic, void(*destroy)(void *))
void ae_v_cmulc(ae_complex *vdst, ae_int_t stride_dst, ae_int_t n, ae_complex alpha)
ae_bool ae_db_malloc(ae_dyn_block *block, ae_int_t size, ae_state *state, ae_bool make_automatic)
double ae_cos(double x, ae_state *state)
ae_bool ae_isneginf(double x, ae_state *state)
ae_bool ae_matrix_init_copy(ae_matrix *dst, ae_matrix *src, ae_state *state, ae_bool make_automatic)
ae_bool ae_matrix_init(ae_matrix *dst, ae_int_t rows, ae_int_t cols, ae_datatype datatype, ae_state *state, ae_bool make_automatic)
void * ae_align(void *ptr, size_t alignment)
void ae_matrix_destroy(ae_matrix *dst)
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
void * aligned_malloc(size_t size, size_t alignment)
void ae_v_add(double *vdst, ae_int_t stride_dst, const double *vsrc, ae_int_t stride_src, ae_int_t n)
ae_complex ae_c_sub_d(ae_complex lhs, double rhs)
ae_bool x_is_hermitian(x_matrix *a)
void ae_serializer_clear(ae_serializer *serializer)
void ae_db_attach(ae_dyn_block *block, ae_state *state)
void ae_smart_ptr_release(ae_smart_ptr *dst)
void vsub(alglib::complex *vdst, const alglib::complex *vsrc, ae_int_t N, alglib::complex alpha)
alglib::complex csqr(const alglib::complex &z)
double ae_asin(double x, ae_state *state)
void ae_serializer_serialize_int(ae_serializer *serializer, ae_int_t v, ae_state *state)
ae_int_t ae_randominteger(ae_int_t maxv, ae_state *state)
void _rcommstate_destroy(rcommstate *p)
alglib_impl::ae_vector * p_vec
bool fp_isfinite(double x)
ae_int_t ae_v_len(ae_int_t a, ae_int_t b)
void ae_vector_destroy(ae_vector *dst)
ae_bool _rcommstate_init(rcommstate *p, ae_state *_state, ae_bool make_automatic)
void vadd(alglib::complex *vdst, const alglib::complex *vsrc, ae_int_t N, alglib::complex alpha)
struct ae_dyn_block *volatile p_next
void ae_shared_pool_set_seed(ae_shared_pool *dst, void *seed_object, ae_int_t size_of_object, ae_bool(*init)(void *dst, ae_state *state, ae_bool make_automatic), ae_bool(*init_copy)(void *dst, void *src, ae_state *state, ae_bool make_automatic), void(*destroy)(void *ptr), ae_state *state)
ae_complex ae_c_sqr(ae_complex lhs, ae_state *state)
void ae_v_csub(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n)
void ae_free_lock(ae_lock *lock)
void vmove(alglib::complex *vdst, const alglib::complex *vsrc, ae_int_t N, alglib::complex alpha)
ae_bool _rcommstate_init_copy(rcommstate *dst, rcommstate *src, ae_state *_state, ae_bool make_automatic)
void ae_v_move(double *vdst, ae_int_t stride_dst, const double *vsrc, ae_int_t stride_src, ae_int_t n)
void ae_swap_matrices(ae_matrix *mat1, ae_matrix *mat2)
ae_complex ae_c_div(ae_complex lhs, ae_complex rhs)
struct alglib_impl::rcommstate rcommstate
ae_bool _ialglib_i_cmatrixgemmf(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_bool ae_isinf(double x, ae_state *state)
ae_complex ae_c_add(ae_complex lhs, ae_complex rhs)
void ae_vector_clear(ae_vector *dst)
void(* ae_deallocator)(void *)
void ae_v_csubc(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n, ae_complex alpha)
ae_bool ae_force_symmetric(ae_matrix *a)
ae_bool ae_isinf_stateless(double x, ae_int_t endianness)
struct alglib_impl::ae_smart_ptr ae_smart_ptr
ae_bool ae_fp_less(double v1, double v2)
ae_bool ae_shared_pool_is_initialized(void *_dst)
void ae_break(ae_state *state, ae_error_type error_type, const char *msg)
ae_int_t ae_iceil(double x, ae_state *state)
void ae_serializer_unserialize_int(ae_serializer *serializer, ae_int_t *v, ae_state *state)
double ae_randomreal(ae_state *state)
void vmoveneg(alglib::complex *vdst, const alglib::complex *vsrc, ae_int_t N)
ae_int_t ae_sizeof(ae_datatype datatype)
ae_bool ae_smart_ptr_init(ae_smart_ptr *dst, void **subscriber, ae_state *state, ae_bool make_automatic)
const double machineepsilon
const alglib::complex operator*(const double &lhs, const alglib::complex &rhs)
ae_complex ae_c_add_d(ae_complex lhs, double rhs)
struct alglib_impl::ae_frame ae_frame
struct alglib_impl::ae_dyn_block ae_dyn_block
const alglib::complex operator/(const alglib::complex &lhs, const double &rhs)
void ae_v_caddc(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n, ae_complex alpha)
void ae_release_lock(ae_lock *lock)
ae_bool ae_shared_pool_init_copy(void *_dst, void *_src, ae_state *state, ae_bool make_automatic)
ae_complex ae_c_sub(ae_complex lhs, ae_complex rhs)
ae_bool ae_fp_neq(double v1, double v2)
ae_bool ae_isnan(double x, ae_state *state)
void ae_v_cmoved(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n, double alpha)
ae_bool x_force_hermitian(x_matrix *a)
void ae_shared_pool_next_recycled(ae_shared_pool *pool, ae_smart_ptr *pptr, ae_state *state)
__host__ __device__ float length(float2 v)
void ae_v_cadd(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n)
const double minrealnumber
void ae_touch_ptr(void *p)
void ae_acquire_lock(ae_lock *lock)
std::vector< T > & operator-=(std::vector< T > &_v1, const std::vector< T > &_v2)
double ae_maxreal(double m1, double m2, ae_state *state)
ae_complex ae_v_cdotproduct(const ae_complex *v0, ae_int_t stride0, const char *conj0, const ae_complex *v1, ae_int_t stride1, const char *conj1, ae_int_t n)
struct alglib_impl::ae_shared_pool_entry ae_shared_pool_entry
void ae_x_attach_to_vector(x_vector *dst, ae_vector *src)
ae_bool ae_vector_set_length(ae_vector *dst, ae_int_t newsize, ae_state *state)
void ae_v_sub(double *vdst, ae_int_t stride_dst, const double *vsrc, ae_int_t stride_src, ae_int_t n)
bool fp_neq(double v1, double v2)
void _ialglib_vcopy_complex(ae_int_t n, const ae_complex *a, ae_int_t stridea, double *b, ae_int_t strideb, const char *conj)
ae_bool _ialglib_i_cmatrixlefttrsmf(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)
void _ialglib_mcopyunblock(ae_int_t m, ae_int_t n, const double *a, ae_int_t op, double *b, ae_int_t stride)
double ae_log(double x, ae_state *state)
void ae_serializer_serialize_double(ae_serializer *serializer, double v, ae_state *state)
void _ialglib_vzero(ae_int_t n, double *p, ae_int_t stride)
struct alglib_impl::ae_vector ae_vector
void ae_smart_ptr_clear(void *_dst)
double ae_minreal(double m1, double m2, ae_state *state)
struct alglib_impl::ae_shared_pool ae_shared_pool
ae_bool _ialglib_i_rmatrixrank1f(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 uoffs, ae_vector *v, ae_int_t voffs)
ae_complex ae_c_d_div(double lhs, ae_complex rhs)
ae_bool ae_c_eq(ae_complex lhs, ae_complex rhs)
bool fp_greater_eq(double v1, double v2)
ae_bool _ialglib_i_rmatrixlefttrsmf(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)
void _ialglib_vcopy(ae_int_t n, const double *a, ae_int_t stridea, double *b, ae_int_t strideb)
struct alglib_impl::ae_matrix ae_matrix
double ae_acos(double x, ae_state *state)
void vmul(alglib::complex *vdst, ae_int_t N, alglib::complex alpha)
double ae_sinh(double x, ae_state *state)
ae_bool _ialglib_i_rmatrixsyrkf(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_int_t ae_ifloor(double x, ae_state *state)
ae_int_t ae_iabs(ae_int_t x, ae_state *state)
ae_shared_pool_entry *volatile recycled_entries
double abscomplex(const alglib::complex &z)
void ae_state_init(ae_state *state)
void ae_shared_pool_destroy(void *_dst)
void ae_matrix_init_from_x(ae_matrix *dst, x_matrix *src, ae_state *state, ae_bool make_automatic)
double ae_sqrt(double x, ae_state *state)
void ae_assert(ae_bool cond, const char *msg, ae_state *state)
signed long long ae_int64_t
ae_bool _ialglib_i_cmatrixrank1f(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 uoffs, ae_vector *v, ae_int_t voffs)
void ae_v_moveneg(double *vdst, ae_int_t stride_dst, const double *vsrc, ae_int_t stride_src, ae_int_t n)
void ae_v_cmoveneg(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n)
ae_bool ae_isfinite_stateless(double x, ae_int_t endianness)
void ae_serializer_alloc_start(ae_serializer *serializer)
const char *volatile error_msg
const double maxrealnumber
ae_bool x_is_symmetric(x_matrix *a)
double ae_atan(double x, ae_state *state)
ae_bool _ialglib_i_cmatrixsyrkf(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)
double ae_cosh(double x, ae_state *state)
double ae_exp(double x, ae_state *state)
void ae_x_attach_to_matrix(x_matrix *dst, ae_matrix *src)
double ae_tanh(double x, ae_state *state)
void ae_v_caddd(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n, double alpha)
ae_complex ae_c_mul(ae_complex lhs, ae_complex rhs)
void *volatile seed_object
ae_dyn_block *volatile p_top_block
void _rcommstate_clear(rcommstate *p)
void ae_v_subd(double *vdst, ae_int_t stride_dst, const double *vsrc, ae_int_t stride_src, ae_int_t n, double alpha)
ae_bool ae_vector_init(ae_vector *dst, ae_int_t size, ae_datatype datatype, ae_state *state, ae_bool make_automatic)
double maxreal(double m1, double m2)
ae_int_t ae_maxint(ae_int_t m1, ae_int_t m2, ae_state *state)
ae_complex ae_c_mul_d(ae_complex lhs, double rhs)
std::vector< T > & operator+=(std::vector< T > &_v1, const std::vector< T > &_v2)
void aligned_free(void *block)
ae_bool ae_isfinite(double x, ae_state *state)
void ae_serializer_unserialize_bool(ae_serializer *serializer, ae_bool *v, ae_state *state)
double ae_sqr(double x, ae_state *state)
void ae_serializer_alloc_entry(ae_serializer *serializer)
void ae_serializer_unserialize_double(ae_serializer *serializer, double *v, ae_state *state)
void ae_init_lock(ae_lock *lock)
ae_complex ae_c_neg(ae_complex lhs)
std::vector< T > & operator/=(std::vector< T > &v1, const std::vector< T > &v2)
const alglib::complex operator-(const double &lhs, const alglib::complex &rhs)
void ae_v_addd(double *vdst, ae_int_t stride_dst, const double *vsrc, ae_int_t stride_src, ae_int_t n, double alpha)
ae_bool ae_vector_init_copy(ae_vector *dst, ae_vector *src, ae_state *state, ae_bool make_automatic)
void ae_v_cmuld(ae_complex *vdst, ae_int_t stride_dst, ae_int_t n, double alpha)
ae_bool ae_is_hermitian(ae_matrix *a)
ae_bool ae_fp_less_eq(double v1, double v2)
ae_int_t ae_round(double x, ae_state *state)
ae_int_t ae_get_endianness()
alglib_impl::ae_int_t ae_int_t
ae_shared_pool_entry *volatile enumeration_counter
void ae_frame_leave(ae_state *state)
void ae_matrix_clear(ae_matrix *dst)
ae_int_t vlen(ae_int_t n1, ae_int_t n2)
void ae_shared_pool_clear_recycled(ae_shared_pool *pool, ae_state *state)
std::vector< T > & operator*=(std::vector< T > &v1, const std::vector< T > &v2)
ae_bool ae_fp_greater(double v1, double v2)
void ae_smart_ptr_destroy(void *_dst)
ae_bool ae_matrix_set_length(ae_matrix *dst, ae_int_t rows, ae_int_t cols, ae_state *state)
void ae_x_set_vector(x_vector *dst, ae_vector *src, ae_state *state)
ae_bool _ialglib_i_rmatrixrighttrsmf(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_int_t ae_minint(ae_int_t m1, ae_int_t m2, ae_state *state)
alglib_impl::ae_vector vec
ae_int_t ae_trunc(double x, ae_state *state)
double ae_atan2(double y, double x, ae_state *state)
double ae_v_dotproduct(const double *v0, ae_int_t stride0, const double *v1, ae_int_t stride1, ae_int_t n)
ae_int_t randominteger(ae_int_t maxv)
bool fp_isposinf(double x)
void _ialglib_mcopyunblock_complex(ae_int_t m, ae_int_t n, const double *a, ae_int_t op, ae_complex *b, ae_int_t stride)
void ae_v_cmovec(ae_complex *vdst, ae_int_t stride_dst, const ae_complex *vsrc, ae_int_t stride_src, const char *conj_src, ae_int_t n, ae_complex alpha)