Density and distribution function for zero_one_inflated_beta

Currently with version 2.13.0 of brms, there are no density or distribution functions for zero_one_inflated_beta.

From the brms output, I can pull mu, phi, zoi and coi. I can derive the a and b parameters for pbeta from mu and phi. Is there a way to calculate the equivalent of pzero_one_inflated_beta with those pieces?

Could you use some of the functions available in brms to instead simulate data from the posterior, rather than feeding this information to a distribution function?