Thanks Guido,
You were correct. I had not spotted this as the declaration of the
SchurComplement
SchurComplement schur_complement (system_matrix, m_inverse);
does nothing until schur_complement is used in the call to cg.solve,
when it attempts to find an inverse.
Many thanks to you and other respondents.
John
On 21/12/11 06:47, Guido Kanschat wrote:
John, these InverseMatrix objects, do they use CG to approximate the
action of the inverse? If so, that's probably where your exception
gets thrown.
Best,
Guido
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
