Update outfct

1. unify P, nz
2. add blockwise cov
3. multinom.simdata cov=cov