<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi,<br>
<br>
Sorry for the mistake. I tried again, changing to "call
VecView(x,0,ierr)"<br>
<br>
but still got the same error:<br>
<br>
[wtay@hpc12:tutorials]$ ./ex22f<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation
Violation, probably memory access out of range<br>
[0]PETSC ERROR: Try option -start_in_debugger or
-on_error_attach_debugger<br>
[0]PETSC ERROR: or see
<a class="moz-txt-link-freetext" href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/faq.html#valgrind">http://www.mcs.anl.gov/petsc/petsc-as/documentation/faq.html#valgrind</a>[0]PETSC
ERROR: or try <a class="moz-txt-link-freetext" href="http://valgrind.org">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to
find memory corruption errors<br>
[0]PETSC ERROR: likely location of problem given in stack below<br>
[0]PETSC ERROR: --------------------- Stack Frames
------------------------------------<br>
[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not
available,<br>
[0]PETSC ERROR: INSTEAD the line number of the start of the
function<br>
[0]PETSC ERROR: is given.<br>
[0]PETSC ERROR: [0] VecView line 735
/home/wtay/Codes/petsc-3.2-p5/src/vec/vec/interface/vector.c<br>
[0]PETSC ERROR: --------------------- Error Message
------------------------------------<br>
[0]PETSC ERROR: Signal received!<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 5, Sat Oct 29
13:45:54 CDT 2011 <br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: ./ex22f on a arch-linu named hpc12 by wtay Sun Apr
22 21:25:01 2012<br>
[0]PETSC ERROR: Libraries linked from
/home/wtay/Lib/petsc-3.2-p5_mumps_debug/lib<br>
[0]PETSC ERROR: Configure run at Sun Nov 27 15:39:26 2011<br>
[0]PETSC ERROR: Configure options --with-mpi-dir=/opt/openmpi-1.5.3/
--with-blas-lapack-dir=/opt/intel_xe_2011/mkl/lib/intel64/
--with-debugging=1 --download-hypre=1
--prefix=/home/wtay/Lib/petsc-3.2-p5_mumps_debug COPTFLAGS=-O0
FOPTFLAGS=-O0 --download-mumps=1 --download-parmetis=1
--download-scalapack=1 --download-blacs=1<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: User provided function() line 0 in unknown directory
unknown file<br>
--------------------------------------------------------------------------<br>
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD <br>
with errorcode 59.<br>
<br>
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.<br>
You may or may not see output from other processes, depending on<br>
exactly when Open MPI kills them.<br>
--------------------------------------------------------------<br>
<pre class="moz-signature" cols="72">Yours sincerely,
TAY wee-beng</pre>
<br>
On 22/4/2012 9:20 PM, Jed Brown wrote:
<blockquote
cite="mid:CAM9tzSmN7r4LaQiViKV8phtY3_NwTmPAr+Q+GtmS2CKqF4-9Mw@mail.gmail.com"
type="cite">
<p>Look at the man page for VecView, it needs a PetscViewer
argument.</p>
<div class="gmail_quote">On Apr 22, 2012 2:19 PM, "TAY wee-beng"
<<a moz-do-not-send="true" href="mailto:zonexo@gmail.com">zonexo@gmail.com</a>>
wrote:<br type="attribution">
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> Hi,<br>
<br>
I have attached the ex22f.F file. The changes I added are
given in bold:<br>
<br>
...<br>
<br>
PetscErrorCode ierr<br>
DM da<br>
KSP ksp<br>
<b> Vec x,b</b><br>
external ComputeRHS,ComputeMatrix<br>
<br>
....<br>
<br>
call KSPSetUp(ksp,ierr)<br>
call
KSPSolve(ksp,PETSC_NULL_OBJECT,PETSC_NULL_OBJECT,ierr)<br>
<b>call KSPGetSolution(ksp,x,ierr)<br>
call VecView(x,ierr)</b><br>
call KSPDestroy(ksp,ierr)<br>
call DMDestroy(da,ierr)<br>
<br>
The error is:<br>
<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation
Violation, probably memory access out of range<br>
[0]PETSC ERROR: Try option -start_in_debugger or
-on_error_attach_debugger<br>
[0]PETSC ERROR: or see <a moz-do-not-send="true"
href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/faq.html#valgrind"
target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/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<br>
[0]PETSC ERROR: likely location of problem given in stack
below<br>
[0]PETSC ERROR: --------------------- Stack Frames
------------------------------------<br>
[0]PETSC ERROR: Note: The EXACT line numbers in the stack
are not available,<br>
[0]PETSC ERROR: INSTEAD the line number of the start
of the function<br>
[0]PETSC ERROR: is given.<br>
[0]PETSC ERROR: [0] VecView line 735
/home/wtay/Codes/petsc-3.2-p5/src/vec/vec/interface/vector.c<br>
[0]PETSC ERROR: --------------------- Error Message
------------------------------------<br>
[0]PETSC ERROR: Signal received!<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 5, Sat
Oct 29 13:45:54 CDT 2011 <br>
[0]PETSC ERROR: See docs/changes/index.html for recent
updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble
shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: ./ex22f on a arch-linu named hpc12 by wtay
Sun Apr 22 21:11:39 2012<br>
[0]PETSC ERROR: Libraries linked from
/home/wtay/Lib/petsc-3.2-p5_mumps_debug/lib<br>
[0]PETSC ERROR: Configure run at Sun Nov 27 15:39:26 2011<br>
[0]PETSC ERROR: Configure options
--with-mpi-dir=/opt/openmpi-1.5.3/
--with-blas-lapack-dir=/opt/intel_xe_2011/mkl/lib/intel64/
--with-debugging=1 --download-hypre=1
--prefix=/home/wtay/Lib/petsc-3.2-p5_mumps_debug
COPTFLAGS=-O0 FOPTFLAGS=-O0 --download-mumps=1
--download-parmetis=1 --download-scalapack=1
--download-blacs=1<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: User provided function() line 0 in unknown
directory unknown file<br>
--------------------------------------------------------------------------<br>
MPI_ABORT was invoked on rank 0 in communicator
MPI_COMM_WORLD <br>
with errorcode 59.<br>
<br>
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI
processes.<br>
You may or may not see output from other processes,
depending on<br>
exactly when Open MPI kills them.<br>
<br>
<pre cols="72">Yours sincerely,
TAY wee-beng</pre>
<br>
On 22/4/2012 9:06 PM, Jed Brown wrote:
<blockquote type="cite">
<p>Run in a debugger and/or use --with-debugging=1 so that
the error trace has more information. You could also
show us the exact code that you used.</p>
<div class="gmail_quote">On Apr 22, 2012 2:03 PM, "TAY
wee-beng" <<a moz-do-not-send="true"
href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</a>>
wrote:<br type="attribution">
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> Hi,<br>
<br>
I added "Vec x,b" after "KSP ksp"<br>
and then "call KSPGetSolution(ksp, x, ierr)"<br>
<br>
I wanted to see the output so I added "call
VecView(x,ierr)" but I got this error:<br>
<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Caught signal number 11 SEGV:
Segmentation Violation, probably memory access out
of range<br>
[0]PETSC ERROR: Try option -start_in_debugger or
-on_error_attach_debugger<br>
[0]PETSC ERROR: or see <a moz-do-not-send="true"
href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/faq.html#valgrind"
target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/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<br>
[0]PETSC ERROR: configure using
--with-debugging=yes, recompile, link, and run <br>
[0]PETSC ERROR: to get more information on the
crash.<br>
[0]PETSC ERROR: --------------------- Error Message
------------------------------------<br>
[0]PETSC ERROR: Signal received!<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch
5, Sat Oct 29 13:45:54 CDT 2011 <br>
[0]PETSC ERROR: See docs/changes/index.html for
recent updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about
trouble shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual
pages.<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: ./ex22f on a arch-linu named hpc12
by wtay Sun Apr 22 21:02:14 2012<br>
[0]PETSC ERROR: Libraries linked from
/home/wtay/Lib/petsc-3.2-p5_mumps_rel/lib<br>
[0]PETSC ERROR: Configure run at Sun Nov 27 15:18:15
2011<br>
[0]PETSC ERROR: Configure options
--with-mpi-dir=/opt/openmpi-1.5.3/
--with-blas-lapack-dir=/opt/intel_xe_2011/mkl/lib/intel64/
--with-debugging=0 --download-hypre=1
--prefix=/home/wtay/Lib/petsc-3.2-p5_mumps_rel
COPTFLAGS=-O3 FOPTFLAGS=-O3 --download-mumps=1
--download-parmetis=1 --download-scalapack=1
--download-blacs=1<br>
[0]PETSC ERROR:
------------------------------------------------------------------------<br>
[0]PETSC ERROR: User provided function() line 0 in
unknown directory unknown file<br>
--------------------------------------------------------------------------<br>
MPI_ABORT was invoked on rank 0 in communicator
MPI_COMM_WORLD <br>
with errorcode 59.<br>
<br>
<br>
<pre cols="72">Yours sincerely,
TAY wee-beng</pre>
<br>
On 22/4/2012 2:53 PM, Matthew Knepley wrote:
<blockquote type="cite">On Sun, Apr 22, 2012 at 3:31
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_quote">
<blockquote class="gmail_quote" style="margin:0
0 0 .8ex;border-left:1px #ccc
solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> Hi,<br>
<br>
I am using petsc-dev 2012-04-20.<br>
<br>
Btw, I'm referring to :<br>
<br>
<a moz-do-not-send="true"
href="http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex22f.F.html"
target="_blank">http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex22f.F.html</a><br>
<br>
Part of the code is :<br>
<br>
<pre>call KSPSetFromOptions(ksp,ierr)
call KSPSetUp(ksp,ierr)
call KSPSolve(ksp,PETSC_NULL_OBJECT,PETSC_NULL_OBJECT,ierr)
call KSPDestroy(ksp,ierr)
call DMDestroy(da,ierr)
call PetscFinalize(ierr)</pre>
<br>
<br>
Unlike other codes like ex29c or ex45c,
there isn't a "call
KSPGetSolution(ksp,x,ierr)"<br>
</div>
</blockquote>
<div><br>
</div>
<div>You need to declare "Vec x", and then you
can call KSPGetSolution(ksp, x, ierr)</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0
0 0 .8ex;border-left:1px #ccc
solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> Also I
want to add "call VecView(x,ierr)" to print
out the results, which is usally added after
the above.<br>
<br>
Thank you<br>
<pre cols="72">Yours sincerely,
TAY wee-beng</pre>
<br>
On 22/4/2012 1:14 AM, Matthew Knepley wrote:
<blockquote type="cite">On Sat, Apr 21, 2012
at 6:31 PM, 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_quote">
<blockquote class="gmail_quote"
style="margin:0 0 0
.8ex;border-left:1px #ccc
solid;padding-left:1ex"> Hi,<br>
<br>
May I know if ex22f is complete? I
can't find :<br>
<br>
call KSPGetSolution(ksp,x,ierr)<br>
<br>
If I entered it, it says x not found.<br>
</blockquote>
<div><br>
</div>
<div>This is correct in petsc-dev. What
version are you using?</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote"
style="margin:0 0 0
.8ex;border-left:1px #ccc
solid;padding-left:1ex"> Thank you!<span><font
color="#888888"><span><font
color="#888888"><br>
<br>
-- <br>
Yours sincerely,<br>
<br>
TAY wee-beng<br>
<br>
</font></span></font></span></blockquote>
<span><font color="#888888"> </font></span></div>
<span><font color="#888888"> <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<br>
</font></span></blockquote>
</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<br>
</blockquote>
</div>
</blockquote>
</div>
</blockquote>
</div>
</blockquote>
</div>
</blockquote>
</body>
</html>