Imported Upstream version 0.0.34

1 files changed, 76 insertions, 0 deletions

diff --git a/CSXCAD/src/CSFunctionParser.cpp b/CSXCAD/src/CSFunctionParser.cpp
new file mode 100644
index 0000000..c8741ff
--- /dev/null
+++ b/CSXCAD/src/CSFunctionParser.cpp
@@ -0,0 +1,76 @@
+/*
+*	Copyright (C) 2010,2011 Thorsten Liebig (Thorsten.Liebig@gmx.de)
+*
+*	This program is free software: you can redistribute it and/or modify
+*	it under the terms of the GNU Lesser General Public License as published
+*	by the Free Software Foundation, either version 3 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 Lesser General Public License for more details.
+*
+*	You should have received a copy of the GNU Lesser General Public License
+*	along with this program.  If not, see <http://www.gnu.org/licenses/>.
+*/
+
+#include "CSFunctionParser.h"
+#include <math.h>
+#include <iostream>
+
+double bessel_first_kind_0(const double* p)
+{
+	return j0(p[0]);
+}
+double bessel_first_kind_1(const double* p)
+{
+	return j1(p[0]);
+}
+double bessel_first_kind_n(const double* p)
+{
+	int n=p[0];
+	if (n<0)
+	{
+		std::cerr << "CSFunctionParser::bessel_first_kind_n (jn): first argument must be integer larger than zero! found: " << p[0] << std::endl;
+		return 0;
+	}
+	return jn(n,p[1]);
+}
+
+double bessel_second_kind_0(const double* p)
+{
+	return y0(p[0]);
+}
+double bessel_second_kind_1(const double* p)
+{
+	return y1(p[0]);
+}
+double bessel_second_kind_n(const double* p)
+{
+	int n=p[0];
+	if (n<0)
+	{
+		std::cerr << "CSFunctionParser::bessel_second_kind_n (yn): first argument must be integer larger than zero! found: " << p[0] << std::endl;
+		return 0;
+	}
+	return yn(n,p[1]);
+}
+
+
+CSFunctionParser::CSFunctionParser()
+{
+	//some usefull constants
+	AddConstant("pi", 3.14159265358979323846);
+	AddConstant("e", 2.71828182845904523536);
+
+	//some bessel functions of first kind
+	AddFunction("j0",bessel_first_kind_0,1);
+	AddFunction("j1",bessel_first_kind_1,1);
+	AddFunction("jn",bessel_first_kind_n,2);
+
+	//some bessel functions of second kind
+	AddFunction("y0",bessel_second_kind_0,1);
+	AddFunction("y1",bessel_second_kind_1,1);
+	AddFunction("yn",bessel_second_kind_n,2);
+}