Mulgen

as per #54 

I made it perform a transpose when given a ConjTrans because DGEMM's docs have `TRANSA = 'C' or 'c',  op( A ) = A'. `.  The conjugate transpose is defined as the element-wise complex conjugate of the transpose, so in the case of real numbers, it is the same as the transpose, so the behavior is implemented as defined.  It sounds like panicing would be a special case where none is needed.

I don't have an issue with putting constants into mat64 to make things easier on users.