compute covariance matrix more efficiently
Merge pull request #1440 from votjakovr/feature/gp_refactoring

Compute covariance matrix more efficiently