new interface:
Instead of -f one uses -A and -B as flags for input files corresponding
to matrices A and B. If both flags are set C=A*B is computed (if dimensions
are correct). If only A is given, C=A*A^T is computed.