
<html>

  <head>


    <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">

  </head>

  <body text="#000000" bgcolor="#FFFFFF">

    <font face="Ubuntu">Hi,<br>

      <br>

      I am solving a variable coefficients Poisson equation with

      periodic BCs.<br>

      The equation is discretized by using the standard 5-points stencil

      finite differencing scheme.<br>

      I managed to solve the system successfully</font><font

      face="Ubuntu"> with the PCG method and now I would like to add<br>

      a preconditioner to speed up the calculation. My idea is to use

      the multigrid preconditioner.<br>

      <br>

      Example ex22f.F implements what I think I need. <br>

      If I understand correctly example ex22f.F, the subroutines

      "ComputeRHS" and "ComputeMatrix" define how the <br>

      matrix and rhs-vector have to be computed at each level.<br>

      In my case tough, both the jacobian and the rhs-vector cannot be

      computed "analytically", that is, they depend on variables<br>

      whose values are available only at the finest grid.<br>

      <br>

      How can I overcome this difficulty? <br>

      <br>

      Thank you,<br>

      <br>

      Michele<br>

      <br>

      <br>

    </font>

  </body>

</html>