X-Git-Url: https://ginac.de/CLN/cln.git//cln.git?a=blobdiff_plain;ds=inline;f=src%2Freal%2Ftranscendental%2Fcl_R_atan2.cc;fp=src%2Freal%2Ftranscendental%2Fcl_R_atan2.cc;h=1e967a99b341e8977ddf2df7d8d68bee6322453e;hb=dd9e0f894eec7e2a8cf85078330ddc0a6639090b;hp=0000000000000000000000000000000000000000;hpb=d57000b193438122b366aa0d450674be9402d06c;p=cln.git

diff --git a/src/real/transcendental/cl_R_atan2.cc b/src/real/transcendental/cl_R_atan2.cc
new file mode 100644
index 0000000..1e967a9
--- /dev/null
+++ b/src/real/transcendental/cl_R_atan2.cc
@@ -0,0 +1,84 @@
+// atan().
+
+// General includes.
+#include "cl_sysdep.h"
+
+// Specification.
+#include "cl_real.h"
+
+
+// Implementation.
+
+#include "cl_float.h"
+#include "cl_F_tran.h"
+#include "cl_N.h"
+#include "cl_R.h"
+
+const cl_R atan (const cl_R& x, const cl_R& y)
+{
+// Methode:
+// y=0 -> bei x>0: 0 als Ergebnis,
+//        bei x<0: pi als Ergebnis.
+//        bei x=0: Error.
+// x=0 -> bei y>0: pi/2 als Ergebnis.
+//        bei y<0: -pi/2 als Ergebnis.
+//        bei y=0: Error.
+// Falls x und y beide rational: beide in Floats umwandeln.
+// 0 <= |y| <= x  ->  atan(y/x)
+// 0 <= |x| <= y  ->  pi/2 - atan(x/y)
+// 0 <= |x| <= -y  ->  -pi/2 - atan(x/y)
+// 0 <= |y| <= -x  ->  für y>=0: pi + atan(y/x), für y<0: -pi + atan(y/x)
+
+	if (eq(y,0)) {
+		// y=0 (exakt)
+		if (zerop(x)) // x=0 -> Error
+			{ cl_error_division_by_0(); }
+		if (minusp(x)) // x<0 -> pi in Default-Float-Genauigkeit
+			{ return cl_pi(); }
+		return 0; // x>0 -> 0
+	}
+	elif (eq(x,0)) {
+		// x=0 (exakt)
+		if (zerop(y)) // y=0 -> Error
+			{ cl_error_division_by_0(); }
+		if (minusp(y)) // y<0 -> -pi/2
+			{ return - scale_float(cl_pi(),-1); }
+		return scale_float(cl_pi(),-1); // y>0 -> pi/2
+	} else {
+		Mutable(cl_R,x);
+		Mutable(cl_R,y);
+		// Check special case of rational numbers:
+		if (rationalp(x))
+			if (rationalp(y)) {
+				// x,y in Floats umwandeln:
+				x = cl_float(The(cl_RA)(x));
+				y = cl_float(The(cl_RA)(y));
+			}
+		// x,y nicht exakt =0, x/y und y/x werden Floats sein.
+		if (abs(x) >= abs(y)) {
+			// |x| >= |y|
+			var cl_F z = atanx(The(cl_F)(y/x));
+			// Division war erfolgreich, also x/=0.
+			if (minusp(x))
+				// x<0 -> pi bzw. -pi addieren:
+				if (!minusp(y))
+					// y>=0 -> atan(y/x) + pi
+					return z + cl_pi(z);
+				else
+					// y<0 -> atan(y/x) - pi
+					return z - cl_pi(z);
+			else
+				return z;
+		} else {
+			// |x| < |y|
+			var cl_F z = atanx(The(cl_F)(x/y));
+			// von pi/2 bzw. -pi/2 subtrahieren:
+			if (!minusp(y))
+				// y>=0 -> pi/2 - atan(x/y)
+				return scale_float(cl_pi(z),-1) - z;
+			else
+				// y<0 -> -pi/2 - atan(x/y)
+				return - scale_float(cl_pi(z),-1) - z;
+		}
+	}
+}