| /petsc/src/mat/tests/ |
| H A D | ex192.c | 16 PetscBool data_provided, herm, symm, use_lu, cuda = PETSC_FALSE; in main() local
29 herm = PETSC_FALSE; in main()
32 PetscCall(PetscOptionsGetBool(NULL, NULL, "-hermitian_solve", &herm, NULL)); in main()
33 if (herm) symm = PETSC_TRUE; in main()
142 … schur %" PetscInt_FMT ", size mat %" PetscInt_FMT "\n", solver, nrhs, symm, herm, size_schur, m)); in main()
148 if (cuda && symm && !herm) use_lu = PETSC_TRUE; in main()
150 if (herm && !use_lu) { /* test also conversion routines inside the solver packages */ in main()
158 if (herm) { in main()
182 if (symm && herm) PetscCall(VecAbs(x)); in main()
200 if (nfact == 1 && (!cuda || (herm && symm))) PetscCall(MatFactorInvertSchurComplement(F)); in main()
|
| H A D | ex215.c | 13 PetscBool data_provided, herm, symm, hpd; in main() local
23 herm = PETSC_FALSE; in main()
27 PetscCall(PetscOptionsGetBool(NULL, NULL, "-hermitian_solve", &herm, NULL)); in main()
105 if (herm) { in main()
119 PetscCall(MatSetOption(F, MAT_HERMITIAN, (PetscBool)(hpd || herm))); in main()
|
| /petsc/src/mat/tests/output/ |
| H A D | ex192_mumps_2.out | 1 Solving with mumps: nrhs 2, sym 1, herm 0, size schur 5, size mat 12
|
| H A D | ex192_mumps.out | 1 Solving with mumps: nrhs 2, sym 0, herm 0, size schur 5, size mat 12
|
| H A D | ex192_mkl_pardiso_1.out | 1 Solving with mkl_pardiso: nrhs 2, sym 1, herm 0, size schur 5, size mat 12
|
| H A D | ex192_mkl_pardiso_3.out | 1 Solving with mkl_pardiso: nrhs 2, sym 1, herm 1, size schur 5, size mat 12
|
| H A D | ex192_mumps_3.out | 1 Solving with mumps: nrhs 2, sym 1, herm 1, size schur 5, size mat 12
|
| H A D | ex192_mkl_pardiso.out | 1 Solving with mkl_pardiso: nrhs 2, sym 0, herm 0, size schur 5, size mat 12
|
| /petsc/src/ksp/ksp/guess/impls/pod/ |
| H A D | pod.c | 241 PetscBool set, herm; in KSPGuessUpdate_POD() local
244 PetscCall(MatIsHermitianKnown(guess->A, &set, &herm)); in KSPGuessUpdate_POD()
246 PetscCall(MatIsSymmetricKnown(guess->A, &set, &herm)); in KSPGuessUpdate_POD()
248 off = (guess->ksp->transpose_solve && (!set || !herm)) ? 2 * pod->n : pod->n; in KSPGuessUpdate_POD()
253 if (!set || !herm) { in KSPGuessUpdate_POD()
|
| /petsc/src/mat/impls/nest/ |
| H A D | matnest.c | 261 static PetscErrorCode MatMultTransposeKernel_Nest(Mat A, Vec x, Vec y, PetscBool herm) in MatMultTransposeKernel_Nest() argument
274 …if (herm) PetscCall(MatMultHermitianTransposeAdd(bA->m[i][j], bx[i], by[j], by[j])); /* y[j] <- y[… in MatMultTransposeKernel_Nest()
297 static PetscErrorCode MatMultTransposeAddKernel_Nest(Mat A, Vec x, Vec y, Vec z, PetscBool herm) in MatMultTransposeAddKernel_Nest() argument
315 …if (herm) PetscCall(MatMultHermitianTransposeAdd(bA->m[i][j], bx[i], bz[j], bz[j])); /* z[j] <- y[… in MatMultTransposeAddKernel_Nest()
|
| /petsc/src/mat/impls/dense/mpi/ |
| H A D | mpidense.c | 582 static PetscErrorCode MatMultTransposeKernel_MPIDense(Mat A, Vec xx, Vec yy, PetscBool herm) in MatMultTransposeKernel_MPIDense() argument
592 if (herm) PetscCall((*a->A->ops->multhermitiantranspose)(a->A, xx, a->lvec)); in MatMultTransposeKernel_MPIDense()
624 …ic PetscErrorCode MatMultTransposeAddKernel_MPIDense(Mat A, Vec xx, Vec yy, Vec zz, PetscBool herm) in MatMultTransposeAddKernel_MPIDense() argument
634 if (herm) PetscCall((*a->A->ops->multhermitiantranspose)(a->A, xx, a->lvec)); in MatMultTransposeAddKernel_MPIDense()
|
| /petsc/src/mat/impls/dense/seq/ |
| H A D | dense.c | 1045 …_SeqDense(Mat A, Vec xx, Vec yy, PetscInt c_start, PetscInt c_end, PetscBool trans, PetscBool herm) in MatMultColumnRangeKernel_SeqDense() argument
1065 …if (herm) PetscCallBLAS("BLASgemv", BLASgemv_("C", &m, &n, &_DOne, v + c_start * mat->lda, &mat->l… in MatMultColumnRangeKernel_SeqDense()
1105 …e(Mat A, Vec xx, Vec zz, Vec yy, PetscInt c_start, PetscInt c_end, PetscBool trans, PetscBool herm) in MatMultAddColumnRangeKernel_SeqDense() argument
1120 …if (herm) PetscCallBLAS("BLASgemv", BLASgemv_("C", &m, &n, &_DOne, v + c_start * mat->lda, &mat->l… in MatMultAddColumnRangeKernel_SeqDense()
|
| /petsc/src/mat/impls/aij/seq/seqhipsparse/ |
| H A D | aijhipsparse.hip.cxx | 3021 …de MatMultAddKernel_SeqAIJHIPSPARSE(Mat A, Vec xx, Vec yy, Vec zz, PetscBool trans, PetscBool herm) in MatMultAddKernel_SeqAIJHIPSPARSE() argument
3032 …PetscCheck(!herm || trans, PetscObjectComm((PetscObject)A), PETSC_ERR_GPU, "Hermitian and not tran… in MatMultAddKernel_SeqAIJHIPSPARSE()
3044 if (herm || !A->form_explicit_transpose) { in MatMultAddKernel_SeqAIJHIPSPARSE()
3045 opA = herm ? HIPSPARSE_OPERATION_CONJUGATE_TRANSPOSE : HIPSPARSE_OPERATION_TRANSPOSE; in MatMultAddKernel_SeqAIJHIPSPARSE()
|
| /petsc/src/mat/impls/aij/seq/seqcusparse/ |
| H A D | aijcusparse.cu | 3555 …ode MatMultAddKernel_SeqAIJCUSPARSE(Mat A, Vec xx, Vec yy, Vec zz, PetscBool trans, PetscBool herm) in MatMultAddKernel_SeqAIJCUSPARSE() argument
3568 …PetscCheck(!herm || trans, PetscObjectComm((PetscObject)A), PETSC_ERR_GPU, "Hermitian and not tran… in MatMultAddKernel_SeqAIJCUSPARSE()
3580 if (herm || !A->form_explicit_transpose) { in MatMultAddKernel_SeqAIJCUSPARSE()
3581 opA = herm ? CUSPARSE_OPERATION_CONJUGATE_TRANSPOSE : CUSPARSE_OPERATION_TRANSPOSE; in MatMultAddKernel_SeqAIJCUSPARSE()
|