In the Worksheet "Walker-Penrose Killing tensor in Kerr spacetime" of the SageManifolds project, the sheet gives the principal null directions for Kerr spacetime in Boyer-Lindquist coordinates. I have difficulties verifying (with Mathematica) that this is indeed the case and would appreciate a source.