On Fri, Dec 9, 2022 at 10:51 AM Hongrui Yu wrote:
> Thank you for your reply! Unfortunately yes.. I’ll need to modify
> stiffness between nodes on the boundary so most of them are going to be in
> completely general location.
>
Hmm, there is usually a better way to do this. This is a mesh and
Thank you for your reply! Unfortunately yes.. I’ll need to modify stiffness
between nodes on the boundary so most of them are going to be in completely
general location.
I can create an IS after distribution using DMPlexCreatePointNumbering() but
they are Global numbering. Is there a way to
On Thu, Dec 8, 2022 at 6:06 PM Hongrui Yu wrote:
> Hello! I’m trying to adapt a serial Finite Element code using PETSc. In
> this code it reads in special stiffness terms between the boundary DoFs
> from an input file, and add them to corresponding locations in the global
> Jacobian matrix.
>
Hello! I'm trying to adapt a serial Finite Element code using PETSc. In this
code it reads in special stiffness terms between the boundary DoFs from an
input file, and add them to corresponding locations in the global Jacobian
matrix. I currently use a DM Plex object to store the mesh information.