Hi,
Sorry for the very late response...
I attached the slepc_variables and slepcconf files, I couldn't find
slepc_common_variables. Where exactly should that be?
Even using Arpack, I still have issues. I have a generalized condensed
eigenvalue problem, and I'm interested in the smallest eigenvalue(s) and
-vector(s). It's a 3D 3-variable problem with n_nodes()=576, even though
I would like to refine the mesh further. The eigensolver is not
converging or converging very slowly, depending on the parameters, for
example:
equation_systems.parameters.set <unsigned int> ("eigenpairs") = 3;
equation_systems.parameters.set <unsigned int> ("basis vectors") = 200;
equation_systems.parameters.set<unsigned int> ("linear solver maximum
iterations") = 400;
On the other hand, if I print_matlab the matrices and solve there the
problem using eigs, I get the correct solution in a few seconds. But
shouldn't the solver be the same?
Thanks for helping
Lorenzo
On 4/30/14 9:15 PM, John Peterson wrote:
On Wed, Apr 30, 2014 at 12:04 PM, Lorenzo Zanon
<[email protected] <mailto:[email protected]>> wrote:
Hello,
Now I should have slepc correctly installed with arpack. But still it
looks like I cannot use arpack in libmesh, since in the file
slepc_eigen_solver.C the arpack lines are commented out:
// case ARPACK:
// ierr = EPSSetType (_eps, (char*) EPSARPACK);
LIBMESH_CHKERRABORT(ierr); return;
Is there maybe something else I should add when libmesh is configured,
other than --enable-slepc?
This kind of looks like something we commented out because we didn't
have a good way of configuring for the situation where slepc has
arpack vs. doesn't have arpack.
Since you have it, you should be safe to uncomment this case and
recompile libmesh in order to start using Arpack with Slepc.
By the way, would you mind sharing
your slepcconf/slepc_common_variables/slepc_variables file with Arpack
enabled? There is probably something we can look for there to see if
slepc has Arpack enabled and handle this case automatically...
--
John
#
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# SLEPc - Scalable Library for Eigenvalue Problem Computations
# Copyright (c) 2002-2013, Universitat Politecnica de Valencia, Spain
#
# This file is part of SLEPc.
#
# SLEPc is free software: you can redistribute it and/or modify it under the
# terms of version 3 of the GNU Lesser General Public License as published by
# the Free Software Foundation.
#
# SLEPc is distributed in the hope that it will be useful, but WITHOUT ANY
# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
# FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
# more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with SLEPc. If not, see <http://www.gnu.org/licenses/>.
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
#
include ${PETSC_DIR}/conf/variables
include ${SLEPC_DIR}/${PETSC_ARCH}/conf/slepcvariables
SLEPC_LIB_DIR = ${SLEPC_DIR}/${PETSC_ARCH}/lib
SLEPC_INCLUDE = -I${SLEPC_DIR} -I${SLEPC_DIR}/${PETSC_ARCH}/include
-I${SLEPC_DIR}/include
CC_INCLUDES = ${PETSC_INCLUDE} ${SLEPC_INCLUDE}
FC_INCLUDES = ${PETSC_INCLUDE} ${SLEPC_INCLUDE}
SLEPC_EXTERNAL_LIB = ${ARPACK_LIB} ${BLZPACK_LIB} ${TRLAN_LIB} ${PRIMME_LIB}
${BLOPEX_LIB} ${FEAST_LIB}
INSTALL_LIB_DIR = ${SLEPC_LIB_DIR}
CCPPFLAGS = ${PETSC_CCPPFLAGS} ${SLEPC_INCLUDE}
FCPPFLAGS = ${PETSC_FCPPFLAGS} ${SLEPC_INCLUDE}
SLEPC_LIB = ${CC_LINKER_SLFLAG}${SLEPC_LIB_DIR} -L${SLEPC_LIB_DIR}
-lslepc ${SLEPC_EXTERNAL_LIB} ${PETSC_KSP_LIB}
# no longer required
SLEPC_FORTRAN_LIB =
SHLIBS = libslepc
LIBNAME = ${INSTALL_LIB_DIR}/${LIBBASE}.${AR_LIB_SUFFIX}
#if !defined(__SLEPCCONF_H)
#define __SLEPCCONF_H
#ifndef SLEPC_LIB_DIR
#define SLEPC_LIB_DIR "/usr/local/macports108/lib/petsc/lib"
#endif
#ifndef SLEPC_HAVE_ARPACK
#define SLEPC_HAVE_ARPACK 1
#define SLEPC_ARPACK_HAVE_UNDERSCORE 1
#endif
#endif
------------------------------------------------------------------------------
"Accelerate Dev Cycles with Automated Cross-Browser Testing - For FREE
Instantly run your Selenium tests across 300+ browser/OS combos.
Get unparalleled scalability from the best Selenium testing platform available
Simple to use. Nothing to install. Get started now for free."
http://p.sf.net/sfu/SauceLabs
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
