implement new lat_cor algorithm: lat_cor_direct()