[petsc-users] Questions abt ex22f

[TAY wee-beng]

Hi,

I added "Vec x,b" after "KSP ksp"
and then "call KSPGetSolution(ksp, x, ierr)"

I wanted to see the output so I added "call VecView(x,ierr)" but I got 
this error:

[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, 
probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see 
http://www.mcs.anl.gov/petsc/petsc-as/documentation/faq.html#valgrind[0]PETSC 
ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to 
find memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, 
and run
[0]PETSC ERROR: to get more information on the crash.
[0]PETSC ERROR: --------------------- Error Message 
------------------------------------
[0]PETSC ERROR: Signal received!
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 5, Sat Oct 29 
13:45:54 CDT 2011
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: ./ex22f on a arch-linu named hpc12 by wtay Sun Apr 22 
21:02:14 2012
[0]PETSC ERROR: Libraries linked from 
/home/wtay/Lib/petsc-3.2-p5_mumps_rel/lib
[0]PETSC ERROR: Configure run at Sun Nov 27 15:18:15 2011
[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
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: User provided function() line 0 in unknown directory 
unknown file
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 59.


Yours sincerely,

TAY wee-beng


On 22/4/2012 2:53 PM, Matthew Knepley wrote:
> On Sun, Apr 22, 2012 at 3:31 AM, TAY wee-beng <zonexo at gmail.com 
> <mailto:zonexo at gmail.com>> wrote:
>
>     Hi,
>
>     I am using petsc-dev 2012-04-20.
>
>     Btw, I'm referring to :
>
>     http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex22f.F.html
>
>     Part of the code is :
>
>     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)
>
>
>
>     Unlike other codes like ex29c or ex45c, there isn't a "call
>     KSPGetSolution(ksp,x,ierr)"
>
>
> You need to declare "Vec x", and then you can call KSPGetSolution(ksp, 
> x, ierr)
>
>    Matt
>
>     Also I want to add "call VecView(x,ierr)" to print out the
>     results, which is usally added after the above.
>
>     Thank you
>
>     Yours sincerely,
>
>     TAY wee-beng
>
>
>     On 22/4/2012 1:14 AM, Matthew Knepley wrote:
>>     On Sat, Apr 21, 2012 at 6:31 PM, TAY wee-beng <zonexo at gmail.com
>>     <mailto:zonexo at gmail.com>> wrote:
>>
>>         Hi,
>>
>>         May I know if ex22f is complete? I can't find :
>>
>>         call KSPGetSolution(ksp,x,ierr)
>>
>>         If I entered it, it says x not found.
>>
>>
>>     This is correct in petsc-dev. What version are you using?
>>
>>       Thanks,
>>
>>         Matt
>>
>>         Thank you!
>>
>>         -- 
>>         Yours sincerely,
>>
>>         TAY wee-beng
>>
>>
>>
>>
>>     -- 
>>     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
>
>
>
>
> -- 
> 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-users/attachments/20120422/dc9d912b/attachment.htm>