def METHOD_NAME(unit_mesh_2d):
    fes = VectorH1(unit_mesh_2d, order=2)
    u, v = fes.TnT()
    a = BilinearForm(fes, symmetric=True)
    a += InnerProduct(Grad(u), Grad(v)).Compile(realcompile=True, wait=True, keep_files=True) * dx
    a.Assemble()