Call MatSetOption(mat,MAT_KEEP_NONZERO_PATTERN,PETSC_TRUE);
Before calling MatZeroRows() the first time, insert the periodic condition
extra column value into the matrix for each periodic row so the matrix "has
room" for that extra value.
With this change, there should never be any mallocs, an
> On Jul 14, 2024, at 6:04 PM, Alex Sp wrote:
>
> This Message Is From an External Sender
> This message came from outside your organization.
> Hello,
>
> I am trying to solve the problem of the 3D fluid flow around a cylinder using
> an in-house finite element Fortran code. I want to apply p
Hello,
I am trying to solve the problem of the 3D fluid flow around a cylinder
using an in-house finite element Fortran code. I want to apply periodic
boundary conditions in the neutral direction (the z direction which
coincides with the axis of the cylinder). The method works fine and it
gives co