@jeffriley says they sometimes do (to be reproduced -- he will check for the right parameters?), in which case
WhiteDwarfs::CalculateRadiusOnPhase_Static() will have a square root of a negative number.
[Also, the comparison with zero in that function should be a direct comparison, without using Compare() ]