I have written a c-code for for generating samples from joint posterior density using Gibbs sampling technique. I have used Adaptive Rejection sampling (Wild & Gilks (1993)) to generate
samples from the univariate conditional density. However, the code is working fine for some cases,
for some cases it is not working. Is anybody aware of the source code for ARS? I can share my code if you can help me to find out the bug in the code.