/* * 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 . */ #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; }