[petsc-users] Looking for PETSc Example: Poisson Equation in 2D using FEM - scale upto 1000 procs

[Matthew Knepley]

On Mon, Sep 18, 2023 at 4:16 PM Khurana, Parv <p.khurana22 at imperial.ac.uk>
wrote:

> Hello PETSc users,
>
>
>
> Thank you for this very active community of users and the mailing list.
>
>
>
> I am looking for a PETSc example which solves the Poisson equation on a 2D
> domain using FEM (or HO-FEM if possible). I would like the following:
>
>    1. The example should be formulated fully on PETSc, and to be solved
>    with KSP objects in PETSc.
>    2. The problem should scale up to a few hundred processors (ideally
>    1000 procs).
>    3. Ideally on a unit square with either Square or Triangle element
>    discretization.
>    4. There should be an option to specify Dirichlet/Neumann-type BCs on
>    the boundaries.
>
>
>
> I was wondering if someone could point me to such an example as I am
> relatively new to PETSc – and I am trying to avoid reinventing the wheel. I
> have had a look at Examples 29,32,50 and 66 in the PETSc tutorials – while
> they are very close to what I need I am not sure if they scale to a few
> hundred processors. Furthermore, I am aware that I can formulate such a
> problem with FENICS/Firedrake with relative ease and these software
> interfaces with PETSc quite well. However, I am just trying to see if such
> an application already exists purely in PETSc.
>
>
>
It looks like Mark has everything working. However, if you have problems
getting that to run, I am here at ICL this week.

  Thanks,

     Matt


> Thanks, and Best
>
> Parv
>
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230918/b5963eec/attachment.html>