I am writing my results as scalar results, using std::vector<std::string> solution_names = std::vector<std::string>{"A_x" , "A_y", "B_x", "B_y", "C_x", "C_y"};
std::vector<DataComponentInterpretation::DataComponentInterpretation> data_component_interpretation(n_components * dim, DataComponentInterpretation::component_is_scalar); LinearAlgebraTrilinos::MPI::Vector distributed_solution(old_solution); data_out.add_data_vector (distributed_solution, solution_names, DataOut <dim>::type_dof_data, data_component_interpretation); //And so on Now I would expect that A_x and A_y will look exactly the same, after I use the same equation. But according to my viewer (ParaView/Visit) A_x is different compared to A_y (and the same goes for B and C). Partially I get wild fluctuations in the result for A_x, while the result for A_y is smooth. Why does that happen? Thanks! -- 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.