Hierarchical GPs with varying lengthscales

[jessegrabowski]

@mike-lawrence

then computing the matrix-multiply more efficiently through knowledge that it’s a lower-tri and using the foward-substitution algorithm

Help wanted on this issue