@@ -1033,24 +1033,15 @@ def test_cumulant_function(self):
1033
1033
cumulant_function_second_order = numeric .calculate_cumulant_function (
1034
1034
pulse , spectrum , omega , second_order = True
1035
1035
)
1036
- # Make sure first and second order are --roughly -- of same order
1037
- # of magnitude. Unlike the frequency shifts themselves, the
1038
- # second order contributions to the cumulant function vanish on the
1039
- # diagonal, whereas the first order contributions dominate. Hence,
1040
- # be quite lenient.
1041
1036
second_order_contribution = (cumulant_function_second_order
1042
1037
- cumulant_function_first_order )
1043
- rel = (np .linalg .norm (cumulant_function_first_order )
1044
- / np .linalg .norm (second_order_contribution ))
1045
1038
1046
1039
# Second order terms should be anti-hermitian
1047
1040
self .assertArrayAlmostEqual (second_order_contribution ,
1048
1041
- second_order_contribution .transpose (0 , 2 , 1 ),
1049
1042
atol = 1e-16 )
1050
1043
self .assertEqual (cumulant_function_first_order .shape ,
1051
1044
cumulant_function_second_order .shape )
1052
- self .assertLessEqual (rel , 200 )
1053
- self .assertGreaterEqual (rel , 1 / 10 )
1054
1045
1055
1046
def test_error_transfer_matrix (self ):
1056
1047
"""Test raises of numeric.error_transfer_matrix."""
0 commit comments