/*
 * Copyright 2011 Steven Gribble
 *
 *  This file is the solution to an exercise problem posed during
 *  one of the UW CSE 333 lectures (333exercises).
 *
 *  333exercises is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  333exercises is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with 333exercises.  If not, see <http://www.gnu.org/licenses/>.
 */

#include "complex.h"

Complex ComplexAdd(Complex a, Complex b) {
  Complex res;

  res.real = a.real + b.real;
  res.imaginary = a.imaginary + b.imaginary;
  return res;
}

Complex ComplexSubtract(Complex a, Complex b) {
  Complex res;

  res.real = a.real - b.real;
  res.imaginary = a.imaginary - b.imaginary;
  return res;
}

Complex ComplexMultiply(Complex a, Complex b) {
  Complex res;

  res.real = (a.real * b.real) - (a.imaginary * b.imaginary);
  res.imaginary = (a.imaginary * b.real) - (a.real * b.imaginary);
  return res;
}

Complex ComplexDivide(Complex a, Complex b) {
  // This is the hard one.  See Wikipedia for the formula. ;)
  Complex res;

  res.real =
    ((a.real * b.real) + (a.imaginary * b.imaginary)) /
    ((b.real * b.real) + (b.imaginary * b.imaginary));

  res.imaginary =
    ((a.imaginary * b.real) - (a.real * b.imaginary)) /
    ((b.real * b.real) + (b.imaginary * b.imaginary));

  return res;
}