addmatrix mit broadcasten

makefile

@@ -57,12 +57,9 @@ imageInputTests: imageInput.o imageInputTests.c $(unityfolder)/unity.c
 # --------------------------
 # Clean
 # --------------------------
-#clean:
-#ifeq ($(OS),Windows_NT)
-#	del /f *.o *.exe
-#else
-#	rm -f *.o mnist runMatrixTests runNeuralNetworkTests runImageInputTests
-#endif
-# clean für windows
 clean:
-	rm -f *.o *.exe 
+ifeq ($(OS),Windows_NT)
+	del /f *.o *.exe
+else
+	rm -f *.o mnist runMatrixTests runNeuralNetworkTests runImageInputTests
+endif

matrix.c

@@ -1,6 +1,8 @@
 #include "matrix.h"
+#include <stdio.h>
 #include <stdlib.h>
 #include <string.h>
+
 // TODO Matrix-Funktionen implementieren
 /*typedef struct {
     unsigned int rows;  //Zeilen
@@ -8,6 +10,12 @@
     MatrixType *buffer; //Zeiger auf Speicherbereich Reihen*Spalten
 } Matrix;*/
 Matrix createMatrix(unsigned int rows, unsigned int cols) {
+
+  Matrix errorMatrix = {0, 0, NULL};
+  if (rows == 0 || cols == 0) {
+
+    return errorMatrix;
+  }
   MatrixType *buffer =
       malloc(rows * cols * sizeof(MatrixType)); // Speicher reservieren, malloc
                                                 // liefert Zeiger auf Speicher
@@ -20,6 +28,7 @@ void clearMatrix(Matrix *matrix) {
   matrix->cols = UNDEFINED_MATRIX_VALUE;
   free((*matrix).buffer); // Speicher freigeben
 }
+
 void setMatrixAt(const MatrixType value, Matrix matrix,
                  const unsigned int rowIdx, // Kopie der Matrix wird übergeben
                  const unsigned int colIdx) {
@@ -40,26 +49,168 @@ MatrixType getMatrixAt(const Matrix matrix,
 
   return value;
 }
+
+Matrix broadCastCols(const Matrix matrix, const unsigned int rows,
+                     const unsigned int cols) {
+
+  Matrix copy = createMatrix(
+      rows, cols); // Matrix 1 Kopie erstellen mit Dimensionen von Matrix2
+  for (int r = 0; r < rows; r++) {
+
+    MatrixType value = getMatrixAt(matrix, r, 0);
+    for (int c = 0; c < cols; c++) {
+
+      setMatrixAt(value, copy, r, c);
+    }
+  }
+
+  return copy;
+}
+Matrix broadCastRows(const Matrix matrix, const unsigned int rows,
+                     const unsigned int cols) {
+
+  Matrix copy = createMatrix(rows, cols);
+
+  for (int c = 0; c < cols; c++) {
+    MatrixType value = getMatrixAt(matrix, 0, c);
+
+    for (int r = 0; r < rows; r++) {
+
+      setMatrixAt(value, copy, r, c);
+    }
+  }
+
+  return copy;
+}
+
 Matrix add(const Matrix matrix1, const Matrix matrix2) {
 
-  // Ergebnismatrix
+  // Broadcasting nur bei Vektor und Matrix, Fehlermeldung bei zwei unpassenden
-  Matrix result;
+  // Matrizen
 
-  // Broadcasting nur bei Vektor und Matrix, Fehlermeldung bei zwei unpassender
+  const unsigned int rows1 = matrix1.rows;
-  // Matrix
+  const unsigned int rows2 = matrix2.rows;
-  if (matrix1.rows != matrix2.rows) {
+  const unsigned int cols1 = matrix1.cols;
+  const unsigned int cols2 = matrix2.cols;
 
-    // check, which one is smaller
+  const int rowsEqual = ((rows1 == rows2) ? 1 : 0);
-    // realloc
+
+  const int colsEqual = ((cols1 == cols2) ? 1 : 0);
+
+  if (rowsEqual && colsEqual) // addieren
+
+  {
+    Matrix result = createMatrix(rows1, cols1); // Speicher reservieren
+
+    for (int i = 0; i < (rows1 * cols1); i++) { // addieren
+
+      result.buffer[i] =
+          (matrix1.buffer[i] +
+           matrix2.buffer[i]); // buffer[i]  ⇔  *(buffer + i) Adresse =
+                               // Startadresse + (i * sizeof(MatrixType))
+    }
+
+    return result; // zurückgeben
   }
 
-  if (matrix1.cols != matrix2.cols) {
+  else if (rowsEqual && !colsEqual) {
+
+    if (cols1 == 1) {
+
+      Matrix result = createMatrix(rows2, cols2);
+
+      Matrix copy1 = broadCastCols(matrix1, rows2, cols2);
+
+      for (int i = 0; i < (rows2 * cols2); i++) { // addieren
+
+        result.buffer[i] =
+            (copy1.buffer[i] +
+             matrix2.buffer[i]); // buffer[i]  ⇔  *(buffer + i) Adresse =
+                                 // Startadresse + (i * sizeof(MatrixType))
+      }
+      return result;
+
+      // add und return
+
+    } else if (cols2 == 1) {
+
+      Matrix result = createMatrix(rows1, cols1);
+
+      Matrix copy2 = broadCastCols(matrix2, rows1, cols1);
+
+      for (int i = 0; i < (rows1 * cols1); i++) { // addieren
+
+        result.buffer[i] =
+            (matrix1.buffer[i] +
+             copy2.buffer[i]); // buffer[i]  ⇔  *(buffer + i) Adresse =
+                               // Startadresse + (i * sizeof(MatrixType))
+      }
+      return result;
+      // add und return
+
+    }
+
+    else {
+
+      printf("Fehlermeldung"); // vielleicht Fehlermeldung ändern zu
+                               // Programmabbruch
+      Matrix error = {0, 0, NULL};
+      return error;
+    }
+
   }
 
-  // Speicher reservieren
+  else if (!rowsEqual && colsEqual) {
 
-  // Matrix addieren
+    if (rows1 == 1) {
 
-  return result;
+      Matrix result = createMatrix(rows2, cols2);
+
+      Matrix copy1 = broadCastRows(matrix1, rows2, cols2);
+
+      for (int i = 0; i < (rows2 * cols2); i++) { // addieren
+
+        result.buffer[i] =
+            (copy1.buffer[i] +
+             matrix2.buffer[i]); // buffer[i]  ⇔  *(buffer + i) Adresse =
+                                 // Startadresse + (i * sizeof(MatrixType))
+      }
+      return result;
+
+      // add und return
+
+    } else if (rows2 == 1) {
+
+      Matrix result = createMatrix(rows1, cols1);
+
+      Matrix copy2 = broadCastCols(matrix2, rows1, cols1);
+      // add und return
+
+      for (int i = 0; i < (rows1 * cols1); i++) { // addieren
+
+        result.buffer[i] =
+            (matrix1.buffer[i] +
+             copy2.buffer[i]); // buffer[i]  ⇔  *(buffer + i) Adresse =
+                               // Startadresse + (i * sizeof(MatrixType))
+      }
+      return result;
+
+    }
+
+    else {
+
+      printf("Fehlermeldung"); // vielleicht Fehlermeldung ändern zu
+                               // Programmabbruch
+      Matrix error = {0, 0, NULL};
+      return error;
+    }
+  }
+
+  else {
+    printf(
+        "Fehlermeldung"); // vielleicht Fehlermeldung ändern zu Programmabbruch
+    Matrix error = {0, 0, NULL};
+    return error;
+  }
 }
 Matrix multiply(const Matrix matrix1, const Matrix matrix2) { return matrix1; }

matrix.h

@@ -19,6 +19,11 @@ void setMatrixAt(MatrixType value, Matrix matrix, unsigned int rowIdx,
                  unsigned int colIdx);
 MatrixType getMatrixAt(const Matrix matrix, unsigned int rowIdx,
                        unsigned int colIdx);
+
+Matrix broadCastCols(const Matrix matrix, const unsigned int rows,
+                     const unsigned int cols);
+Matrix broadCastRows(const Matrix matrix, const unsigned int rows,
+                     const unsigned int cols);
 Matrix add(const Matrix matrix1, const Matrix matrix2);
 Matrix multiply(const Matrix matrix1, const Matrix matrix2);