-
Notifications
You must be signed in to change notification settings - Fork 15
Open
Description
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).
Metadata
Metadata
Assignees
Labels
No labels