Initial commitv1.0

1 files changed, 81 insertions, 0 deletions

diff --git a/src/integral_double.h b/src/integral_double.h
new file mode 100644
index 0000000..de179fb
--- /dev/null
+++ b/src/integral_double.h
@@ -0,0 +1,81 @@
+/*
+Copyright 2016 Ian Jauslin
+
+Licensed under the Apache License, Version 2.0 (the "License");
+you may not use this file except in compliance with the License.
+You may obtain a copy of the License at
+
+    http://www.apache.org/licenses/LICENSE-2.0
+
+Unless required by applicable law or agreed to in writing, software
+distributed under the License is distributed on an "AS IS" BASIS,
+WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+See the License for the specific language governing permissions and
+limitations under the License.
+*/
+
+/*
+  Preprocessor macros for integration using with doubles
+*/
+
+
+// reset CPP macros
+#undef INTEGRAL_FUNC
+#undef INTEGRAL_FLOAT_TYPE
+#undef INTEGRAL_FLOAT_INIT
+#undef INTEGRAL_FLOAT_FREE
+#undef INTEGRAL_FLOAT_SET_D
+#undef INTEGRAL_FLOAT_SET_UI
+#undef INTEGRAL_FLOAT_ADD
+#undef INTEGRAL_FLOAT_SUB
+#undef INTEGRAL_FLOAT_UI_SUB
+#undef INTEGRAL_FLOAT_MUL
+#undef INTEGRAL_FLOAT_POW_UI
+#undef INTEGRAL_FLOAT_DIV
+#undef INTEGRAL_FLOAT_DIV_UI
+#undef INTEGRAL_FLOAT_ISNUMBER
+#undef INTEGRAL_FLOATARRAY_TYPE
+#undef INTEGRAL_FLOATARRAY_FUNC
+#undef INTEGRAL_POLYNOMIAL_FUNC
+#undef INTEGRAL_POLYNOMIALARRAY_TYPE
+#undef INTEGRAL_POLYNOMIALARRAY_FUNC
+
+// suffix of function names
+#define INTEGRAL_FUNC(NAME) NAME ## _double
+
+// type of floats
+#define INTEGRAL_FLOAT_TYPE double
+// set float from double
+#define INTEGRAL_FLOAT_SET_D(FLOAT, VAL) FLOAT=VAL
+// set float from unsigned int
+#define INTEGRAL_FLOAT_SET_UI(FLOAT, VAL) FLOAT=(double)VAL
+// add floats
+#define INTEGRAL_FLOAT_ADD(FLOAT, VAL1, VAL2) FLOAT=VAL1+VAL2
+// subtract floats
+#define INTEGRAL_FLOAT_SUB(FLOAT, VAL1, VAL2) FLOAT=VAL1-VAL2
+// subtract floats in which the first is specified as an unsigned int
+#define INTEGRAL_FLOAT_UI_SUB(FLOAT, VAL1, VAL2) FLOAT=(double)VAL1-VAL2
+// multiply floats
+#define INTEGRAL_FLOAT_MUL(FLOAT, VAL1, VAL2) FLOAT=VAL1*VAL2
+// power of a float, specified as an unsigned int
+#define INTEGRAL_FLOAT_POW_UI(FLOAT, VAL1, VAL2) FLOAT=pow(VAL1, (double)VAL2)
+// divide floats
+#define INTEGRAL_FLOAT_DIV(FLOAT, VAL1, VAL2) FLOAT=VAL1/VAL2
+// divide floats, one of which is specified as an unsigned int
+#define INTEGRAL_FLOAT_DIV_UI(FLOAT, VAL1, VAL2) FLOAT=VAL1/(double)VAL2
+// check whether a float is a regular number
+#define INTEGRAL_FLOAT_ISNUMBER(FLOAT) (fpclassify(FLOAT)>=FP_ZERO)
+
+// type of float arrays
+#define INTEGRAL_FLOATARRAY_TYPE array_double
+// prefix of float array function names
+#define INTEGRAL_FLOATARRAY_FUNC(NAME) array_double_ ## NAME
+
+// prefix of polynomial function names
+#define INTEGRAL_POLYNOMIAL_FUNC(NAME) polynomial_double_ ## NAME
+
+// type of polynomial arrays
+#define INTEGRAL_POLYNOMIALARRAY_TYPE array_polynomial_double
+// prefix of polynomial array function names
+#define INTEGRAL_POLYNOMIALARRAY_FUNC(NAME) array_polynomial_double_ ## NAME
+