From 49d9026592b8aa6426181b578f110b390cf07a49 Mon Sep 17 00:00:00 2001
From: sunwen <sunwen@realcan.cn>
Date: Mon, 24 Apr 2023 11:21:31 +0800
Subject: [PATCH] Add acosd and acosd's unittest.

---
 src/AuroraDefs.h         |  3 +++
 src/Function1D.cpp       | 18 ++++++++++++++++++
 src/Function1D.h         |  2 ++
 test/Function1D_Test.cpp | 10 ++++++++++
 4 files changed, 33 insertions(+)

diff --git a/src/AuroraDefs.h b/src/AuroraDefs.h
index fe80a1e..2c4cad4 100644
--- a/src/AuroraDefs.h
+++ b/src/AuroraDefs.h
@@ -4,4 +4,7 @@
 //必须在mkl.h和Eigen的头之前，<complex>之后
 #define MKL_Complex16 std::complex<double>
 #include "mkl.h"
+
+#define PI 3.141592653589793238462
+
 #endif //AURORA_AURORADEFS_H
diff --git a/src/Function1D.cpp b/src/Function1D.cpp
index 709968f..2602f4e 100644
--- a/src/Function1D.cpp
+++ b/src/Function1D.cpp
@@ -385,3 +385,21 @@ Matrix Aurora::acos(const Matrix& aMatrix)
     vdAcos(size, matrixData, resultData);
     return Matrix::New(resultData, aMatrix);
 }
+
+Matrix Aurora::acosd(const Matrix& aMatrix)
+{
+    if(aMatrix.isComplex() || aMatrix.isNull())
+    {
+        return Matrix();
+    }
+
+    size_t size = aMatrix.getDataSize();
+    double* matrixData = aMatrix.getData();
+    double* resultData = Aurora::malloc(size);
+    vdAcos(size, matrixData, resultData);
+    for(size_t i=0; i<size; ++i)
+    {
+        resultData[i] = resultData[i] * 180 / PI;
+    }
+    return Matrix::New(resultData, aMatrix);
+}
diff --git a/src/Function1D.h b/src/Function1D.h
index 456b5cf..4fcfdbd 100644
--- a/src/Function1D.h
+++ b/src/Function1D.h
@@ -55,6 +55,8 @@ namespace Aurora {
 
     Matrix acos(const Matrix& aMatrix);
 
+    Matrix acosd(const Matrix& aMatrix);
+
     /**
      * 多项式计算
      * @brief 例如p[1 0 1],x[3 2 5],代表对多项式 y = x^2 + 1 求(x=3, x=2, x=5)时所有的y
diff --git a/test/Function1D_Test.cpp b/test/Function1D_Test.cpp
index ecf70ec..d8cb1ea 100644
--- a/test/Function1D_Test.cpp
+++ b/test/Function1D_Test.cpp
@@ -255,3 +255,13 @@ TEST_F(Function1D_Test, acos) {
     EXPECT_DOUBLE_AE(result.getData()[1],0.5917);
     EXPECT_DOUBLE_AE(result.getData()[2],1.1593);
 }
+
+TEST_F(Function1D_Test, acosd) {
+    double *data = new double[3]{0.02,0.83,0.4,};
+    auto matrix = Aurora::Matrix::fromRawData(data, 3);
+    auto result = Aurora::acosd(matrix);
+    result.printf();
+    EXPECT_DOUBLE_AE(result.getData()[0],88.854);
+    EXPECT_DOUBLE_AE(result.getData()[1],33.9013);
+    EXPECT_DOUBLE_AE(result.getData()[2],66.4218);
+}