> One option is to write into a temporary vector without ghost elements > and skipping values that are not locally_owned. Before writing out the > solution, one assigns the temporary vector into (the ghosted) > joint_solution. I have no time to do work on that over the next days, > though.
I've just implemented this and it should now work. There was another bug downstream: we were taking the linfty norm of the old temperature at one point, but we have disallowed that recently for ghosted vectors. But I think that this particular operation is actually safe for ghosted operations, unlike the other norms, so I've replaced the assertion by a comment. If I'm missing something obvious here, please let me know. With this, the program should now run just fine. Best W. -- ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
