seq/spqr/aijspqr.c

114   cholmod_dense *Y_handle = NULL, *QTB_handle = NULL, *Z_handle = NULL;
120     Y_handle = SuiteSparseQR_C_solve(SPQR_RETX_EQUALS_B, chol->spqrfact, QTB_handle, chol->common);
121     PetscCheck(Y_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_solve failed");
125     Y_handle = SuiteSparseQR_C_solve(SPQR_RETX_EQUALS_B, chol->spqrfact, Z_handle, chol->common);
126     PetscCheck(Y_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_solve failed");
129   *_Y_handle = Y_handle;
137   cholmod_dense cholB, *Y_handle = NULL;
143   PetscCall(MatSolve_SPQR_Internal(F, &cholB, &Y_handle));
146   PetscCall(PetscArraycpy(v, (PetscScalar *)Y_handle->x, n));
148   PetscCallExternal(!cholmod_l_free_dense, &Y_handle, chol->common);
157   cholmod_dense cholB, *Y_handle = NULL;
163   PetscCall(MatSolve_SPQR_Internal(F, &cholB, &Y_handle));
166   if ((size_t)lda == Y_handle->d) {
167     PetscCall(PetscArraycpy(v, (PetscScalar *)Y_handle->x, lda * Y_handle->ncol));
169     for (size_t j = 0; j < Y_handle->ncol; j++) PetscCall(PetscArraycpy(&v[j * lda], &(((PetscScalar *)Y_handle->x)[j * Y_handle->d]), Y_handle->nrow));
172   PetscCallExternal(!cholmod_l_free_dense, &Y_handle, chol->common);
181   cholmod_dense *Y_handle = NULL, *RTB_handle = NULL;
186   Y_handle = SuiteSparseQR_C_qmult(SPQR_QX, chol->spqrfact, RTB_handle, chol->common);
187   PetscCheck(Y_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_qmult failed");
188   *_Y_handle = Y_handle;
196   cholmod_dense cholB, *Y_handle = NULL;
202   PetscCall(MatSolveTranspose_SPQR_Internal(F, &cholB, &Y_handle));
205   PetscCall(PetscArraycpy(v, (PetscScalar *)Y_handle->x, n));
207   PetscCallExternal(!cholmod_l_free_dense, &Y_handle, chol->common);
215   cholmod_dense cholB, *Y_handle = NULL;
221   PetscCall(MatSolveTranspose_SPQR_Internal(F, &cholB, &Y_handle));
224   if ((size_t)lda == Y_handle->d) {
225     PetscCall(PetscArraycpy(v, (PetscScalar *)Y_handle->x, lda * Y_handle->ncol));
227     for (size_t j = 0; j < Y_handle->ncol; j++) PetscCall(PetscArraycpy(&v[j * lda], &(((PetscScalar *)Y_handle->x)[j * Y_handle->d]), Y_handle->nrow));
230   PetscCallExternal(!cholmod_l_free_dense, &Y_handle, chol->common);