Dr, Bangerth,
>
> u_h(x) = \sum_i U_i \varphi_i(x)
>
>
First of all, I'm really sorry to bother you continuously...
Could I have please one more question....?
I also wanted to calculate l2 norm of |u_h(x)-u_h_exact(x)| (where
u_h_exact(x) is exact solution of u_h(x)).
So what I did is following...
1. ExactSolution<dim>() is defined like this...
template <int dim>
double ExactSolution<dim>::value (const Point<dim> &p,
const unsigned int /*component*/) const
{
double exact;
exact =
psi0*p(0)*p(0)/(R0*R0*R0*R0)*(2.*R0*R0-p(0)*p(0)-4.*alp*alp*p(1)*p(1));
return exact;
}
2.I used project library
VectorTools::project(dof_handler,
constraints,
QGauss<dim>(2),
ExactSolution<dim>(),
exact_solution);
3. I calculated l2 norm
Difference = solution;
Difference -= exact_solution;
norm = Difference.l2norm();
But, according to what you mentioned, I calculated l2norm of U_i-U_i_exact
not l2norm of u_h(x)-u_h_exact(x)
Am I right....?
Thank you.
Kyusik.
--
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.
