def METHOD_NAME():
    filename = os.path.join(mesh_dir, 'collela_2d.h5')
    solution = (cos(pi*x/2)*cos(pi*y/2))**2
    f        = laplace(laplace(solution))
    dir_zero_boundary    = get_boundaries(1, 2, 3, 4)
    dir_nonzero_boundary = get_boundaries()
    l2_error, h1_error, h2_error = run_biharmonic_2d_dir(filename, solution, f,
            dir_zero_boundary, dir_nonzero_boundary)
    expected_l2_error = 0.10977627980052021
    expected_h1_error = 0.32254511059711766
    expected_h2_error = 1.87205519824758
    assert( abs(l2_error - expected_l2_error) < 1.e-7)
    assert( abs(h1_error - expected_h1_error) < 1.e-7)
    assert( abs(h2_error - expected_h2_error) < 1.e-7)