
<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 need some help to use the PETSc library inside my Fortran 95

      code.<br>

      My goal is to solve a symmetric linear system that derives from

      the finite difference discretization<br>

      of the Poisson equation. I will use the preconditioned conjugate

      method.<br>

      <br>

      I am mostly interested in how to decompose my data among the

      different processes. <br>

      In particular:<br>

      <br>

      1) Since my code already implements the 2D-decomposition, would it

      be best to build the matrix with the  DMDA object type, DA object

      type<br>

           or the regular Mat type? <br>

      <br>

      2) After inserting values into a vector/matrix, PETSc performs any

      needed message passing of nonlocal components, thus the values

      locally contained on a process  may be communicated to another

      process.  How can I revert this at the end of the computation,

      that is, how can I be sure that the local solution vector contains

      the values associated to the grid nodes contained into the hosting

      process?<br>

      <br>

            <br>

      Thank you,<br>

      <br>

      Michele</font>

  </body>

</html>