numpy.cc

Go to the documentation of this file.

#include "numpy.hh"

double interp(double v, double* x, double* y, int n)
{  int i;
   for(i=0; i<n; ++i)if(v<x[i])break;
   if(i==0)return y[0];
   if(i==n)return y[n-1];
   return y[i-1] + (v-x[i-1])/(x[i] - x[i-1]) *(y[i] - y[i-1]);
}

double interp(double v, double* x, double* y, int n, double left, double right)
{  int i;
   for(i=0; i<n; ++i)if(v<x[i])break;
   if(i==0)return left;
   if(i==n)return right;
   return y[i-1] + (v-x[i-1])/(x[i] - x[i-1]) *(y[i] - y[i-1]);
}

double* interp(double* v, int Nv, double* x, double* y, int n, double left, double right)
{  double* res = new double[Nv];
   for(int i=0; i<Nv; ++i)res[i] = interp(v[i], x, y, n, left, right);
   return res;
}

Complex operator+(const Complex& x, const Complex& y)
{  return Complex(x.Re() + y.Re(), x.Im() + y.Im());
}

Complex operator*(const Complex& x, const Complex& y)
{  return Complex(x.Re()*y.Re() - x.Im()*y.Im(), x.Re()*y.Im() + x.Im()*y.Re());
}

Complex Exp(double phase)
{  return Complex(cos(phase), sin(phase));
}

double fabs(const Complex& x)
{  return sqrt(x.Re()*x.Re() + x.Im()*x.Im());
}