Re: [Libmesh-users] Performance of EquationSystems::reinit() with ParallelMesh

2008-09-04 Thread Tim Kroeger
Dear Ben, On Thu, 4 Sep 2008, Benjamin Kirk wrote: > I noticed something very reminiscent of this just two days ago. > [...] > I can't quantify it at the moment but this took a lot longer than expected. > I just "seemed really long." It was somewhat acceptable in my case since I > then did a hun

Re: [Libmesh-users] Performance of EquationSystems::reinit() with ParallelMesh

2008-09-04 Thread Benjamin Kirk
I have just caught up to speed with the mailing list after a few distractions. I noticed something very reminiscent of this just two days ago. In my case I run a transient solution to steady-state and then stop the simulation. I then re-read this result, refine the mesh, project the solution,

Re: [Libmesh-users] Accessing DiffSolver::_solve_result

2008-09-04 Thread John Peterson
On Thu, Sep 4, 2008 at 10:46 AM, Roy Stogner <[EMAIL PROTECTED]> wrote: > > On Thu, 4 Sep 2008, David Knezevic wrote: > >> The patch is attached (I must've checked out libmesh anonomously cos I can't >> seem to commit myself, I'll look into that). > > There've been some server side SVN problems jus

Re: [Libmesh-users] Accessing DiffSolver::_solve_result

2008-09-04 Thread David Knezevic
Hi Roy, So far IIRC we've only been using that value internally to let the transient solvers and continuation solvers which call DiffSolver::solve back off their time step or continuation step size on failure. My thinking was that whether such reduction succeeded in making the problem well-beha

Re: [Libmesh-users] Accessing DiffSolver::_solve_result

2008-09-04 Thread Roy Stogner
On Thu, 4 Sep 2008, David Knezevic wrote: > The patch is attached (I must've checked out libmesh anonomously cos I can't > seem to commit myself, I'll look into that). There've been some server side SVN problems just recently. But try it again now; I just double-checked and got Tim's logging t

Re: [Libmesh-users] Performance of EquationSystems::reinit() with ParallelMesh

2008-09-04 Thread Roy Stogner
On Thu, 4 Sep 2008, Tim Kroeger wrote: > On Wed, 3 Sep 2008, Roy Stogner wrote: > > What about a better documentation of this? I have attached a patch for this. > >>> but I can't find a corresponding "serialize()" method. >> >> There is no reverse method to produce a parallel vector from a glob

Re: [Libmesh-users] Performance of EquationSystems::reinit() with ParallelMesh

2008-09-04 Thread Tim Kroeger
Dear Roy, On Wed, 3 Sep 2008, Roy Stogner wrote: As far as I understand, currently the NumericVector interface class (with its implementations such as PetscVector) is used as a parallel as well as a serial vector. Yes. I notice that there is a method localize() that seems to transform a se

Re: [Libmesh-users] Accessing DiffSolver::_solve_result

2008-09-04 Thread Roy Stogner
On Thu, 4 Sep 2008, David Knezevic wrote: > I'm using (a subclass) of FEMSystem, with NewtonSolver, and I'd like to > detect whether or not a quasi-Newton solve converged, i.e. I'd like to > detect whether NewtonSolver::solve returns the enum > DiffSolver::DIVERGED_MAX_NONLINEAR_ITERATIONS. Maybe

[Libmesh-users] Accessing DiffSolver::_solve_result

2008-09-04 Thread David Knezevic
I'm using (a subclass) of FEMSystem, with NewtonSolver, and I'd like to detect whether or not a quasi-Newton solve converged, i.e. I'd like to detect whether NewtonSolver::solve returns the enum DiffSolver::DIVERGED_MAX_NONLINEAR_ITERATIONS. Maybe I'm missing something, but I didn't see any eas