IMPL: C64 type (#35)

Axect · Jan 17, 2021 · a1fb540 · a1fb540
1 parent 97fcd46
commit a1fb540
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 0 deletions.
diff --git a/src/fuga/mod.rs b/src/fuga/mod.rs
@@ -125,6 +125,7 @@ pub use crate::structure::{
  vector::*, 
  dataframe::*,
  ad::*,
+ complex::C64,
 };
 
 pub use crate::util::{api::*, low_level::*, non_macro::*, print::*, useful::*, wrapper::*};

diff --git a/src/prelude/mod.rs b/src/prelude/mod.rs
@@ -128,6 +128,7 @@ pub use crate::structure::{
  polynomial::*,
  vector::*,
  dataframe::*,
+ complex::C64,
 };
 
 pub use simpler::{solve, SimplerLinearAlgebra};

diff --git a/src/structure/complex.rs b/src/structure/complex.rs
@@ -3,6 +3,8 @@ use crate::traits::fp::FPVector;
 use crate::traits::math::{Vector, Normed, Norm, InnerProduct};
 use crate::traits::sugar::VecOps;
 
+pub type C64 = Complex<f64>;
+
 impl Vector for Complex<f64> {
  type Scalar = Self;
 
@@ -121,4 +123,10 @@ impl Normed for Vec<Complex<f64>> {
  }
 }
 
+impl InnerProduct for Vec<Complex<f64>> {
+ fn dot(&self, rhs: &Self) -> Self::Scalar {
+ self.zip_with(|x, y| x.conj() * y, rhs).sum()
+ }
+}
+
 impl VecOps for Vec<Complex<f64>> {}