August 2012 - petsc-dev - mailman3.cels.anl.gov

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Peter Brune 14 Aug '12

14 Aug '12

On Tue, Aug 14, 2012 at 4:47 PM, Jed Brown <jedbrown(a)mcs.anl.gov> wrote: > On Tue, Aug 14, 2012 at 3:44 PM, Matthew Knepley <knepley(a)gmail.com>wrote: > >> You can, but the setter routines interpret NULL as "ignore", i.e. >>> >>> if (func) ctx->func = func; >>> >> >> Does this make sense to anyone? Can we remove that if? > > > Then we would have to revisit lots of locations and get the old value > before resetting it. Passing NULL was a convenience. > Correct; the sheer amount of code using this convention makes it exceedingly painful to fix. It would be much easier to alter our internals. I'm *this* close to fugueing on TS especially and removing all PetscObjects from ctx arguments internal to PETSc. We can pass the TS around by PetscObjectCompose using the coarsen/restrict hooks. - Peter

1 0

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Matthew Knepley 14 Aug '12

14 Aug '12

On Tue, Aug 14, 2012 at 4:47 PM, Jed Brown <jedbrown(a)mcs.anl.gov> wrote: > On Tue, Aug 14, 2012 at 3:44 PM, Matthew Knepley <knepley(a)gmail.com>wrote: > >> You can, but the setter routines interpret NULL as "ignore", i.e. >>> >>> if (func) ctx->func = func; >>> >> >> Does this make sense to anyone? Can we remove that if? > > > Then we would have to revisit lots of locations and get the old value > before resetting it. Passing NULL was a convenience. > Ah, I see. This is crappy design. Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

1 0

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Jed Brown 14 Aug '12

14 Aug '12

On Tue, Aug 14, 2012 at 3:44 PM, Matthew Knepley <knepley(a)gmail.com> wrote: > You can, but the setter routines interpret NULL as "ignore", i.e. >> >> if (func) ctx->func = func; >> > > Does this make sense to anyone? Can we remove that if? Then we would have to revisit lots of locations and get the old value before resetting it. Passing NULL was a convenience.

1 0

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Matthew Knepley 14 Aug '12

14 Aug '12

On Tue, Aug 14, 2012 at 4:40 PM, Jed Brown <jedbrown(a)mcs.anl.gov> wrote: > On Tue, Aug 14, 2012 at 3:35 PM, Matthew Knepley <knepley(a)gmail.com>wrote: > >> I guess I do not understand. What is the context set API? Why can't you >> pass NULL? > > > You can, but the setter routines interpret NULL as "ignore", i.e. > > if (func) ctx->func = func; > Does this make sense to anyone? Can we remove that if? Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

1 0

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Peter Brune 14 Aug '12

14 Aug '12

