From fb6edbff0c674508c212e1ce127f5cf6ae451e3e Mon Sep 17 00:00:00 2001
From: NunoSempere <nuno.sempere@protonmail.com>
Date: Fri, 11 Aug 2023 14:01:02 +0200
Subject: [PATCH] add sample_laplace

---
 squiggle.c | 4 ++++
 squiggle.h | 1 +
 2 files changed, 5 insertions(+)

diff --git a/squiggle.c b/squiggle.c
index 51076c4..0150b6a 100644
--- a/squiggle.c
+++ b/squiggle.c
@@ -134,6 +134,10 @@ double sample_beta(double a, double b, uint64_t* seed)
     return gamma_a / (gamma_a + gamma_b);
 }
 
+double sample_laplace(double successes, double failures, uint64_t* seed){
+	return sample_beta(successes + 1, failures + 1, seed);
+}
+
 // Array helpers
 double array_sum(double* array, int length)
 {
diff --git a/squiggle.h b/squiggle.h
index 8d8bd82..cd1eaad 100644
--- a/squiggle.h
+++ b/squiggle.h
@@ -19,6 +19,7 @@ double sample_to(double low, double high, uint64_t* seed);
 
 double sample_gamma(double alpha, uint64_t* seed);
 double sample_beta(double a, double b, uint64_t* seed);
+double sample_laplace(double successes, double failures, uint64_t* seed);
 
 // Array helpers
 double array_sum(double* array, int length);