summaryrefslogtreecommitdiff
path: root/src/operations/subtraction/testSubtraction.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/operations/subtraction/testSubtraction.c')
-rw-r--r--src/operations/subtraction/testSubtraction.c132
1 files changed, 0 insertions, 132 deletions
diff --git a/src/operations/subtraction/testSubtraction.c b/src/operations/subtraction/testSubtraction.c
deleted file mode 100644
index ce29e2a3..00000000
--- a/src/operations/subtraction/testSubtraction.c
+++ /dev/null
@@ -1,132 +0,0 @@
-/*
- * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
- * Copyright (C) 2008-2008 - INRIA - Bruno JOFRET
- *
- * 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 <stdlib.h>
-#include <stdio.h>
-#include <assert.h>
-#include <math.h>
-#include "subtraction.h"
-
-#define SIZE 10000
-
-/* #define LOCAL_DEBUG */
-
-static void zdiffaTest(void) {
- doubleComplex M1[SIZE];
- doubleComplex M2[SIZE];
- doubleComplex M1_but_M2[SIZE];
-
- int i = 0;
-
- srand(1);
- for (i = 0; i < SIZE; ++i) {
- M1[i] = DoubleComplex((double) rand(), (double) rand());
- M2[i] = DoubleComplex((double) rand(), (double) rand());
- }
- zdiffa(M1, SIZE, M2, SIZE, M1_but_M2);
- for (i = 0; i < SIZE; ++i) {
-#ifdef LOCAL_DEBUG
- printf("M1_but_M2 = %e + %e i\n", zreals(M1_but_M2[i]), zimags(M1_but_M2[i]));
-#endif
- assert(zreals(M1_but_M2[i]) - (zreals(M1[i]) - zreals(M2[i])) == 0);
- assert(zimags(M1_but_M2[i]) - (zimags(M1[i]) - zimags(M2[i])) == 0);
- }
-}
-
-static void cdiffaTest(void) {
- floatComplex M1[SIZE];
- floatComplex M2[SIZE];
- floatComplex M1_but_M2[SIZE];
-
- int i = 0;
-
- srand(1);
- for (i = 0; i < SIZE; ++i) {
- M1[i] = FloatComplex((float) rand(), (float) rand());
- M2[i] = FloatComplex((float) rand(), (float) rand());
- }
- cdiffa(M1, SIZE, M2, SIZE, M1_but_M2);
- for (i = 0; i < SIZE; ++i) {
-#ifdef LOCAL_DEBUG
- printf("M1_but_M2 = %e + %e i\n", creals(M1_but_M2[i]), cimags(M1_but_M2[i]));
-#endif
- assert(fabsf(creals(M1_but_M2[i]) - (creals(M1[i]) - creals(M2[i]))) / creals(M1_but_M2[i]) < 1e-07);
- assert(fabsf(cimags(M1_but_M2[i]) - (cimags(M1[i]) - cimags(M2[i]))) / cimags(M1_but_M2[i]) < 1e-07);
- }
-}
-
-static void sdiffaTest(void) {
- float M1[SIZE];
- float M2[SIZE];
- float M1_but_M2[SIZE];
-
- int i = 0;
-
- srand(1);
- for (i = 0; i < SIZE; ++i) {
- M1[i] = (float) rand();
- M2[i] = (float) rand();
- }
- sdiffa(M1, SIZE, M2, SIZE, M1_but_M2);
- for (i = 0; i < SIZE; ++i) {
-#ifdef LOCAL_DEBUG
- printf("M1_but_M2 = %e\n", M1_but_M2[i]);
- printf("M1[i] - M2[i] = %e\n", (M1[i] - M2[i]));
- printf("Error = %e\n",(float) fabsf(M1_but_M2[i] - (M1[i] - M2[i])) / M1_but_M2[i]);
-#endif
- assert(fabsf(M1_but_M2[i] - (M1[i] - M2[i])) / M1_but_M2[i] < 1e-07);
- }
-}
-
-static void ddiffaTest(void) {
- double M1[SIZE];
- double M2[SIZE];
- double M1_but_M2[SIZE];
-
- int i = 0;
-
- srand(1);
- for (i = 0; i < SIZE; ++i) {
- M1[i] = (double) rand();
- M2[i] = (double) rand();
- }
- ddiffa(M1, SIZE, M2, SIZE, M1_but_M2);
- for (i = 0; i < SIZE; ++i) {
-#ifdef LOCAL_DEBUG
- printf("M1_but_M2 = %e\n", M1_but_M2[i]);
- printf("M1[i] = %e\n", M1[i]);
- printf("M2[i] = %e\n", M2[i]);
- printf("M1[i] - M2[i] = %e\n", (M1[i] - M2[i]));
- printf("Error = %e\n", fabs(M1_but_M2[i] - (M1[i] - M2[i])));
-#endif
- assert(fabs(M1_but_M2[i] - (M1[i] - M2[i])) == 0);
- }
-
-}
-
-static int testSubtraction(void) {
-
- printf("\n>>>> Subtraction Tests\n");
- sdiffaTest();
- ddiffaTest();
- cdiffaTest();
- zdiffaTest();
-
- return 0;
-}
-
-
-
-int main(void) {
- assert(testSubtraction() == 0);
- return 0;
-}