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