Add dependent double observer sampling as built-in option for gmultmix and multinomPois, with type = 'depDouble'.
I also did a little reorganization to split the piFuns out into their own R file (and their own cpp file), to make it more clear they are shared between functions, and to make adding future piFuns easier. Updated the docs to reflect the changes and added relevant tests.
Finally, I corrected a related bug when dealing with missing values in the gmultmix C engine.
Add dependent double observer sampling as built-in option for
gmultmixandmultinomPois, withtype = 'depDouble'.I also did a little reorganization to split the piFuns out into their own R file (and their own cpp file), to make it more clear they are shared between functions, and to make adding future piFuns easier. Updated the docs to reflect the changes and added relevant tests.
Finally, I corrected a related bug when dealing with missing values in the gmultmix C engine.