polarized tanh interface: discrepancy between code and doc
In the DevRef, section "Interface with tanh profile", we have diagonal matrices with elements like h_a^_+ etc. The "+" elements are on top, the "-" elements on bottom. In the implementation, in function Compute::refractionMatrixBlocksTanh, it is just the other way round.
If I swap the elements in the implementation, then several tests break; one with unphysical T>1.