Hi Wolfgang, >> For me it runs okay until >> output_results(): we are setting values in the ghosted vector >> joint_solution there, which we are no longer allowed to do >> (ExcGhostsPresent). I am too tired to fix this now, though. :-) > > Aw, indeed. I admit not even having an idea how to work around that. What do > you think?
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. -- Timo Heister http://num.math.uni-goettingen.de/~heister _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
