We present the reasonableness of the extension of a two-level domain decomposition method to a multilevel method as a preconditioner for interpolation with radial basis functions (RBF) on distributed memory systems. The arising subproblems are efficiently solved using the FGP algorithm, a method that is well-suited for shared memory settings.