def test_quantum_number_S(self, name, rho, expected_S):
        # Setup
        Job = InitJob("test_data/JobDef_scase.json")
        # Spin 0 density matrix
        Job.Electron.rho = np.matrix(rho)

        # Action
        S = Job.Electron.quantum_number_S(Job)

        # Result
        assert S == expected_S