Re: [petsc-users] PETSc initialization error

2020-06-19 Thread Sam Guo
To clarify, I call PETSc initialize and PETSc finalize everytime I call SLEPc: PetscInitializeNoPointers(argc,args,nullptr,nullptr); SlepcInitialize(,,static_cast(nullptr),help); //calling slepc SlepcFinalize(); PetscFinalize(); On Fri, Jun 19, 2020 at 10:32 PM Sam Guo wrote:

[petsc-users] PETSc initialization error

2020-06-19 Thread Sam Guo
Dear PETSc team, When I called SLEPc multiple time, I eventually got following error: MPI operation not supported by PETSc's sequential MPI wrappers [0]PETSC ERROR: #1 PetscInitialize() line 967 in ../../../petsc/src/sys/objects/pinit.c [0]PETSC ERROR: #2 SlepcInitialize() line 262 in

[petsc-users] SNES Line search: fixing the non zero solution values of dirichlet nodes during line search

2020-06-19 Thread Ashish Patel
Dear PETSc users, We use PETSc as part of a finite element method program and we are trying to properly implement Dirichlet boundary conditions for non-linear, transient problems. We find that when we use a line search method it also changes the non-zero solution value of Dirichlet nodes as it

Re: [petsc-users] Implementing periodicity using DMPlex

2020-06-19 Thread Matthew Knepley
On Thu, Jun 18, 2020 at 1:20 PM Shashwat Tiwari wrote: > Thank you for the help. I followed your suggestion and removed the > "DMPlexConstructGhostCells" function and now I am using a gmsh generated > mesh with periodicity in both directions, and distributing the mesh with 1 > level of overlap.

Re: [petsc-users] Cohesive Element Support

2020-06-19 Thread Jacob Faibussowitsch
Hello, Thank you both for your comprehensive replies. Matt: Thanks for the rundown, I will take a look at PyLith and its cohesive element impls. Blaise: The project is still in very early stages as you can tell, so everything is up in the air. Currently we are weighing using a cohesive

Re: [petsc-users] MatFindOffBlockDiagonalEntries error

2020-06-19 Thread Barry Smith
Looks like a completely different location now. Are you calling MatAXPY ? If so perhaps the arguments are wrong. One is likely a parallel matrix and one sequential. > On Jun 19, 2020, at 11:46 AM, Eda Oktay wrote: > > Dear Mark, > > I updated Petsc and again, here's the error: > >

Re: [petsc-users] MatFindOffBlockDiagonalEntries error

2020-06-19 Thread Mark Adams
This is more useful and you are getting a new error. Maybe the old error is fixed. It looks like you are trying to copy a matrix into another but they have different communicators. For example, one is created with MPI_COMM_WOLD and the other with PETSC_COMM_WORLD. On Fri, Jun 19, 2020 at 12:48

Re: [petsc-users] MatFindOffBlockDiagonalEntries error

2020-06-19 Thread Junchao Zhang
Do you have a small test example? --Junchao Zhang On Fri, Jun 19, 2020 at 11:47 AM Eda Oktay wrote: > Dear Mark, > > I updated Petsc and again, here's the error: > > [0]PETSC ERROR: - Error Message > -- > [0]PETSC

Re: [petsc-users] MatFindOffBlockDiagonalEntries error

2020-06-19 Thread Eda Oktay
Dear Mark, I updated Petsc and again, here's the error: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Arguments must have same communicators [0]PETSC ERROR: [1]PETSC ERROR: - Error Message

Re: [petsc-users] MatFindOffBlockDiagonalEntries error

2020-06-19 Thread Mark Adams
I don't know what is going on here. There was an update to this function about a year ago, so that might fix your problem. We would need you to test with a current version. Mark On Fri, Jun 19, 2020 at 11:23 AM Eda Oktay wrote: > Hi all, > > I am trying to find off block diagonal entries of a

[petsc-users] MatFindOffBlockDiagonalEntries error

2020-06-19 Thread Eda Oktay
Hi all, I am trying to find off block diagonal entries of a matrix and I am trying to use MatFindOffBlockDiagonalEntries. However, although my matrix is not NULL, I am getting an error message as follows: 1]PETSC ERROR: - Error Message

Re: [petsc-users] Cohesive Element Support

2020-06-19 Thread Blaise A Bourdin
On Jun 18, 2020, at 5:28 AM, Matthew Knepley mailto:knep...@gmail.com>> wrote: On Wed, Jun 17, 2020 at 4:05 PM Jacob Faibussowitsch mailto:jacob@gmail.com>> wrote: Hello, I am looking to perform large scale fracture and crack propagation simulations and have a few questions regarding