added missing file OperatorUtils.hh

Utils/OperatorUtils.hh

+#pragma once
+
+#include <memory>
+#include <complex>
+#include <codi.hpp>
+#include "ErrLogger/ErrLogger.hh"
+
+std::complex<codi::RealReverse> operator*(const std::complex<codi::RealReverse>& realRev, const double inputDouble)
+{
+  std::complex<codi::RealReverse> result(realRev.real()*inputDouble, realRev.imag()*inputDouble);
+  return result;
+}
+
+std::complex<codi::RealReverse> operator*(const double inputDouble, const std::complex<codi::RealReverse>& realRev)
+{
+  std::complex<codi::RealReverse> result(realRev.real()*inputDouble, realRev.imag()*inputDouble);
+  return result;
+}
+
+std::complex<codi::RealReverse> operator*(const std::complex<codi::RealReverse>& realRev, const std::complex<double>& inputDoubleCompl)
+{
+  std::complex<codi::RealReverse> result(realRev.real()*inputDoubleCompl.real()-realRev.imag()*inputDoubleCompl.imag(), realRev.real()*inputDoubleCompl.imag()+realRev.imag()*inputDoubleCompl.real());
+  return result;
+}
+std::complex<codi::RealReverse> operator*(const std::complex<double>& inputDoubleCompl, const std::complex<codi::RealReverse>& realRev)
+{
+  std::complex<codi::RealReverse> result(realRev.real()*inputDoubleCompl.real()-realRev.imag()*inputDoubleCompl.imag(), realRev.real()*inputDoubleCompl.imag()+realRev.imag()*inputDoubleCompl.real());
+  return result;
+}
+
+std::complex<codi::RealReverse> operator/(const std::complex<codi::RealReverse>& realRev, const double inputDouble)
+{
+  std::complex<codi::RealReverse> result(realRev.real()/inputDouble, realRev.imag()/inputDouble);
+  return result;
+}
+
+std::complex<codi::RealReverse> operator/(const std::complex<codi::RealReverse>& realRev, const std::complex<double> inputDoubleCompl)
+{
+  double theDenom=norm(inputDoubleCompl);
+  if(theDenom < 1e-20){
+    Alert << "Division by 0; devision of the complex number " << inputDoubleCompl << " is not allowed!!!" << endmsg;
+    exit(1);
+  }
+  std::complex<codi::RealReverse> result((realRev.real()+inputDoubleCompl.real()+realRev.imag()+inputDoubleCompl.imag())/theDenom, (realRev.real()+inputDoubleCompl.imag()-realRev.imag()+inputDoubleCompl.real())/theDenom);
+
+  return result;
+}
+
+std::complex<codi::RealReverse> operator+(const std::complex<codi::RealReverse>& realRev, const std::complex<double>& inputDoubleCompl)
+{
+std::complex<codi::RealReverse> result(realRev.real()+inputDoubleCompl.real(), realRev.imag()+inputDoubleCompl.imag());
+  return result;
+  }
+
+std::complex<codi::RealReverse> operator+(const std::complex<double>& inputDoubleCompl, const std::complex<codi::RealReverse>& realRev)
+{
+  std::complex<codi::RealReverse> result(realRev.real()+inputDoubleCompl.real(), realRev.imag()+inputDoubleCompl.imag());
+  return result;
+}
+
+std::complex<codi::RealReverse> operator-(const std::complex<codi::RealReverse>& realRev, const std::complex<double>& inputDoubleCompl)
+{
+  std::complex<codi::RealReverse> result(realRev.real()-inputDoubleCompl.real(), realRev.imag()-inputDoubleCompl.imag());
+  return result;
+}
+
+std::complex<codi::RealReverse> operator-(const std::complex<double>& inputDoubleCompl, const std::complex<codi::RealReverse>& realRev)
+{
+  std::complex<codi::RealReverse> result(realRev.real()-inputDoubleCompl.real(), realRev.imag()-inputDoubleCompl.imag());
+  return result;
+}
+
+