C++ code for stochastic spacial neutral/hierarchical models