[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Barry Smith
On Dec 13, 2012, at 7:03 PM, Jed Brown wrote: > On Thu, Dec 13, 2012 at 5:00 PM, Barry Smith wrote: > I suggest > > -snes_check_jacobian_view > > then MatView is just called on the appropriate viewer instead of ugly variant > options like snes_check_jacobian_draw etc > > In fact

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Barry Smith
On Dec 13, 2012, at 5:34 PM, Jed Brown wrote: > Options for drawing the result, using coloring instead of dense differencing, > and thresholding using rtol and atol to avoid false positives. IIRC, the > reason I added all of those variants was because I was given some Jacobian > code that was

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Peter Brune
Run with -snes_check_jacobian_display to show difference of >> hand-coded and finite difference Jacobian. >> > 1.32338e-09 = ||J - Jfd||//J|| 5.27058e-07 = ||J - Jfd|| >> > 1 SNES Function norm 4.625448863612e-05 >> > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >> hand-coded Jacobian is probably correct. >> > Run with -snes_check_jacobian_display to show difference of >> hand-coded and finite difference Jacobian. >> > 1.57582e-09 = ||J - Jfd||//J|| 6.27588e-07 = ||J - Jfd|| >> > 2 SNES Function norm 4.391784647745e-14 >> > 0 SNES Function norm 6.761591873515e-01 >> > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >> hand-coded Jacobian is probably correct. >> > >> > >> > >> > >> >> > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/0827afed/attachment-0001.html>

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Peter Brune
heck_jacobian_display to show difference of > hand-coded and finite difference Jacobian. > > 1.57582e-09 = ||J - Jfd||//J|| 6.27588e-07 = ||J - Jfd|| > > 2 SNES Function norm 4.391784647745e-14 > > 0 SNES Function norm 6.761591873515e-01 > > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the > hand-coded Jacobian is probably correct. > > > > > > > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/d649e8fb/attachment.html>

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Barry Smith
On Dec 13, 2012, at 5:04 PM, Peter Brune wrote: > Dangit Barry! I already started a patch moving all Jed's -snes_compare_* > garbage in SNESComputeJacobian() over to a SNESJacobianMonitor registration > type deal where it has functions you can register to run after computing the > Jacobian.

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Barry Smith
On Dec 13, 2012, at 4:59 PM, Jed Brown wrote: > Surely you were aware of the -snes_compare_* options that can be activated in > SNESComputeJacobian(). Absolutely not. And why in hell would someone put them in a computational routine SNESComputeJacobian() (and not tell anyone about them?).

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Peter Brune
io is O(1.e-8), the >> hand-coded Jacobian is probably correct. >> Run with -snes_check_jacobian_display to show difference of >> hand-coded and finite difference Jacobian. >> 1.57582e-09 = ||J - Jfd||//J|| 6.27588e-07 = ||J - Jfd|| >> 2 SNES Function norm 4.391784647745e-14 >> 0 SNES Function norm 6.761591873515e-01 >> Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >> hand-coded Jacobian is probably correct. >> >> >> > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/0fd37597/attachment.html>

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Jed Brown
lename.dat"? > >Barry > > We'll write a utility routine that when given a viewer type returns the > appropriate canonical form. > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/4a0983f3/attachment.html>

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Barry Smith
I got sick of -snes_type test only checking the Jacobian at the initial solution and not being able to test later Jacobians so I added -snes_check_jacobian that allows the regular running of SNESSolve() (also inside TS) and compares every Jacobian computed with a finite difference version.

[petsc-dev] example incorrectly labeled please fix for nightly tests

2012-12-13 Thread Jed Brown
6 > < timestep 5: time 0.05, solution norm 0.437213 > /Users/barrysmith/Src/petsc-dev/src/ts/examples/tutorials > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/53e8554e/attachment.html>

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Satish Balay
Just a note: about the invoking only 'make -C {PETSC_DIR}/${PETSC_ARCH}' It will miss redoing the dsymutil stuff - that 'cmake:' currently does. Also I avoid -C in target 'cmake' [and use 'cd' instead] becuase -C is gnu-make extension - and not portable. Satish On Thu, 13 Dec 2012, Satish Balay

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Satish Balay
you would get that behavior with: make cmake VERBOSE=1 Satish On Thu, 13 Dec 2012, Barry Smith wrote: > > Satish, > >No need to do this change. I would like the current make cmake to continue > to behave exactly as it does now. > >Barry > > On Dec 13, 2012, at 3:37 PM, Satish Bal

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Barry Smith
Satish, No need to do this change. I would like the current make cmake to continue to behave exactly as it does now. Barry On Dec 13, 2012, at 3:37 PM, Satish Balay wrote: > On Thu, 13 Dec 2012, Barry Smith wrote: > >> >> >> On Dec 13, 2012, at 3:18 PM, Jed Brown wrote: >> >>> I

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Satish Balay
On Thu, 13 Dec 2012, Barry Smith wrote: > > > On Dec 13, 2012, at 3:18 PM, Jed Brown wrote: > > > It is displaying all the debugging information. You'll only get the useful > > information if you run this: > > > > make -j5 -C arch-gnu > >Yes, your suggestion is not verbose because it do

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Jed Brown
> > 1.10818e-10 = ||J - Jfd||//J|| 4.40468e-07 = ||J - Jfd|| >>> > 2 SNES Function norm 7.677760150017e-14 >>> > 0 SNES Function norm 2.179048757949e+00 >>> > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >>> hand-coded Jacobian is probably correct. >>> > Run with -snes_check_jacobian_display to show difference of >>> hand-coded and finite difference Jacobian. >>> > 1.32338e-09 = ||J - Jfd||//J|| 5.27058e-07 = ||J - Jfd|| >>> > 1 SNES Function norm 4.625448863612e-05 >>> > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >>> hand-coded Jacobian is probably correct. >>> > Run with -snes_check_jacobian_display to show difference of >>> hand-coded and finite difference Jacobian. >>> > 1.57582e-09 = ||J - Jfd||//J|| 6.27588e-07 = ||J - Jfd|| >>> > 2 SNES Function norm 4.391784647745e-14 >>> > 0 SNES Function norm 6.761591873515e-01 >>> > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >>> hand-coded Jacobian is probably correct. >>> > >>> > >>> > >>> > >>> >>> >> > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/1ca7718e/attachment.html>

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Barry Smith
I don't know the meaning of -C for make and was too lazy to do man make until after receiving you alls mail :-( I'm pretty sure make didn't support -C in 1986 On Dec 13, 2012, at 3:27 PM, Sean Farley wrote: > On Thu, Dec 13, 2012 at 3:23 PM, Barry Smith wrote: >> >> >> On Dec 13,

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Sean Farley
On Thu, Dec 13, 2012 at 3:23 PM, Barry Smith wrote: > > > On Dec 13, 2012, at 3:18 PM, Jed Brown wrote: > >> It is displaying all the debugging information. You'll only get the useful >> information if you run this: >> >> make -j5 -C arch-gnu > >Yes, your suggestion is not verbose because it

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Barry Smith
On Dec 13, 2012, at 3:18 PM, Jed Brown wrote: > It is displaying all the debugging information. You'll only get the useful > information if you run this: > > make -j5 -C arch-gnu Yes, your suggestion is not verbose because it doesn't do anything Barrys-MacBook-Pro:tutorials barrysmith$ m

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Jed Brown
3612e-05 > > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the > hand-coded Jacobian is probably correct. > > Run with -snes_check_jacobian_display to show difference of > hand-coded and finite difference Jacobian. > > 1.57582e-09 = ||J - Jfd||//J|| 6.27588e-07 = ||J - Jfd|| > > 2 SNES Function norm 4.391784647745e-14 > > 0 SNES Function norm 6.761591873515e-01 > > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the > hand-coded Jacobian is probably correct. > > > > > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/9c84c69e/attachment.html>

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Jed Brown
- Jfd|| >>> 1 SNES Function norm 4.625448863612e-05 >>> Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >>> hand-coded Jacobian is probably correct. >>> Run with -snes_check_jacobian_display to show difference of >>> hand-coded and finite difference Jacobian. >>> 1.57582e-09 = ||J - Jfd||//J|| 6.27588e-07 = ||J - Jfd|| >>> 2 SNES Function norm 4.391784647745e-14 >>> 0 SNES Function norm 6.761591873515e-01 >>> Testing hand-coded Jacobian, if the ratio is O(1.e-8), the >>> hand-coded Jacobian is probably correct. >>> >>> >>> >> > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/37b44a0a/attachment-0001.html>

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Barry Smith
Jed, When I have edited PETSc source and want to update the libraries I use make cmake. It produces so much clutter that it is hard to see the relevant bits. Could there be an alternative make cmakeclear that only displays the files that have changed (and hence will be compiled), the comp

[petsc-dev] -snes_check_jacobian

2012-12-13 Thread Jed Brown
> 2 SNES Function norm 4.391784647745e-14 > 0 SNES Function norm 6.761591873515e-01 > Testing hand-coded Jacobian, if the ratio is O(1.e-8), the > hand-coded Jacobian is probably correct. > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/6784addf/attachment.html>

[petsc-dev] a TS feature

2012-12-13 Thread Shri
s is something Shri needs. > > > > Is this general enough to add to the library itself? He can write an > > adaptive controller that hits any times he wants. > > > > As a library feature, I'd be more interested in having a controller that > > attempted to hit the zeros of a user-provided functional (because that > > could be used for other forms of model nonsmoothness). > > > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/3ed91320/attachment.html>

[petsc-dev] "Libraries don't have to suck"

2012-12-13 Thread Jed Brown
-built-on.html> > > > > On 11/26/2012 06:08 PM, Jed Brown wrote: > >> Point in favor of evolutionary libraries over the Matryoshka dolls that >> arise when interfaces are frozen forever. >> >> http://akkartik.name/blog/**libraries2<http://akkartik.name/blog/libraries2> >> > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/bfac2cdc/attachment-0001.html>

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Jed Brown
Just do make -C $PETSC_DIR/$PETSC_ARCH On Thu, Dec 13, 2012 at 1:23 PM, Barry Smith wrote: > > > On Dec 13, 2012, at 3:18 PM, Jed Brown wrote: > > > It is displaying all the debugging information. You'll only get the > useful information if you run this: > > > > make -j5 -C arch-gnu > >Yes

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Matthew Knepley
On Thu, Dec 13, 2012 at 1:23 PM, Barry Smith wrote: > > > On Dec 13, 2012, at 3:18 PM, Jed Brown wrote: > >> It is displaying all the debugging information. You'll only get the useful >> information if you run this: >> >> make -j5 -C arch-gnu > >Yes, your suggestion is not verbose because it

[petsc-dev] making make cmake less verbose

2012-12-13 Thread Jed Brown
It is displaying all the debugging information. You'll only get the useful information if you run this: make -j5 -C arch-gnu We can make a top-level target that does this without the verbose info, but make doesn't have a nice way to send verbose information to a file while displaying only terse i

[petsc-dev] a TS feature

2012-12-13 Thread Jed Brown
ntroller that hits any times he wants. > > > > As a library feature, I'd be more interested in having a controller that > attempted to hit the zeros of a user-provided functional (because that > could be used for other forms of model nonsmoothness). > > > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/f592c42a/attachment.html>

[petsc-dev] Fwd: [Seniors] Talk on Peer Education For Scientific Computing

2012-12-13 Thread Barry Smith
for > Research Scientists," September 6, 2006, http://nanohub.org/resources/1811. > > ___ > Seniors mailing list > Seniors at ci.uchicago.edu > https://lists.ci.uchicago.edu/cgi-bin/mailman/listinfo/seniors -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/0e3ab22a/attachment.html>

[petsc-dev] Fwd: Is there any example that allows time-integration provided by users

2012-12-13 Thread Peter Brune
niversity of Science and Technology > >> Tel??966?0544701599 > >> > >> > >> This message and its contents, including attachments are intended > solely for the original recipient. If you are not the intended recipient or > have received this message in error, please notify me immediately and > delete this message from your computer system. Any unauthorized use or > distribution is prohibited. Please consider the environment before printing > this email. > >> > > > > > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/8fa34fb7/attachment.html>

[petsc-dev] Fwd: Is there any example that allows time-integration provided by users

2012-12-13 Thread Barry Smith
On Dec 13, 2012, at 11:27 AM, Peter Brune wrote: > An abandoned attempt at this lies dormant in src/snes/impls/multiblock. We > could try to revive it. Or at least give it a reasonable name :-( like SNESFIELDSPLIT > > - Peter > > > On Thu, Dec 13, 2012 at 11:25 AM, Jed Brown wrote: >

[petsc-dev] Fwd: Is there any example that allows time-integration provided by users

2012-12-13 Thread Peter Brune
ng Abdullah University of Science and Technology >>> Tel??966?0544701599 >>> >>> >>> -- >>> This message and its contents, including attachments are intended solely >>> for the original recipient. If you are not the i

[petsc-dev] Fwd: Is there any example that allows time-integration provided by users

2012-12-13 Thread Barry Smith
s email. > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/65986587/attachment.html>

[petsc-dev] Fwd: Is there any example that allows time-integration provided by users

2012-12-13 Thread Matthew Knepley
On Thu, Dec 13, 2012 at 9:30 AM, Peter Brune wrote: > > > > On Thu, Dec 13, 2012 at 11:29 AM, Barry Smith wrote: >> >> >> On Dec 13, 2012, at 11:27 AM, Peter Brune wrote: >> >> > An abandoned attempt at this lies dormant in src/snes/impls/multiblock. >> > We could try to revive it. >> >>Or a

[petsc-dev] Fwd: [Seniors] Talk on Peer Education For Scientific Computing

2012-12-13 Thread Jed Brown
mailing list > Seniors at ci.uchicago.edu > https://lists.ci.uchicago.edu/cgi-bin/mailman/listinfo/seniors > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/17e6c092/attachment-0001.html>

[petsc-dev] Fwd: Is there any example that allows time-integration provided by users

2012-12-13 Thread Jed Brown
ge from your computer system. Any unauthorized use or >> distribution is prohibited. Please consider the environment before printing >> this email. > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121213/23eafd6f/attachment.html>