Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

I have an equation -$u''=1,~~x\in(0,1)$. I solved it numerically by Finite element method. and find an approximate solution. $u_{h}$. As you know to do this I defined bilinear and linear functionals for weak formulation s.t

a(u,v)=l(v) for every v in $H_{0}^{1}$ where

a(u,v)=$\int_{0}^{1}u'v'dx$ and l(v)=$\int_{0}^{1}vdx$.

But to find a functinal estimates, |J($u-u_{h}$ )|, I defined a functinal J(v)=$\int_{0}^{1}vdx$ and try to find z in dual problem a(v,z)=J(v) for every v in $H_{0}^{1}$

and I solved the dual problem then find approximation for z.

But J($u-u_{h}$)=a($u-u_{h}$,z)=a(u,z)-a($u_{h}$,z) =l(z)-a($u_{h}$,z)$~~~~~~~$ (**)

(**) is called the error indicator and my questions is that how can I find it ?

share|improve this question

1 Answer 1

You compute $z$ approximately as you described. Then plug this $z$ into the formula $J(u-u_h)=\ell(z)-a(u_h,z)$.

share|improve this answer
    
my quuestion is that. i ve just values of z at grid points not a explicit formula for z. so how can I put them into formula ? do I need to use gauss quadrature ? –  Brhn Aug 15 '12 at 17:34
    
Actually you can do it analytically. The coefficients of the expansion of $z$ in terms of the hat functions are given. So you just need to analytically compute a few integrals involving hat functions. –  timur Aug 15 '12 at 18:57

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.