diff options
Diffstat (limited to 'src/matrixOperations/hilbert/shilba.c')
-rw-r--r-- | src/matrixOperations/hilbert/shilba.c | 44 |
1 files changed, 0 insertions, 44 deletions
diff --git a/src/matrixOperations/hilbert/shilba.c b/src/matrixOperations/hilbert/shilba.c deleted file mode 100644 index bd3b5ee5..00000000 --- a/src/matrixOperations/hilbert/shilba.c +++ /dev/null @@ -1,44 +0,0 @@ -/* - * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab - * Copyright (C) 2008-2008 - INRIA - Allan SIMON - * - * 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 "hilbert.h" - -void shilba ( float* out, int size ){ - - int iIndex1 = 0; - int iIndex2 = 0; - float dblVal = (float) size; - float dblTemp = 0; - - for(iIndex1 = 0 ; iIndex1 < size ; iIndex1++) - { - if(iIndex1 != 0) - dblVal = ((float)(size - iIndex1) * dblVal * (float)(size + iIndex1)) /(float) pow((double)iIndex1,2); - - dblTemp = dblVal * dblVal; - - out[iIndex1 * size + iIndex1] =(float) dblTemp /(float) ( 2 * iIndex1 + 1); - - if(iIndex1 == size - 1) - break; - - for(iIndex2 = iIndex1 + 1 ; iIndex2 < size ; iIndex2++) - { - dblTemp = -((float)(size - iIndex2) * dblTemp *(float) (size + iIndex2)) /(float) pow((double)iIndex2,2); - out[iIndex1 * size + iIndex2] = dblTemp /(float) (iIndex1 + iIndex2 + 1); - out[iIndex2 * size + iIndex1] = out[iIndex1 * size + iIndex2]; - } - } - - -} |