Lines Matching refs:iparm
41 void pardiso_64init(void *pt, INT_TYPE *mtype, INT_TYPE iparm[]) in pardiso_64init() argument
47 for (PetscInt i = 0; i < IPARM_SIZE; i++) iparm[i] = iparm_copy[i]; in pardiso_64init()
63 INT_TYPE iparm[IPARM_SIZE]; member
254 …PetscCheck(mpardiso->iparm[12 - 1] != 1, PetscObjectComm((PetscObject)F), PETSC_ERR_SUP, "Hermitia… in MatMKLPardisoSolveSchur_Private()
259 if (!mpardiso->iparm[12 - 1]) { in MatMKLPardisoSolveSchur_Private()
267 if (!mpardiso->iparm[12 - 1]) { in MatMKLPardisoSolveSchur_Private()
321 mpardiso->iparm[36 - 1] = 2; in MatFactorSetSchurIS_MKL_PARDISO()
334 …phase, &mat_mkl_pardiso->n, NULL, NULL, NULL, NULL, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, in MatDestroy_MKL_PARDISO()
342 if (mat_mkl_pardiso->iparm[34] == 1) PetscCall(PetscFree(mat_mkl_pardiso->a)); in MatDestroy_MKL_PARDISO()
397 mat_mkl_pardiso->iparm[6 - 1] = 1; in MatSolve_MKL_PARDISO()
399 …&mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (void *)xarray, (void *)… in MatSolve_MKL_PARDISO()
402 mat_mkl_pardiso->iparm[6 - 1] = 0; in MatSolve_MKL_PARDISO()
404 …mat_mkl_pardiso->perm, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (… in MatSolve_MKL_PARDISO()
428 mat_mkl_pardiso->iparm[6 - 1] = 1; in MatSolve_MKL_PARDISO()
431 …mat_mkl_pardiso->perm, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (… in MatSolve_MKL_PARDISO()
434 mat_mkl_pardiso->iparm[6 - 1] = 0; in MatSolve_MKL_PARDISO()
457 …&mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (void *)barray, (void *)… in MatForwardSolve_MKL_PARDISO()
482 …&mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (void *)barray, (void *)… in MatBackwardSolve_MKL_PARDISO()
497 oiparm12 = mat_mkl_pardiso->iparm[12 - 1]; in MatSolveTranspose_MKL_PARDISO()
498 mat_mkl_pardiso->iparm[12 - 1] = 2; in MatSolveTranspose_MKL_PARDISO()
500 mat_mkl_pardiso->iparm[12 - 1] = oiparm12; in MatSolveTranspose_MKL_PARDISO()
530 …mat_mkl_pardiso->perm, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (… in MatMatSolve_MKL_PARDISO()
565 mat_mkl_pardiso->iparm[6 - 1] = 1; in MatMatSolve_MKL_PARDISO()
568 …mat_mkl_pardiso->perm, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (… in MatMatSolve_MKL_PARDISO()
575 mat_mkl_pardiso->iparm[6 - 1] = 0; in MatMatSolve_MKL_PARDISO()
593 …&mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, NULL, (void *)mat_mkl_pa… in MatFactorNumeric_MKL_PARDISO()
597 …if (mat_mkl_pardiso->iparm[18] > 0) PetscCall(PetscLogFlops(PetscPowRealInt(10., 6) * mat_mkl_pard… in MatFactorNumeric_MKL_PARDISO()
636 icntl = mat_mkl_pardiso->iparm[34]; in MatSetFromOptions_MKL_PARDISO()
637 bs = mat_mkl_pardiso->iparm[36]; in MatSetFromOptions_MKL_PARDISO()
638 MKL_PARDISO_INIT(pt, &mat_mkl_pardiso->mtype, mat_mkl_pardiso->iparm); in MatSetFromOptions_MKL_PARDISO()
640 mat_mkl_pardiso->iparm[27] = 1; in MatSetFromOptions_MKL_PARDISO()
642 mat_mkl_pardiso->iparm[27] = 0; in MatSetFromOptions_MKL_PARDISO()
644 mat_mkl_pardiso->iparm[34] = icntl; in MatSetFromOptions_MKL_PARDISO()
645 mat_mkl_pardiso->iparm[36] = bs; in MatSetFromOptions_MKL_PARDISO()
648 …mat_mkl_pardiso_1", "Use default values (if 0)", "None", mat_mkl_pardiso->iparm[0], &icntl, &flg)); in MatSetFromOptions_MKL_PARDISO()
649 if (flg) mat_mkl_pardiso->iparm[0] = icntl; in MatSetFromOptions_MKL_PARDISO()
651 …Fill-in reducing ordering for the input matrix", "None", mat_mkl_pardiso->iparm[1], &icntl, &flg)); in MatSetFromOptions_MKL_PARDISO()
652 if (flg) mat_mkl_pardiso->iparm[1] = icntl; in MatSetFromOptions_MKL_PARDISO()
654 …t("-mat_mkl_pardiso_4", "Preconditioned CGS/CG", "None", mat_mkl_pardiso->iparm[3], &icntl, &flg)); in MatSetFromOptions_MKL_PARDISO()
655 if (flg) mat_mkl_pardiso->iparm[3] = icntl; in MatSetFromOptions_MKL_PARDISO()
657 …onsInt("-mat_mkl_pardiso_5", "User permutation", "None", mat_mkl_pardiso->iparm[4], &icntl, &flg)); in MatSetFromOptions_MKL_PARDISO()
658 if (flg) mat_mkl_pardiso->iparm[4] = icntl; in MatSetFromOptions_MKL_PARDISO()
660 …Int("-mat_mkl_pardiso_6", "Write solution on x", "None", mat_mkl_pardiso->iparm[5], &icntl, &flg)); in MatSetFromOptions_MKL_PARDISO()
661 if (flg) mat_mkl_pardiso->iparm[5] = icntl; in MatSetFromOptions_MKL_PARDISO()
663 …mat_mkl_pardiso_8", "Iterative refinement step", "None", mat_mkl_pardiso->iparm[7], &icntl, &flg)); in MatSetFromOptions_MKL_PARDISO()
664 if (flg) mat_mkl_pardiso->iparm[7] = icntl; in MatSetFromOptions_MKL_PARDISO()
666 …("-mat_mkl_pardiso_10", "Pivoting perturbation", "None", mat_mkl_pardiso->iparm[9], &icntl, &flg)); in MatSetFromOptions_MKL_PARDISO()
667 if (flg) mat_mkl_pardiso->iparm[9] = icntl; in MatSetFromOptions_MKL_PARDISO()
669 …onsInt("-mat_mkl_pardiso_11", "Scaling vectors", "None", mat_mkl_pardiso->iparm[10], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
670 if (flg) mat_mkl_pardiso->iparm[10] = icntl; in MatSetFromOptions_MKL_PARDISO()
672 …th transposed or conjugate transposed matrix A", "None", mat_mkl_pardiso->iparm[11], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
673 if (flg) mat_mkl_pardiso->iparm[11] = icntl; in MatSetFromOptions_MKL_PARDISO()
675 …uracy using (non-) symmetric weighted matching", "None", mat_mkl_pardiso->iparm[12], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
676 if (flg) mat_mkl_pardiso->iparm[12] = icntl; in MatSetFromOptions_MKL_PARDISO()
678 …mkl_pardiso_18", "Numbers of non-zero elements", "None", mat_mkl_pardiso->iparm[17], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
679 if (flg) mat_mkl_pardiso->iparm[17] = icntl; in MatSetFromOptions_MKL_PARDISO()
681 …er of floating point operations (0 to disable)", "None", mat_mkl_pardiso->iparm[18], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
682 if (flg) mat_mkl_pardiso->iparm[18] = icntl; in MatSetFromOptions_MKL_PARDISO()
684 …", "Pivoting for symmetric indefinite matrices", "None", mat_mkl_pardiso->iparm[20], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
685 if (flg) mat_mkl_pardiso->iparm[20] = icntl; in MatSetFromOptions_MKL_PARDISO()
687 …l_pardiso_24", "Parallel factorization control", "None", mat_mkl_pardiso->iparm[23], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
688 if (flg) mat_mkl_pardiso->iparm[23] = icntl; in MatSetFromOptions_MKL_PARDISO()
690 …_25", "Parallel forward/backward solve control", "None", mat_mkl_pardiso->iparm[24], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
691 if (flg) mat_mkl_pardiso->iparm[24] = icntl; in MatSetFromOptions_MKL_PARDISO()
693 …PetscCall(PetscOptionsInt("-mat_mkl_pardiso_27", "Matrix checker", "None", mat_mkl_pardiso->iparm[… in MatSetFromOptions_MKL_PARDISO()
694 if (flg) mat_mkl_pardiso->iparm[26] = icntl; in MatSetFromOptions_MKL_PARDISO()
696 …ng selected components of the solution vectors", "None", mat_mkl_pardiso->iparm[30], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
697 if (flg) mat_mkl_pardiso->iparm[30] = icntl; in MatSetFromOptions_MKL_PARDISO()
699 …nditional numerical reproducibility (CNR) mode", "None", mat_mkl_pardiso->iparm[33], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
700 if (flg) mat_mkl_pardiso->iparm[33] = icntl; in MatSetFromOptions_MKL_PARDISO()
702 …"-mat_mkl_pardiso_60", "Intel MKL PARDISO mode", "None", mat_mkl_pardiso->iparm[59], &icntl, &flg)… in MatSetFromOptions_MKL_PARDISO()
703 if (flg) mat_mkl_pardiso->iparm[59] = icntl; in MatSetFromOptions_MKL_PARDISO()
714 for (i = 0; i < IPARM_SIZE; i++) mat_mkl_pardiso->iparm[i] = 0; in MatFactorMKL_PARDISOInitialize_Private()
717 mat_mkl_pardiso->iparm[27] = 1; in MatFactorMKL_PARDISOInitialize_Private()
719 mat_mkl_pardiso->iparm[27] = 0; in MatFactorMKL_PARDISOInitialize_Private()
722 mat_mkl_pardiso->iparm[0] = 1; /* Solver default parameters overridden with provided by iparm */ in MatFactorMKL_PARDISOInitialize_Private()
723 mat_mkl_pardiso->iparm[1] = 2; /* Metis reordering */ in MatFactorMKL_PARDISOInitialize_Private()
724 mat_mkl_pardiso->iparm[5] = 0; /* Write solution into x */ in MatFactorMKL_PARDISOInitialize_Private()
725 mat_mkl_pardiso->iparm[7] = 0; /* Max number of iterative refinement steps */ in MatFactorMKL_PARDISOInitialize_Private()
726 mat_mkl_pardiso->iparm[17] = -1; /* Output: Number of nonzeros in the factor LU */ in MatFactorMKL_PARDISOInitialize_Private()
727 mat_mkl_pardiso->iparm[18] = -1; /* Output: Mflops for LU factorization */ in MatFactorMKL_PARDISOInitialize_Private()
729 mat_mkl_pardiso->iparm[23] = 1; /* Parallel factorization control*/ in MatFactorMKL_PARDISOInitialize_Private()
734 …mat_mkl_pardiso->iparm[34] = 1; /* Cluster Sparse Solver use C-style indexing for ia and ja arrays… in MatFactorMKL_PARDISOInitialize_Private()
737 …mat_mkl_pardiso->iparm[34] = 0; /* Cluster Sparse Solver use Fortran-style indexing for ia and ja … in MatFactorMKL_PARDISOInitialize_Private()
738 mat_mkl_pardiso->iparm[36] = bs; in MatFactorMKL_PARDISOInitialize_Private()
741 mat_mkl_pardiso->iparm[39] = 0; /* Input: matrix/rhs/solution stored on rank-0 */ in MatFactorMKL_PARDISOInitialize_Private()
755 mat_mkl_pardiso->iparm[9] = 13; /* Perturb the pivot elements with 1E-13 */ in MatFactorMKL_PARDISOInitialize_Private()
756 mat_mkl_pardiso->iparm[10] = 1; /* Use nonsymmetric permutation and scaling MPS */ in MatFactorMKL_PARDISOInitialize_Private()
757 …mat_mkl_pardiso->iparm[12] = 1; /* Switch on Maximum Weighted Matching algorithm (default for non… in MatFactorMKL_PARDISOInitialize_Private()
759 mat_mkl_pardiso->iparm[9] = 8; /* Perturb the pivot elements with 1E-8 */ in MatFactorMKL_PARDISOInitialize_Private()
760 mat_mkl_pardiso->iparm[10] = 0; /* Use nonsymmetric permutation and scaling MPS */ in MatFactorMKL_PARDISOInitialize_Private()
761 …mat_mkl_pardiso->iparm[12] = 1; /* Switch on Maximum Weighted Matching algorithm (default for non-… in MatFactorMKL_PARDISOInitialize_Private()
763 mat_mkl_pardiso->iparm[26] = 1; /* Matrix checker */ in MatFactorMKL_PARDISOInitialize_Private()
781 if (mat_mkl_pardiso->iparm[34] == 1) PetscCall(PetscFree(mat_mkl_pardiso->a)); in MatFactorSymbolic_AIJMKL_PARDISO_Private()
784 if (mat_mkl_pardiso->iparm[34] == 1) mat_mkl_pardiso->n = A->rmap->N; in MatFactorSymbolic_AIJMKL_PARDISO_Private()
790 if (mat_mkl_pardiso->iparm[18]) mat_mkl_pardiso->iparm[18] = -1; in MatFactorSymbolic_AIJMKL_PARDISO_Private()
793 …&mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, NULL, NULL, &mat_mkl_par… in MatFactorSymbolic_AIJMKL_PARDISO_Private()
824 if (nneg) *nneg = mat_mkl_pardiso->iparm[22]; in MatGetInertia_MKL_PARDISO()
825 if (npos) *npos = mat_mkl_pardiso->iparm[21]; in MatGetInertia_MKL_PARDISO()
826 if (nzero) *nzero = F->rmap->N - (mat_mkl_pardiso->iparm[22] + mat_mkl_pardiso->iparm[21]); in MatGetInertia_MKL_PARDISO()
858 …O iparm[%" PetscInt_FMT "]: %" PetscInt_FMT "\n", i, (PetscInt)mat_mkl_pardiso->iparm[i - 1])); in MatView_MKL_PARDISO()
876 info->nz_used = mat_mkl_pardiso->iparm[17]; in MatGetInfo_MKL_PARDISO()
877 info->nz_allocated = mat_mkl_pardiso->iparm[17]; in MatGetInfo_MKL_PARDISO()
895 mat_mkl_pardiso->iparm[icntl - 1] = ival; in MatMkl_PardisoSetCntl_MKL_PARDISO()
903 backup = mat_mkl_pardiso->iparm[34]; in MatMkl_PardisoSetCntl_MKL_PARDISO()
904 bs = mat_mkl_pardiso->iparm[36]; in MatMkl_PardisoSetCntl_MKL_PARDISO()
906 MKL_PARDISO_INIT(pt, &mat_mkl_pardiso->mtype, mat_mkl_pardiso->iparm); in MatMkl_PardisoSetCntl_MKL_PARDISO()
908 mat_mkl_pardiso->iparm[27] = 1; in MatMkl_PardisoSetCntl_MKL_PARDISO()
910 mat_mkl_pardiso->iparm[27] = 0; in MatMkl_PardisoSetCntl_MKL_PARDISO()
912 mat_mkl_pardiso->iparm[34] = backup; in MatMkl_PardisoSetCntl_MKL_PARDISO()
913 mat_mkl_pardiso->iparm[36] = bs; in MatMkl_PardisoSetCntl_MKL_PARDISO()