1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
use core::ops::{Add, Mul};
use num_traits::Float;
#[derive(Debug)]
pub struct GaussianDistribution<F: Float> {
pub mean: F,
pub var: F,
}
impl<F: Float> GaussianDistribution<F> {
pub fn new(mean: F, var: F) -> Self {
GaussianDistribution { mean, var }
}
}
impl<F: Float> Add for GaussianDistribution<F> {
type Output = GaussianDistribution<F>;
fn add(self, other: GaussianDistribution<F>) -> GaussianDistribution<F> {
GaussianDistribution {
mean: self.mean + other.mean,
var: self.var + other.var,
}
}
}
impl<F: Float> Mul for GaussianDistribution<F> {
type Output = GaussianDistribution<F>;
fn mul(self, other: GaussianDistribution<F>) -> GaussianDistribution<F> {
let mean = (self.var * other.mean + other.var * self.mean) / (self.var + other.var);
let var = F::one() / (F::one() / self.var + F::one() / other.var);
GaussianDistribution { mean, var }
}
}