/* $Id: fakecomplex.c,v 1.1 2004/04/24 20:38:30 geuzaine Exp $ */

#include "fakecomplex.h"

void cotofl(complex a, double b , double c)
{
  b = a.r;
  c = a.i;
}

complex fltoco(double a, double b)
{
  complex s;

  s.r = a;
  s.i = b;
  return(s);
}

complex Csum(complex a, complex b)
{
  complex s;

  s.r = a.r + b.r;
  s.i = a.i + b.i;
  return(s);
}

complex Csub(complex a, complex b)
{
  complex s;

  s.r = a.r - b.r;
  s.i = a.i - b.i;
  return(s);
}

complex Cprod(complex a, complex b)
{
  complex s;

  s.r = a.r * b.r - a.i * b.i;
  s.i = a.r * b.i + a.i * b.r;
  return(s);
}

complex Cdiv(complex a, complex b)
{
  complex s;
  double den;

  den = b.r * b.r + b.i * b.i;
  s.r = (a.r * b.r + a.i * b.i) / den;
  s.i = (a.i * b.r - a.r * b.i) / den;
  return(s);
}

complex Cconj(complex a)
{
  complex s;

  s.r = a.r;
  s.i = -a.i;
  return(s);
}

complex Cneg(complex a)
{
  complex s;

  s.r = -a.r;
  s.i = -a.i;
  return(s);
}

double Cmodu(complex a)
{
  return(sqrt(a.r * a.r + a.i * a.i));
}

complex Cpow(complex a, double b)
{
  complex s;
  double mod, arg;

  mod = a.r * a.r + a.i * a.i;
  arg = atan2(a.i,a.r);
  mod = pow(mod,0.5*b);
  arg *= b;
  s.r = mod * cos(arg);
  s.i = mod * sin(arg);

  return(s);
}

complex Cprodr(double a, complex b)
{
  complex s;

  s.r = a * b.r;
  s.i = a * b.i;
  return(s);
}

void Cisum(complex *a, complex b)
{
  (*a).r += b.r;
  (*a).i += b.i;
}

void Cisub(complex *a, complex b)
{
  (*a).r -= b.r;
  (*a).i -= b.i;
}

void Ciprod(complex *a, complex b)
{
  complex aa,d;

  aa.r=(*a).r;
  aa.i=(*a).i;

  d.r = aa.r * b.r - aa.i * b.i;
  d.i = aa.r * b.i + aa.i * b.r;

  (*a).r = d.r;
  (*a).i = d.i;
}