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

Re: Laplacian solver



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




--
"Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness