powered by
get penalty function
get_penalty(y, X, locs, covfun_name)
response
design matrix
locations
string name of covariance function