You have one call to set the TS and SNES jacobian nested as: TSSetRHSJacobian(ts ..); SNESSetJacobian(snes, .. SNESTSComputeJacobian, (void *)ts); And then a subsequent and seemingly unrelated call to SNESSetJacobian(snes, .. .. SNESDefaultComputeJacobianColor, PETSC_NULL); (Where the last argument is the context in both calls, of course). SNESDefaultComputeJacobianColor takes an optional context argument as a user-defined coloring. This allows people to set their own coloring if need be; otherwise it automatically pulls one off the DM or Mat. However, in this case, we have that there is no possibility of removing that TS from the context spot, as PETSC_NULL is ignored. - Peter On Tue, Aug 14, 2012 at 4:35 PM, Matthew Knepley <knepley(a)gmail.com> wrote: > On Tue, Aug 14, 2012 at 4:28 PM, Peter Brune <prbrune(a)gmail.com> wrote: > >> This is a problem I've been harping on but should harp on more; It is >> IMPOSSIBLE to reset the context on a lot of things to null after it's been >> set. What's happening is that the context is being set as the TS object in >> the previous setting of the SNES jacobian in TSSetRHSJacobian. Then, this >> call is expecting it to have a MatFDColoring. I could just make >> SNESDefaultComputeJacobianColor ignore its context argument, but I left it >> in there so that a user could pass their own coloring in. >> > > I guess I do not understand. What is the context set API? Why can't you > pass NULL? > > Matt > > >> >> It's my opinion that sticking PetscObjects in the context field like the >> SNES wrappers for TS function and jacobian evaluations do is a bad habit. >> However, just composing it with, say, the DM won't give us the flexibility >> we need if, for instance, we want multilevel methods to work on this. >> >> - Peter >> >> On Tue, Aug 14, 2012 at 4:12 PM, Barry Smith <bsmith(a)mcs.anl.gov> wrote: >> >>> >>> Things are still mighty fragile ---- >>> >>> > [0]PETSC ERROR: --------------------- Error Message >>> ------------------------------------ >>> > [0]PETSC ERROR: Invalid argument! >>> > [0]PETSC ERROR: Wrong type of object: Parameter # 6! >>> > [0]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [0]PETSC ERROR: Petsc Development HG revision: >>> bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 >>> -0500 >>> > [0]PETSC ERROR: See docs/changes/index.html for recent updates. >>> > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. >>> > [0]PETSC ERROR: See docs/index.html for manual pages. >>> > [1]PETSC ERROR: --------------------- Error Message >>> ------------------------------------ >>> > [1]PETSC ERROR: Invalid argument! >>> > [1]PETSC ERROR: Wrong type of object: Parameter # 6! >>> > [1]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [1]PETSC ERROR: Petsc Development HG revision: >>> bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 >>> -0500 >>> > [1]PETSC ERROR: See docs/changes/index.html for recent updates. >>> > [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting. >>> > [1]PETSC ERROR: See docs/index.html for manual pages. >>> > [1]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [1]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by >>> barrysmith Mon Aug 13 21:36:58 2012 >>> > [2]PETSC ERROR: --------------------- Error Message >>> ------------------------------------ >>> > [2]PETSC ERROR: Invalid argument! >>> > [2]PETSC ERROR: Wrong type of object: Parameter # 6! >>> > [2]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [2]PETSC ERROR: Petsc Development HG revision: >>> bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 >>> -0500 >>> > [2]PETSC ERROR: See docs/changes/index.html for recent updates. >>> > [2]PETSC ERROR: See docs/faq.html for hints about trouble shooting. >>> > [2]PETSC ERROR: See docs/index.html for manual pages. >>> > [2]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [2]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by >>> barrysmith Mon Aug 13 21:36:58 2012 >>> > [2]PETSC ERROR: Libraries linked from >>> /Users/barrysmith/Src/petsc-dev/arch-gnu/lib >>> > [2]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 >>> > [2]PETSC ERROR: Configure options --download-blacs --download-fftw >>> --download-hypre --download-libyaml --download-metis --download-ml >>> --download-mpich --download-mumps --download-parmetis --download-ptscotch >>> --download-scalapack --download-superlu --download-superlu_dist >>> --download-triangle --download-yaml >>> --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp >>> --with-server --with-shared-libraries PETSC_ARCH=arch-gnu >>> > [2]PETSC ERROR: [0]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [0]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by >>> barrysmith Mon Aug 13 21:36:58 2012 >>> > [0]PETSC ERROR: Libraries linked from >>> /Users/barrysmith/Src/petsc-dev/arch-gnu/lib >>> > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 >>> > [0]PETSC ERROR: Configure options --download-blacs --download-fftw >>> --download-hypre --download-libyaml --download-metis --download-ml >>> --download-mpich --download-mumps --download-parmetis --download-ptscotch >>> --download-scalapack --download-superlu --download-superlu_dist >>> --download-triangle --download-yaml >>> --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp >>> --with-server --with-shared-libraries PETSC_ARCH=arch-gnu >>> > [0]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [0]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c >>> > [0]PETSC ERROR: SNESComputeJacobian() line 2079 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >>> > [0]PETSC ERROR: SNESSolve_LS() line 198 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c >>> > [0]PETSC ERROR: SNESSolve() line 3570 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >>> > [1]PETSC ERROR: Libraries linked from >>> /Users/barrysmith/Src/petsc-dev/arch-gnu/lib >>> > [1]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 >>> > [1]PETSC ERROR: Configure options --download-blacs --download-fftw >>> --download-hypre --download-libyaml --download-metis --download-ml >>> --download-mpich --download-mumps --download-parmetis --download-ptscotch >>> --download-scalapack --download-superlu --download-superlu_dist >>> --download-triangle --download-yaml >>> --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp >>> --with-server --with-shared-libraries PETSC_ARCH=arch-gnu >>> > [1]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [1]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c >>> > [1]PETSC ERROR: SNESComputeJacobian() line 2079 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >>> > [1]PETSC ERROR: SNESSolve_LS() line 198 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c >>> > [1]PETSC ERROR: SNESSolve() line 3570 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >>> > [1]PETSC ERROR: TSStep_Theta() line 116 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c >>> > [1]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> > [2]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c >>> > [2]PETSC ERROR: SNESComputeJacobian() line 2079 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >>> > [2]PETSC ERROR: SNESSolve_LS() line 198 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c >>> > [2]PETSC ERROR: SNESSolve() line 3570 in >>> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >>> > [2]PETSC ERROR: TSStep_Theta() line 116 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c >>> > [2]PETSC ERROR: TSStep() line 1977 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >>> > [2]PETSC ERROR: TSSolve() line 2100 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >>> > [0]PETSC ERROR: TSStep_Theta() line 116 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c >>> > [0]PETSC ERROR: TSStep() line 1977 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >>> > [0]PETSC ERROR: TSSolve() line 2100 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >>> > [0]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c >>> > TSStep() line 1977 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >>> > [1]PETSC ERROR: TSSolve() line 2100 in >>> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >>> > [1]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c >>> > [2]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c >>> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 >>> > [cli_0]: aborting job: >>> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 >>> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 >>> > [cli_1]: aborting job: >>> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 >>> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 >>> > [cli_2]: aborting job: >>> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 >>> /Users/barrysmith/Src/petsc-dev/src/ts/examples/tutorials >>> Possible problem with ex13_3, diffs above >>> ========================================= >>> >>> >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener >

