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
[EMAIL PROTECTED], [EMAIL PROTECTED] Fax +49-421-218-4236
MeVis Research GmbH, Universitaetsallee 29, 28359 Bremen, Germany
Amtsgericht Bremen HRB 16222
Geschaeftsfuehrer: Prof. Dr. H.-O. Peitgen
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
