use do concurrent in covariance calc
use cshift in cross_vec