Matmul for large (“tall”) data

This is diagv = true*diagv + true*X*Xi' (using true==1)