summaryrefslogtreecommitdiff
path: root/src/c/matrixOperations/hilb/dhilba.c
diff options
context:
space:
mode:
authorSiddhesh Wani2015-05-25 14:46:31 +0530
committerSiddhesh Wani2015-05-25 14:46:31 +0530
commitdb464f35f5a10b58d9ed1085e0b462689adee583 (patch)
treede5cdbc71a54765d9fec33414630ae2c8904c9b8 /src/c/matrixOperations/hilb/dhilba.c
downloadScilab2C_fossee_old-db464f35f5a10b58d9ed1085e0b462689adee583.tar.gz
Scilab2C_fossee_old-db464f35f5a10b58d9ed1085e0b462689adee583.tar.bz2
Scilab2C_fossee_old-db464f35f5a10b58d9ed1085e0b462689adee583.zip
Original Version
Diffstat (limited to 'src/c/matrixOperations/hilb/dhilba.c')
-rw-r--r--src/c/matrixOperations/hilb/dhilba.c43
1 files changed, 43 insertions, 0 deletions
diff --git a/src/c/matrixOperations/hilb/dhilba.c b/src/c/matrixOperations/hilb/dhilba.c
new file mode 100644
index 0000000..3a03c0f
--- /dev/null
+++ b/src/c/matrixOperations/hilb/dhilba.c
@@ -0,0 +1,43 @@
+/*
+ * 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 "hilb.h"
+
+void dhilba ( double* out, int size ){
+
+ int iIndex1 = 0;
+ int iIndex2 = 0;
+ double dblVal = size;
+ double dblTemp = 0;
+
+ for(iIndex1 = 0 ; iIndex1 < size ; iIndex1++)
+ {
+ if(iIndex1 != 0)
+ dblVal = ((size - iIndex1) * dblVal * (size + iIndex1)) / pow((double)iIndex1,2);
+
+ dblTemp = dblVal * dblVal;
+
+ out[iIndex1 * size + iIndex1] = dblTemp / ( 2 * iIndex1 + 1);
+
+ if(iIndex1 == size - 1)
+ break;
+
+ for(iIndex2 = iIndex1 + 1 ; iIndex2 < size ; iIndex2++)
+ {
+ dblTemp = -((size - iIndex2) * dblTemp * (size + iIndex2)) / pow((double)iIndex2,2);
+ out[iIndex1 * size + iIndex2] = dblTemp / (iIndex1 + iIndex2 + 1);
+ out[iIndex2 * size + iIndex1] = out[iIndex1 * size + iIndex2];
+ }
+ }
+
+
+}