Dear developers, Let me point out a small inconsistency in the way siesta outputs stress, that confused me for a while: The off-diagonal elements of "Stress-tensor-Voigt" are written in the order 12,23,13 (write_subs.F). This is at variance with the standard convention, 23,13,12, which by the way is followed in the GeometryConstraints block (fixed.F).
Best regards, Roberto
