[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