From 6a320264c2de3d6dd8cc1d1327b3c30df4c8cb26 Mon Sep 17 00:00:00 2001 From: Siddhesh Wani Date: Mon, 25 May 2015 14:46:31 +0530 Subject: Original Version --- 2.3-1/src/c/operations/division/cldivs.c | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) create mode 100644 2.3-1/src/c/operations/division/cldivs.c (limited to '2.3-1/src/c/operations/division/cldivs.c') diff --git a/2.3-1/src/c/operations/division/cldivs.c b/2.3-1/src/c/operations/division/cldivs.c new file mode 100644 index 00000000..8ada5da4 --- /dev/null +++ b/2.3-1/src/c/operations/division/cldivs.c @@ -0,0 +1,27 @@ +/* + * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab + * Copyright (C) 2008-2008 - INRIA - Arnaud TORSET + * + * This file must be used under the terms of the CeCILL. + * This source file is licensed as described in the file COPYING, which + * you should have received as part of this distribution. The terms + * are also available at + * http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt + * + */ + + +#include "division.h" +#include "conj.h" +#include "multiplication.h" +/* + A\B = B*A^-1 = B*conj(A)/|A|^2 + |A|^2 = Re(A)^2 + Im(A)^2 + */ + +floatComplex cldivs (floatComplex in1, floatComplex in2){ + floatComplex fC1, fC2; + fC1 = cmuls(in2,cconjs(in1)); + fC2 = FloatComplex(creals(in1)*creals(in1)+cimags(in1)*cimags(in1),0); + return crdivs(fC1,fC2); +} -- cgit