1 0

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Jed Brown 14 Aug '12

14 Aug '12

On Tue, Aug 14, 2012 at 3:35 PM, Matthew Knepley <knepley(a)gmail.com> wrote: > I guess I do not understand. What is the context set API? Why can't you > pass NULL? You can, but the setter routines interpret NULL as "ignore", i.e. if (func) ctx->func = func;

1 0

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Matthew Knepley 14 Aug '12

14 Aug '12

On Tue, Aug 14, 2012 at 4:28 PM, Peter Brune <prbrune(a)gmail.com> wrote: > This is a problem I've been harping on but should harp on more; It is > IMPOSSIBLE to reset the context on a lot of things to null after it's been > set. What's happening is that the context is being set as the TS object in > the previous setting of the SNES jacobian in TSSetRHSJacobian. Then, this > call is expecting it to have a MatFDColoring. I could just make > SNESDefaultComputeJacobianColor ignore its context argument, but I left it > in there so that a user could pass their own coloring in. > I guess I do not understand. What is the context set API? Why can't you pass NULL? Matt > > It's my opinion that sticking PetscObjects in the context field like the > SNES wrappers for TS function and jacobian evaluations do is a bad habit. > However, just composing it with, say, the DM won't give us the flexibility > we need if, for instance, we want multilevel methods to work on this. > > - Peter > > On Tue, Aug 14, 2012 at 4:12 PM, Barry Smith <bsmith(a)mcs.anl.gov> wrote: > >> >> Things are still mighty fragile ---- >> >> > [0]PETSC ERROR: --------------------- Error Message >> ------------------------------------ >> > [0]PETSC ERROR: Invalid argument! >> > [0]PETSC ERROR: Wrong type of object: Parameter # 6! >> > [0]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [0]PETSC ERROR: Petsc Development HG revision: >> bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 >> -0500 >> > [0]PETSC ERROR: See docs/changes/index.html for recent updates. >> > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. >> > [0]PETSC ERROR: See docs/index.html for manual pages. >> > [1]PETSC ERROR: --------------------- Error Message >> ------------------------------------ >> > [1]PETSC ERROR: Invalid argument! >> > [1]PETSC ERROR: Wrong type of object: Parameter # 6! >> > [1]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [1]PETSC ERROR: Petsc Development HG revision: >> bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 >> -0500 >> > [1]PETSC ERROR: See docs/changes/index.html for recent updates. >> > [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting. >> > [1]PETSC ERROR: See docs/index.html for manual pages. >> > [1]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [1]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by >> barrysmith Mon Aug 13 21:36:58 2012 >> > [2]PETSC ERROR: --------------------- Error Message >> ------------------------------------ >> > [2]PETSC ERROR: Invalid argument! >> > [2]PETSC ERROR: Wrong type of object: Parameter # 6! >> > [2]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [2]PETSC ERROR: Petsc Development HG revision: >> bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 >> -0500 >> > [2]PETSC ERROR: See docs/changes/index.html for recent updates. >> > [2]PETSC ERROR: See docs/faq.html for hints about trouble shooting. >> > [2]PETSC ERROR: See docs/index.html for manual pages. >> > [2]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [2]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by >> barrysmith Mon Aug 13 21:36:58 2012 >> > [2]PETSC ERROR: Libraries linked from >> /Users/barrysmith/Src/petsc-dev/arch-gnu/lib >> > [2]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 >> > [2]PETSC ERROR: Configure options --download-blacs --download-fftw >> --download-hypre --download-libyaml --download-metis --download-ml >> --download-mpich --download-mumps --download-parmetis --download-ptscotch >> --download-scalapack --download-superlu --download-superlu_dist >> --download-triangle --download-yaml >> --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp >> --with-server --with-shared-libraries PETSC_ARCH=arch-gnu >> > [2]PETSC ERROR: [0]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [0]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by >> barrysmith Mon Aug 13 21:36:58 2012 >> > [0]PETSC ERROR: Libraries linked from >> /Users/barrysmith/Src/petsc-dev/arch-gnu/lib >> > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 >> > [0]PETSC ERROR: Configure options --download-blacs --download-fftw >> --download-hypre --download-libyaml --download-metis --download-ml >> --download-mpich --download-mumps --download-parmetis --download-ptscotch >> --download-scalapack --download-superlu --download-superlu_dist >> --download-triangle --download-yaml >> --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp >> --with-server --with-shared-libraries PETSC_ARCH=arch-gnu >> > [0]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [0]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c >> > [0]PETSC ERROR: SNESComputeJacobian() line 2079 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >> > [0]PETSC ERROR: SNESSolve_LS() line 198 in >> /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c >> > [0]PETSC ERROR: SNESSolve() line 3570 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >> > [1]PETSC ERROR: Libraries linked from >> /Users/barrysmith/Src/petsc-dev/arch-gnu/lib >> > [1]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 >> > [1]PETSC ERROR: Configure options --download-blacs --download-fftw >> --download-hypre --download-libyaml --download-metis --download-ml >> --download-mpich --download-mumps --download-parmetis --download-ptscotch >> --download-scalapack --download-superlu --download-superlu_dist >> --download-triangle --download-yaml >> --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp >> --with-server --with-shared-libraries PETSC_ARCH=arch-gnu >> > [1]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [1]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c >> > [1]PETSC ERROR: SNESComputeJacobian() line 2079 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >> > [1]PETSC ERROR: SNESSolve_LS() line 198 in >> /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c >> > [1]PETSC ERROR: SNESSolve() line 3570 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >> > [1]PETSC ERROR: TSStep_Theta() line 116 in >> /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c >> > [1]PETSC ERROR: >> ------------------------------------------------------------------------ >> > [2]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c >> > [2]PETSC ERROR: SNESComputeJacobian() line 2079 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >> > [2]PETSC ERROR: SNESSolve_LS() line 198 in >> /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c >> > [2]PETSC ERROR: SNESSolve() line 3570 in >> /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c >> > [2]PETSC ERROR: TSStep_Theta() line 116 in >> /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c >> > [2]PETSC ERROR: TSStep() line 1977 in >> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >> > [2]PETSC ERROR: TSSolve() line 2100 in >> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >> > [0]PETSC ERROR: TSStep_Theta() line 116 in >> /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c >> > [0]PETSC ERROR: TSStep() line 1977 in >> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >> > [0]PETSC ERROR: TSSolve() line 2100 in >> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >> > [0]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c >> > TSStep() line 1977 in >> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >> > [1]PETSC ERROR: TSSolve() line 2100 in >> /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c >> > [1]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c >> > [2]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c >> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 >> > [cli_0]: aborting job: >> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 >> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 >> > [cli_1]: aborting job: >> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 >> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 >> > [cli_2]: aborting job: >> > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 >> /Users/barrysmith/Src/petsc-dev/src/ts/examples/tutorials >> Possible problem with ex13_3, diffs above >> ========================================= >> >> > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

