feat: Default and overridable solve and root computation methods.

[daniel-dodd]

Post #370, akin to kernel / likelihood quadrature computation, would be good to have on the top level default methods for these defined on the Prior that say what solver and root method we'd like to have, defaulting to Cholesky-based methods.