p0 does calcs too.. may or may not be working
Still don't work
Should work for np > dim
works!
Tidied up if branches for np==0