Re: [petsc-users] MatSetNullSpace results in nan with MUMPS

2022-04-04 Thread Bhargav Subramanya
Dear Matt, Thank you for the reply. On Sun, Apr 3, 2022 at 4:00 PM Matthew Knepley wrote: > On Sat, Apr 2, 2022 at 8:59 PM Bhargav Subramanya < > bhargav.subrama...@kaust.edu.sa> wrote: > >> Dear All, >> >> I am trying to solve Ax = b in parallel using MUMPS, w

[petsc-users] MatSetNullSpace results in nan with MUMPS

2022-04-02 Thread Bhargav Subramanya
Dear All, I am trying to solve Ax = b in parallel using MUMPS, where x is composed of velocity, pressure, and temperature. There is a null space due to the homogeneous Neumann pressure boundary conditions. 1. Without explicitly removing the null space, the solution is bounded and robust, and I

Re: [petsc-users] Reuse MUMPS factorization

2022-03-13 Thread Bhargav Subramanya
Thanks a lot, Matt. On Sat, Mar 12, 2022 at 9:50 PM Matthew Knepley wrote: > On Sat, Mar 12, 2022 at 1:01 PM Bhargav Subramanya < > bhargav.subrama...@kaust.edu.sa> wrote: > >> Dear All, >> >> I have the following two queries: >> >> 1. I am

[petsc-users] Reuse MUMPS factorization

2022-03-12 Thread Bhargav Subramanya
Dear All, I have the following two queries: 1. I am running simulations using MUMPS through a job submission system. Since the job run time is limited, I need to restart the simulations periodically. Since I am restarting the simulation with exactly the same parameters, for instance, number of

Re: [petsc-users] dmplex extruded layers

2021-11-24 Thread Bhargav Subramanya
On Wed, Nov 24, 2021 at 8:59 PM Matthew Knepley wrote: > On Wed, Nov 24, 2021 at 12:27 PM Bhargav Subramanya < > bhargav.subrama...@kaust.edu.sa> wrote: > >> Dear Matt and Mark, >> >> Thank you very much for your reply. Your inputs are very useful to me. >

[petsc-users] dmplex extruded layers

2021-11-22 Thread Bhargav Subramanya
Dear All, I have a prismatic mesh generated by extruding a base icosahedron spherical surface mesh. The mesh is distributed on 2 processes after extrusion. The dof for my problem are defined only on the simplices of the spherical layers. 1. For a given spherical layer, I want to get points,

Re: [petsc-users] dmplex extrude in parallel

2021-11-11 Thread Bhargav Subramanya
if (dmDist) { DMDestroy(); dm = dmDist; } } ierr = DMViewFromOptions(dm, NULL, "-dm_view");CHKERRQ(ierr); --- On Thu, Nov 11, 2021 at 9:41 PM Matthew Knepley wrote: > On Thu, Nov 11, 2021 at 1:17 PM Bhargav Subramanya < > bhargav.s

Re: [petsc-users] dmplex extrude in parallel

2021-11-11 Thread Bhargav Subramanya
10, 2021 at 8:26 AM Bhargav Subramanya < > bhargav.subrama...@kaust.edu.sa> wrote: > >> Dear Matt, >> >> Thanks a lot for the reply. I am now able to generate the prismatic mesh >> properly. >> > > Cool. > > >> In the case of mpiexec -n 2 ./me

Re: [petsc-users] dmplex extrude in parallel

2021-11-10 Thread Bhargav Subramanya
AM Bhargav Subramanya < > bhargav.subrama...@kaust.edu.sa> wrote: > >> Dear All, >> >> I want to generate a prismatic mesh from a base spherical surface mesh. I >> first refine the base icosahedron sequential mesh to some extent, >> distribute the mesh,