<div dir="ltr"><div>Thanks for your reply.</div><div>            
<br></div><div dir="ltr">I feel sorry that my last email was not very clear. Recently, I use Petsc to write a thermal program (including temperature, pressure, velocity field). 
I use the Matrix-Free method with user-defined precinditioning, in which the PC_TYPE is ILU (0). 
<br></div><div dir="ltr"><br></div><div dir="ltr">Due to the complexity of the direct relationship between the fields, I can only construct an approximate Jacobian matrix 
(which only contains the Jacobian matrix inside a single physical field, similar to block Jacobi preconditioning), 
and use this approximate Jacobian matrix as the preconditioning matrix.  </div><div dir="ltr">My preconditioning matrix construction method is to use SNESSetJacobian to construct the preconditioning matrix <b>B</b> (MatSetValues),
but not to operate on the Jacobian matrix <b>J</b> (due to the use of Matrix Free method).  <br></div><div dir="ltr">At the same time, the elements of the preconditioning matrix are transformed, which changes with the change of the solution vector u, 
but the position and number of non-zero elements remain unchanged. 
<br></div><div dir="ltr"><br></div><div dir="ltr">So I wonder in particular whether the preprocessing matrix will be updated in the linear step (Gmres iteration step)? 
Because the solution vector will change in each linear step, the preconditioning matrix will change.            
<br></div><div dir="ltr"><br></div><div dir="ltr">I read an article about fgmres. This method applies to variable preconditioning matrices. Is it applicable to my problem? <br></div><div dir="ltr"><br></div><div>Thanks,</div><div>Yingjie<br></div></div><br><div class="gmail_quote"><div dir="ltr">Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> 于2018年12月19日周三 下午1:51写道：<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
<br>
> On Dec 19, 2018, at 9:07 AM, Yingjie Wu via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br>
> <br>
> Respected Petsc developers:             <br>
> Hi,             <br>
> Recently, I developed a two-dimensional non-linear equations solver (SNES) using Petsc to solve temperature, velocity and pressure fields.  Since the coupling characteristics of the problem are complex, it is difficult to construct Jacobian matrix explicitly, so Matrix-Free method with Preconditioner is adopted.  I have the following questions about the Petsc's Preconditioner. <br>
> <br>
> 1. Is the preconditioning matrix updated in linear or non-linear steps?  <br>
<br>
   At each nonlinear iteration the routine set with SNESSetJacobian() is called for you to provide the new Jacobian. <br>
<br>
> Because the elements in my preconditioning matrix are functions of variables (such as density is a function of temperature and pressure, and preconditioning matrix elements often contain density).  <br>
> <br>
> 2.For some problems the preconditioning matrix element are constant, if it possible that the preconditioning matrix is constructed and decomposed only once at the beginning of the program and used for subsequent calculation? <br>
<br>
    Yes. PETSc only rebuilds the preconditioner after YOU change entries in the matrix with MatSetValues() etc. it knows if the matrix has not been changed (for example if you solve a linear problem with SNES) and doesn't rebuild the preconditioner.<br>
<br>
> <br>
> 3. What is the difference between using Matrix-Free method and using normal Newton method for preconditioning matrix?  <br>
<br>
   With Matrix-free you do not provide the elements to the matrix with MatSetValues() etc, with "normal" Newton's method you do. Really the only difference between the two is how the linear system is solved; Newton's method is the same in both cases.<br>
<br>
> As far as I know, when using the Matrix-Free method, the type of preconditioning matrix needs to be defined as MATMFFD. Are there any other differences? <br>
> <br>
> Thanks,<br>
> Yingjie<br>
<br>
</blockquote></div>