First of all, Thank you very much for your reply And, I'm sorry, I should have been more specific about my problem....
In the lecture 31.65, the professor change the step-6.cc a little to apply Picard iteration. Doing this, he uses previous_solution_gradients[q_index]. So, like this, I just used previous_solution[q_index] But, he also say that previous_solution_gradients is tensor(it is declared as std::vector<Tensor<1,dim>>) But previous_solution is not tensor it is declared as vector<double> So, my question is that it is OK, If I use previous_solution as like previous_solution_gradients is used, even if two types are not same? And, In step-6.cc, There is no comments about Lagrange finite elements in this code. So, according to what you mention, the way I did was wrong, What should I have to do to get maximum values of previous_solution? Best, Kyusik -- 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.