dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

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
>