I suggest looking at the examples. For instance, this is the FD Laplacian:<br>
<br>
<a href="http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/ksp/ksp/examples/tutorials/ex2.c.html">http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/ksp/ksp/examples/tutorials/ex2.c.html
</a><br>
<br>
If you have a logically Cartesian grid, I advise you to use the DA object:<br>
<br>
<a href="http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/snes/examples/tutorials/ex5.c.html">http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/snes/examples/tutorials/ex5.c.html
</a><br>
<br>
Matt<br><br><div><span class="gmail_quote">On 5/22/06, <b class="gmail_sendername"><a href="mailto:mep4gk01@ucy.ac.cy">mep4gk01@ucy.ac.cy</a></b> <<a href="mailto:mep4gk01@ucy.ac.cy">mep4gk01@ucy.ac.cy</a>> wrote:
</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear Sirs,<br><br>My
name is George Katsambas and I am an MSc student at the University of
Cyprus (mechanical engineering department). I am trying to parallelize
a Laplacian sovler using BiCGStab method (in Fortran). After a lot of
effort I create, using MPI, the system Ax=b but it was very difficult
for me, to parallelize the serial solver so I download and installed
PETSc. I read the manual and I run the tutorials but I
cannot understand how I can pass my data into the PETSc routines in
order to prepare the data for the KSP solver. Can you please give some
quitelines how I can do this and if it is possible to be done. Your
help will be very valuable for me.<br><br>Best regards,<br>George Katsambas<br><br></blockquote></div><br><br clear="all"><br>-- <br>"Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness