c++ - Indexing into CHOLMOD dense vector array -
i have cholmod_dense data structure:
cholmod_dense* ex = cholmod_l_solve(cholmod_a, l, b, &com);
and want extract values , copy them variable. means need index double array , copy values over.
for (int k=0; k<ncols; k++) t_x[k]=((double*)ex->x)[k]; which compiler ok segmenation fault. or think should able do:
double* e_x =(double*)ex->x; (int k=0; k<ncols; k++) t_x[k]=*e_x[k]; but compiler dislikes this:
error: invalid type argument of unary ‘*’ (have ‘double’) (int k=0; k<ncols; k++) t_x[k]= *e_x[k]; according cholmod userguide:
- cholmod dense: dense matrix, either real, complex or zomplex, in column-major order. differs row-major convention used in c. dense matrix x contains • x->x, double array of size x->nzmax or twice complex case. • x->z, double array of size x->nzmax if x zomplex.
so should able grab ex->x , index double array, cannot without getting segmentation fault. can me out?
the cholmod library written in c , code linking cholmod library (the code snippet shown above) c++.
ok, looks likes made couple mistakes.
first of all, running segmentation fault because using cholmod_l_zeros(); assumes long integers. instead, should using cholmod_zeros(); since using doubles.
after fixing ran error cholmod error: invalid xtype right after cholmod_solve(cholmod_a, l, b, &com); statement. because cholmod_factor* l defined out of scope. after fixing both of these issues, code copies values cholmod_dense ex->x double array on t_x double vector:
cholmod_dense* ex = cholmod_solve(cholmod_a, l, b, &com); double* e_x = (double*)ex->x; (int k=0; k<ncols; k++) t_x[k] = e_x[k]; i unaware [] operator automatically dereferences pointers. know!
Comments
Post a Comment