Re: [petsc-users] MPI Iterative solver crash on HPC

2019-01-17 Thread Zhang, Junchao via petsc-users
Usually when I meet a SEGV error, I will run it again with a parallel debugger like DDT and wait for it to segfault, and then examine the stack trace to see what is wrong. Can you compress your matrix and upload it to google drive, so we can try to reproduce the error. --Junchao Zhang On Thu,

Re: [petsc-users] I suggest to create a demo code for the problem of "Get a value from global vector to a local processor"

2019-01-17 Thread leejearl via petsc-users
Hi Matt,Thanks again for your helps. You have help me a lot. Petsc is a great library. I need to learn more about it, thanks. LeejearlOn Thu, 2019-01-17 at 21:41 -0500, Matthew Knepley wrote: > On Thu, Jan 17, 2019 at 9:39 PM leejearl > wrote: > > Hi, Matt, > > On Thu, 2019-01-17 at

Re: [petsc-users] I suggest to create a demo code for the problem of "Get a value from global vector to a local processor"

2019-01-17 Thread leejearl via petsc-users
Hi, Matt, On Thu, 2019-01-17 at 21:20 -0500, Matthew Knepley wrote: > On Thu, Jan 17, 2019 at 9:04 PM leejearl > wrote: > > Hi Matt, > > Thanks for your reply. We have got to a consensus. For such a > > problem, I have sought the mail list for helps and found a great > > many misleading

Re: [petsc-users] I suggest to create a demo code for the problem of "Get a value from global vector to a local processor"

2019-01-17 Thread leejearl via petsc-users
Hi Matt,Thanks for your reply. We have got to a consensus. For such a problem, I have sought the mail list for helps and found a great many misleading answers. I think some of answers may work for the DMDA. Now I have overcome this by creating a PetscSF and PetscSFBcast, and the

Re: [petsc-users] VecCopy fails after VecDuplicate

2019-01-17 Thread Sajid Ali via petsc-users
Nevermind, it was my fault in thinking the error was with u_abs and not u. I switched from local array based value setting for initial conditions to VecSetValues when converting the uniprocessor example to an MPI program. While I removed VecRestoreArray and swapped u_local[*ptr] assignments with

Re: [petsc-users] MPI Iterative solver crash on HPC

2019-01-17 Thread Smith, Barry F. via petsc-users
0]PETSC ERROR: Caught signal number 15 Terminate: Some process (or the batch system) has told this process to end This happens because you have run out of time requested on your batch system so the batch system is terminating the processes or you are running low on memory (or some other

Re: [petsc-users] VecCopy fails after VecDuplicate

2019-01-17 Thread Sajid Ali via petsc-users
As requested : [sajid@xrm free_space]$ ./ex_modify Solving a linear TS problem on 1 processor m : 256, slices : 1000.00, lambda : 1.239800e-10 [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Object is in wrong

Re: [petsc-users] VecCopy fails after VecDuplicate

2019-01-17 Thread Smith, Barry F. via petsc-users
It is likely the u vector that is not assembled not the u_abs; the complete error message has this information. If the monitor is called from TS then I cannot imagine how the u could possibly be unassembled. Please send the complete error message. Barry > On Jan 17, 2019, at 1:18 PM,

[petsc-users] VecCopy fails after VecDuplicate

2019-01-17 Thread Sajid Ali via petsc-users
Hi, I have the following 2 lines in a function in my code : ierr = VecDuplicate(u,_abs);CHKERRQ(ierr); ierr = VecCopy(u,u_abs);CHKERRQ(ierr); The VecCopy fails with the error message : [0]PETSC ERROR: Object is in wrong state [0]PETSC ERROR: Not for unassembled vector adding the following

[petsc-users] I suggest to create a demo code for the problem of "Get a value from global vector to a local processor"

2019-01-17 Thread leejearl via petsc-users
Hi all Petscer, I have ask helps for some questions. Thanks for the replies from developer. I have known more about the petsc. I have also sought helps in the mail lists, and I find that there are many subjects focused on such a problem. Some subjects are listed as follow