* add function la_matrix_info that writes various information about a matrix to
  a stream
* use scaled quantities
* use SI units for input and output and convert them to scaled quantities
  during calculation
* enable scaling
* disable debugging output