Re: [deal.II] Re: hp convergence for non-linear solver.

Thanks for replying me ! I moved to squared domain to avoid other possible error when I imply boundary condition, but I still have same manufactured solution e^(-x^2-y^2) I have to say that I have some changed error curve after moving to new mesh domain. Which is slightly better, but still uns

Re: [deal.II] Re: hp convergence for non-linear solver.

On 11/02/2016 07:17 PM, Jaekwang Kim wrote: I'd like check the visualized error per cell. Is there any module in deal.ii that enables this...? Yes, you just need to create a vector with as many entries as there are cell (e.g., the output of VectorTools::integrate_difference) and attach that

Re: [deal.II] Re: hp convergence for non-linear solver.

Thank you for replying me! I'd like check the visualized error per cell. Is there any module in deal.ii that enables this...? 2016년 11월 2일 수요일 오전 7시 27분 56초 UTC-5, Wolfgang Bangerth 님의 말: > > On 11/01/2016 11:06 PM, Jaekwang Kim wrote: > > > > do you have any other idea? > > > > Output the

Re: [deal.II] Re: hp convergence for non-linear solver.

On 11/01/2016 11:06 PM, Jaekwang Kim wrote: do you have any other idea? Output the errors per cell and visualize. Is there a pattern? Is the error localized on the boundary, for example? Best W. -- Wolfgang Bangerth

[deal.II] Re: hp convergence for non-linear solver.

Thank for the reply! I just fixed my code according to your suggestion. Now, I am evaluating my error with following code lines. template void msurface::evaluate_error() { Vector difference_per_cell (triangulation.n_active_cells()); *const MappingQ mapping (degree);*

[deal.II] Re: hp convergence for non-linear solver.

Thanks for replying! Could you explain more? I considered mapping degree only when I assemble system... for example. template void msurface::assemble_system () { const MappingQ mapping (degree); const QGauss quadrature_formula(degree+2); FEValues fe_values (mapping, fe, quadratu

[deal.II] Re: hp convergence for non-linear solver.

Jaekwang, [...] > > VectorTools::integrate_difference > (dof_handler,solution,Solution(), > >difference_per_cell > *,QGauss(degree+**2**)*,VectorTools::L2_norm); > somegthing that immediately comes to mind is that you don't use a non-default Mapping i

[deal.II] Re: hp convergence for non-linear solver.

Jaekwang, On Tuesday, November 1, 2016 at 5:43:36 PM UTC-4, Jaekwang Kim wrote: > > I think that I didn't make any mistake on my code lines because at least > my numerical solution is converging for the case p=1. > > but I couldn't get such convergence for the case p=2, and 3. I have tried > to

[deal.II] Re: hp convergence for non-linear solver.

Thanks for replying! Here I attached solution for p=2, and it seems that it has right share. When you said 'right quadrature' dose that meanquadrature for