58 s = x * sqrt (1.0 + r*r) ;
70 s = y * sqrt (1.0 + r*r) ;
100 double *cr,
double *ci
103 double tr, ti, r, den ;
104 if (fabs (br) >= fabs (bi))
108 tr = (ar + ai * r) / den ;
109 ti = (ai - ar * r) / den ;
115 tr = (ar * r + ai) / den ;
116 ti = (ai * r - ar) / den ;
159 double *Xold, *Zold, *Xnew, *Znew ;
163 if (xtype_out < xtype1 || xtype_out > xtype2)
196 for (k = 0 ; k < nz ; k++)
214 for (k = 0 ; k < nz ; k++)
236 for (k = 0 ; k < nz ; k++)
276 for (k = 0 ; k < nz ; k++)
278 Xnew [2*k ] = Xold [k] ;
296 for (k = 0 ; k < nz ; k++)
329 for (k = 0 ; k < nz ; k++)
331 Xold [k] = Xold [2*k] ;
335 *XX =
CHOLMOD(realloc) (nz,
sizeof (double), *XX, &nz2,
353 for (k = 0 ; k < nz ; k++)
355 Xnew [k] = Xold [2*k ] ;
356 Znew [k] = Xold [2*k+1] ;
403 for (k = 0 ; k < nz ; k++)
405 Xnew [2*k ] = Xold [k] ;
406 Xnew [2*k+1] = Zold [k] ;
448 A->xtype = to_xtype ;
479 T->xtype = to_xtype ;
510 X->xtype = to_xtype ;
547 L->xtype = to_xtype ;
double CHOLMOD() hypot(double x, double y)
int CHOLMOD() sparse_xtype(int to_xtype, cholmod_sparse *A, cholmod_common *Common)
int CHOLMOD() triplet_xtype(int to_xtype, cholmod_triplet *T, cholmod_common *Common)
int CHOLMOD() factor_xtype(int to_xtype, cholmod_factor *L, cholmod_common *Common)
#define RETURN_IF_NULL_COMMON(result)
static int change_complexity(Int nz, int xtype_in, int xtype_out, int xtype1, int xtype2, void **XX, void **ZZ, cholmod_common *Common)
int CHOLMOD() divcomplex(double ar, double ai, double br, double bi, double *cr, double *ci)
#define RETURN_IF_NULL(A, result)
#define ERROR(status, msg)
#define RETURN_IF_XTYPE_INVALID(A, xtype1, xtype2, result)
int CHOLMOD() dense_xtype(int to_xtype, cholmod_dense *X, cholmod_common *Common)