Dear Matt,
Thank you for your help last time.
I want to get more detail about the Petsc-MUMPS factorisation;
so I go to look the code "/src/mat/impls/aij/mpi/mumps/mumps.c".
And I found the following functions are quite important to
the question:
PetscErrorCode MatCholeskyFactorSymbolic_MUMPS(Mat F,Mat A,IS r,const
MatFactorInfo *info);
PetscErrorCode MatFactorNumeric_MUMPS(Mat F,Mat A,const MatFactorInfo
*info);
PetscErrorCode MatSolve_MUMPS(Mat A,Vec b,Vec x);
I print some sentence to trace when these functions are called.
Then I test my code; the values in the matrix is changing but the
structure stays the same. Below is the output.
We can see that at 0th step, all the symbolic, numeric and solve
are called; in the subsequent steps only the solve stage is called,
the numeric step is not called.
Iteration 0 Step 0.0005 Time 0.0005
[INFO]: Direct Solver setup
MatCholeskyFactorSymbolic_MUMPS
finish MatCholeskyFactorSymbolic_MUMPS
MatFactorNumeric_MUMPS
finish MatFactorNumeric_MUMPS
MatSolve_MUMPS
Iteration 1 Step 0.0005 Time 0.0005
MatSolve_MUMPS
Iteration 2 Step 0.0005 Time 0.001
MatSolve_MUMPS
[INFO]: End of program!!!
I am wondering if there is any possibility to split the numeric
and solve stage (as you mentioned using KSPSolve).
Thank you very much indeed.
Kind Regards,
Shidi
On 2018-05-04 21:10, Y. Shidi wrote:
Thank you very much for your reply.
That is really clear.
Kind Regards,
Shidi
On 2018-05-04 21:05, Matthew Knepley wrote:
On Fri, May 4, 2018 at 3:54 PM, Y. Shidi <[email protected]> wrote:
Dear Matt,
Thank you very much for your reply!
So what you mean is that I can just do the KSPSolve() every
iteration
once the MUMPS is set?
Yes.
That means inside the KSPSolve() the numerical factorization is
performed. If that is the case, it seems that the ksp object is
not changed when the values in the matrix are changed.
Yes.
Or do I need to call both KSPSetOperators() and KSPSolve()?
If you do SetOperators, it will redo the factorization. If you do not,
it will look
at the Mat object, determine that the structure has not changed, and
just redo
the numerical factorization.
Thanks,
Matt
On 2018-05-04 14:44, Matthew Knepley wrote:
On Fri, May 4, 2018 at 9:40 AM, Y. Shidi <[email protected]> wrote:
Dear PETSc users,
I am currently using MUMPS to solve linear systems directly.
Generally, we use ICNTL(7) or ICNTL(29) to do the preprocessing
step and then solve the system.
In my code, the values in the matrix is changed in each iteration,
but the structure of the matrix stays the same, which means the
performance can be improved if symbolic factorisation is only
performed once. Hence, it is necessary to split the symbolic
and numeric factorisation. However, I cannot find a specific step
(control parameter) to perform the numeric factorisation.
I have used ICNTL(3) and ICNTL(4) to print the MUMPS information,
it seems that the symbolic and numeric factorisation always perform
together.
If you use KSPSolve instead, it will automatically preserve the
symbolic
factorization.
Thanks,
Matt
So I am wondering if anyone has an idea about it.
Below is how I set up MUMPS solver:
PC pc;
PetscBool flg_mumps, flg_mumps_ch;
flg_mumps = PETSC_FALSE;
flg_mumps_ch = PETSC_FALSE;
PetscOptionsGetBool(NULL, NULL, "-use_mumps_lu", &flg_mumps,
NULL);
PetscOptionsGetBool(NULL, NULL, "-use_mumps_ch", &flg_mumps_ch,
NULL);
if(flg_mumps ||flg_mumps_ch)
{
KSPSetType(_ksp, KSPPREONLY);
PetscInt ival,icntl;
PetscReal val;
KSPGetPC(_ksp, &pc);
/// Set preconditioner type
if(flg_mumps)
{
PCSetType(pc, PCLU);
}
else if(flg_mumps_ch)
{
MatSetOption(A, MAT_SPD, PETSC_TRUE);
PCSetType(pc, PCCHOLESKY);
}
PCFactorSetMatSolverPackage(pc, MATSOLVERMUMPS);
PCFactorSetUpMatSolverPackage(pc);
PCFactorGetMatrix(pc, &_F);
icntl = 7; ival = 0;
MatMumpsSetIcntl( _F, icntl, ival );
MatMumpsSetIcntl(_F, 3, 6);
MatMumpsSetIcntl(_F, 4, 2);
}
KSPSetUp(_ksp);
Kind Regards,
Shidi
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to
which
their experiments lead.
-- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/ [1] [1]
Links:
------
[1] http://www.caam.rice.edu/~mk51/ [2]
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/ [2]
Links:
------
[1] https://www.cse.buffalo.edu/~knepley/
[2] http://www.caam.rice.edu/~mk51/
