I tried the procedure exactly like shown in step-40 and these are the results:
*DEAL.II:* <https://lh3.googleusercontent.com/-ZwWP5Xfoe6U/WJdmEMU4v8I/AAAAAAAAAFs/5hR3YA0RkeoYiSGZCDbEiHWH_I_pJJYWwCLcB/s1600/dealii_result.png> *C++ CODE:* <https://lh3.googleusercontent.com/-MnonNmQiobQ/WJdlxPUSPDI/AAAAAAAAAFk/aKBykIr6vmk98aCzEuUPzpa9IGybtLswQCLcB/s1600/cpp_result.png> I already tried some boundary interpolate_boundary_values approach, but I don't think it is something related to boundaries since I just have boundary conditions at the bottom and upper boundary. We have still bad results at the left and right boundary though. Any ideas what could cause such a behavior? I am really not an expert yet to dive deeply inside the deal.II parallel structure. Is it even correct to initialize an MPI vector within setup_system using the information before assembly and solution of the system, then fill this vector with data and pass it to the output_results() function? Kind regards, Seyed Ali -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.