On 11/01/2011 03:28 PM, Roy Stogner wrote: > We're testing SLEPc + complex now, and seeing the same IPNorm error. > > The problem seems to be some issue with the default iterative solvers > we tell SLEPc to use; it looks like for some reason they're getting a > "norm" of their initial random starting vector that isn't a > non-negative real number. I'm not sure how this would happen, though > a Google search suggests that we're not the first people to run into > this IPNorm error when doing a generalized eigenproblem.
Interesting. So it's not a bug per se, good to know. > Running with LIBMESH_OPTIONS="-eps_type lapack" comes up with results; > I presume that's switching to a direct solver. The lapack solver uses the QR algorithm (or something similar) to find all the eigenvalues I believe. Its very robust, but only works for small problems since it treats the matrix as dense. Dave >> What options did you use to configure libMesh? >> >> I'll see if I can reproduce this. We currently have --enable-complex >> BuildBot tests and --enable-slepc tests but I don't think we're >> testing the combination of the two... >> --- >> Roy >> >> On Fri, 28 Oct 2011, Kyunghoon Lee wrote: >> >>> Hi all, >>> >>> I compiled petsc with complex variable support using the following options: >>> >>> ./configure --prefix=/Users/aeronova/ >>> Development/local/lib64/petsc/petsc-3.2-p4 --download-mpich=1 >>> --download-blacs=1 --download-parmetis=1 --download-scalapack=1 >>> --download-mumps=1 --with-scalar-type=complex --with-clanguage=cxx >>> >>> Then, rebuild slepc& libmesh (svn version). As I test libmesh with >>> examples, I got an error for ex17 (plz see the below). I'd appreciate if >>> someone could help me with this problem. (ex7, complex variable example, >>> seemed ok, though there were a couple of warnings such as WARNING! There >>> are >>> options you set that were not used!, WARNING! could be spelling mistake, >>> etc!, Option left: name:-f value: .5) >>> >>> >>> [0]PETSC ERROR: --------------------- Error Message >>> ------------------------------------ >>> [0]PETSC ERROR: IPNorm: The inner product is not well defined! >>> [0]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17 13:37:48 >>> 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: ./ex17-dbg on a darwin10. named spike.local by aeronova Fri >>> Oct 28 05:13:59 2011 >>> [0]PETSC ERROR: Libraries linked from >>> /Users/aeronova/Development/local/lib64/petsc/petsc-3.1-p8/lib >>> [0]PETSC ERROR: Configure run at Fri Oct 28 04:12:19 2011 >>> [0]PETSC ERROR: Configure options >>> --prefix=/Users/aeronova/Development/local/lib64/petsc/petsc-3.1-p8 >>> --download-mpich=1 --download-blacs=1 --download-parmetis=1 >>> --download-scalapack=1 --download-mumps=1 --with-scalar-type=complex >>> --with-clanguage=cxx >>> [0]PETSC ERROR: >>> ------------------------------------------------------------------------ >>> [0]PETSC ERROR: IPNorm() line 70 in src/ip/ipdot.c >>> [0]PETSC ERROR: IPOrthogonalize() line 338 in src/ip/iporthog.c >>> [0]PETSC ERROR: EPSGetStartVector() line 1487 in src/eps/interface/solve.c >>> [0]PETSC ERROR: EPSSolve_ARNOLDI() line 387 in >>> src/eps/impls/krylov/arnoldi/arnoldi.c >>> [0]PETSC ERROR: EPSSolve() line 146 in src/eps/interface/solve.c >>> [0]PETSC ERROR: User provided function() line 488 in >>> "unknowndirectory/"src/solvers/slepc_eigen_solver.C >>> application called MPI_Abort(comm=0x84000000, 1) - process 0[unset]: >>> aborting job: >>> application called MPI_Abort(comm=0x84000000, 1) - process 0 >>> make[2]: *** [run] Error 1 >>> make[1]: *** [run] Error 1 >>> make: *** [run_examples] Error 2 >>> >>> ------------------------------------------------------------------------------ >>> The demand for IT networking professionals continues to grow, and the >>> demand for specialized networking skills is growing even more rapidly. >>> Take a complimentary Learning@Cisco Self-Assessment and learn >>> about Cisco certifications, training, and career opportunities. >>> http://p.sf.net/sfu/cisco-dev2dev >>> _______________________________________________ >>> Libmesh-users mailing list >>> [email protected] >>> https://lists.sourceforge.net/lists/listinfo/libmesh-users >>> >>> > ------------------------------------------------------------------------------ > RSA® Conference 2012 > Save $700 by Nov 18 > Register now! > http://p.sf.net/sfu/rsa-sfdev2dev1 > _______________________________________________ > Libmesh-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------------ RSA® Conference 2012 Save $700 by Nov 18 Register now! http://p.sf.net/sfu/rsa-sfdev2dev1 _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
