This function constructs basis functions. It assumes coordinates form a metric.
MakeBasis(q = q, model.input)A list with
K spatial basis functions
KtQK which is literally the matrix operation transpose(K) times Q times K, and is the variance of the multivariate #' random effect gamma.
The desired number of basis functions. Must be an integer greater than or equal to 1.
A list of output created by running the create.data.R function