One should use the primitiveVectors(), not conventionalVectors(), when you compute the fractional coordinates.
If you change that single line in the routine findFractionalCoordinates() I believe it will work as intended, that is line 12 in stress_optimizers.py should be:
uc = configuration.bravaisLattice().primitiveVectors()