Produce a thin-plane spline matrix based on a given location matrix
thinPlateSplineMatrix(location)
A thin-plane spline matrix
pesudo_sequence <- seq(-5, 5, length = 5)
two_dim_location <- as.matrix(expand.grid(x = pesudo_sequence, y = pesudo_sequence))
thin_plate_matrix <- thinPlateSplineMatrix(two_dim_location)