Overall, the algorithm described makes sense, but crucial details are missing and no attempts to replicate this paper succeeded so far. For instance, paragraph  says "The numerical implementation involves a series of intrinsic concurrent interpolations and solutions of tridiagonal systems of linear equations" with no further information given. What kind of interpolations? How is this tridiagonal system constructed in 3D? Such kind of instruction, as well as data, is missing.