On 19. okt. 2010 14:43, Anders Logg wrote: > Confirmed. Can you add a bug report for this and I will try to fix it. > >
Yes and great! -- Marie > -- > Anders > > > On Thu, Oct 07, 2010 at 10:54:41AM +0200, Marie Rognes wrote: > >> >> Something slightly peculiar seems to happen with 3d refinement. Refining >> a 1 x 1 x 1 Unit Cube 3 times gives largest cell diameter that is >> increasing (but smallest cell diameter that is decreasing as expected). >> It does not seem to (only) be a bug in hmax()/hmin(), odd convergence >> rates are also observed. Any ideas? >> >> Example: >> >> from dolfin import * >> >> mesh = UnitCube(1, 1, 1) >> >> for i in range(3): >> mesh = refine(mesh) >> >> print "Refinement level = ", i >> print "h_max = ", mesh.hmax() >> print "h_min = ", mesh.hmin() >> >> gives: >> >> Refinement level = 0 >> h_max = 1.65831239518 >> h_min = 0.866025403784 >> Refinement level = 1 >> h_max = 1.47901994577 >> h_min = 0.433012701892 >> Refinement level = 2 >> h_max = 1.7809758561 >> h_min = 0.216506350946 >> >> > -- > Anders > _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp

