A = size 2 by 2 Matrix of std::complex<double>:
[ (1,2) (5,6) ]
[ (3,4) (7,8) ]

B = size 2 by 2 Matrix of std::complex<double>:
[ (9,10) (13,14) ]
[ (11,12) (15,16) ]

A*B = size 2 by 2 Matrix of std::complex<double>:
[ (-28,154) (-36,210) ]
[ (-32,238) (-40,326) ]

transpose(A)*B = size 2 by 2 Matrix of std::complex<double>:
[ (-26,108) (-34,148) ]
[ (-34,276) (-42,380) ]

hermitian(A)*B = size 2 by 2 Matrix of std::complex<double>:
[ (110,-16) (150,-24) ]
[ (278,-8) (382,-16) ]

conjugate(A)*B = size 2 by 2 Matrix of std::complex<double>:
[ (156,-14) (212,-22) ]
[ (240,-10) (328,-18) ]

