aCos is inversed ?

.gitignore

@@ -31,4 +31,4 @@
 *.bin
 
 # Temp
-*~
+*~

go.sh

@@ -1,6 +1,21 @@
 #!/bin/bash
 echo Run test of complex libs
 cd ./tests/dev
+echo **Debug
+echo ***Clang
 rm -rf *.bin
-clang++ test-complex-lib.cc ../../ucomplex.cc -o test-complex-lib.cc.bin
+clang++ -std=c++11 test-complex-lib.cc ../../ucomplex.cc -o test-complex-lib.cc.bin
+./test-complex-lib.cc.bin
+echo ***Gcc
+rm -rf *.bin
+g++ -std=c++11 test-complex-lib.cc ../../ucomplex.cc -o test-complex-lib.cc.bin
+./test-complex-lib.cc.bin
+echo **Build
+echo ***Clang
+rm -rf *.bin
+clang++ -std=c++11 -O2 test-complex-lib.cc ../../ucomplex.cc -o test-complex-lib.cc.bin
+./test-complex-lib.cc.bin
+echo ***Gcc
+rm -rf *.bin
+g++ -std=c++11 -O2  test-complex-lib.cc ../../ucomplex.cc -o test-complex-lib.cc.bin
 ./test-complex-lib.cc.bin

tests/dev/test-complex-lib.cc

@@ -11,7 +11,41 @@ void output(std::string operation, complex c1, std::complex<double> c2) {
 	   c1.r, c1.i,
 	   c2.real(), c2.imag());
     // Epsilon for double is 1e-16, check it 1.0 + 1e-16 = 1.0
-    if (diff_r > 1e-16||diff_i > 1e-16) printf("WARNING!! Non-zero diff!!!\n");
+    //if (diff_r > 1e-16||diff_i > 1e-16) 
+    if (diff_r > 1e-15||diff_i > 1e-15) 
+	printf("\n********\n\t\tWARNING!! Non-zero diff!!!\n********\n");
+}
+/********************************************************************/
+/********************************************************************/
+/********************************************************************/
+void output(std::string operation, complex c1, complex c2) {
+    double diff_r = c1.r - c2.r;
+    double diff_i = c1.i - c2.i;
+    printf("%s: \tdiff=(%g,%g), \tc1 = (%g,%g),  \tc2 = (%g,%g)\n",
+	   operation.c_str(),
+	   diff_r, diff_i,
+	   c1.r, c1.i,
+	   c2.r, c2.i);
+    // Epsilon for double is 1e-16, check it 1.0 + 1e-16 = 1.0
+    if (diff_r > 1e-16||diff_i > 1e-16) 
+	//if (diff_r > 1e-15||diff_i > 1e-15) 
+	printf("\n********\n\t\tWARNING!! Non-zero diff!!!\n********\n");
+}
+/********************************************************************/
+/********************************************************************/
+/********************************************************************/
+void output(std::string operation, std::complex<double> c1, std::complex<double> c2) {
+    double diff_r = c1.real() - c2.real();
+    double diff_i = c1.imag() - c2.imag();
+    printf("%s: \tdiff=(%g,%g), \tc1 = (%g,%g),  \tc2 = (%g,%g)\n",
+	   operation.c_str(),
+	   diff_r, diff_i,
+	   c1.real(), c1.imag(),
+	   c2.real(), c2.imag());
+    // Epsilon for double is 1e-16, check it 1.0 + 1e-16 = 1.0
+    if (diff_r > 1e-16||diff_i > 1e-16) 
+	//if (diff_r > 1e-15||diff_i > 1e-15) 
+	printf("\n********\n\t\tWARNING!! Non-zero diff!!!\n********\n");
 }
 /********************************************************************/
 /********************************************************************/
@@ -22,7 +56,8 @@ void output_double(std::string operation, double c1, double c2) {
 	   operation.c_str(),
 	   diff, c1, c2);
     // Epsilon for double is 1e-16, check it 1.0 + 1e-16 = 1.0
-    if (diff > 1e-16) printf("WARNING!! Non-zero diff!!!\n");
+    if (diff > 1e-16)
+	printf("\n********\n\t\tWARNING!! Non-zero diff!!!\n********\n");
 }
 /********************************************************************/
 /********************************************************************/
@@ -72,10 +107,10 @@ int main() {
     output("sDiv", c1, c2);
 
     c1 = Conjg(a1);
-    c2 = conj(a2);
+    c2 = std::conj(a2);
     output("Conj", c1, c2);
     c1 = Conjg(b1);
-    c2 = conj(b2);
+    c2 = std::conj(b2);
     output("sConj", c1, c2);
 
     c1 = Cinv(a1);
@@ -90,9 +125,74 @@ int main() {
     output_double("Arg", Carc(a1), std::arg(a2));
     output_double("sArg", Carc(b1), std::arg(b2));
 
-    c1 = Cinv(a1);
-    c2 = 1.0/a2;
-    output("Inv", c1, c2);
+    c1 = Cexp(a1);
+    c2 = std::exp(a2);
+    output("Exp", c1, c2);
+    c1 = Cexp(b1);
+    c2 = std::exp(b2);
+    output("sExp", c1, c2);
+
+    c1 = Clog(a1);
+    c2 = std::log(a2);
+    output("Log", c1, c2);
+    c1 = Clog(b1);
+    c2 = std::log(b2);
+    output("sLog", c1, c2);
+
+    c1 = Csqrt(a1);
+    c2 = std::sqrt(a2);
+    output("Sqrt", c1, c2);
+    c1 = Csqrt(b1);
+    c2 = std::sqrt(b2);
+    output("sSqrt", c1, c2);
+
+    c1 = Ccos(a1);
+    c2 = std::cos(a2);
+    output("Cos", c1, c2);
+    c1 = Ccos(b1);
+    c2 = std::cos(b2);
+    output("sCos", c1, c2);
+
+    c1 = Csin(a1);
+    c2 = std::sin(a2);
+    output("Sin", c1, c2);
+    c1 = Csin(b1);
+    c2 = std::sin(b2);
+    output("sSin", c1, c2);
+
+    c1 = Ctan(a1);
+    c2 = std::tan(a2);
+    output("Tan", c1, c2);
+    c1 = Ctan(b1);
+    c2 = std::tan(b2);
+    output("sTan", c1, c2);
+
+    c1 = Carc_cos(a1);
+    c2 = std::acos(a2);
+    output("aCos", c1, c2);
+    c1 = Carc_cos(b1);
+    c2 = std::acos(b2);
+    output("saCos", c1, c2);
+
+    c1 = Ccos(Carc_cos(b1));
+    output("saCos_cos_b1", c1, b1);
+    c2 = std::cos(std::acos(b2));
+    output("saCos_cos_b2", c2, b2);
+
+
+    c1 = Carc_sin(a1);
+    c2 = std::asin(a2);
+    output("aSin", c1, c2);
+    c1 = Carc_sin(b1);
+    c2 = std::asin(b2);
+    output("saSin", c1, c2);
+
+    c1 = Carc_tan(a1);
+    c2 = std::atan(a2);
+    output("aTan", c1, c2);
+    c1 = Carc_tan(b1);
+    c2 = std::atan(b2);
+    output("saTan", c1, c2);
 
     
     return 0;