Thin plate splines are only visually smooth if the order of the
wiggliness penalty, \(m\), satisfies \(2m > d+1\). If \(2m<d+1\) then this routine finds the smallest
\(m\) giving visual smoothness
for the given \(d\), otherwise the supplied \(m\) is used. The null space dimension is given by:
\(M=(m+d-1)!/(d!(m-1)!)\)
which is the value returned.