<div dir="ltr">Please respond to all, to get the petsc mailing list.<div><br></div><div>The out file here looks fine.  It would help of you run in debug mode (and get -g) so that gdb can give line numbers.</div><div><br></div><div>If the error is in the dot product it looks like you are getting an Inf or Nan someplace.  If the first evaluation of the norm was the problem then run the bad code and add a VecView and MatView, and add a VecNorm and MatNorm so that the error happens in your code, just to check.  If you had 1e200 number in the vector then the 2 norm would give you this trap for overflow.  You could test with the inf norm and it works then it is an overflow problem.</div><div><div><br></div><div><br><div><br></div><div><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Greg Miller</b> <span dir="ltr"><<a href="mailto:grgmiller@ucdavis.edu" target="_blank">grgmiller@ucdavis.edu</a>></span><br>Date: Mon, Jul 13, 2015 at 3:54 PM<br>Subject: Re: Fwd: same petsc problem<br>To: Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>><br>Cc: "<a href="mailto:treb@lbl.gov" target="_blank">treb@lbl.gov</a>" <<a href="mailto:treb@lbl.gov" target="_blank">treb@lbl.gov</a>><br><br><br>thanks.  here it is.<br>
<br>
I'm not getting the nan problem with this code - not sure why.  NaN  was showing up KSP_solve in the first evaluation of the vector norm.<br>
In going from my original example to the minimal code I mailed today the matrix coefficients got rounded through asci printing.<br>
That's the only difference I can think of.  The gdb stack from the old code is shown below.<br>
G<br>
<br>
<br>
<br>
Program received signal SIGFPE, Arithmetic exception.<br>
0x00007ffff6c4ed15 in ddot_ () from /usr/lib/libblas.so.3<br>
(gdb) where<br>
#0  0x00007ffff6c4ed15 in ddot_ () from /usr/lib/libblas.so.3<br>
#1  0x00000000008f87da in VecNorm_Seq (xin=0x22e7ca0, type=NORM_2, z=0x7fffffffba30)<br>
    at /home/usr/local/src/petsc-3.5.3/src/vec/vec/impls/seq/bvec2.c:614<br>
#2  0x00000000008c9eea in VecNorm (x=0x22e7ca0, type=NORM_2, val=0x7fffffffba30)<br>
    at /home/usr/local/src/petsc-3.5.3/src/vec/vec/interface/rvector.c:242<br>
#3  0x00000000008cab40 in VecNormalize (x=0x22e7ca0, val=0x7fffffffbaa0) at /home/usr/local/src/petsc-3.5.3/src/vec/vec/interface/rvector.c:337<br>
#4  0x0000000000d44a9d in KSPGMRESCycle (itcount=0x7fffffffbb08, ksp=0x2288ff0)<br>
    at /home/usr/local/src/petsc-3.5.3/src/ksp/ksp/impls/gmres/gmres.c:161<br>
#5  0x0000000000d453e3 in KSPSolve_GMRES (ksp=0x2288ff0) at /home/usr/local/src/petsc-3.5.3/src/ksp/ksp/impls/gmres/gmres.c:235<br>
<span><br>
<br>
On 07/13/2015 11:54 AM, Mark Adams wrote:<br>
> Greg, I am forwarding this to the PETSc mailing list.<br>
><br>
> Please send the entire output from this run.  As I recall you were getting a message that all values were not the same on all processors in GMRES.  I have seen this when I get NaNs in the system.<br>
><br>
> While you are doing this you should use a simple solver like change:<br>
><br>
> -pressure_pc_type gamg<br>
><br>
> to<br>
><br>
> -pressure_pc_type jacobi<br>
><br>
><br>
> And add:<br>
><br>
</span>> *-*pressure_*ksp_monitor_true_residual*<br>
<span>><br>
> Mark<br>
><br>
><br>
> ---------- Forwarded message ----------<br>
</span><span>> From: *Greg Miller* <<a href="mailto:grgmiller@ucdavis.edu" target="_blank">grgmiller@ucdavis.edu</a> <mailto:<a href="mailto:grgmiller@ucdavis.edu" target="_blank">grgmiller@ucdavis.edu</a>>><br>
> Date: Mon, Jul 13, 2015 at 2:08 PM<br>
> Subject: same petsc problem<br>
</span><span>> To: Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a> <mailto:<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>>><br>
> Cc: David Trebotich <<a href="mailto:treb@hpcrd.lbl.gov" target="_blank">treb@hpcrd.lbl.gov</a> <mailto:<a href="mailto:treb@hpcrd.lbl.gov" target="_blank">treb@hpcrd.lbl.gov</a>>><br>
><br>
><br>
> Hi Mark.  I'm still stuck on the same petsc problem.   Would you please try the attached minimal example and advise me?<br>
><br>
> I'm running this without MPI:<br>
> make DIM=2 DEBUG=TRUE MPI=FALSE USE_PETSC=TRUE test<br>
><br>
> There is no input file.<br>
><br>
> Thank you,<br>
> Greg<br>
><br>
> --<br>
> Greg Miller<br>
> Department of Chemical Engineering and Materials Science<br>
> University of California, Davis<br>
> One Shields Avenue<br>
> Davis, CA 95616<br>
</span>> <a href="mailto:grgmiller@ucdavis.edu" target="_blank">grgmiller@ucdavis.edu</a> <mailto:<a href="mailto:grgmiller@ucdavis.edu" target="_blank">grgmiller@ucdavis.edu</a>><br>
<div><div>><br>
<br>
--<br>
Greg Miller<br>
Department of Chemical Engineering and Materials Science<br>
University of California, Davis<br>
One Shields Avenue<br>
Davis, CA 95616<br>
<a href="mailto:grgmiller@ucdavis.edu" target="_blank">grgmiller@ucdavis.edu</a><br>
</div></div></div><br></div></div></div></div>