[petsc-dev] broken TS example related to Jacobian and coloring etc.

Tue Aug 14 16:35:35 CDT 2012

On Tue, Aug 14, 2012 at 4:28 PM, Peter Brune <prbrune at 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 at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/d44b585a/attachment.html>