NonlinearDiscreteProblem should use system assembler if symmetric==True

Issue #225 resolved
Jan Blechta created an issue

Result of the long discussion is the revelation of this issue:

NonlinearVariationalSolver::NonlinearDiscreteProblem::J/F uses normal Assembler breaking a symmetry even if parameter NonlinearVariationalSolver.parameters['symmetric'] == true. In fact, this parameter is not used anywhere.

Comments (9)

  1. Jan Blechta reporter

    Option 'symmetric' yet to be removed.

    Not necessarily. Cholesky/MINRES can be invoked rather than LU/GMRES.

  2. Log in to comment