def METHOD_NAME(output_dir):
    mesh = Mesh.from_file(data_dir + '/meshes/elements/2_4_1.mesh')
    domain = refine(FEDomain('domain', mesh), output_dir, 1)
    ok = compare_mesh('2_4', domain.mesh.coors, domain.mesh.get_conn('2_4'))
    assert ok