make covariance mapping to data vector idiot proof

[matroxel]

Store covariance as list of elements with theta1, theta2, z1, z2, d1, d2 with function to map to 2d array matching optionally with an enforced order of d.