matrixMult.c

double** multiplyMatrices(double** A, int Am, int An, double** B, int Bm, int Bn) {
	int i, j, k, l;

	double** retAr = create2DDoubleArray(Am, Bn);

	if(An != Bm) {
		printf("OH FUCK YO\n");
		return NULL;
	}

	for(i = 0; i < Am; i++) {
		for(j = 0; j < Bn; j++) {
			for(k = 0; k < An; k++) {
				//printf("%d, %d, %d\n", i,j,k);
				retAr[i][j] += A[i][k] * B[k][j];
			}
		}
	}

	return retAr;
}