<div dir="ltr"><div>Hi, </div><div><br></div><div>I want to know how to configure the FAS so that it solves a problem on a coarse grid of size 4h, interpolate the solution, and then stop. </div><div><br></div><div>Here is the context: I am using Newton LS to solve a problem on square domain discretized with DMDA meshed with step size h. I have a subroutine to compute the initial guess. I want this subroutine to first do a Newton solve on a coarse grid of size 4h. Then it interpolates the solution to the main mesh. I think this is achievable by using one iteration of FAS. Below is the gist of what my subroutine looks like: </div><div><br></div><div><font face="monospace">PetscErrorCode InitialState(DM da, Vec u) { // </font></div><div><font face="monospace">SNES snes; </font></div><font face="monospace">SNESCreate(PETSC_COMM_WORLD,&snes);  <br>SNESSetDM(snes,da);<br>SNESSetType(snes,SNESFAS); // solve with multigrid</font><div><font face="monospace">SNESSetTolerances(snes,PETSC_DEFAULT,PETSC_DEFAULT,PETSC_DEFAULT,1,PETSC_DEFAULT); // just one iteration<br>VecSet(u,0.0);            // start with zeros <br>SNESSolve(snes,NULL,u);   // cheap solve <br></font><div><font face="monospace">SNESGetSolution(snes,&u); // extract solution </font></div><div><font face="monospace">}</font></div></div><div><font face="monospace"><br></font></div><div><font face="arial, sans-serif">For some reason, my initial guess is too accurate. The initial guess produced by this subroutine looks exactly like the final solution. My guess is that it's doing more than what I want it to do. I'm still new to FAS. How can I tell FAS to do just one crude solve and an interpolation?</font></div></div>