The error message when `to_vec` of `y` and `y\bar` are not the same length is quite cryptic ```julia DimensionMismatch("second dimension of A, 4, does not match length of x, 7") ``` we could catch it and give a better error