Lines Matching refs:ICNTL

507 #define ICNTL(I) icntl[(I) - 1]  macro
658 if (mumps->id.internal_id) mumps->id.ICNTL(19) = 0; // sometimes, the inner id is yet built in MatMumpsResetSchur_Private()
688 if (!mumps->id.ICNTL(9)) { /* transpose solve */ in MatMumpsSolveSchur_Private()
707 if (!mumps->id.ICNTL(9)) { /* transpose solve */ in MatMumpsSolveSchur_Private()
736 if (!mumps->id.ICNTL(19)) { /* do nothing when Schur complement has not been computed */ in MatMumpsHandleSchur_Private()
752 mumps->id.ICNTL(26) = 2; /* expansion phase */ in MatMumpsHandleSchur_Private()
756 mumps->id.ICNTL(26) = -1; in MatMumpsHandleSchur_Private()
1689 if (mumps->id.ICNTL(15) == 1) { in MatConvertToTriples_nest_xaij()
1716 } else mumps->id.ICNTL(15) = 0; in MatConvertToTriples_nest_xaij()
1949 mumps->id.ICNTL(20) = 10; /* dense distributed RHS, need to set rhs_loc[], irhs_loc[] */ in MatSolve_MUMPS()
1953 mumps->id.ICNTL(20) = 0; /* dense centralized RHS; Scatter b into a sequential b_seq vector*/ in MatSolve_MUMPS()
1962 mumps->id.ICNTL(20) = 0; in MatSolve_MUMPS()
1977 if (mumps->id.ICNTL(26) < 0 || mumps->id.ICNTL(26) > 2) { in MatSolve_MUMPS()
1980 mumps->id.ICNTL(26) = 1; /* condensation phase */ in MatSolve_MUMPS()
1981 } else if (mumps->id.ICNTL(26) == 1) PetscCall(MatMumpsHandleSchur_Private(A, PETSC_FALSE)); in MatSolve_MUMPS()
1990 else if (mumps->id.ICNTL(26) == 1) { // condense the right hand side in MatSolve_MUMPS()
1996 if (mumps->scat_sol && mumps->ICNTL9_pre != mumps->id.ICNTL(9)) { in MatSolve_MUMPS()
2009 mumps->ICNTL9_pre = mumps->id.ICNTL(9); /* save current value of id.ICNTL(9) */ in MatSolve_MUMPS()
2037 const PetscMUMPSInt value = mumps->id.ICNTL(9); in MatSolveTranspose_MUMPS()
2040 mumps->id.ICNTL(9) = 0; in MatSolveTranspose_MUMPS()
2042 mumps->id.ICNTL(9) = value; in MatSolveTranspose_MUMPS()
2078 mumps->id.ICNTL(20) = 0; /* dense RHS */ in MatMatSolve_MUMPS()
2087 mumps->id.ICNTL(20) = 1; /* sparse RHS */ in MatMatSolve_MUMPS()
2117 if (mumps->id.ICNTL(26) < 0 || mumps->id.ICNTL(26) > 2) { in MatMatSolve_MUMPS()
2120 …mumps->id.ICNTL(26) = 1; /* condensation phase, i.e, to solve id.re… in MatMatSolve_MUMPS()
2121 } else if (mumps->id.ICNTL(26) == 1) PetscCall(MatMumpsHandleSchur_Private(A, PETSC_FALSE)); in MatMatSolve_MUMPS()
2130 else if (mumps->id.ICNTL(26) == 1) { // condense the right hand side in MatMatSolve_MUMPS()
2149 …PetscCheck(!mumps->id.ICNTL(19), PetscObjectComm((PetscObject)A), PETSC_ERR_SUP, "Parallel Schur c… in MatMatSolve_MUMPS()
2169 mumps->id.ICNTL(20) = 10; /* dense distributed RHS */ in MatMatSolve_MUMPS()
2176 mumps->id.ICNTL(20) = 0; /* dense centralized RHS */ in MatMatSolve_MUMPS()
2322 const PetscMUMPSInt value = mumps->id.ICNTL(9); in MatMatSolveTranspose_MUMPS()
2325 mumps->id.ICNTL(9) = 0; in MatMatSolveTranspose_MUMPS()
2327 mumps->id.ICNTL(9) = value; in MatMatSolveTranspose_MUMPS()
2365 …PetscCheck(size <= 1 || mumps->id.ICNTL(13) == 1, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "ICNTL(13)… in MatGetInertia_SBAIJMUMPS()
2369 …PetscCheck(mumps->id.ICNTL(24) == 1, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "-mat_mumps_icntl_24 mu… in MatGetInertia_SBAIJMUMPS()
2513 if (!mumps->id.ICNTL(18)) { /* A is centralized */ in MatFactorNumeric_MUMPS()
2527 …if (mumps->id.ICNTL(22)) PetscCall(PetscStrncpy(mumps->id.ooc_prefix, ((PetscObject)F)->prefix, si… in MatFactorNumeric_MUMPS()
2546 …PetscCheck(mumps->myid || mumps->id.ICNTL(16) <= 0, PETSC_COMM_SELF, PETSC_ERR_LIB, "MUMPS error i… in MatFactorNumeric_MUMPS()
2558 if (mumps->id.ICNTL(19) == 1) { /* stored by rows */ in MatFactorNumeric_MUMPS()
2559 mumps->id.ICNTL(19) = 2; in MatFactorNumeric_MUMPS()
2566 if (!mumps->sym && mumps->id.ICNTL(19) && mumps->id.ICNTL(19) != 1) mumps->id.ICNTL(19) = 3; in MatFactorNumeric_MUMPS()
2606 …PetscBool schur = mumps->id.icntl ? (PetscBool)(mumps->id.ICNTL(26) == -1) : (PetscBool)(mumps… in MatSetFromOptions_MUMPS()
2679 mumps->id.ICNTL(3) = 0; in MatSetFromOptions_MUMPS()
2680 mumps->id.ICNTL(4) = 0; in MatSetFromOptions_MUMPS()
2682 mumps->id.ICNTL(18) = 0; /* centralized assembled matrix input */ in MatSetFromOptions_MUMPS()
2683 mumps->id.ICNTL(7) = 7; /* automatic choice of ordering done by the package */ in MatSetFromOptions_MUMPS()
2685 mumps->id.ICNTL(18) = 3; /* distributed assembled matrix input */ in MatSetFromOptions_MUMPS()
2686 mumps->id.ICNTL(21) = 1; /* distributed solution */ in MatSetFromOptions_MUMPS()
2689 mumps->id.ICNTL(15) = 1; in MatSetFromOptions_MUMPS()
2693 } else mumps->id.ICNTL(15) = 0; in MatSetFromOptions_MUMPS()
2696 …for (icntl = 0; icntl < nICNTL_pre; ++icntl) mumps->id.ICNTL(mumps->ICNTL_pre[1 + 2 * icntl]) = mu… in MatSetFromOptions_MUMPS()
2710 …mumps->id.ICNTL(19) = 1; … in MatSetFromOptions_MUMPS()
2716 mumps->id.ICNTL(19) = 3; /* MUMPS returns full matrix */ in MatSetFromOptions_MUMPS()
2718 mumps->id.ICNTL(19) = 2; /* MUMPS returns lower triangular part */ in MatSetFromOptions_MUMPS()
2721 mumps->id.ICNTL(26) = -1; in MatSetFromOptions_MUMPS()
2733 …icntl_1", "ICNTL(1): output stream for error messages", "None", mumps->id.ICNTL(1), &icntl, &flg)); in MatSetFromOptions_MUMPS()
2734 if (flg) mumps->id.ICNTL(1) = icntl; in MatSetFromOptions_MUMPS()
2735 …ream for diagnostic printing, statistics, and warning", "None", mumps->id.ICNTL(2), &icntl, &flg)); in MatSetFromOptions_MUMPS()
2736 if (flg) mumps->id.ICNTL(2) = icntl; in MatSetFromOptions_MUMPS()
2737 … stream for global information, collected on the host", "None", mumps->id.ICNTL(3), &icntl, &flg)); in MatSetFromOptions_MUMPS()
2738 if (flg) mumps->id.ICNTL(3) = icntl; in MatSetFromOptions_MUMPS()
2740 …mumps_icntl_4", "ICNTL(4): level of printing (0 to 4)", "None", mumps->id.ICNTL(4), &icntl, &flg)); in MatSetFromOptions_MUMPS()
2741 if (flg) mumps->id.ICNTL(4) = icntl; in MatSetFromOptions_MUMPS()
2742 …if (mumps->id.ICNTL(4) || PetscLogPrintInfo) mumps->id.ICNTL(3) = 6; /* resume MUMPS default id.IC… in MatSetFromOptions_MUMPS()
2744 …a zero-free diagonal and/or scale the matrix (0 to 7)", "None", mumps->id.ICNTL(6), &icntl, &flg)); in MatSetFromOptions_MUMPS()
2745 if (flg) mumps->id.ICNTL(6) = icntl; in MatSetFromOptions_MUMPS()
2747 …=Scotch, 4=PORD, 5=Metis, 6=QAMD, and 7=auto(default)", "None", mumps->id.ICNTL(7), &icntl, &flg)); in MatSetFromOptions_MUMPS()
2750 mumps->id.ICNTL(7) = icntl; in MatSetFromOptions_MUMPS()
2753 …ICNTL(8): scaling strategy (-2 to 8 or 77)", "None", mumps->id.ICNTL(8), &mumps->id.ICNTL(8), NULL… in MatSetFromOptions_MUMPS()
2755 …tl_10", "ICNTL(10): max num of refinements", "None", mumps->id.ICNTL(10), &mumps->id.ICNTL(10), NU… in MatSetFromOptions_MUMPS()
2756 …lated to an error analysis (via -ksp_view)", "None", mumps->id.ICNTL(11), &mumps->id.ICNTL(11), NU… in MatSetFromOptions_MUMPS()
2757 …g strategy for symmetric matrices (0 to 3)", "None", mumps->id.ICNTL(12), &mumps->id.ICNTL(12), NU… in MatSetFromOptions_MUMPS()
2758 … node (enable ScaLAPACK) and its splitting", "None", mumps->id.ICNTL(13), &mumps->id.ICNTL(13), NU… in MatSetFromOptions_MUMPS()
2759 …ge increase in the estimated working space", "None", mumps->id.ICNTL(14), &mumps->id.ICNTL(14), NU… in MatSetFromOptions_MUMPS()
2761 if (rbs == cbs && rbs > 1) mumps->id.ICNTL(15) = (PetscMUMPSInt)-rbs; in MatSetFromOptions_MUMPS()
2762 …input matrix resulting from a block format", "None", mumps->id.ICNTL(15), &mumps->id.ICNTL(15), &f… in MatSetFromOptions_MUMPS()
2764 …if (mumps->id.ICNTL(15) < 0) PetscCheck((-mumps->id.ICNTL(15) % cbs == 0) && (-mumps->id.ICNTL(15)… in MatSetFromOptions_MUMPS()
2765 else if (mumps->id.ICNTL(15) > 0) { in MatSetFromOptions_MUMPS()
2798 … "ICNTL(19): computes the Schur complement", "None", mumps->id.ICNTL(19), &mumps->id.ICNTL(19), NU… in MatSetFromOptions_MUMPS()
2799 if (mumps->id.ICNTL(19) <= 0 || mumps->id.ICNTL(19) > 3) { /* reset any Schur data (if any) */ in MatSetFromOptions_MUMPS()
2822 …t-of-core factorization and solve (0 or 1)", "None", mumps->id.ICNTL(22), &mumps->id.ICNTL(22), &f… in MatSetFromOptions_MUMPS()
2823 …if (flg && mumps->id.ICNTL(22) != 1) mumps->id.ICNTL(22) = 0; // MUMPS treats values other than 1 … in MatSetFromOptions_MUMPS()
2824 if (mumps->id.ICNTL(22)) { in MatSetFromOptions_MUMPS()
2829 …emory (MB) that can allocate per processor", "None", mumps->id.ICNTL(23), &mumps->id.ICNTL(23), NU… in MatSetFromOptions_MUMPS()
2830 …24): detection of null pivot rows (0 or 1)", "None", mumps->id.ICNTL(24), &mumps->id.ICNTL(24), NU… in MatSetFromOptions_MUMPS()
2831 …if (mumps->id.ICNTL(24)) mumps->id.ICNTL(13) = 1; /* turn-off ScaLAPACK to help with the correct d… in MatSetFromOptions_MUMPS()
2833 … a deficient matrix and a null space basis", "None", mumps->id.ICNTL(25), &mumps->id.ICNTL(25), NU… in MatSetFromOptions_MUMPS()
2834 …olution phase if a Schur complement matrix", "None", mumps->id.ICNTL(26), &mumps->id.ICNTL(26), NU… in MatSetFromOptions_MUMPS()
2835 …locking size for multiple right-hand sides", "None", mumps->id.ICNTL(27), &mumps->id.ICNTL(27), NU… in MatSetFromOptions_MUMPS()
2836 …r parallel analysis and ICNTL(29) ordering", "None", mumps->id.ICNTL(28), &mumps->id.ICNTL(28), NU… in MatSetFromOptions_MUMPS()
2837 …rallel ordering 1 = ptscotch, 2 = parmetis", "None", mumps->id.ICNTL(29), &mumps->id.ICNTL(29), NU… in MatSetFromOptions_MUMPS()
2839 …tors may be discarded during factorization", "None", mumps->id.ICNTL(31), &mumps->id.ICNTL(31), NU… in MatSetFromOptions_MUMPS()
2841 …icntl_33", "ICNTL(33): compute determinant", "None", mumps->id.ICNTL(33), &mumps->id.ICNTL(33), NU… in MatSetFromOptions_MUMPS()
2842 …s Block Low Rank (BLR) based factorization", "None", mumps->id.ICNTL(35), &mumps->id.ICNTL(35), NU… in MatSetFromOptions_MUMPS()
2843 …L(36): choice of BLR factorization variant", "None", mumps->id.ICNTL(36), &mumps->id.ICNTL(36), NU… in MatSetFromOptions_MUMPS()
2844 …ompression of the contribution blocks (CB)", "None", mumps->id.ICNTL(37), &mumps->id.ICNTL(37), NU… in MatSetFromOptions_MUMPS()
2845 …ed compression rate of LU factors with BLR", "None", mumps->id.ICNTL(38), &mumps->id.ICNTL(38), NU… in MatSetFromOptions_MUMPS()
2846 …(48): multithreading with tree parallelism", "None", mumps->id.ICNTL(48), &mumps->id.ICNTL(48), NU… in MatSetFromOptions_MUMPS()
2847 …orkarray at the end of factorization phase", "None", mumps->id.ICNTL(49), &mumps->id.ICNTL(49), NU… in MatSetFromOptions_MUMPS()
2848 …ostponing and rank-revealing factorization", "None", mumps->id.ICNTL(56), &mumps->id.ICNTL(56), NU… in MatSetFromOptions_MUMPS()
2849 …defines options for symbolic factorization", "None", mumps->id.ICNTL(58), &mumps->id.ICNTL(58), NU… in MatSetFromOptions_MUMPS()
2919 switch (mumps->id.ICNTL(18)) { in MatLUFactorSymbolic_AIJMUMPS()
2925 …if (1 < mumps->id.ICNTL(6) && mumps->id.ICNTL(6) < 7) PetscCall(MatMumpsMakeMumpsScalarArray(PETSC… in MatLUFactorSymbolic_AIJMUMPS()
2926 if (r && mumps->id.ICNTL(7) == 7) { in MatLUFactorSymbolic_AIJMUMPS()
2927 mumps->id.ICNTL(7) = 1; in MatLUFactorSymbolic_AIJMUMPS()
2944 …if (1 < mumps->id.ICNTL(6) && mumps->id.ICNTL(6) < 7) PetscCall(MatMumpsMakeMumpsScalarArray(PETSC… in MatLUFactorSymbolic_AIJMUMPS()
2988 switch (mumps->id.ICNTL(18)) { in MatLUFactorSymbolic_BAIJMUMPS()
2994 …if (1 < mumps->id.ICNTL(6) && mumps->id.ICNTL(6) < 7) PetscCall(MatMumpsMakeMumpsScalarArray(PETSC… in MatLUFactorSymbolic_BAIJMUMPS()
3001 …if (1 < mumps->id.ICNTL(6) && mumps->id.ICNTL(6) < 7) PetscCall(MatMumpsMakeMumpsScalarArray(PETSC… in MatLUFactorSymbolic_BAIJMUMPS()
3043 switch (mumps->id.ICNTL(18)) { in MatCholeskyFactorSymbolic_MUMPS()
3049 …if (1 < mumps->id.ICNTL(6) && mumps->id.ICNTL(6) < 7) PetscCall(MatMumpsMakeMumpsScalarArray(PETSC… in MatCholeskyFactorSymbolic_MUMPS()
3056 …if (1 < mumps->id.ICNTL(6) && mumps->id.ICNTL(6) < 7) PetscCall(MatMumpsMakeMumpsScalarArray(PETSC… in MatCholeskyFactorSymbolic_MUMPS()
3101 …PetscViewerASCIIPrintf(viewer, " ICNTL(1) (output for error): %d\n", mumps->id.ICNTL(1))); in MatView_MUMPS()
3102 …PetscViewerASCIIPrintf(viewer, " ICNTL(2) (output of diagnostic msg): %d\n", mumps->id.ICNTL(2))); in MatView_MUMPS()
3103 …PetscViewerASCIIPrintf(viewer, " ICNTL(3) (output for global info): %d\n", mumps->id.ICNTL(3))); in MatView_MUMPS()
3104 …PetscViewerASCIIPrintf(viewer, " ICNTL(4) (level of printing): %d\n", mumps->id.ICNTL(4))); in MatView_MUMPS()
3105 …PetscViewerASCIIPrintf(viewer, " ICNTL(5) (input mat struct): %d\n", mumps->id.ICNTL(5))); in MatView_MUMPS()
3106 …PetscViewerASCIIPrintf(viewer, " ICNTL(6) (matrix prescaling): %d\n", mumps->id.ICNTL(6))); in MatView_MUMPS()
3107 …etscViewerASCIIPrintf(viewer, " ICNTL(7) (sequential matrix ordering):%d\n", mumps->id.ICNTL(7))); in MatView_MUMPS()
3108 …PetscViewerASCIIPrintf(viewer, " ICNTL(8) (scaling strategy): %d\n", mumps->id.ICNTL(8))); in MatView_MUMPS()
3109 …etscViewerASCIIPrintf(viewer, " ICNTL(10) (max num of refinements): %d\n", mumps->id.ICNTL(10))); in MatView_MUMPS()
3110 …etscViewerASCIIPrintf(viewer, " ICNTL(11) (error analysis): %d\n", mumps->id.ICNTL(11))); in MatView_MUMPS()
3111 if (mumps->id.ICNTL(11) > 0) { in MatView_MUMPS()
3119 …tf(viewer, " ICNTL(12) (efficiency control): %d\n", mumps->id.ICNTL(12))); in MatView_MUMPS()
3120 …tf(viewer, " ICNTL(13) (sequential factorization of the root node): %d\n", mumps->id.ICNTL(13))); in MatView_MUMPS()
3121 …tf(viewer, " ICNTL(14) (percentage of estimated workspace increase): %d\n", mumps->id.ICNTL(14))); in MatView_MUMPS()
3122 …tf(viewer, " ICNTL(15) (compression of the input matrix): %d\n", mumps->id.ICNTL(15))); in MatView_MUMPS()
3124 …tf(viewer, " ICNTL(18) (input mat struct): %d\n", mumps->id.ICNTL(18))); in MatView_MUMPS()
3125 …tf(viewer, " ICNTL(19) (Schur complement info): %d\n", mumps->id.ICNTL(19))); in MatView_MUMPS()
3126 …tf(viewer, " ICNTL(20) (RHS sparse pattern): %d\n", mumps->id.ICNTL(20))); in MatView_MUMPS()
3127 …tf(viewer, " ICNTL(21) (solution struct): %d\n", mumps->id.ICNTL(21))); in MatView_MUMPS()
3128 …tf(viewer, " ICNTL(22) (in-core/out-of-core facility): %d\n", mumps->id.ICNTL(22))); in MatView_MUMPS()
3129 …tf(viewer, " ICNTL(23) (max size of memory can be allocated locally):%d\n", mumps->id.ICNTL(23))); in MatView_MUMPS()
3131 …tf(viewer, " ICNTL(24) (detection of null pivot rows): %d\n", mumps->id.ICNTL(24))); in MatView_MUMPS()
3132 …tf(viewer, " ICNTL(25) (computation of a null space basis): %d\n", mumps->id.ICNTL(25))); in MatView_MUMPS()
3133 …tf(viewer, " ICNTL(26) (Schur options for RHS or solution): %d\n", mumps->id.ICNTL(26))); in MatView_MUMPS()
3134 …tf(viewer, " ICNTL(27) (blocking size for multiple RHS): %d\n", mumps->id.ICNTL(27))); in MatView_MUMPS()
3135 …tf(viewer, " ICNTL(28) (use parallel or sequential ordering): %d\n", mumps->id.ICNTL(28))); in MatView_MUMPS()
3136 …tf(viewer, " ICNTL(29) (parallel ordering): %d\n", mumps->id.ICNTL(29))); in MatView_MUMPS()
3138 …tf(viewer, " ICNTL(30) (user-specified set of entries in inv(A)): %d\n", mumps->id.ICNTL(30))); in MatView_MUMPS()
3139 …tf(viewer, " ICNTL(31) (factors is discarded in the solve phase): %d\n", mumps->id.ICNTL(31))); in MatView_MUMPS()
3140 …tf(viewer, " ICNTL(33) (compute determinant): %d\n", mumps->id.ICNTL(33))); in MatView_MUMPS()
3141 …tf(viewer, " ICNTL(35) (activate BLR based factorization): %d\n", mumps->id.ICNTL(35))); in MatView_MUMPS()
3142 …tf(viewer, " ICNTL(36) (choice of BLR factorization variant): %d\n", mumps->id.ICNTL(36))); in MatView_MUMPS()
3143 …tf(viewer, " ICNTL(37) (compression of the contribution blocks): %d\n", mumps->id.ICNTL(37))); in MatView_MUMPS()
3144 …tf(viewer, " ICNTL(38) (estimated compression rate of LU factors): %d\n", mumps->id.ICNTL(38))); in MatView_MUMPS()
3145 …tf(viewer, " ICNTL(48) (multithreading with tree parallelism): %d\n", mumps->id.ICNTL(48))); in MatView_MUMPS()
3146 …r, " ICNTL(49) (compact workarray at the end of factorization phase):%d\n", mumps->id.ICNTL(49))); in MatView_MUMPS()
3147 …tf(viewer, " ICNTL(56) (postponing and rank-revealing factorization):%d\n", mumps->id.ICNTL(56))); in MatView_MUMPS()
3148 …tf(viewer, " ICNTL(58) (options for symbolic factorization): %d\n", mumps->id.ICNTL(58))); in MatView_MUMPS()
3279 if (mumps->id.icntl) mumps->id.ICNTL(26) = -1; in MatFactorSetSchurIS_MUMPS()
3292 …PetscCheck(mumps->id.ICNTL(19), PetscObjectComm((PetscObject)F), PETSC_ERR_ORDER, "Schur complemen… in MatFactorCreateSchurComplement_MUMPS()
3299 if (mumps->id.ICNTL(19) == 1) { /* stored by rows */ in MatFactorCreateSchurComplement_MUMPS()
3307 if (mumps->id.ICNTL(19) == 2) { /* lower triangular stored by columns */ in MatFactorCreateSchurComplement_MUMPS()
3311 } else if (mumps->id.ICNTL(19) == 3) { /* full matrix */ in MatFactorCreateSchurComplement_MUMPS()
3340 } else PetscCall(PetscMUMPSIntCast(ival, &mumps->id.ICNTL(icntl))); in MatMumpsSetIcntl_MUMPS()
3355 } else *ival = mumps->id.ICNTL(icntl); in MatMumpsGetIcntl_MUMPS()
3612 …PetscCheck(mumps->id.ICNTL(24) == 1, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "-mat_mumps_icntl_24 mu… in MatMumpsGetNullPivots_MUMPS()
3667 mumps->id.ICNTL(20) = 1; /* rhs is sparse */ in MatMumpsGetInverse_MUMPS()
3668 mumps->id.ICNTL(21) = 0; /* solution is in assembled centralized format */ in MatMumpsGetInverse_MUMPS()