[petsc-users] Error in Lapack xTREXC 1

[Jose E. Roman]

El 31/05/2013, a las 21:47, Nachiket Gokhale escribió:

> 
> 
> I am trying various options to solve a quadratic eigenvalue problem using SlepC.   I get the following error message 
> 
> Calling QEPSolve
> [0]PETSC ERROR: --------------------- Error Message ------------------------------------
> [0]PETSC ERROR: Error in external library!
> [0]PETSC ERROR: Error in Lapack xTREXC 1!
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Development GIT revision: ff331c249a8bbf694711b310f25fec0e839b33db  GIT Date: 2013-05-19 22:14:01 -0500
> [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: /home/gokhale/hyser.clone/src/forward on a linux-gcc-gpp-mumps named asd1.wai.com by gokhale Fri May 31 15:09:49 2013
> [0]PETSC ERROR: Libraries linked from /opt/petsc/petsc-git/linux-gcc-gpp-mumps/lib
> [0]PETSC ERROR: Configure run at Mon May 20 17:02:54 2013
> [0]PETSC ERROR: Configure options --download-blacs=1 --download-f-blas-lapack=1 --download-metis --download-mpich=yes --download-mumps=yes --download-parmetis=1 --download-scalapack=1 --with-cc=/usr/local/bin/gcc --with-clanguage=C++ --with-cmake=/usr/bin/cmake28 --with-cxx=/usr/local/bin/g++ --with-debugging=0 --with-fc=/usr/local/bin/gfortran --with-mpi=1 --with-shared-libraries=1 --with-x11=0 --with-x=0 PETSC_ARCH=linux-gcc-gpp-mumps
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: DSSort_NHEP_Total() line 459 in /opt/SlepC/slepc-git/src/ds/impls/nhep/dsnhep.c
> [0]PETSC ERROR: DSSort_NHEP() line 493 in /opt/SlepC/slepc-git/src/ds/impls/nhep/dsnhep.c
> [0]PETSC ERROR: DSSort() line 510 in /opt/SlepC/slepc-git/src/ds/interface/dsops.c
> [0]PETSC ERROR: QEPSolve_QArnoldi() line 259 in /opt/SlepC/slepc-git/src/qep/impls/qarnoldi/qarnoldi.c
> [0]PETSC ERROR: QEPSolve() line 107 in /opt/SlepC/slepc-git/src/qep/interface/qepsolve.c
> [0]PETSC ERROR: solvecomplex() line 161 in "unknowndirectory/"/home/gokhale/hyser.clone/lib/eqnlib/slepc_quad_system.cpp
> 
> 
> I am not sure what it means; Google/mailing list search didn't seem to  help.  My command line is 
> 
> ~/hyser.clone/src/forward -qep_nev 10 -qep_ncv 32 -qep_monitor_conv -qep_max_it 10000 -qep_smallest_imaginary -qep_type qarnoldi -qep_general -qep_linear_explicitmatrix -qep_eps_type krylovschur   -qep_eps_shift -0.1
> 
> If it helps, the matrices are small  - 60x60. 
> 
> -Nachiket
> 
> 

This is quite strange. Can you try with --download-f2cblaslapack instead of --download-f-blas-lapack? If it does not help, send the matrices to slepc-maint.

By the way, the -qep_eps_shift option does not exist.

Jose