fortran77 	: row-major	= contiguous row, separate column
c		: column-major

a(m,n) = m rows by n columns

a11 a12 a13 ... a1n
a21 a22 a23 ... a2n
a31 a32 a33 ... a3n
... ... ... ... ...
... ... ... ... ...
... ... ... ... ...
am1 am2 am3 ... amn


ds = B dx

Br1x1 Br1x2 Br1x3 ..... Br1xn
Br2x1 Br2x2 Br2x3 ..... Br2xn

B    ( NR X NC ) LDB=NR
B^-1 ( NC X NR ) LDB=NC
