use SLD = real + i*imag rather than N b_C - i * N b''