Package: libdeal.ii6.2.0 Version: 6.2.0-3 Something in the Trilinos Sacado linkage isn't working. Although step-31 works fine with (block) sparse matrices and the AMG preconditioner, step-33 fails as follows:
% make run
============================ Remaking Makefile.dep
==============optimized===== step-33.cc
============================ Linking step-33
mpicxx -o step-33 step-33.o /usr/lib/libdeal_II_2d.so /usr/lib/liblac.so
/usr/lib/libpetscall.so /usr/lib/slepc/linux-gnu-c-opt/lib/libslepc.so
/usr/lib/libmetis.so /usr/lib/libbase.so
/usr/lib/libtrilinos_stratimikosamesos.so
/usr/lib/libtrilinos_stratimikosaztecoo.so
/usr/lib/libtrilinos_stratimikosifpack.so /usr/lib/libtrilinos_stratimikosml.so
/usr/lib/libtrilinos_stratimikos.so /usr/lib/libtrilinos_ml.so
/usr/lib/libtrilinos_amesos.so /usr/lib/libtrilinos_aztecoo.so
/usr/lib/libtrilinos_belos.so /usr/lib/libtrilinos_ifpack.so
/usr/lib/libtrilinos_rtop.so /usr/lib/libtrilinos_sacado.so
/usr/lib/libtrilinos_thyra.so /usr/lib/libtrilinos_thyraepetra.so
/usr/lib/libtrilinos_thyraepetraext.so /usr/lib/libtrilinos_epetraext.so
/usr/lib/libtrilinos_epetra.so /usr/lib/libtrilinos_teuchos.so
/usr/lib/libtrilinos_triutils.so -llapack -lnetcdf_c++ -lnetcdf -lz -lblas
-lgfortran -lamd -lumfpack -L/usr/lib -lpthread -pthread -rdynamic
-L/usr/lib -L/usr/lib -lmetis -lparpack -lexpat -ltrilinos_galeri
============================ Running step-33
T=0
Number of active cells: 1633
Number of degrees of freedom: 6924
NonLin Res Lin Iter Lin Res
_____________________________________
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
Aztec status AZ_breakdown: numerical breakdown
nan 0000 nan
----------------------------------------------------
Exception on processing:
--------------------------------------------------------
An error occurred in line <3042> of file <step-33.cc> in function
void ConservationLaw<dim>::run() [with int dim = 2]
The violated condition was:
nonlin_iter <= 10
The name and call sequence of the exception was:
ExcMessage ("No convergence in nonlinear solver")
Additional Information:
No convergence in nonlinear solver
--------------------------------------------------------
Aborting!
----------------------------------------------------
make: *** [run] Error 1
Something seems to be wrong in the finite differencing, but it will take
some effort to track the error down. Upstream reports that stock
deal.II and Trilinos works fine for them.
-Adam
--
GPG fingerprint: D54D 1AEE B11C CE9B A02B C5DD 526F 01E8 564E E4B6
Engineering consulting with open source tools
http://www.opennovation.com/
signature.asc
Description: This is a digitally signed message part
