Re: [Libmesh-users] solution and old_local_solution

[Tim K. <tim...@ce...>]

Dear myself,

On Thu, 13 Nov 2008, Tim Kroeger wrote:

> *system.solution += *system.old_local_solution;
> system.update();
>
> Unfortunately, this doesn't work in parallel since solution is a
> parallel vector and old_local_solution is a serial vector.

I have now found the following solution that seems to work fine for 
me:


   for(unsigned int i=system.solution->first_local_index(); i<system.solution->last_local_index(); i++)
     {
       system.solution->add(i,(*system.old_local_solution)(i));
     }
   system.solution->close();
   system.update();


I'm not quite sure whether this is efficient, but it doesn't seem to 
be a bottleneck at the moment.

Best Regards,

Tim

-- 
Dr. Tim Kroeger                                        Phone +49-421-218-7710
tim...@me..., tim...@ce...  Fax   +49-421-218-4236

MeVis Research GmbH, Universitaetsallee 29, 28359 Bremen, Germany

Amtsgericht Bremen HRB 16222
Geschaeftsfuehrer: Prof. Dr. H.-O. Peitgen