Laplacian solver

Barry Smith bsmith at mcs.anl.gov
Mon May 22 17:11:36 CDT 2006


   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 at ucy.ac.cy <mep4gk01 at ucy.ac.cy> 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
>> 
>> 
>
>
>




More information about the petsc-users mailing list