Kyusik, 2016-09-23 9:11 GMT-04:00 <hanks0...@gmail.com>: > 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? A gradient is a tensor because it has dim components. In 2D, you have the derivatives with respect to x and y -> you have two components. The solution is only a scalar so you get only a double. So if you need to use the solution instead of the gradient, the type will be different but that's not a problem. The logic of the code should be the same.
> And, In step-6.cc, There is no comments about Lagrange finite elements in > this code. step-6 uses FE_Q which are Lagrange finite elements http://dealii.org/developer/doxygen/deal.II/classFE__Q.html So what you did is correct. Best, Bruno -- 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.
