<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix"><br>
<pre class="moz-signature" cols="72">
</pre>
On 20/7/2013 4:08 PM, Matthew Knepley wrote:<br>
</div>
<blockquote
cite="mid:CAMYG4G=93aMhMv45eCNBkYy7BARvghCGF_vD0zsrn4KkoEFRYg@mail.gmail.com"
type="cite">
<div dir="ltr">On Sat, Jul 20, 2013 at 8:39 AM, TAY wee-beng <span
dir="ltr"><<a moz-do-not-send="true"
href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</a>></span>
wrote:<br>
<div class="gmail_extra">
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF"> Hi,<br>
<br>
I'm trying to use GAMG to speed up solving of Poisson
eqn. I used:<br>
<br>
call KSPSetOptionsPrefix(ksp,"poisson_",ierr)<br>
<br>
-poisson_pc_gamg_agg_nsmooths 1 -poisson_pc_type gamg<br>
<br>
I remember it used to work in some problems but now it
can't work, with error:<br>
</div>
</blockquote>
<div><br>
</div>
<div>Upgrade to the latest release and we will help you
debug this.</div>
<div><br>
</div>
<div> Matt</div>
</div>
</div>
</div>
</blockquote>
Hi Matt,<br>
<br>
I've used the latest 3.4.2. What else do I need to provide?<br>
<br>
Thanks!<br>
<blockquote
cite="mid:CAMYG4G=93aMhMv45eCNBkYy7BARvghCGF_vD0zsrn4KkoEFRYg@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF"> <b><i>[0]PETSC
ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Caught signal number 8
FPE: Floating Point Exception,probably divide by
zero</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Try option
-start_in_debugger or -on_error_attach_debugger</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: or see <a
moz-do-not-send="true"
href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind"
target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind</a>[0]PETSC
ERROR: or try <a moz-do-not-send="true"
href="http://valgrind.org" target="_blank">http://valgrind.org</a>
on GNU/linux and Apple Mac OS X to find memory
corruption errors</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: likely location of
problem given in stack below</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: ---------------------
Stack Frames ------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Note: The EXACT line
numbers in the stack are not available,</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: INSTEAD the line
number of the start of the function</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: is given.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0]
KSPComputeExtremeSingularValues_GMRES line 24
src/ksp/ksp/impls/gmres/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\impls\gmres\gmreig.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0]
KSPComputeExtremeSingularValues line 40
src/ksp/ksp/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\INTERF~1\itfunc.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] PCGAMGOptprol_AGG
line 1295
src/ksp/pc/impls/gamg/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\pc\impls\gamg\agg.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] PCSetUp_GAMG line
564
src/ksp/pc/impls/gamg/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\pc\impls\gamg\gamg.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] PCSetUp line 810
src/ksp/pc/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\pc\INTERF~1\precon.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] KSPSetUp line 182
src/ksp/ksp/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\INTERF~1\itfunc.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] KSPSolve line 351
src/ksp/ksp/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\INTERF~1\itfunc.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: ---------------------
Error Message ------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Signal received!</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Petsc Development HG
revision: 9850aeb5d33f0b33bc931843c4b3b3b4f8df6a3b
HG Date: Tue Oct 02 22:18:53 2012 -0500</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: See
docs/changes/index.html for recent updates.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: See docs/faq.html for
hints about trouble shooting.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: See docs/index.html for
manual pages.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
C:\Obj_tmp\ibm3d_high_Re_staggered_AB2\Debug\ibm3d_high_Re_staggered_AB2.exe
on a petsc-3.3 named USER-PC by User Sat Jul 20
15:37:40 2013</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Libraries linked from
/cygdrive/d/wtay/Lib/petsc-3.3-dev_win32_vs2008/lib</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Configure run at Thu
Oct 4 10:01:13 2012</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Configure options
--with-cc="win32fe cl" --with-fc="win32fe ifort"
--with-cxx="win32fe cl"
--with-mpi-dir=/cygdrive/c/MPICH2/
--download-f-blas-lapack=1
--prefix=/cygdrive/d/wtay/Lib/petsc-3.3-dev_win32_vs2008
--with-debugging=1 --useThreads=0</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: User provided function()
line 0 in unknown directory unknown file</i></b><b><i><br>
</i></b><b><i>application called
MPI_Abort(MPI_COMM_WORLD, 59) - process 0</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i>job aborted:</i></b><b><i><br>
</i></b><b><i>rank: node: exit code[: error message]</i></b><b><i><br>
</i></b><b><i>0: User-PC: 59: process 0 exited without
calling finalize</i></b><i><br>
</i><br>
<br>
I read in one of the threads that I can use:<br>
<br>
-pc_type gamg -pc_gamg_agg_nsmooths 1
-mg_levels_ksp_type richardson -mg_levels_pc_type sor<br>
<br>
It worked but I got the msg:<br>
<br>
WARNING! There are options you set that were not used!<br>
WARNING! could be spelling mistake, etc!<br>
Option left: name:-mg_levels_ksp_type value: richardson<br>
Option left: name:-mg_levels_pc_type value: sor<br>
Option left: name:-pc_gamg_agg_nsmooths value: 1<br>
Option left: name:-pc_type value: gamg<br>
Press any key to continue . . .<br>
<br>
If I used this:<br>
<br>
-poisson_pc_type gamg -poisson_pc_gamg_agg_nsmooths 1
-poisson_mg_levels_ksp_type richardson
-poisson_mg_levels_pc_type sor<br>
<br>
It aborts with error:<br>
<br>
<b><i>[0]PETSC ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Caught signal number 8
FPE: Floating Point Exception,probably divide by
zero</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Try option
-start_in_debugger or -on_error_attach_debugger</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: or see <a
moz-do-not-send="true"
href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind"
target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind</a>[0]PETSC
ERROR: or try <a moz-do-not-send="true"
href="http://valgrind.org" target="_blank">http://valgrind.org</a>
on GNU/linux and Apple Mac OS X to find memory
corruption errors</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: likely location of
problem given in stack below</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: ---------------------
Stack Frames ------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Note: The EXACT line
numbers in the stack are not available,</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: INSTEAD the line
number of the start of the function</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: is given.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0]
KSPComputeExtremeSingularValues_GMRES line 24
src/ksp/ksp/impls/gmres/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\impls\gmres\gmreig.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0]
KSPComputeExtremeSingularValues line 40
src/ksp/ksp/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\INTERF~1\itfunc.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] PCGAMGOptprol_AGG
line 1295
src/ksp/pc/impls/gamg/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\pc\impls\gamg\agg.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] PCSetUp_GAMG line
564
src/ksp/pc/impls/gamg/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\pc\impls\gamg\gamg.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] PCSetUp line 810
src/ksp/pc/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\pc\INTERF~1\precon.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] KSPSetUp line 182
src/ksp/ksp/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\INTERF~1\itfunc.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: [0] KSPSolve line 351
src/ksp/ksp/interface/C:\wtay\DOWNLO~1\Codes\PETSC-~1\src\ksp\ksp\INTERF~1\itfunc.c</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: ---------------------
Error Message ------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Signal received!</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Petsc Development HG
revision: 9850aeb5d33f0b33bc931843c4b3b3b4f8df6a3b
HG Date: Tue Oct 02 22:18:53 2012 -0500</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: See
docs/changes/index.html for recent updates.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: See docs/faq.html for
hints about trouble shooting.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: See docs/index.html for
manual pages.</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
C:\Obj_tmp\ibm3d_high_Re_staggered_AB2\Debug\ibm3d_high_Re_staggered_AB2.exe
on a petsc-3.3 named USER-PC by User Sat Jul 20
15:36:24 2013</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Libraries linked from
/cygdrive/d/wtay/Lib/petsc-3.3-dev_win32_vs2008/lib</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Configure run at Thu
Oct 4 10:01:13 2012</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: Configure options
--with-cc="win32fe cl" --with-fc="win32fe ifort"
--with-cxx="win32fe cl"
--with-mpi-dir=/cygdrive/c/MPICH2/
--download-f-blas-lapack=1
--prefix=/cygdrive/d/wtay/Lib/petsc-3.3-dev_win32_vs2008
--with-debugging=1 --useThreads=0</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR:
------------------------------------------------------------------------</i></b><b><i><br>
</i></b><b><i>[0]PETSC ERROR: User provided function()
line 0 in unknown directory unknown file</i></b><b><i><br>
</i></b><b><i>application called
MPI_Abort(MPI_COMM_WORLD, 59) - process 0</i></b><b><i><br>
</i></b><b><i><br>
</i></b><b><i>job aborted:</i></b><b><i><br>
</i></b><b><i>rank: node: exit code[: error message]</i></b><b><i><br>
</i></b><b><i>0: User-PC: 59: process 0 exited without
calling finalize</i></b><b><i><br>
</i></b><br>
<br>
<br>
So is there a recommended test command or method?<br>
<br>
Thank!<span class="HOEnZb"><font color="#888888"><br>
<pre cols="72">--
Yours sincerely,
TAY wee-beng</pre>
</font></span></div>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
What most experimenters take for granted before they begin
their experiments is infinitely more interesting than any
results to which their experiments lead.<br>
-- Norbert Wiener
</div>
</div>
</blockquote>
<br>
</body>
</html>