1 0

Re: [petsc-dev] broken TS example related to Jacobian and coloring etc.
by Peter Brune 14 Aug '12

14 Aug '12

This is a problem I've been harping on but should harp on more; It is IMPOSSIBLE to reset the context on a lot of things to null after it's been set. What's happening is that the context is being set as the TS object in the previous setting of the SNES jacobian in TSSetRHSJacobian. Then, this call is expecting it to have a MatFDColoring. I could just make SNESDefaultComputeJacobianColor ignore its context argument, but I left it in there so that a user could pass their own coloring in. It's my opinion that sticking PetscObjects in the context field like the SNES wrappers for TS function and jacobian evaluations do is a bad habit. However, just composing it with, say, the DM won't give us the flexibility we need if, for instance, we want multilevel methods to work on this. - Peter On Tue, Aug 14, 2012 at 4:12 PM, Barry Smith <bsmith(a)mcs.anl.gov> wrote: > > Things are still mighty fragile ---- > > > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > > [0]PETSC ERROR: Invalid argument! > > [0]PETSC ERROR: Wrong type of object: Parameter # 6! > > [0]PETSC ERROR: > ------------------------------------------------------------------------ > > [0]PETSC ERROR: Petsc Development HG revision: > bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 > -0500 > > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [0]PETSC ERROR: See docs/index.html for manual pages. > > [1]PETSC ERROR: --------------------- Error Message > ------------------------------------ > > [1]PETSC ERROR: Invalid argument! > > [1]PETSC ERROR: Wrong type of object: Parameter # 6! > > [1]PETSC ERROR: > ------------------------------------------------------------------------ > > [1]PETSC ERROR: Petsc Development HG revision: > bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 > -0500 > > [1]PETSC ERROR: See docs/changes/index.html for recent updates. > > [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [1]PETSC ERROR: See docs/index.html for manual pages. > > [1]PETSC ERROR: > ------------------------------------------------------------------------ > > [1]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by > barrysmith Mon Aug 13 21:36:58 2012 > > [2]PETSC ERROR: --------------------- Error Message > ------------------------------------ > > [2]PETSC ERROR: Invalid argument! > > [2]PETSC ERROR: Wrong type of object: Parameter # 6! > > [2]PETSC ERROR: > ------------------------------------------------------------------------ > > [2]PETSC ERROR: Petsc Development HG revision: > bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 > -0500 > > [2]PETSC ERROR: See docs/changes/index.html for recent updates. > > [2]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [2]PETSC ERROR: See docs/index.html for manual pages. > > [2]PETSC ERROR: > ------------------------------------------------------------------------ > > [2]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by > barrysmith Mon Aug 13 21:36:58 2012 > > [2]PETSC ERROR: Libraries linked from > /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > > [2]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > > [2]PETSC ERROR: Configure options --download-blacs --download-fftw > --download-hypre --download-libyaml --download-metis --download-ml > --download-mpich --download-mumps --download-parmetis --download-ptscotch > --download-scalapack --download-superlu --download-superlu_dist > --download-triangle --download-yaml > --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp > --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > > [2]PETSC ERROR: [0]PETSC ERROR: > ------------------------------------------------------------------------ > > [0]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by > barrysmith Mon Aug 13 21:36:58 2012 > > [0]PETSC ERROR: Libraries linked from > /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > > [0]PETSC ERROR: Configure options --download-blacs --download-fftw > --download-hypre --download-libyaml --download-metis --download-ml > --download-mpich --download-mumps --download-parmetis --download-ptscotch > --download-scalapack --download-superlu --download-superlu_dist > --download-triangle --download-yaml > --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp > --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > > [0]PETSC ERROR: > ------------------------------------------------------------------------ > > [0]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c > > [0]PETSC ERROR: SNESComputeJacobian() line 2079 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESSolve_LS() line 198 in > /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c > > [0]PETSC ERROR: SNESSolve() line 3570 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [1]PETSC ERROR: Libraries linked from > /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > > [1]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > > [1]PETSC ERROR: Configure options --download-blacs --download-fftw > --download-hypre --download-libyaml --download-metis --download-ml > --download-mpich --download-mumps --download-parmetis --download-ptscotch > --download-scalapack --download-superlu --download-superlu_dist > --download-triangle --download-yaml > --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp > --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > > [1]PETSC ERROR: > ------------------------------------------------------------------------ > > [1]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c > > [1]PETSC ERROR: SNESComputeJacobian() line 2079 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [1]PETSC ERROR: SNESSolve_LS() line 198 in > /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c > > [1]PETSC ERROR: SNESSolve() line 3570 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [1]PETSC ERROR: TSStep_Theta() line 116 in > /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c > > [1]PETSC ERROR: > ------------------------------------------------------------------------ > > [2]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c > > [2]PETSC ERROR: SNESComputeJacobian() line 2079 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [2]PETSC ERROR: SNESSolve_LS() line 198 in > /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c > > [2]PETSC ERROR: SNESSolve() line 3570 in > /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [2]PETSC ERROR: TSStep_Theta() line 116 in > /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c > > [2]PETSC ERROR: TSStep() line 1977 in > /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > > [2]PETSC ERROR: TSSolve() line 2100 in > /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > > [0]PETSC ERROR: TSStep_Theta() line 116 in > /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c > > [0]PETSC ERROR: TSStep() line 1977 in > /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > > [0]PETSC ERROR: TSSolve() line 2100 in > /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > > [0]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c > > TSStep() line 1977 in > /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > > [1]PETSC ERROR: TSSolve() line 2100 in > /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > > [1]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c > > [2]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c > > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 > > [cli_0]: aborting job: > > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 > > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 > > [cli_1]: aborting job: > > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 > > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 > > [cli_2]: aborting job: > > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 > /Users/barrysmith/Src/petsc-dev/src/ts/examples/tutorials > Possible problem with ex13_3, diffs above > ========================================= > >

