metromc 0.2.0

Markov chain Monte Carlo sampling using the Independence Metropolis-Hastings algorithm.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

use statrs::distribution::Continuous;
use crate::Pdf;

pub struct Gamma {
    dist: statrs::distribution::Gamma,
}
impl Gamma {
    #[inline(always)]
    pub fn new(shape: f64, rate: f64) -> Self {
        Self {
            dist: statrs::distribution::Gamma::new(shape, rate).unwrap(),
        }
    }
}

impl Pdf for Gamma {
    #[inline(always)]
    fn prob(&self, x: f64) -> f64 {
        self.dist.pdf(x)
    }
}