Due to !176 (merged) an assertion was introduced in the creation of the symbolic relaxation matrix. The assertion tests that all entries of the symbolic relaxation matrix are sympy symbols. However, when a relaxation_rates_modifier this might not be the case.
This MR fixes the problem and introduces a test case with the psm method