On 1/30/20 8:45 AM, Andrew Davis wrote: > I thought the same thing---that 'old_solution' would be a zero vector---but > when I add the line `std::cout << old_solution << std::endl;` in the function > IntegrateCell<dim>::operator() (just after the `get_function_values' call in > the original post) prints a non-zero vector that I expect.
That does not sound credible to me. Can you construct a small example that does exactly that? I.e., print out, for example, the l2-norm of old_solution right before the call to fe_values.get_function_values(), and the l2-norm of old_solution_values right after. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/34bf95a7-7fca-c04f-6f7c-b97deb97e1bf%40colostate.edu.