def matrix_element(self, j, m, jp, mp):
        result = hbar*sqrt(mp*(mp - S.One) - j*(j + S.One))
        result *= KroneckerDelta(m, mp - 1)
        result *= KroneckerDelta(j, jp)
        return result
dataset/ETHPy150Open sympy/sympy/sympy/physics/quantum/spin.py/JminusOp.matrix_element/swapped operands of `-`