Re: [Libmesh-devel] Matrix Free Memory Scaling with ParallelMesh

[Derek G. <fri...@gm...>]

On Jun 24, 2008, at 1:19 PM, Roy Stogner wrote:

> I've got a good guess as to where the memory spikes are occuring:
> check out __libmesh_petsc_snes_residual() in petsc_nonlinear_solver.C.
> Is that X_local a serial vector?

Good catch Roy... that's along the lines of what I was thinking.

> I think I noticed this problem and tried to avoid it in
> petsc_diff_solver.C; you might cut and paste some of that code into
> Ben's solver to see if using my "swap vectors, System::update"
> localization works any better.

I will try that.  Thanks!

Derek