1 0

broken TS example related to Jacobian and coloring etc.
by Barry Smith 14 Aug '12

14 Aug '12

Things are still mighty fragile ---- > [0]PETSC ERROR: --------------------- Error Message ------------------------------------ > [0]PETSC ERROR: Invalid argument! > [0]PETSC ERROR: Wrong type of object: Parameter # 6! > [0]PETSC ERROR: ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Development HG revision: bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 -0500 > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [0]PETSC ERROR: See docs/index.html for manual pages. > [1]PETSC ERROR: --------------------- Error Message ------------------------------------ > [1]PETSC ERROR: Invalid argument! > [1]PETSC ERROR: Wrong type of object: Parameter # 6! > [1]PETSC ERROR: ------------------------------------------------------------------------ > [1]PETSC ERROR: Petsc Development HG revision: bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 -0500 > [1]PETSC ERROR: See docs/changes/index.html for recent updates. > [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [1]PETSC ERROR: See docs/index.html for manual pages. > [1]PETSC ERROR: ------------------------------------------------------------------------ > [1]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by barrysmith Mon Aug 13 21:36:58 2012 > [2]PETSC ERROR: --------------------- Error Message ------------------------------------ > [2]PETSC ERROR: Invalid argument! > [2]PETSC ERROR: Wrong type of object: Parameter # 6! > [2]PETSC ERROR: ------------------------------------------------------------------------ > [2]PETSC ERROR: Petsc Development HG revision: bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 -0500 > [2]PETSC ERROR: See docs/changes/index.html for recent updates. > [2]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [2]PETSC ERROR: See docs/index.html for manual pages. > [2]PETSC ERROR: ------------------------------------------------------------------------ > [2]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by barrysmith Mon Aug 13 21:36:58 2012 > [2]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > [2]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > [2]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-triangle --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > [2]PETSC ERROR: [0]PETSC ERROR: ------------------------------------------------------------------------ > [0]PETSC ERROR: ./ex13 on a arch-gnu named Barrys-MacBook-Pro.local by barrysmith Mon Aug 13 21:36:58 2012 > [0]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > [0]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-triangle --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > [0]PETSC ERROR: ------------------------------------------------------------------------ > [0]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c > [0]PETSC ERROR: SNESComputeJacobian() line 2079 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > [0]PETSC ERROR: SNESSolve_LS() line 198 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > [1]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > [1]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > [1]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-triangle --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > [1]PETSC ERROR: ------------------------------------------------------------------------ > [1]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c > [1]PETSC ERROR: SNESComputeJacobian() line 2079 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > [1]PETSC ERROR: SNESSolve_LS() line 198 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c > [1]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > [1]PETSC ERROR: TSStep_Theta() line 116 in /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c > [1]PETSC ERROR: ------------------------------------------------------------------------ > [2]PETSC ERROR: SNESDefaultComputeJacobianColor() line 50 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snesj2.c > [2]PETSC ERROR: SNESComputeJacobian() line 2079 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > [2]PETSC ERROR: SNESSolve_LS() line 198 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/ls/ls.c > [2]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > [2]PETSC ERROR: TSStep_Theta() line 116 in /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c > [2]PETSC ERROR: TSStep() line 1977 in /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > [2]PETSC ERROR: TSSolve() line 2100 in /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > [0]PETSC ERROR: TSStep_Theta() line 116 in /Users/barrysmith/Src/petsc-dev/src/ts/impls/implicit/theta/theta.c > [0]PETSC ERROR: TSStep() line 1977 in /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > [0]PETSC ERROR: TSSolve() line 2100 in /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > [0]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c > TSStep() line 1977 in /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > [1]PETSC ERROR: TSSolve() line 2100 in /Users/barrysmith/Src/petsc-dev/src/ts/interface/ts.c > [1]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c > [2]PETSC ERROR: main() line 108 in src/ts/examples/tutorials/ex13.c > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 > [cli_0]: aborting job: > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0 > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 > [cli_1]: aborting job: > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 1 > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 > [cli_2]: aborting job: > application called MPI_Abort(MPI_COMM_WORLD, 62) - process 2 /Users/barrysmith/Src/petsc-dev/src/ts/examples/tutorials Possible problem with ex13_3, diffs above =========================================

