On Wed, 23 Feb 2011, Roman Vetter wrote:
> - The PETSc linear solver provides a print_converged_reason() member,
> which includes the KSPConvergedReason enums, but not all of them! It
> seems that _ksp is not destroyed right after the solve(), which allows
> that function to fetch the convergenc
On Wed, 23 Feb 2011, Roman Vetter wrote:
> 1. Replace the contents of PetscLinearSolver::print_converged_reason() by
>
> KSPConvergedReason reason;
> KSPGetConvergedReason(_ksp, &reason);
> libMesh::out << "Linear solver convergence/divergence reason: " <<
> KSPConvergedReasons[reason] << std::e
On Mon, 28 Feb 2011, Johannes Huber wrote:
> sorry, I tried to submit a patch file as descriped on the libMesh-Homepage,
> but it didn't work. What exactly do I have to type in? svn ci didn't do the
> job, I guess, I don't have the permission to add a file (which makes sense).
Right - we give an
Copying this to libmesh-devel and libmesh-users; I don't think I see
the answer (or even any productive way to debug such a problem)
myself, I'm afraid...
On Mon, 28 Feb 2011, Andreas Günther wrote:
the Libmesh problem I observe reminds me of the problem described in
http://www.mail-archive.c
On Thu, 24 Feb 2011, Roman Vetter wrote:
>>> There is a third possibility, which is a compromise:
>>> In the case of higher order variables, one could write each component
>>> separately as scalar variables, perhaps with a suffix in the name
>>> denoting the component. This way, all information
