Conjugation of the right argument in tensor contraction

[DmitryLyakh]

Turns out this is due to wrong handling of conjugation of the second argument in a contraction. This only happens when your own implementation of BLAS is used, when the code calls a standard BLAS everything is fine. For conjugation on the first argument the result is OK, so this gives the following work-around (red = old code that shows the error, green = work-around that gives correct results). Hope this gives a clue for you to fix.

• ierr=talsh_tensor_contract("B(p,q)+=A(i,q)*M+(i,p)",mo_tensor,ao_tensor,mocoef_alpha)
• ierr=talsh_tensor_contract("B(p,q)+=M+(i,p)*A(i,q)",mo_tensor,mocoef_alpha,ao_tensor)

I can update you on the details on your reguar bi-weekly meeting if you can make that one (in 30 minutes from now).

[LucasVisscher]

Made a work-around for exacorr now. But would be nice to have this fixed.