[petsc-users] PETSC ERROR: KSPComputeEigenvalues()

[Oo]

Hi,


I met a problem when I tried to use "KSPComputeEigenvalues()".
 
The following is the solving part of my code.
When this programme run to the line 
" KSPComputeEigenvalues(ksp, numberEigen, r, c,neig);"
Then, PETSc Errors appear.


Do you have any suggestions?


////////////////////////////////////////////
    //Solving the Linear System Part======================================================     cout<<"solve the Linear system"<<endl;     t0=(double)clock();     KSP ksp;     PC pc;     KSPCreate(PETSC_COMM_WORLD, &ksp);    // KSPSetType(ksp, KSPCG);     KSPSetType(ksp, KSPGMRES);     KSPSetOperators(ksp, coeff_matrix, coeff_matrix, DIFFERENT_NONZERO_PATTERN); 
     PetscReal rtol=(1.e-4)*pow(h,4.0);     PetscReal stol=(1.e-3)*pow(h,4.0); 
     KSPSetTolerances(ksp, rtol, stol, PETSC_DEFAULT, PETSC_DEFAULT); 
     KSPGetPC(ksp, &pc);     PCSetType(pc, PCGAMG);     KSPSetFromOptions(ksp); 
 
     //====KSPSetComputeEigenvalues()     KSPSetComputeEigenvalues(ksp, PETSC_TRUE);     //KSPSetComputeEigenvalues(ksp, PETSC_FALSE);     //============================= 
 
     KSPSetUp(ksp); 
     KSPSolve(ksp, load_vect, m_solut); 
 
      cout<<"Before KSPComputeEignvalues()"<<endl;     //======KSPComputeEigenvalues()     numberEigen=1;//4*P_FEMMesh->P_Vertices->size();     KSPComputeEigenvalues(ksp,numberEigen,r,c,neig);     cout<<"After KSPComputeEignvalues()"<<endl; 
     cout<<"EignValues:"<<endl;     for(unsigned int j=0; j<numberEigen; j++)     {         cout<<r[numberEigen]<<" + "<<c[numberEigen]<<"i"<<endl;     }     //========================


Here is Petsc errors.when the programme runs to "KSPComputeEigenvalues(ksp,numberEigen,r,c,neig);"These errors appear.
PETSC ERRORs
[0]PETSC ERROR: --------------------- Error Message ------------------------------------
 
[0]PETSC ERROR: Null argument, when expecting valid pointer!
 
[0]PETSC ERROR: Null Pointer: Parameter # 3!
 
[0]PETSC ERROR: ------------------------------------------------------------------------
 
[0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 2013 
 
[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: /Users/wumeng/MyWork/BuildMSplineTools/bin/msplinePDE_PFEM_2 on a arch-darwin-c-debug named vis032b.sophia.inria.fr by wumeng Thu Apr 24 11:50:16 2014
 
[0]PETSC ERROR: Libraries linked from /Users/wumeng/MyWork/PETSc/arch-darwin-c-debug/lib
 
[0]PETSC ERROR: Configure run at Wed Dec 18 15:31:31 2013
 
[0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=gfortran --download-f-blas-lapack --download-mpich
 
[0]PETSC ERROR: ------------------------------------------------------------------------
 
[0]PETSC ERROR: KSPComputeEigenvalues() line 118 in /Users/wumeng/MyWork/PETSc/src/ksp/ksp/interface/itfunc.c
 
After KSPComputeEignvalues()
 
EignValues:
 
[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/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: likely location of problem given in stack below
 
[0]PETSC ERROR: ---------------------  Stack Frames ------------------------------------
 
[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,
 
[0]PETSC ERROR:       INSTEAD the line number of the start of the function
 
[0]PETSC ERROR:       is given.
 
[0]PETSC ERROR: [0] KSPComputeEigenvalues line 116 /Users/wumeng/MyWork/PETSc/src/ksp/ksp/interface/itfunc.c
 
[0]PETSC ERROR: --------------------- Error Message ------------------------------------
 
[0]PETSC ERROR: Signal received!
 
[0]PETSC ERROR: ------------------------------------------------------------------------
 
[0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 2013 
 
[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: /Users/wumeng/MyWork/BuildMSplineTools/bin/msplinePDE_PFEM_2 on a arch-darwin-c-debug named vis032b.sophia.inria.fr by wumeng Thu Apr 24 11:50:16 2014
 
[0]PETSC ERROR: Libraries linked from /Users/wumeng/MyWork/PETSc/arch-darwin-c-debug/lib
 
[0]PETSC ERROR: Configure run at Wed Dec 18 15:31:31 2013
 
[0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=gfortran --download-f-blas-lapack --download-mpich
 
[0]PETSC ERROR: ------------------------------------------------------------------------
 
[0]PETSC ERROR: User provided function() line 0 in unknown directory unknown file
 
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0
 
[unset]: aborting job:
 
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0










Thanks,




M.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140424/980c85ec/attachment-0001.html>