[Trilinos-Users] aztecOO thyra trouble
Marco Morandini
morandini at aero.polimi.it
Tue Apr 3 01:42:58 MDT 2007
I'm having trouble running
test_single_aztecoo_thyra_solver_driver.exe
with more than 1 process (on a single cpu machine,
but this should not make any difference, right?).
Is it supposed to work in parallel?
I'm the only one seeing this?
I've partially tracked down the problem:
after the call
// Here we will be careful not to reset the forward operator in
fears that
// it will blow out the internally created stuff.
aztecFwdSolver->SetUserOperator(const_cast<Epetra_Operator*>(&*aztec_epetra_epetraFwdOp));
in Thyra_AztecOOLinearOpWithSolveFactory.cpp
the code enters
int AztecOO::SetUserOperator(Epetra_Operator * UserOperator) {
in AztecOO.cpp
and hangs inside the call to UserOperator->NormInf():
// Aztec needs upper bound for matrix norm if doing polynomial
preconditioning
if (UserOperator->HasNormInf()) {
AZ_set_MATFREE_matrix_norm(Amat_, UserOperator->NormInf());
I'm lost after that.
I was not able to get a meaningful backtrace.
The command line is
mpirun -np 2 test_single_aztecoo_thyra_solver_driver.exe
--matrix-file=KheadSorted.mtx
The mpi implementation is openmpi-1.2; the same happens
with mpich.
Operating system: Suse 10.2,
gcc version 4.1.2 20061115 (prerelease) (SUSE Linux)
Trilinos 7.0.5, configured with
../configure --with-gnumake --enable-mpi --with-mpi-compilers
--enable-teuchos --enable-teuchos-extended --enable-teuchos-complex
--enable-teuchos-expat --enable-thyra --enable-epetra --enable-galeri
--enable-moertel --enable-stratimikos --enable-amesos
--enable-amesos-umfpack --enable-ifpack --enable-aztecoo
--enable-aztecoo-azlu --enable-triutils --enable-nox --enable-nox-epetra
--enable-nox-lapack --enable-new_package
--prefix=/home/marco/local/trilinos-7.0.5 --enable-expat
--enable-epetraext --enable-ml --enable-anasazi --enable-loca-anasazi
--disable-thyra-tests --disable-thyra-examples --enable-nox-examples
--enable-rtop --with-trilinos
Thanks in advance,
Marco
More information about the Trilinos-Users
mailing list