Hi,
I was wondering what to do to run a convergence analysis when you don't
know the exact solution to your problem.
'VectorTools::integrate_difference' requires a class Function<dim> for
the exact solution.
So what I did, so far, is to define a nested class derived from
Function<dim>, in the class for my problem, with method 'value'
returning the solution of the current problem. I export that function
for my finest grid (convergence analysis in space) that I consider as my
exact solution, then when I run on coarser grids, I use that function in
'VectorTools::integrate_difference' as my exact solution.
Results are a bit weird, so far. Do you see something wrong in
principle? Any idea of what I should be doing instead?
Best,
Ben
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
