Code with 1 matrix and 1 vector runs (there were alignment problems with icol_t since I am
reading 16 elements at a time. Alignment must be 64 and no longer 16 floats.  64 bytes.
Need a transpose to generate solution and check it for 1/1 case. Results do not quite match yet
so there is a bug.
still not working 100%
Problems with method_8a_base. Can't get agreement.
I fixed the code, which now runs with vector instructions for the 1 matrix/1 vector case.
I get 14 Gflops/sec in the Reverse Cuthill-McKee case on 64^3.