C/ZGEMMT changes the inputs in the UNC case

[grisuthedragon]

The Z/CGEMMT routine changes the second input when calling

cgemmt_("U", "N", "C", &n, &k, &alpha, A, &lda, B, &ldb, &beta, C, &ldc);

or

zgemmt_("U", "N", "C", &n, &k, &alpha, A, &lda, B, &ldb, &beta, C, &ldc);

After the call B is its conjugate.

In the ZGEMMT case this leads:

Matrix A:
(1.000000, 2.000000) (2.000000, 3.000000) (3.000000, 4.000000) 
(4.000000, 5.000000) (5.000000, 6.000000) (6.000000, 7.000000) 
(7.000000, 8.000000) (8.000000, 9.000000) (9.000000, 10.000000) 
Matrix B:
(1.000000, 1.000000) (2.000000, 2.000000) (3.000000, 3.000000) 
(4.000000, 4.000000) (5.000000, 5.000000) (6.000000, 6.000000) 
(7.000000, 7.000000) (8.000000, 8.000000) (9.000000, 9.000000) 
Matrix C:
(144.000000, 12.000000) (0.000000, 0.000000) (0.000000, 0.000000) 
(171.000000, 15.000000) (201.000000, 15.000000) (0.000000, 0.000000) 
(198.000000, 18.000000) (234.000000, 18.000000) (270.000000, 18.000000) 
After cgemmtMatrix A:
(1.000000, 2.000000) (2.000000, 3.000000) (3.000000, 4.000000) 
(4.000000, 5.000000) (5.000000, 6.000000) (6.000000, 7.000000) 
(7.000000, 8.000000) (8.000000, 9.000000) (9.000000, 10.000000) 
Matrix B:
(1.000000, -1.000000) (2.000000, -2.000000) (3.000000, -3.000000) 
(4.000000, -4.000000) (5.000000, -5.000000) (6.000000, -6.000000) 
(7.000000, -7.000000) (8.000000, -8.000000) (9.000000, -9.000000) 

cgemmt_example.c.txt

zgemmt_example.c.txt