raise exception if power function gives P<0 for any k values that we use
(#305)
unit tests of tabulated P(k)
(#305)