summaryrefslogtreecommitdiff
path: root/src/libs/math
diff options
context:
space:
mode:
Diffstat (limited to 'src/libs/math')
-rw-r--r--src/libs/math/Makeinfo2
-rw-r--r--src/libs/math/mthintrn.h25
-rw-r--r--src/libs/math/random.c101
-rw-r--r--src/libs/math/random.h56
-rw-r--r--src/libs/math/random2.c89
-rw-r--r--src/libs/math/sqrt.c97
6 files changed, 370 insertions, 0 deletions
diff --git a/src/libs/math/Makeinfo b/src/libs/math/Makeinfo
new file mode 100644
index 0000000..da10867
--- /dev/null
+++ b/src/libs/math/Makeinfo
@@ -0,0 +1,2 @@
+uqm_CFILES="random.c random2.c sqrt.c"
+uqm_HFILES="mthintrn.h random.h"
diff --git a/src/libs/math/mthintrn.h b/src/libs/math/mthintrn.h
new file mode 100644
index 0000000..cfda767
--- /dev/null
+++ b/src/libs/math/mthintrn.h
@@ -0,0 +1,25 @@
+//Copyright Paul Reiche, Fred Ford. 1992-2002
+
+/*
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
+ */
+
+#ifndef LIBS_MATH_MTHINTRN_H_
+#define LIBS_MATH_MTHINTRN_H_
+
+#include "libs/mathlib.h"
+
+#endif /* LIBS_MATH_MTHINTRN_H_ */
+
diff --git a/src/libs/math/random.c b/src/libs/math/random.c
new file mode 100644
index 0000000..83fcff8
--- /dev/null
+++ b/src/libs/math/random.c
@@ -0,0 +1,101 @@
+///Copyright Paul Reiche, Fred Ford. 1992-2002
+
+/*
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
+ */
+
+ /****************************************************************************
+* FILE: random.c
+* DESC: a library of random number generators for general purpose use.
+*
+* References:
+* "Random Number Generators: Good ones are hard to find" S.K.Park & K.W.Miller
+* Communications of the ACM, Vol31 Number 10, October 1988, Pp 1192-1201
+*
+* HISTORY: Created 1/23/1989
+* LAST CHANGED:
+*
+* Copyright (c) 1989, Robert Leyland and Fred Ford
+****************************************************************************/
+
+/* ----------------------------INCLUDES----------------------------------- */
+#include "mthintrn.h" /* get the externs for error checking */
+
+/* ----------------------------DEFINES------------------------------------ */
+/* constants for licongruential random number generator from CACM article
+ referenced above */
+#define A 16807 /* a relatively prime number -- also M div Q */
+#define M 2147483647L /* 0xFFFFFFFF / 2 */
+#define Q 127773L /* M div A */
+#define R 2836 /* M mod A */
+
+/* ----------------------------STATIC DATA-------------------------------- */
+
+static DWORD seed = 12345L; /* random number seed */
+
+/* ----------------------------CODE--------------------------------------- */
+
+/*****************************************************************************
+* FUNC: DWORD TFB_Random()
+*
+* DESC: random number generator
+*
+* NOTES:
+*
+* HISTORY: Created By Robert leyland
+*
+*****************************************************************************/
+
+DWORD
+TFB_Random (void)
+{
+ seed = A * (seed % Q) - R * (seed / Q);
+ if (seed > M)
+ return (seed -= M);
+ else if (seed)
+ return (seed);
+ else
+ return (seed = 1L);
+}
+
+/*****************************************************************************
+* FUNC: DWORD TFB_SeedRandom(DWORD l)
+*
+* DESC: set the seed for the random number generator to parameter "l", and
+* return the value of the previously active seed, to allow for multiple
+* random number streams.
+*
+* NOTES: if the seed is not valid it will be coerced into a valid range
+*
+* HISTORY: Created By Robert leyland
+*
+*****************************************************************************/
+
+DWORD
+TFB_SeedRandom (DWORD new_seed)
+{
+ DWORD old_seed;
+
+ /* coerce the seed to be in the range 1..M */
+ if (new_seed == 0L) /* 0 becomes 1 */
+ new_seed = 1;
+ else if (new_seed > M) /* and less than M */
+ new_seed -= M;
+
+ old_seed = seed;
+ seed = new_seed;
+ return (old_seed);
+}
+
diff --git a/src/libs/math/random.h b/src/libs/math/random.h
new file mode 100644
index 0000000..0e8d602
--- /dev/null
+++ b/src/libs/math/random.h
@@ -0,0 +1,56 @@
+//Copyright Paul Reiche, Fred Ford. 1992-2002
+
+/*
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
+ */
+
+/****************************************************************************
+* FILE: random.h
+* DESC: definitions and externs for random number generators
+*
+* HISTORY: Created 6/ 6/1989
+* LAST CHANGED:
+*
+* Copyright (c) 1989, Robert Leyland and Scott Anderson
+****************************************************************************/
+
+#ifndef LIBS_MATH_RANDOM_H_
+#define LIBS_MATH_RANDOM_H_
+
+/* ----------------------------GLOBALS/EXTERNS---------------------------- */
+
+DWORD TFB_SeedRandom (DWORD seed);
+DWORD TFB_Random (void);
+
+
+typedef struct RandomContext RandomContext;
+
+#ifdef RANDOM2_INTERNAL
+struct RandomContext {
+ DWORD seed;
+};
+#endif
+
+RandomContext *RandomContext_New (void);
+void RandomContext_Delete (RandomContext *context);
+RandomContext *RandomContext_Copy (const RandomContext *source);
+DWORD RandomContext_Random (RandomContext *context);
+DWORD RandomContext_SeedRandom (RandomContext *context, DWORD new_seed);
+DWORD RandomContext_GetSeed (RandomContext *context);
+
+
+#endif /* LIBS_MATH_RANDOM_H_ */
+
+
diff --git a/src/libs/math/random2.c b/src/libs/math/random2.c
new file mode 100644
index 0000000..9d354b4
--- /dev/null
+++ b/src/libs/math/random2.c
@@ -0,0 +1,89 @@
+/*
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
+ */
+
+// This file contains variants of the random functions in random.c
+// that store the state of the RNG in a context, allowing for multiple
+// independant RNGs to be used simultaneously.
+// The RNG behavior itself is the same.
+
+#include "libs/compiler.h"
+
+#define RANDOM2_INTERNAL
+#include "random.h"
+
+#include "libs/memlib.h"
+
+
+#define A 16807 /* a relatively prime number -- also M div Q */
+#define M 2147483647 /* 0xFFFFFFFF / 2 */
+#define Q 127773 /* M div A */
+#define R 2836 /* M mod A */
+
+RandomContext *
+RandomContext_New (void)
+{
+ RandomContext *result = (RandomContext *) HMalloc (sizeof (RandomContext));
+ result->seed = 12345;
+ return result;
+}
+
+void
+RandomContext_Delete (RandomContext *context)
+{
+ HFree ((void *) context);
+}
+
+RandomContext *
+RandomContext_Copy (const RandomContext *source)
+{
+ RandomContext *result = (RandomContext *) HMalloc (sizeof (RandomContext));
+ *result = *source;
+ return result;
+}
+
+DWORD
+RandomContext_Random (RandomContext *context)
+{
+ context->seed = A * (context->seed % Q) - R * (context->seed / Q);
+ if (context->seed > M) {
+ context->seed -= M;
+ } else if (context->seed == 0)
+ context->seed = 1;
+
+ return context->seed;
+}
+
+DWORD
+RandomContext_SeedRandom (RandomContext *context, DWORD new_seed)
+{
+ DWORD old_seed;
+
+ /* coerce the seed to be in the range 1..M */
+ if (new_seed == 0) /* 0 becomes 1 */
+ new_seed = 1;
+ else if (new_seed > M) /* and less than M */
+ new_seed -= M;
+
+ old_seed = context->seed;
+ context->seed = new_seed;
+ return old_seed;
+}
+
+DWORD
+RandomContext_GetSeed (RandomContext *context)
+{
+ return context->seed;
+}
diff --git a/src/libs/math/sqrt.c b/src/libs/math/sqrt.c
new file mode 100644
index 0000000..1f02cac
--- /dev/null
+++ b/src/libs/math/sqrt.c
@@ -0,0 +1,97 @@
+//Copyright Paul Reiche, Fred Ford. 1992-2002
+
+/*
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
+ */
+
+#include "mthintrn.h"
+
+COUNT
+square_root (DWORD value)
+{
+ UWORD sig_word, mask;
+ COUNT result, shift;
+
+ if ((sig_word = HIWORD (value)) > 0)
+ {
+ DWORD mask_squared, result_shift;
+
+ for (mask = 1 << 15, shift = 31;
+ !(mask & sig_word); mask >>= 1, --shift)
+ ;
+ shift >>= 1;
+ mask = 1 << shift;
+
+ result = mask;
+ mask_squared = result_shift = (DWORD)mask << shift;
+ value -= mask_squared;
+ while (mask >>= 1)
+ {
+ DWORD remainder;
+
+ mask_squared >>= 1;
+ mask_squared >>= 1;
+ if ((remainder = result_shift + mask_squared) > value)
+ result_shift >>= 1;
+ else
+ {
+ value -= remainder;
+
+ result_shift = (result_shift >> 1) + mask_squared;
+
+ result |= mask;
+ }
+ }
+
+ return (result);
+ }
+ else if ((sig_word = LOWORD (value)) > 0)
+ {
+ UWORD mask_squared, result_shift;
+
+ for (mask = 1 << 15, shift = 15;
+ !(mask & sig_word); mask >>= 1, --shift)
+ ;
+ shift >>= 1;
+ mask = 1 << shift;
+
+ result = mask;
+ mask_squared = result_shift = mask << shift;
+ sig_word -= mask_squared;
+ while (mask >>= 1)
+ {
+ UWORD remainder;
+
+ mask_squared >>= 1;
+ mask_squared >>= 1;
+ if ((remainder = result_shift + mask_squared) > sig_word)
+ result_shift >>= 1;
+ else
+ {
+ sig_word -= remainder;
+
+ result_shift = (result_shift >> 1) + mask_squared;
+
+ result |= mask;
+ }
+ }
+
+ return (result);
+ }
+
+ return (0);
+}
+
+
generated by cgit v1.2.3 (git 2.25.1) at 2024-05-19 05:51:36 +0000