make use of Cholesky functions in package dense
use Cholesky on cov matrix instead of general Solve