Compare commits

..

No commits in common. "dd2a621f11644921ef8b1cacafb904483c8f827e" and "648cbb1777bf5877cdf5793778454d7bb5cd8e06" have entirely different histories.

157
matrix.c
View File

@ -1,6 +1,5 @@
#include <stdlib.h> #include <stdlib.h>
#include <string.h> #include <string.h>
#include <stdio.h>
#include "matrix.h" #include "matrix.h"
// TODO Matrix-Funktionen implementieren // TODO Matrix-Funktionen implementieren
@ -8,13 +7,6 @@
Matrix createMatrix(unsigned int rows, unsigned int cols) Matrix createMatrix(unsigned int rows, unsigned int cols)
{ {
Matrix matrix; Matrix matrix;
Matrix empty = {0, 0, NULL};
if(rows == 0 || cols == 0)
{
//print("Fehler: Dimensionen muessen >= 1 sein!");
return empty;
}
matrix.rows = rows; matrix.rows = rows;
matrix.cols = cols; matrix.cols = cols;
@ -22,7 +14,7 @@ Matrix createMatrix(unsigned int rows, unsigned int cols)
if(matrix.buffer == NULL) if(matrix.buffer == NULL)
{ {
//printf("Fehler bei der Speicherreservierung! Keine Matrix erstellt!"); printf("Fehler bei der Speicherreservierung! Keine Matrix erstellt!");
matrix.rows = 0; matrix.rows = 0;
matrix.cols = 0; matrix.cols = 0;
} }
@ -32,7 +24,7 @@ Matrix createMatrix(unsigned int rows, unsigned int cols)
void clearMatrix(Matrix *matrix) void clearMatrix(Matrix *matrix)
{ {
free(matrix->buffer); free(matrix.buffer);
matrix->buffer = NULL; matrix->buffer = NULL;
matrix->rows = 0; matrix->rows = 0;
@ -41,32 +33,32 @@ void clearMatrix(Matrix *matrix)
void setMatrixAt(MatrixType value, Matrix *matrix, unsigned int rowIdx, unsigned int colIdx) void setMatrixAt(MatrixType value, Matrix *matrix, unsigned int rowIdx, unsigned int colIdx)
{ {
if(matrix.buffer == NULL) if(matrix->buffer == NULL)
{ {
//printf("Fehler beim Setzen! Matrix nicht initialisiert"); printf("Fehler beim Setzen! Matrix nicht initialisiert");
return; return;
} }
if(rowIdx >= matrix.rows || colIdx >= matrix.cols) if(rowIdx >= matrix->rows || colIdx >= matrix->cols)
{ {
//printf("Ungueltige Indizes beim Setzen!\n"); printf("Ungueltige Indizes beim Setzen!\n");
return; return;
} }
matrix.buffer[rowIdx * matrix.cols + colIdx] = value; matrix->buffer[rowIdx * matrix->cols + colIdx] = value;
} }
MatrixType getMatrixAt(const Matrix matrix, unsigned int rowIdx, unsigned int colIdx) MatrixType getMatrixAt(const Matrix matrix, unsigned int rowIdx, unsigned int colIdx)
{ {
if(matrix.buffer == NULL) if(matrix.buffer == NULL)
{ {
//printf("Fehler beim Lesen! Matrix nicht initialisiert"); printf("Fehler beim Lesen! Matrix nicht initialisiert");
return 0; return 0;
} }
if(rowIdx >= matrix.rows || colIdx >= matrix.cols) if(rowIdx >= matrix.rows || colIdx >= matrix.cols)
{ {
//printf("Ungueltige Indizes beim Lesen!\n"); printf("Ungueltige Indizes beim Lesen!\n");
return 0; return 0;
} }
@ -75,134 +67,31 @@ MatrixType getMatrixAt(const Matrix matrix, unsigned int rowIdx, unsigned int co
Matrix add(const Matrix matrix1, const Matrix matrix2) Matrix add(const Matrix matrix1, const Matrix matrix2)
{ {
if((matrix1.rows != matrix2.rows) || (matrix1.cols != matrix2.cols))
{
printf("Fehler bei Addition: Matrix Dimensionen passen nicht ueberein!\n");
Matrix empty = {0, 0, NULL};
return empty;
}
float matrix1Wert = 0; float matrix1Wert = 0;
float matrix2Wert = 0; float matrix2Wert = 0;
float summe = 0; float summe = 0;
int broadcasting = 0;
Matrix broadcastedmatrix1 = {0,0,NULL}; Matrix ergebnisMatrix = createMatrix(matrix1.rows, matrix1.cols);
Matrix broadcastedmatrix2 = {0,0,NULL};
if((matrix1.rows != matrix2.rows) || (matrix1.cols != matrix2.cols))
{
if((matrix1.rows != matrix2.rows) && (matrix1.cols == matrix2.cols))
{
if(matrix1.rows == 1)
{
broadcastedmatrix1 = createMatrix(matrix2.rows, matrix2.cols);
for(int row = 0; row<matrix2.rows; row++)
{
for(int col = 0; col<matrix2.cols; col++)
{
matrix1Wert = getMatrixAt(matrix1, 0, col);
setMatrixAt(matrix1Wert, broadcastedmatrix1, row, col);
}
}
broadcasting = 1;
}
else if(matrix2.rows == 1)
{
broadcastedmatrix2 = createMatrix(matrix1.rows, matrix1.cols);
for(int row = 0; row<matrix1.rows; row++)
{
for(int col = 0; col<matrix1.cols; col++)
{
matrix2Wert = getMatrixAt(matrix2, 0, col);
setMatrixAt(matrix2Wert, broadcastedmatrix2, row, col);
}
}
broadcasting = 2;
}
}
else if((matrix1.rows == matrix2.rows) && (matrix1.cols != matrix2.cols))
{
if(matrix1.cols == 1)
{
broadcastedmatrix1 = createMatrix(matrix2.rows, matrix2.cols);
for(int row = 0; row<matrix2.rows; row++)
{
for(int col = 0; col<matrix2.cols; col++)
{
matrix1Wert = getMatrixAt(matrix1, row, 0);
setMatrixAt(matrix1Wert, broadcastedmatrix1, row, col);
}
}
broadcasting = 1;
}
else if(matrix2.cols == 1)
{
broadcastedmatrix2 = createMatrix(matrix1.rows, matrix1.cols);
for(int row = 0; row<matrix1.rows; row++)
{
for(int col = 0; col<matrix1.cols; col++)
{
matrix2Wert = getMatrixAt(matrix2, row, 0);
setMatrixAt(matrix2Wert, broadcastedmatrix2, row, col);
}
}
broadcasting = 2;
}
}
else
{
//printf("Fehler bei Addition: Matrix Dimensionen stimmen nicht ueberein!\n");
Matrix empty = {0, 0, NULL};
return empty;
}
}
Matrix ergebnisMatrix;
switch (broadcasting)
{
case 0:
case 2:
ergebnisMatrix = createMatrix(matrix1.rows, matrix1.cols);
break;
case 1:
ergebnisMatrix = createMatrix(matrix2.rows, matrix2.cols);
break;
}
for(int row = 0; row < ergebnisMatrix.rows; row++) for(int row = 0; row < ergebnisMatrix.rows; row++)
{ {
for(int col = 0; col < ergebnisMatrix.cols; col++) for(int col = 0; col < ergebnisMatrix.cols; col++)
{ {
if(broadcasting == 0) matrix1Wert = getMatrixAt(matrix1, row, col);
{ matrix2Wert = getMatrixAt(matrix2, row, col);
matrix1Wert = getMatrixAt(matrix1, row, col); summe = matrix1Wert + matrix2Wert;
matrix2Wert = getMatrixAt(matrix2, row, col);
summe = matrix1Wert + matrix2Wert;
}
else if(broadcasting == 1)
{
matrix1Wert = getMatrixAt(broadcastedmatrix1, row, col);
matrix2Wert = getMatrixAt(matrix2, row, col);
summe = matrix1Wert + matrix2Wert;
}
else if(broadcasting == 2)
{
matrix1Wert = getMatrixAt(matrix1, row, col);
matrix2Wert = getMatrixAt(broadcastedmatrix2, row, col);
summe = matrix1Wert + matrix2Wert;
}
setMatrixAt(summe, ergebnisMatrix, row, col); setMatrixAt(summe, &ergebnisMatrix, row, col);
} }
} }
clearMatrix(&broadcastedmatrix1);
clearMatrix(&broadcastedmatrix2);
return ergebnisMatrix; return ergebnisMatrix;
} }
@ -210,7 +99,7 @@ Matrix multiply(const Matrix matrix1, const Matrix matrix2)
{ {
if(matrix1.cols != matrix2.rows) if(matrix1.cols != matrix2.rows)
{ {
//printf("Fehler bei Multiplikation: Matrix Dimensionen passen nicht ueberein!\n"); printf("Fehler bei Multiplikation: Matrix Dimensionen passen nicht ueberein!\n");
Matrix empty = {0, 0, NULL}; Matrix empty = {0, 0, NULL};
return empty; return empty;
} }
@ -230,7 +119,7 @@ Matrix multiply(const Matrix matrix1, const Matrix matrix2)
erg += getMatrixAt(matrix1, row, k) * getMatrixAt(matrix2, k, col); erg += getMatrixAt(matrix1, row, k) * getMatrixAt(matrix2, k, col);
} }
setMatrixAt(erg, ergebnisMatrix, row, col); setMatrixAt(erg, &ergebnisMatrix, row, col);
} }
} }