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.
Running with LIBMESH_OPTIONS="-eps_type lapack" comes up with results; I presume that's switching to a direct solver. This probably isn't SLEPc's fault, but if you're doing eigen problems with complex-valued PDEs then you might want to ask them about this issue anyway; I don't know enough about our SLEPc interface or about the various methods they have available to advise you on which to use. For now I'm going to put "-eps_type lapack" in our ex17 Makefile just to get it passing tests; however if you learn more about what we might be doing wrong that's confusing their Arnoldi solver then I'd appreciate hearing about it. Thanks, --- Roy On Thu, 27 Oct 2011, Roy Stogner wrote: > 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
