> "Exception on processing:
> Iterative method reported convergence failure in step 1 with residual nan"
I haven't looked at your code, but it strikes me that you get a NaN after only
one iteration. Have you done something like
std::cout << system_matrix.frobenius_norm() << " "
<< system_rhs.l2_norm();
to verify that your matrix/rhs aren't already containing NaNs before you even
hand them over to the solver?
Best
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
