tmp

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;