Abstract: A method of data recovery of a remaining set of original data from a partial set of the original data. The original data is organized by a geometry having a locality function. Each original datum occupies a point of the geometry and has a corresponding distance to each other point of the geometry as defined by the locality function. The original datum further has corresponding values. The method includes, for each remaining datum of the remaining set: identifying a family of points of the geometry corresponding to members of the partial set local to the remaining datum by using the locality function; estimating the corresponding values of the remaining datum by using a peridynamic differential operator on the corresponding values of each of the members of the partial set in the family of points; and storing the estimated values of the remaining datum on a nonvolatile storage device.