def METHOD_NAME(mesh, fs, degree, dumpfile):
    V = FunctionSpace(mesh, fs, degree)
    f = Function(V, name="f")
    x = SpatialCoordinate(mesh)
    f.interpolate(x[0]*x[1])
    f2 = Function(V, name="f")
    dumpfile = mesh.comm.bcast(dumpfile, root=0)
    with DumbCheckpoint(dumpfile, mode=FILE_CREATE) as chk:
        chk.store(f)
        chk.load(f2)
    assert np.allclose(f.dat.data_ro, f2.dat.data_ro)