[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Laplacian solver




Please do NOT look at this example. It is too simple and doesn't make clear at all how to properly parallelize the data for the matrix.

   Barry


On Mon, 22 May 2006, Matthew Knepley wrote:

I suggest looking at the examples. For instance, this is the FD Laplacian:


http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/ksp/ksp/examples/tutorials/ex2.c.html

If you have a logically Cartesian grid, I advise you to use the DA object:


http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/snes/examples/tutorials/ex5.c.html

  Matt

On 5/22/06, mep4gk01@xxxxxxxxx <mep4gk01@xxxxxxxxx> wrote:

Dear Sirs,

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.


Best regards,
George Katsambas