seq/spqr/aijspqr.c

111 static PetscErrorCode MatSolve_SPQR_Internal(Mat F, cholmod_dense *cholB, cholmod_dense **_Y_handle)
118     QTB_handle = SuiteSparseQR_C_qmult(SPQR_QTX, chol->spqrfact, cholB, chol->common);
123     Z_handle = SuiteSparseQR_C_solve(SPQR_RTX_EQUALS_ETB, chol->spqrfact, cholB, chol->common);
137   cholmod_dense cholB, *Y_handle = NULL;
142   PetscCall(VecWrapCholmod(B, GET_ARRAY_READ, &cholB));
143   PetscCall(MatSolve_SPQR_Internal(F, &cholB, &Y_handle));
149   PetscCall(VecUnWrapCholmod(B, GET_ARRAY_READ, &cholB));
157   cholmod_dense cholB, *Y_handle = NULL;
162   PetscCall(MatDenseWrapCholmod(B, GET_ARRAY_READ, &cholB));
163   PetscCall(MatSolve_SPQR_Internal(F, &cholB, &Y_handle));
173   PetscCall(MatDenseUnWrapCholmod(B, GET_ARRAY_READ, &cholB));
178 static PetscErrorCode MatSolveTranspose_SPQR_Internal(Mat F, cholmod_dense *cholB, cholmod_dense **_Y_handle)
184   RTB_handle = SuiteSparseQR_C_solve(SPQR_RTX_EQUALS_ETB, chol->spqrfact, cholB, chol->common);
196   cholmod_dense cholB, *Y_handle = NULL;
201   PetscCall(VecWrapCholmod(B, GET_ARRAY_READ, &cholB));
202   PetscCall(MatSolveTranspose_SPQR_Internal(F, &cholB, &Y_handle));
208   PetscCall(VecUnWrapCholmod(B, GET_ARRAY_READ, &cholB));
215   cholmod_dense cholB, *Y_handle = NULL;
220   PetscCall(MatDenseWrapCholmod(B, GET_ARRAY_READ, &cholB));
221   PetscCall(MatSolveTranspose_SPQR_Internal(F, &cholB, &Y_handle));
231   PetscCall(MatDenseUnWrapCholmod(B, GET_ARRAY_READ, &cholB));