Re: [petsc-users] How to choose mat_mffd_err in JFNK

Thanks, it works. snes_mf_jorge works for me. It appears to compute h in every ksp. Without -snes_mf_jorge, it is not working. For some reason, it only computes h once, but that h is bad. My gmres residual is not decaying. Indeed, the noise in my function becomes larger when I refine the mesh.

Re: [petsc-users] How to choose mat_mffd_err in JFNK

> On Aug 14, 2019, at 9:41 PM, Tang, Qi wrote: > > Thanks for the help, Barry. I tired both ds and wp, and again it depends on > if I could find the correct parameter set. It is getting harder as I refine > the mesh. > > So I try to use SNESDefaultMatrixFreeCreate2, SNESMatrixFreeMult2_Priva

Re: [petsc-users] How to choose mat_mffd_err in JFNK

Thanks for the help, Barry. I tired both ds and wp, and again it depends on if I could find the correct parameter set. It is getting harder as I refine the mesh. So I try to use SNESDefaultMatrixFreeCreate2, SNESMatrixFreeMult2_Private or SNESDiffParameterCompute_More in mfem. But it looks like

Re: [petsc-users] Creating a 3D dmplex mesh with cell list and distributing it

On Wed, Aug 14, 2019 at 10:23 PM Swarnava Ghosh wrote: > Hi Matthew, > > I added DMView(pCgdft->dmplex,PETSC_VIEWER_STDOUT_WORLD); before and after > distribution, and I get the following: > It looks like you are running things with the wrong 'mpirun' Thanks, Matt > dmplex before distr