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.

Reply via email to