[petsc-dev] annoying Nan in tests

2012-08-14 Thread Barry Smith
Peter, Thanks On Aug 13, 2012, at 10:43 PM, Peter Brune 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 reas

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

2012-08-14 Thread Barry Smith
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: >

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

2012-08-14 Thread Peter Brune
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 > = > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/1eacf6df/attachment.html>

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

2012-08-14 Thread Matthew Knepley
t;> > 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 &

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

2012-08-14 Thread Jed Brown
t -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/04771e47/attachment.html>

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

2012-08-14 Thread Peter Brune
olve() 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/6528dd2a/attachment.html>

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

2012-08-14 Thread Matthew Knepley
bert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/d2266f72/attachment.html>

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

2012-08-14 Thread Jed Brown
lots of locations and get the old value before resetting it. Passing NULL was a convenience. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/91ff0024/attachment-0001.html>

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

2012-08-14 Thread Matthew Knepley
egin 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/8d26d03c/attachment.html>

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

2012-08-14 Thread Peter Brune
HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/cd0a38be/attachment.html>

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

2012-08-14 Thread Barry Smith
On Aug 14, 2012, at 4:49 PM, Matthew Knepley wrote: > On Tue, Aug 14, 2012 at 4:47 PM, Jed Brown wrote: > On Tue, Aug 14, 2012 at 3:44 PM, Matthew Knepley wrote: > You can, but the setter routines interpret NULL as "ignore", i.e. > > if (func) ctx->func = func; > > Does this make sense to an

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

2012-08-14 Thread Peter Brune
t > > > > -- > > 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/dd814f41/attachment.html>

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

2012-08-14 Thread Barry Smith
On Aug 14, 2012, at 5:06 PM, Peter Brune wrote: > On a related note, it seems that we have this wonky notion of different code > paths for the same thing, in which we have SNESSetApplicationContext, > DMSetApplicationContext, and then the ctx arguments to > SNESSetFunction/Jacobian/GS/Picard.

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

2012-08-14 Thread Jed Brown
vinced to change my mind. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/5843cfbc/attachment.html>

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

2012-08-14 Thread Barry Smith
On Aug 14, 2012, at 5:16 PM, Jed Brown wrote: > On Tue, Aug 14, 2012 at 3:56 PM, Peter Brune wrote: > 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 especiall

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

2012-08-14 Thread Peter Brune
o set a coloring in an appropriate spot. - Peter -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/d8fccee3/attachment.html>

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

2012-08-14 Thread Barry Smith
On Aug 14, 2012, at 5:26 PM, Peter Brune wrote: > > On Tue, Aug 14, 2012 at 5:16 PM, Jed Brown wrote: > On Tue, Aug 14, 2012 at 3:56 PM, Peter Brune wrote: > Correct; the sheer amount of code using this convention makes it exceedingly > painful to fix. It would be much easier to alter our i

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

2012-08-14 Thread Jed Brown
SNES (optimization, UQ, etc), they should be able to compose using standard interfaces. This is why I think that SNES coloring is the misbehaving component. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/3cd662bf/attachment.html>

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

2012-08-14 Thread Jed Brown
On Tue, Aug 14, 2012 at 4:51 PM, Barry Smith wrote: > How about just having a SNESSetTS()? Please no. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120814/d94f520a/attachment.html>

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

2012-08-14 Thread Barry Smith
On Aug 14, 2012, at 5:58 PM, Jed Brown wrote: > On Tue, Aug 14, 2012 at 4:19 PM, Barry Smith wrote: > Ok, but then that context cannot be used for any other purpose because of > the possible conflict. Which means SNEScoloring thing must be changed. > >Why not use SNESSetApplicationContex

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

2012-08-14 Thread Barry Smith
On Aug 14, 2012, at 5:59 PM, Jed Brown wrote: > On Tue, Aug 14, 2012 at 4:51 PM, Barry Smith wrote: > How about just having a SNESSetTS()? > > Please no. But you've just told me that ONLY ts is allowed to use the SNESSetJacobian() context slot (coloring and presumably users) aren't allowe