def fom(reco, true_image):
    """Sobolev type FoM enforcing both gradient and absolute similarity."""
    gradient = odl.Gradient(reco.space)
    return (gradient(reco - true_image).norm() +
            reco.space.dist(reco, true_image))