seq/spqr/aijspqr.c

111 static PetscErrorCode MatSolve_SPQR_Internal(Mat F, cholmod_dense *cholB, cholmod_dense **_Y_handle)  in MatSolve_SPQR_Internal()  argument
113   Mat_CHOLMOD   *chol     = (Mat_CHOLMOD *)F->data;  in MatSolve_SPQR_Internal()
119 …PetscCheck(QTB_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_qmult fail…  in MatSolve_SPQR_Internal()
121 …PetscCheck(Y_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_solve failed…  in MatSolve_SPQR_Internal()
124 …PetscCheck(Z_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_solve failed…  in MatSolve_SPQR_Internal()
126 …PetscCheck(Y_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_solve failed…  in MatSolve_SPQR_Internal()
134 static PetscErrorCode MatSolve_SPQR(Mat F, Vec B, Vec X)  in MatSolve_SPQR()  argument
136   Mat_CHOLMOD  *chol = (Mat_CHOLMOD *)F->data;  in MatSolve_SPQR()
143   PetscCall(MatSolve_SPQR_Internal(F, &cholB, &Y_handle));  in MatSolve_SPQR()
154 static PetscErrorCode MatMatSolve_SPQR(Mat F, Mat B, Mat X)  in MatMatSolve_SPQR()  argument
156   Mat_CHOLMOD  *chol = (Mat_CHOLMOD *)F->data;  in MatMatSolve_SPQR()
163   PetscCall(MatSolve_SPQR_Internal(F, &cholB, &Y_handle));  in MatMatSolve_SPQR()
178 static PetscErrorCode MatSolveTranspose_SPQR_Internal(Mat F, cholmod_dense *cholB, cholmod_dense **…  in MatSolveTranspose_SPQR_Internal()  argument
180   Mat_CHOLMOD   *chol     = (Mat_CHOLMOD *)F->data;  in MatSolveTranspose_SPQR_Internal()
185 …PetscCheck(RTB_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_solve fail…  in MatSolveTranspose_SPQR_Internal()
187 …PetscCheck(Y_handle, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SuiteSparseQR_C_qmult failed…  in MatSolveTranspose_SPQR_Internal()
193 static PetscErrorCode MatSolveTranspose_SPQR(Mat F, Vec B, Vec X)  in MatSolveTranspose_SPQR()  argument
195   Mat_CHOLMOD  *chol = (Mat_CHOLMOD *)F->data;  in MatSolveTranspose_SPQR()
202   PetscCall(MatSolveTranspose_SPQR_Internal(F, &cholB, &Y_handle));  in MatSolveTranspose_SPQR()
212 static PetscErrorCode MatMatSolveTranspose_SPQR(Mat F, Mat B, Mat X)  in MatMatSolveTranspose_SPQR()  argument
214   Mat_CHOLMOD  *chol = (Mat_CHOLMOD *)F->data;  in MatMatSolveTranspose_SPQR()
221   PetscCall(MatSolveTranspose_SPQR_Internal(F, &cholB, &Y_handle));  in MatMatSolveTranspose_SPQR()
235 static PetscErrorCode MatQRFactorNumeric_SPQR(Mat F, Mat A, const MatFactorInfo *info)  in MatQRFactorNumeric_SPQR()  argument
237   Mat_CHOLMOD   *chol = (Mat_CHOLMOD *)F->data;  in MatQRFactorNumeric_SPQR()
247 …PetscCheck(!err, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "SPQR factorization failed with s…  in MatQRFactorNumeric_SPQR()
252   F->ops->solve    = MatSolve_SPQR;  in MatQRFactorNumeric_SPQR()
253   F->ops->matsolve = MatMatSolve_SPQR;  in MatQRFactorNumeric_SPQR()
258     F->ops->solvetranspose    = MatSolve_SPQR;  in MatQRFactorNumeric_SPQR()
259     F->ops->matsolvetranspose = MatMatSolve_SPQR;  in MatQRFactorNumeric_SPQR()
262     F->ops->solvetranspose    = MatSolveTranspose_SPQR;  in MatQRFactorNumeric_SPQR()
263     F->ops->matsolvetranspose = MatMatSolveTranspose_SPQR;  in MatQRFactorNumeric_SPQR()
269 PETSC_INTERN PetscErrorCode MatQRFactorSymbolic_SPQR(Mat F, Mat A, IS perm, const MatFactorInfo *in…  in MatQRFactorSymbolic_SPQR()  argument
271   Mat_CHOLMOD   *chol = (Mat_CHOLMOD *)F->data;  in MatQRFactorSymbolic_SPQR()
282 …PetscCheck(chol->spqrfact, PetscObjectComm((PetscObject)F), PETSC_ERR_LIB, "CHOLMOD analysis faile…  in MatQRFactorSymbolic_SPQR()
287 …PetscCall(PetscObjectComposeFunction((PetscObject)F, "MatQRFactorNumeric_C", MatQRFactorNumeric_SP…  in MatQRFactorSymbolic_SPQR()
307 PETSC_INTERN PetscErrorCode MatGetFactor_seqaij_spqr(Mat A, MatFactorType ftype, Mat *F)  in MatGetFactor_seqaij_spqr()  argument
343   *F                  = B;  in MatGetFactor_seqaij_spqr()