mpi/mkl_cpardiso/mkl_cpardiso.c

508   PetscInt          icntl, threads;
516   PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_66", "Maximum number of factors with identical sparsity structure that must be kept in memory at the same time", "None", mat_mkl_cpardiso->maxfct, &icntl, &flg));
517   if (flg) mat_mkl_cpardiso->maxfct = icntl;
519   PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_67", "Indicates the actual matrix for the solution phase", "None", mat_mkl_cpardiso->mnum, &icntl, &flg));
520   if (flg) mat_mkl_cpardiso->mnum = icntl;
522   PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_68", "Message level information", "None", mat_mkl_cpardiso->msglvl, &icntl, &flg));
523   if (flg) mat_mkl_cpardiso->msglvl = icntl;
525   PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_69", "Defines the matrix type", "None", mat_mkl_cpardiso->mtype, &icntl, &flg));
526   if (flg) mat_mkl_cpardiso->mtype = icntl;
527   PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_1", "Use default values", "None", mat_mkl_cpardiso->iparm[0], &icntl, &flg));
529   if (flg && icntl != 0) {
530     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_2", "Fill-in reducing ordering for the input matrix", "None", mat_mkl_cpardiso->iparm[1], &icntl, &flg));
531     if (flg) mat_mkl_cpardiso->iparm[1] = icntl;
533     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_4", "Preconditioned CGS/CG", "None", mat_mkl_cpardiso->iparm[3], &icntl, &flg));
534     if (flg) mat_mkl_cpardiso->iparm[3] = icntl;
536     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_5", "User permutation", "None", mat_mkl_cpardiso->iparm[4], &icntl, &flg));
537     if (flg) mat_mkl_cpardiso->iparm[4] = icntl;
539     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_6", "Write solution on x", "None", mat_mkl_cpardiso->iparm[5], &icntl, &flg));
540     if (flg) mat_mkl_cpardiso->iparm[5] = icntl;
542     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_8", "Iterative refinement step", "None", mat_mkl_cpardiso->iparm[7], &icntl, &flg));
543     if (flg) mat_mkl_cpardiso->iparm[7] = icntl;
545     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_10", "Pivoting perturbation", "None", mat_mkl_cpardiso->iparm[9], &icntl, &flg));
546     if (flg) mat_mkl_cpardiso->iparm[9] = icntl;
548     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_11", "Scaling vectors", "None", mat_mkl_cpardiso->iparm[10], &icntl, &flg));
549     if (flg) mat_mkl_cpardiso->iparm[10] = icntl;
551     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_12", "Solve with transposed or conjugate transposed matrix A", "None", mat_mkl_cpardiso->iparm[11], &icntl, &flg));
552     if (flg) mat_mkl_cpardiso->iparm[11] = icntl;
554     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_13", "Improved accuracy using (non-) symmetric weighted matching", "None", mat_mkl_cpardiso->iparm[12], &icntl, &flg));
555     if (flg) mat_mkl_cpardiso->iparm[12] = icntl;
557     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_18", "Numbers of non-zero elements", "None", mat_mkl_cpardiso->iparm[17], &icntl, &flg));
558     if (flg) mat_mkl_cpardiso->iparm[17] = icntl;
560     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_19", "Report number of floating point operations", "None", mat_mkl_cpardiso->iparm[18], &icntl, &flg));
561     if (flg) mat_mkl_cpardiso->iparm[18] = icntl;
563     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_21", "Pivoting for symmetric indefinite matrices", "None", mat_mkl_cpardiso->iparm[20], &icntl, &flg));
564     if (flg) mat_mkl_cpardiso->iparm[20] = icntl;
566     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_24", "Parallel factorization control", "None", mat_mkl_cpardiso->iparm[23], &icntl, &flg));
567     if (flg) mat_mkl_cpardiso->iparm[23] = icntl;
569     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_25", "Parallel forward/backward solve control", "None", mat_mkl_cpardiso->iparm[24], &icntl, &flg));
570     if (flg) mat_mkl_cpardiso->iparm[24] = icntl;
572     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_27", "Matrix checker", "None", mat_mkl_cpardiso->iparm[26], &icntl, &flg));
573     if (flg) mat_mkl_cpardiso->iparm[26] = icntl;
575     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_31", "Partial solve and computing selected components of the solution vectors", "None", mat_mkl_cpardiso->iparm[30], &icntl, &flg));
576     if (flg) mat_mkl_cpardiso->iparm[30] = icntl;
578     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_34", "Optimal number of threads for conditional numerical reproducibility (CNR) mode", "None", mat_mkl_cpardiso->iparm[33], &icntl, &flg));
579     if (flg) mat_mkl_cpardiso->iparm[33] = icntl;
581     PetscCall(PetscOptionsInt("-mat_mkl_cpardiso_60", "Intel MKL Cluster PARDISO mode", "None", mat_mkl_cpardiso->iparm[59], &icntl, &flg));
582     if (flg) mat_mkl_cpardiso->iparm[59] = icntl;
764 static PetscErrorCode MatMkl_CPardisoSetCntl_MKL_CPARDISO(Mat F, PetscInt icntl, PetscInt ival)
769   if (icntl <= 64) {
770     mat_mkl_cpardiso->iparm[icntl - 1] = ival;
772     if (icntl == 65) mkl_set_num_threads((int)ival);
773     else if (icntl == 66) mat_mkl_cpardiso->maxfct = ival;
774     else if (icntl == 67) mat_mkl_cpardiso->mnum = ival;
775     else if (icntl == 68) mat_mkl_cpardiso->msglvl = ival;
776     else if (icntl == 69) mat_mkl_cpardiso->mtype = ival;
789 . icntl - index of MKL Cluster PARDISO parameter
793 . -mat_mkl_cpardiso_<icntl> <ival> - set the option numbered icntl to ival
803 PetscErrorCode MatMkl_CPardisoSetCntl(Mat F, PetscInt icntl, PetscInt ival)
806   PetscTryMethod(F, "MatMkl_CPardisoSetCntl_C", (Mat, PetscInt, PetscInt), (F, icntl, ival));