2
2
/* create a Householder reflection [v,beta,s]=house(x), overwrite x with v,
3
3
* where (I-beta*v*v')*x = s*e1 and e1 = [1 0 ... 0]'.
4
4
* Note that this CXSparse version is different than CSparse. */
9
9
if (!x || !beta) return (-1) ; /* check inputs */
10
11
for (i = 0 ; i < n ; i++) s += x [i] * CS_CONJ (x [i]) ;