From d9c2e782c53ce63c9d579c78ac95ae6cd3d9de5f Mon Sep 17 00:00:00 2001
From: Veronique Legrand <vlegrand@pasteur.fr>
Date: Thu, 10 Oct 2024 16:56:31 +0200
Subject: [PATCH] fixed bug in test
---
src/unit_test_math_utils.cpp | 30 +++++++++++++-----------------
1 file changed, 13 insertions(+), 17 deletions(-)
diff --git a/src/unit_test_math_utils.cpp b/src/unit_test_math_utils.cpp
index 392ce89..eeed7f5 100644
--- a/src/unit_test_math_utils.cpp
+++ b/src/unit_test_math_utils.cpp
@@ -33,9 +33,8 @@ void test_getCollisionProba() {
float p=getCollisionProba(2,5000000000,UINT_MAX,2);
assert(round(p*10000)==127);
p=getCollisionProba(2,5000000000,UINT_MAX,1);
- assert(p=0.1128);
+ assert(round(p*10000)==1128);
p =getCollisionProba(1,2,UINT_MAX,3);
- //cout<<p<<endl;
assert(p==0.0);
p =getCollisionProba(1,1000000000,UINT_MAX,1);
assert(floor(p*1000+0.5)==23);
@@ -46,30 +45,27 @@ void test_getCollisionProba() {
}
void test_gammln() {
- float a;
+ float a,tmp;
a=gammln(1);
- assert(exp(a)==1);
+ tmp=exp(a);
+ assert((round(tmp*10000))/10000==1);
a=gammln(2);
- assert(exp(a)==1);
+ tmp=exp(a);
+ assert((round(tmp*10000))/10000==1);
a=gammln(3);
- assert(exp(a)==2);
+ tmp=exp(a);
+ assert((round(tmp*10000))/10000==2);
a=gammln(4);
- assert(exp(a)==6);
+ tmp=exp(a);
+ assert((round(tmp*10000))/10000==6);
a=gammln(5);
- assert(exp(a)==24);
+ tmp=exp(a);
+ assert((round(tmp*10000))/10000==24);
a=gammln(6);
- float tmp=exp(a);
+ tmp=exp(a);
float tmp2=tmp*10000;
float tmp3=round(tmp2)/10000;
- //printf("%s \n",tmp3);
assert(tmp3==120);
- /* unsigned long n=5000;
- unsigned long m=n+1;
- a=gammln(m);
- float b=gammln(n);
- double truc=exp(a-b);
- assert(truc==n);
- */
}
void test_pmf() {
--
GitLab