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.
On Thursday, January 30, 2020 at 10:39:27 AM UTC-5, Wolfgang Bangerth wrote: > > On 1/29/20 1:14 PM, Andrew Davis wrote: > > * > > * > > *For some reason the feValues.get_function_values(old_solution, > > old_solution_values); call in the CellIntegrator function always sets > the > > old_solution_values to zero. Does anyone see what I'm doing wrong?* > > So the output of this function is wrong. Is the input correct? My guess > would > be that if you looked at 'old_solution', you'd find that it's a zero > vector. > Why that is so I don't know, but it's something you should be able to > trace > back to its source. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > <javascript:> > 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/ae51cb34-3209-482a-afbd-c33922b2d78c%40googlegroups.com.
