summaryrefslogtreecommitdiff
path: root/src/matrixOperations/matrixMultiplication.c
diff options
context:
space:
mode:
authorjofret2008-06-19 14:45:46 +0000
committerjofret2008-06-19 14:45:46 +0000
commit6142e72914eb2d0b94a0cf9fef554cd3a6760c3f (patch)
tree3182b744582fa37b02afc7df5d0c9b28870109fa /src/matrixOperations/matrixMultiplication.c
parent7c9f4fd3f6cfa720f9f9bb6d47ac1641f4bc4f77 (diff)
downloadscilab2c-6142e72914eb2d0b94a0cf9fef554cd3a6760c3f.tar.gz
scilab2c-6142e72914eb2d0b94a0cf9fef554cd3a6760c3f.tar.bz2
scilab2c-6142e72914eb2d0b94a0cf9fef554cd3a6760c3f.zip
Reorganize and separate operations
Diffstat (limited to 'src/matrixOperations/matrixMultiplication.c')
-rw-r--r--src/matrixOperations/matrixMultiplication.c221
1 files changed, 0 insertions, 221 deletions
diff --git a/src/matrixOperations/matrixMultiplication.c b/src/matrixOperations/matrixMultiplication.c
deleted file mode 100644
index dbed4330..00000000
--- a/src/matrixOperations/matrixMultiplication.c
+++ /dev/null
@@ -1,221 +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
- *
- */
-
-#ifndef WITHOUT_BLAS
-#include "blas.h"
-#endif
-#include "matrixMultiplication.h"
-
-/*
-** \brief Compute a multiplication for floats matrixes.
-** \param in1 : input matrix.
-** \param lines1 : lines of in1 matrix.
-** \param columns1 : columns of in1 matrix.
-** \param in2 : input arry.
-** \param lines2 : lines of in2 matrix.
-** \param columns2 : columns of in2 matrix.
-** \param out : Matrix that contains the multiplication in1 * in2.
-*/
-void smulma(float *in1, int lines1, int columns1,
- float *in2, int lines2, int columns2,
- float *out)
-{
- int i = 0;
- int k = 0;
- float accu = 0;
-
- /*
- ** How to convert 2 index matrixes to one.
- ** #define in1(a, b) in1[a+b*lines1]
- ** #define in2(c, d) in2[c+d*lines2]
- */
-
- for (i = 0 ; i < lines1 * columns2 ; ++i)
- {
- accu = 0;
- for (k = 0; k < columns1 ; ++k)
- {
- accu += in1[i % lines1 + k * lines1]
- * in2[k + (i / lines1) * lines2];
- }
- out[i] = accu;
- }
-
-}
-
-/*
-** \brief Compute a multiplication for doubles matrixes.
-** \param in1 : input matrix.
-** \param lines1 : lines of in1 matrix.
-** \param columns1 : columns of in1 matrix.
-** \param in2 : input arry.
-** \param lines2 : lines of in2 matrix.
-** \param columns2 : columns of in2 matrix.
-** \param out : Matrix that contains the multiplication in1 * in2.
-*/
-void dmulma(double *in1, int lines1, int columns1,
- double *in2, int lines2, int columns2,
- double *out)
-{
-#ifndef WITHOUT_BLAS
- /*
- ** USES BLAS DGEMM FUNCTION.
- */
- double One = 1;
- double Zero = 0;
-
- /* Cr <- 1*Ar*Br + 0*Cr */
- dgemm_("N","N", &columns2, &columns2, &columns1, &One,
- in1 , &lines1, in2, &lines2, &Zero, out, &columns2);
-#else
- /*
- ** DO NOT USE ANY BLAS FUNCTION.
- */
- int i = 0;
- int k = 0;
- double accu = 0;
-
- /*
- ** How to convert 2 index matrixes to one.
- ** #define in1(a, b) in1[a+b*lines1]
- ** #define in2(c, d) in2[c+d*lines2]
- */
-
- for (i = 0 ; i < lines1 * columns2 ; ++i)
- {
- accu = 0;
- for (k = 0; k < columns1 ; ++k)
- {
- accu += in1[i % lines1 + k * lines1]
- * in2[k + (i / lines1) * lines2];
- }
- out[i] = accu;
- }
-#endif
-}
-
-/*
-** \brief Compute a multiplication for floats complex matrixes.
-** \param in1 : input matrix.
-** \param lines1 : lines of in1 matrix.
-** \param columns1 : columns of in1 matrix.
-** \param in2 : input arry.
-** \param lines2 : lines of in2 matrix.
-** \param columns2 : columns of in2 matrix.
-** \param out : Matrix that contains the multiplication in1 * in2.
-*/
-void cmulma(floatComplex *in1, int lines1, int columns1,
- floatComplex *in2, int lines2, int columns2,
- floatComplex *out)
-{
- int i = 0;
- int k = 0;
- floatComplex accu = FloatComplex(0, 0);
-
- for (i = 0 ; i < lines1 * columns2 ; ++i)
- {
- accu = FloatComplex(0,0);
- for (k = 0; k < columns1 ; ++k)
- {
- accu = cadds(accu,
- ctimess(in1[i % lines1 + k *lines1] ,
- in2[k + (i / lines1) *lines2] ));
- }
- out[i] = accu;
- }
-}
-
-/*
-** \brief Compute a multiplication for doubles matrixes.
-** \param in1 : input matrix.
-** \param lines1 : lines of in1 matrix.
-** \param columns1 : columns of in1 matrix.
-** \param in2 : input arry.
-** \param lines2 : lines of in2 matrix.
-** \param columns2 : columns of in2 matrix.
-** \param out : Matrix that contains the multiplication in1 * in2.
-*/
-void zmulma(doubleComplex *in1, int lines1, int columns1,
- doubleComplex *in2, int lines2, int columns2,
- doubleComplex *out)
-{
-#ifndef WITHOUT_BLAS
- /*
- ** USES BLAS DGEMM FUNCTION.
- */
- int i = 0;
- double One = 1;
- double MinusOne = -1;
- double Zero = 0;
-
- double *in1Real = malloc((uint) lines1 * (uint) columns1 * sizeof(double));
- double *in1Imag = malloc((uint) lines1 * (uint) columns1 * sizeof(double));
- double *in2Real = malloc((uint) lines2 * (uint) columns2 * sizeof(double));
- double *in2Imag = malloc((uint) lines2 * (uint) columns2 * sizeof(double));
-
- double *RealOut = malloc((uint) lines1 * (uint) columns2 * sizeof(double));
- double *ImagOut = malloc((uint) lines1 * (uint) columns2 * sizeof(double));
-
- zreala(in1, lines1 * columns1, in1Real);
- zreala(in2, lines2 * columns2, in2Real);
- zimaga(in1, lines1 * columns1, in1Imag);
- zimaga(in2, lines2 * columns2, in2Imag);
-
- /* Cr <- 1*Ar*Br + 0*Cr */
- dgemm_("N","N", &lines1, &columns2, &columns1, &One,
- in1Real, &lines1, in2Real, &lines2, &Zero, RealOut, &lines1);
-
- /* Cr <- -1*Ai*Bi + 1*Cr */
- dgemm_("N","N", &lines1, &columns2, &columns1, &MinusOne,
- in1Imag, &lines1, in2Imag, &lines2, &One, RealOut, &lines1);
-
- /* Ci <- 1*Ar*Bi + 0*Ci */
- dgemm_("N","N", &lines1, &columns2, &columns1, &One,
- in1Real, &lines1, in2Imag, &lines2, &Zero, ImagOut, &lines1);
-
- /* Ci <- 1*Ai*Br + 1*Ci */
- dgemm_("N","N", &lines1, &columns2, &columns1, &One,
- in1Imag, &lines1, in2Real, &lines2, &One, ImagOut, &lines1);
-
- /* Now fill output matrix */
- for(i = 0 ; i < lines1 * columns2 ; ++i)
- {
- out[i] = DoubleComplex(RealOut[i], ImagOut[i]);
- }
- /* FREE allocated variables */
- free(in1Real);
- free(in2Real);
- free(in1Imag);
- free(in2Imag);
- free(RealOut);
- free(ImagOut);
-#else
- /*
- ** DO NOT USE ANY BLAS FUNCTION.
- */
- int i = 0;
- int k = 0;
- doubleComplex accu = DoubleComplex(0, 0);
-
- for (i = 0 ; i < lines1 * columns2 ; ++i)
- {
- accu = DoubleComplex(0,0);
- for (k = 0; k < columns1 ; ++k)
- {
- accu = zadds(accu,
- ztimess(in1[i % lines1 + k *lines1] ,
- in2[k + (i / lines1) *lines2] ));
- }
- out[i] = accu;
- }
-#endif
-}
generated by cgit (git 2.34.1) at 2025-03-21 16:04:53 +0000