fix normalP for zero and unit ideal input, singular trac 596

fixed `normalP()` for zero ideal input.
`normalP()` should now also behave similar to `normal()` for the unit ideal case.
(please review/check!)

open task: update documentation for normal() and normalP() for the unit ideal case 
 will be done in a different pull request.