Data imputation with mtsdi package R

I found the package mtsdi and the corresponding function mnimput: R Documentation This function gives very good results for my data sets. However, I would like to understand the mathematical background of the function. I am familiar with the EM algorithm but how exactly does the function work? How are splines used here and is there a connection with the algorithm derived by Schafer?