Dear Joaquín,
>From your previous posts I assume that all of PeFT, Pe[qs] and F[qs]
represent rank 2 tensors. If this is correct then no, the contraction
indices that you've specified are not correct. They're out of the range
[0,1] that are valid for rank-2 tensors, and this would presumably trigger
an assertion in debug mode (at the very least).
For simplicity, lets call these three rank-2 tensors "A,B,C". Then we could
compute the two possible single index contraction results, expressed in
index notation, as follows:
(1a) A_{ij} = B_{ik} C_{kj} --> A = contract<1,0>(B,C); // The second
index of the first tensor contracts with the first index of the second
tensor
or, more simply,
(1b) A = B*C;
(2)A_{ij} = B_{ik} C_{kj} --> A = contract<1,1>(B,C); // The second index
of the first tensor contracts with the second index of the second tensor
Does this explain a bit better how the function works?
Regards,
J-P
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