1 0

Re: [petsc-dev] annoying Nan in tests
by Barry Smith 14 Aug '12

14 Aug '12

Peter, Thanks On Aug 13, 2012, at 10:43 PM, Peter Brune <prbrune(a)gmail.com> wrote: > I fixed it; the prefix fix from last week for PFLOTRAN broke it for composed solvers' default options for the line search type, making the richardson smoothers run with the basic linesearch and diverge ASAP. The reason it took so long to show up in dev is due to it being pushed to 3.3 originally, so noticing it was even harder. > > - Peter > > On Mon, Aug 13, 2012 at 9:24 PM, Peter Brune <prbrune(a)gmail.com> wrote: > Something broke; I think I see what it is from your logs. I'll fix it. > > - Peter > > > On Mon, Aug 13, 2012 at 9:17 PM, Barry Smith <bsmith(a)mcs.anl.gov> wrote: > > What's up with this stuff? Are they bugs? Bad algorithms? What? It is hard to make a major change to PETSc (and use alltests to check that you haven't broken anything) when this crap comes up constantly. > > Is there a plan to deal with this or are we just screwed forever? > > > Barry > > > 1 SNES Function norm 59.1937 > > [0]PETSC ERROR: --------------------- Error Message ------------------------------------ > > [0]PETSC ERROR: Floating point exception! > > [0]PETSC ERROR: Vec entry at local location 83 is not-a-number or infinite at end of function: Parameter number 3! > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: Petsc Development HG revision: bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 -0500 > > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [0]PETSC ERROR: See docs/index.html for manual pages. > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: ./ex5 on a arch-gnu named Barrys-MacBook-Pro.local by barrysmith Mon Aug 13 21:12:18 2012 > > [0]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > > [0]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-triangle --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: SNESComputeFunction() line 1927 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESLineSearchApply_Basic() line 41 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: SNESLineSearchApply() line 506 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: SNESSolve_NRichardson() line 181 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/richardson/snesrichardson.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESSolve_NGMRES() line 259 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/ngmres/snesngmres.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: main() line 152 in src/snes/examples/tutorials/ex5.c > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > > [cli_0]: aborting job: > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > /Users/barrysmith/Src/petsc-dev/src/snes/examples/tutorials > Possible problem with with ex5_5_ngmres_richardson, diffs above > ========================================= > 2,21c2,44 > < 1 SNES Function norm 0.18208 > < 2 SNES Function norm 0.0272102 > < 3 SNES Function norm 0.00491305 > < 4 SNES Function norm 0.000996383 > < 5 SNES Function norm 0.000243372 > < 6 SNES Function norm 6.53869e-05 > < 7 SNES Function norm 1.84573e-05 > < 8 SNES Function norm 5.30014e-06 > < 9 SNES Function norm 1.5228e-06 > < 10 SNES Function norm 4.34313e-07 > < 11 SNES Function norm 1.22592e-07 > < 12 SNES Function norm 3.42365e-08 > < 13 SNES Function norm 9.46937e-09 > < 14 SNES Function norm 2.59769e-09 > < 15 SNES Function norm 7.078e-10 > < 16 SNES Function norm 1.918e-10 > < 17 SNES Function norm 5.173e-11 > < 18 SNES Function norm 1.390e-11 > < 19 SNES Function norm < 1.e-11 > < 20 SNES Function norm < 1.e-11 > --- > > [0]PETSC ERROR: --------------------- Error Message ------------------------------------ > > [0]PETSC ERROR: Floating point exception! > > [0]PETSC ERROR: Vec entry at local location 14 is not-a-number or infinite at end of function: Parameter number 3! > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: Petsc Development HG revision: bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 -0500 > > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [0]PETSC ERROR: See docs/index.html for manual pages. > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: ./ex5 on a arch-gnu named Barrys-MacBook-Pro.local by barrysmith Mon Aug 13 21:12:19 2012 > > [0]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > > [0]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-triangle --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: SNESComputeFunction() line 1927 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESLineSearchApply_Basic() line 41 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: SNESLineSearchApply() line 506 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: SNESSolve_NRichardson() line 181 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/richardson/snesrichardson.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASUpSmooth_Private() line 447 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 695 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: main() line 152 in src/snes/examples/tutorials/ex5.c > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > > [cli_0]: aborting job: > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > /Users/barrysmith/Src/petsc-dev/src/snes/examples/tutorials > Possible problem with with ex5_5_fas, diffs above > ========================================= > 2,11c2,46 > < 1 SNES Function norm 0.156469 > < 2 SNES Function norm 0.0115623 > < 3 SNES Function norm 0.00145395 > < 4 SNES Function norm 0.000151214 > < 5 SNES Function norm 1.97033e-05 > < 6 SNES Function norm 2.60149e-06 > < 7 SNES Function norm 3.26556e-07 > < 8 SNES Function norm 4.31521e-08 > < 9 SNES Function norm 5.41224e-09 > < 10 SNES Function norm 1.18135e-09 > --- > > [0]PETSC ERROR: --------------------- Error Message ------------------------------------ > > [0]PETSC ERROR: Floating point exception! > > [0]PETSC ERROR: Vec entry at local location 14 is not-a-number or infinite at end of function: Parameter number 3! > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: Petsc Development HG revision: bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 -0500 > > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [0]PETSC ERROR: See docs/index.html for manual pages. > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: ./ex5 on a arch-gnu named Barrys-MacBook-Pro.local by barrysmith Mon Aug 13 21:12:20 2012 > > [0]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > > [0]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-triangle --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: SNESComputeFunction() line 1927 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESLineSearchApply_Basic() line 41 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: SNESLineSearchApply() line 506 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: SNESSolve_NRichardson() line 181 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/richardson/snesrichardson.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASUpSmooth_Private() line 447 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 695 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCoarseCorrection() line 567 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Multiplicative() line 694 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 774 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESSolve_NGMRES() line 259 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/ngmres/snesngmres.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: main() line 152 in src/snes/examples/tutorials/ex5.c > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > > [cli_0]: aborting job: > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > /Users/barrysmith/Src/petsc-dev/src/snes/examples/tutorials > Possible problem with with ex5_5_ngmres_fas, diffs above > ========================================= > 2,51c2,46 > < 1 SNES Function norm 0.648993 > < 2 SNES Function norm 0.518893 > < 3 SNES Function norm 0.439196 > < 4 SNES Function norm 0.386902 > < 5 SNES Function norm 0.350239 > < 6 SNES Function norm 0.322173 > < 7 SNES Function norm 0.298913 > < 8 SNES Function norm 0.279868 > < 9 SNES Function norm 0.26322 > < 10 SNES Function norm 0.249106 > < 11 SNES Function norm 0.236341 > < 12 SNES Function norm 0.225242 > < 13 SNES Function norm 0.214976 > < 14 SNES Function norm 0.205853 > < 15 SNES Function norm 0.168755 > < 16 SNES Function norm 0.150005 > < 17 SNES Function norm 0.146453 > < 18 SNES Function norm 0.122221 > < 19 SNES Function norm 0.109397 > < 20 SNES Function norm 0.107737 > < 21 SNES Function norm 0.106311 > < 22 SNES Function norm 0.0907226 > < 23 SNES Function norm 0.0813505 > < 24 SNES Function norm 0.0796973 > < 25 SNES Function norm 0.0782834 > < 26 SNES Function norm 0.0770559 > < 27 SNES Function norm 0.0758003 > < 28 SNES Function norm 0.0746744 > < 29 SNES Function norm 0.0740446 > < 30 SNES Function norm 0.061927 > < 31 SNES Function norm 0.0549197 > < 32 SNES Function norm 0.054583 > < 33 SNES Function norm 0.0519673 > < 34 SNES Function norm 0.0498744 > < 35 SNES Function norm 0.0491159 > < 36 SNES Function norm 0.0482318 > < 37 SNES Function norm 0.047538 > < 38 SNES Function norm 0.0468002 > < 39 SNES Function norm 0.0461691 > < 40 SNES Function norm 0.0459809 > < 41 SNES Function norm 0.0420109 > < 42 SNES Function norm 0.0389469 > < 43 SNES Function norm 0.0383224 > < 44 SNES Function norm 0.0378171 > < 45 SNES Function norm 0.0375595 > < 46 SNES Function norm 0.0373267 > < 47 SNES Function norm 0.0370906 > < 48 SNES Function norm 0.0368703 > < 49 SNES Function norm 0.0366439 > < 50 SNES Function norm 0.0364308 > --- > > 1 SNES Function norm 1.73715 > > 2 SNES Function norm 44.1655 > > 3 SNES Function norm 2.72388e+10 > > [0]PETSC ERROR: --------------------- Error Message ------------------------------------ > > [0]PETSC ERROR: Floating point exception! > > [0]PETSC ERROR: Vec entry at local location 8 is not-a-number or infinite at end of function: Parameter number 3! > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: Petsc Development HG revision: bd03fdcd3cf87c5b29ea28a5f4737ae936d9df04 HG Date: Sat Aug 11 20:22:42 2012 -0500 > > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > > [0]PETSC ERROR: See docs/index.html for manual pages. > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: ./ex5 on a arch-gnu named Barrys-MacBook-Pro.local by barrysmith Mon Aug 13 21:12:21 2012 > > [0]PETSC ERROR: Libraries linked from /Users/barrysmith/Src/petsc-dev/arch-gnu/lib > > [0]PETSC ERROR: Configure run at Sat Aug 11 20:25:02 2012 > > [0]PETSC ERROR: Configure options --download-blacs --download-fftw --download-hypre --download-libyaml --download-metis --download-ml --download-mpich --download-mumps --download-parmetis --download-ptscotch --download-scalapack --download-superlu --download-superlu_dist --download-triangle --download-yaml --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java --with-openmp --with-server --with-shared-libraries PETSC_ARCH=arch-gnu > > [0]PETSC ERROR: ------------------------------------------------------------------------ > > [0]PETSC ERROR: SNESComputeFunction() line 1927 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESLineSearchApply_Basic() line 41 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: SNESLineSearchApply() line 506 in /Users/barrysmith/Src/petsc-dev/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: SNESSolve_NRichardson() line 181 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/richardson/snesrichardson.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASDownSmooth_Private() line 420 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESFASCycle_Additive() line 638 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 776 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCycle_Additive() line 635 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 776 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCycle_Additive() line 635 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 776 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCycle_Additive() line 635 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 776 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCycle_Additive() line 635 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 776 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: SNESFASCycle_Additive() line 635 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve_FAS() line 776 in /Users/barrysmith/Src/petsc-dev/src/snes/impls/fas/fas.c > > [0]PETSC ERROR: SNESSolve() line 3570 in /Users/barrysmith/Src/petsc-dev/src/snes/interface/snes.c > > [0]PETSC ERROR: main() line 152 in src/snes/examples/tutorials/ex5.c > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > > [cli_0]: aborting job: > > application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0 > /Users/barrysmith/Src/petsc-dev/src/snes/examples/tutorials > Possible problem with with ex5_5_fas_additive, diffs above > > > >

1 0