From 618a5de3349f8def5bdb909fea42fb447b6aad4a Mon Sep 17 00:00:00 2001
From: "bernard.parisse" <bernard.parisse@23ce0884-8a58-47d3-bc5c-ddf1cd5b9f9e>
Date: Tue, 22 Aug 2023 09:35:21 +0000
Subject: [PATCH] [giac] random: randmax2+1.0 fixes move fgamma definition from
 first.h to usual.cc, insures that config.h defines are defined

git-svn-id: https://dev.geogebra.org/svn/trunk/geogebra/giac@70150 23ce0884-8a58-47d3-bc5c-ddf1cd5b9f9e
---
 src/giac/cpp/prog.cc     |  8 ++++----
 src/giac/cpp/usual.cc    | 10 ++++++++++
 src/giac/headers/first.h |  9 ---------
 3 files changed, 14 insertions(+), 13 deletions(-)

diff --git a/src/giac/cpp/prog.cc b/src/giac/cpp/prog.cc
index d9f1ffde..2ac8162a 100644
--- a/src/giac/cpp/prog.cc
+++ b/src/giac/cpp/prog.cc
@@ -4020,7 +4020,7 @@ namespace giac {
     for (int i=0;i<n-1;++i){
       // j ← random integer such that i ≤ j < n
       // exchange a[i] and a[j]
-      int j=int(i+(giac_rand(contextptr)/double(rand_max2))*(n-i));
+      int j=int(i+(giac_rand(contextptr)/double(rand_max2+1.0))*(n-i));
       std::swap(temp[i],temp[j]);
     }
   }
@@ -4034,7 +4034,7 @@ namespace giac {
       for (int essai=20;essai>=0;--essai){
 	int i;
 	for (i=0;i<k;++i)
-	  ts[i]=t[i]=int(giac_rand(contextptr)/double(rand_max2)*n);
+	  ts[i]=t[i]=int(giac_rand(contextptr)/double(rand_max2+1.0)*n);
 	sort(ts.begin(),ts.end());
 	for (i=1;i<k;++i){
 	  if (ts[i]==ts[i-1])
@@ -4048,7 +4048,7 @@ namespace giac {
       vector<int> t; t.reserve(k);
       // (algorithm suggested by O. Garet)
       while (n>0){
-	int r=int(giac_rand(contextptr)/double(rand_max2)*n);
+	int r=int(giac_rand(contextptr)/double(rand_max2+1.0)*n);
 	if (r<n-k) // (n-k)/n=proba that the current n is not in the list
 	  --n;
 	else {
@@ -4068,7 +4068,7 @@ namespace giac {
     for (int j=0;j<k;++j){
       int r=-1;
       for (;;){
-	r=int(giac_rand(contextptr)/double(rand_max2)*n);
+	r=int(giac_rand(contextptr)/double(rand_max2+1.0)*n);
 	if (tab[r]){ tab[r]=false;  break; }
       }
       v[j]=r;
diff --git a/src/giac/cpp/usual.cc b/src/giac/cpp/usual.cc
index 0647cae4..1a64aca4 100644
--- a/src/giac/cpp/usual.cc
+++ b/src/giac/cpp/usual.cc
@@ -81,6 +81,16 @@ using namespace std;
 #define sprintf256 sprintf
 #endif
 
+#if defined HAVE_TGAMMAF || defined __APPLE__ || defined EMCC || defined EMCC2 || defined NO_BSD 
+inline float fgamma(float f1){ return tgammaf(f1); }
+#else
+#if defined(__MINGW_H) || defined(VISUALC) || defined(FXCG)// FIXME gamma, not used
+inline float fgamma(float f1){ return f1; }
+#else
+inline float fgamma(float f1){ return gammaf(f1); } // or tgammaf(f1) on some versions of emscripten
+#endif
+#endif
+
 #ifndef NO_NAMESPACE_GIAC
 namespace giac {
 #endif // ndef NO_NAMESPACE_GIAC
diff --git a/src/giac/headers/first.h b/src/giac/headers/first.h
index 28d3c563..ad0dd4aa 100644
--- a/src/giac/headers/first.h
+++ b/src/giac/headers/first.h
@@ -578,15 +578,6 @@ inline float ffloor(float f1){
 #endif
 }
 inline float finv(float f1){ return 1/f1; }
-#if defined HAVE_TGAMMAF || defined __APPLE__ || defined EMCC || defined EMCC2 || defined NO_BSD 
-inline float fgamma(float f1){ return tgammaf(f1); }
-#else
-#if defined(__MINGW_H) || defined(VISUALC) || defined(FXCG)// FIXME gamma, not used
-inline float fgamma(float f1){ return f1; }
-#else
-inline float fgamma(float f1){ return gammaf(f1); } // or tgammaf(f1) on some versions of emscripten
-#endif
-#endif
 #ifdef FXCG
 inline float atan2f(float f1,float f2,int rad){ if (rad) return std::atan2(f1,f2); else return std::atan2(f1,f2)*180/3.14159265358979323846;}
 #else