From b39a4d593115f01e8d5815b21c4c9e26acbedcb1 Mon Sep 17 00:00:00 2001
From: Yi Zhang <zhangyiss@icloud.com>
Date: Wed, 23 Jul 2025 11:35:38 +0800
Subject: [PATCH] tmp

---
 src/core/linear_algebra_ex.cpp | 29 +++++++++++++++++++++++++++++
 1 file changed, 29 insertions(+)

diff --git a/src/core/linear_algebra_ex.cpp b/src/core/linear_algebra_ex.cpp
index 22b92e9..9b167cb 100644
--- a/src/core/linear_algebra_ex.cpp
+++ b/src/core/linear_algebra_ex.cpp
@@ -31,6 +31,34 @@ using namespace gctl;
 
 int main(int argc, char const *argv[])
 {
+    _1cd_array a(10, std::complex<double>(1.0, 2.0));
+    std::cout << "a = "; a.show();
+
+    _1cd_array b(10);
+    veccpy(b, a, std::complex<double>(-0.5, 0.0));
+    std::cout << "b = "; b.show();
+
+    _1cd_array c(10);
+    vecadd(c, a, b, std::complex<double>(1.0, 0.0), std::complex<double>(1.0, 0.0));
+    std::cout << "a + b = "; c.show();
+
+    vecdiff(c, a, b, std::complex<double>(1.0, 0.0), std::complex<double>(1.0, 0.0));
+    std::cout << "a - b = "; c.show();
+
+    veccpy(c, a, std::complex<double>(1.0, 0.0));
+    vecapp(c, a, std::complex<double>(0.5, 0.0));
+    std::cout << "a += 0.5*a "; c.show();
+
+    veccpy(c, a, std::complex<double>(1.0, 0.0));
+    vecsub(c, a, std::complex<double>(0.5, 0.0));
+    std::cout << "a -= 0.5*a "; c.show();
+
+    std::complex<double> d = vecdot(a, b);
+    std::complex<double> e = vecinner(a, b);
+    std::cout << "d = " << d << "\n";
+    std::cout << "e = " << e << "\n";
+
+    /*
     _1d_array a(10);
     _1d_array b(10, 1.0);
     a.sequence(1.0, 1.0);
@@ -59,6 +87,7 @@ int main(int argc, char const *argv[])
     _2d_matrix m(10, 10, 1.0);
     matvec(c, m, a);
     c.show();
+    */
 
     /*
     int rows = 20000;