method_0: 0.5 Gflops, method_4: 23 Gflops and Method_5 22 Gflops produce correct serial results
for compact case. .
Implementing method_7: 4 matrices + 4 vectors with tensor product.
method_7 produces a segmentation fault. Do not know why.