From 574ddf08c208a2d1b8c27fe29525f631816c32d5 Mon Sep 17 00:00:00 2001 From: Abhinav Dronamraju Date: Mon, 17 Jul 2017 21:57:08 +0530 Subject: Median added --- src/c/statisticsFunctions/median/zmediana.c | 64 +++++++++++++++++++++++++++++ 1 file changed, 64 insertions(+) create mode 100644 src/c/statisticsFunctions/median/zmediana.c (limited to 'src/c/statisticsFunctions/median/zmediana.c') diff --git a/src/c/statisticsFunctions/median/zmediana.c b/src/c/statisticsFunctions/median/zmediana.c new file mode 100644 index 00000000..32726e77 --- /dev/null +++ b/src/c/statisticsFunctions/median/zmediana.c @@ -0,0 +1,64 @@ +/* Copyright (C) 2016 - IIT Bombay - FOSSEE + + 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 + Author: Abhinav Dronamraju + Organization: FOSSEE, IIT Bombay + Email: toolbox@scilab.in +*/ + + +#include "median.h" +#include "types.h" +#include "uint16.h" +#include "doubleComplex.h" +#include "addition.h" +#include "division.h" +#include "abs.h" + +doubleComplex zmediana(doubleComplex *in, int size) +{ + doubleComplex a; doubleComplex fin; doubleComplex middle; + + + for (int i = 0; i < size; ++i) + + { + + for (int j = i + 1; j < size; ++j) + + { + + if (zabss(in[i]) > zabss(in[j])) + + { + + a = in[i]; + + in[i] = in[j]; + + in[j] = a; + + } + + } + + } + + + + if(size%2 ==0) + { + middle= zadds(in[size/2], in[(size/2)-1]); + fin= zrdivs(middle, DoubleComplex(2,0)); + } + else + { + fin= in[(size-1)/2]; + } + + return fin; +} -- cgit From 6422b5c0a45161d6571fd6d7077b6d551cbb5e26 Mon Sep 17 00:00:00 2001 From: Abhinav Dronamraju Date: Thu, 10 Aug 2017 20:30:45 +0530 Subject: rank and gsort added --- src/c/statisticsFunctions/median/zmediana.c | 16 +++++++++------- 1 file changed, 9 insertions(+), 7 deletions(-) (limited to 'src/c/statisticsFunctions/median/zmediana.c') diff --git a/src/c/statisticsFunctions/median/zmediana.c b/src/c/statisticsFunctions/median/zmediana.c index 32726e77..36804562 100644 --- a/src/c/statisticsFunctions/median/zmediana.c +++ b/src/c/statisticsFunctions/median/zmediana.c @@ -21,8 +21,10 @@ doubleComplex zmediana(doubleComplex *in, int size) { - doubleComplex a; doubleComplex fin; doubleComplex middle; + doubleComplex a; doubleComplex fin; doubleComplex middle; doubleComplex in_copy[size]; +for(int i=0; i< size; i++) + in_copy[i]= in[i]; for (int i = 0; i < size; ++i) @@ -32,15 +34,15 @@ doubleComplex zmediana(doubleComplex *in, int size) { - if (zabss(in[i]) > zabss(in[j])) + if (zabss(in_copy[i]) > zabss(in_copy[j])) { - a = in[i]; + a = in_copy[i]; - in[i] = in[j]; + in_copy[i] = in_copy[j]; - in[j] = a; + in_copy[j] = a; } @@ -52,12 +54,12 @@ doubleComplex zmediana(doubleComplex *in, int size) if(size%2 ==0) { - middle= zadds(in[size/2], in[(size/2)-1]); + middle= zadds(in_copy[size/2], in_copy[(size/2)-1]); fin= zrdivs(middle, DoubleComplex(2,0)); } else { - fin= in[(size-1)/2]; + fin= in_copy[(size-1)/2]; } return fin; -- cgit