Dear all,
For a transient linear system in my application, the equation is quite
complicted, and it seemed easier to me to solve for (u_{n+1}-u_n) than
for u_{n+1} (that is, I solve for the update rather than for the
solution). To get the new solution, I then do:
*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. Of course,
*system.local_solution += *system.old_local_solution;
system.update();
would run, but do the wrong thing since update() copies things in the
wrong direction.
Is there any nice solution to this problem?
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
