[petsc-users] petscprint problem

2017-02-09 Thread Sharp Stone
Hi, I use "PetscPrintf" and "PetscFPrintf" to make outputs in my code. In Linux, it works fine. But now I run the code in Cray cluster, I found every process will make a output, which make me go crazy and confused. Any ideas? Thanks! -- Best regards, Feng

Re: [petsc-users] ksp solve error with nested matrix

2017-02-09 Thread Lawrence Mitchell
On 09/02/17 17:29, Barry Smith wrote: > >> On Feb 8, 2017, at 11:39 PM, Dave May wrote: >> >> Any time you modify one of the submats, you need to call assembly begin/end >> on that sub matrix AND on the outer matnest object. > >Weird, and prone to errors it seems

Re: [petsc-users] ksp solve error with nested matrix

2017-02-09 Thread Barry Smith
> On Feb 8, 2017, at 11:39 PM, Dave May wrote: > > Any time you modify one of the submats, you need to call assembly begin/end > on that sub matrix AND on the outer matnest object. Weird, and prone to errors it seems to me. Perhaps this needs to be rethought > >

Re: [petsc-users] projection methods in TS

2017-02-09 Thread Gideon Simpson
Ok, thanks for the information. I am looking at a Hamiltonian system, and I wanted to try out different schemes for preserving the energy. I was curious about textbook projected RK, so I think I’ll try Barry’s suggestion about TSSetPostStep. -gideon > On Feb 4, 2017, at 10:26 PM, Emil