#include <stdio.h>
/* Varataan matriiseille tilaa N*N alkion verran. 
 * Laskuissa käytettävä koko kysytään jälleen käyttäjältä. */
#define N 100

/* Matriisien summa ja tulo eivät välttämättä ole määriteltyjä
 * (jos matriisin väärän kokoiset), joten laskutoimituksen tekeville
 * funktioille määritellään palautusarvo: 1 jos lasku onnistui, 0 jos ei. 
 * Vakio kertaa matriisi ja transpoosi ovat aina määriteltyjä. */
int matr_sum(double A[N][N], double B[N][N], double C[N][N], int nA, int mA, int nB, int mB);
void scal_mult(double A[N][N], double a, double C[N][N], int nA, int mA);
int matr_prod(double A[N][N], double B[N][N], double C[N][N], int nA, int mA, int nB, int mB);
void transpose(double A[N][N], double C[N][N], int nA, int mA);

void prn_matr(double A[N][N], int nA, int mA);

int main()
{
	int mA,nA,mB,nB;
	int i,j;

	double A[N][N], B[N][N];

	printf("Enter dimensions of A (separated by a space): "); scanf("%i %i",&mA,&nA);

	printf("Enter components of A:\n");
	for(i=0;i<mA;i++)
	{
		for(j=0;j<nA;j++)
		{
			printf("A[%2i][%2i] = ",i,j); scanf("%lf",&A[i][j]);
		}
	}

	printf("Enter dimensions of B (separated by a space): "); scanf("%i %i",&mB,&nB);

	printf("Enter components of B:\n");
	for(i=0;i<mB;i++)
	{
		for(j=0;j<nB;j++)
		{
			printf("B[%2i][%2i] = ",i,j
				); scanf("%lf",&B[i][j]);
		}
	}

	prn_matr(A,nA,mA);

	prn_matr(A,nB,mB);
}

int matr_sum(double A[N][N], double B[N][N], double C[N][N], int nA, int mA, int nB, int mB)
{
	int i,j;

	if(nA != nB || mA != mB)
	{
		printf("Illegal dimensions! Matrix sum not defined!\n");
		return 0;
	}

	for(i=0;i<mA;i++)
	{
		for(j=0;j<nA;j++)
		{
			C[i][j]=A[i][j]+B[i][j];
		}
	}

	return 1;
}

void scal_mult(double A[N][N], double a, double C[N][N], int nA, int mA)
{
	int i,j;
	for(i=0;i<mA;i++)
	{
		for(j=0;j<nA;j++)
		{
			C[i][j]=a*A[i][j];
		}
	}
}

int matr_prod(double A[N][N], double B[N][N], double C[N][N], int nA, int mA, int nB, int mB)
{
	int i,j,k;

	if(nA != mB)
	{
		printf("Illegal dimensions! Matrix product not defined!\n");
		return 0;
	}

	for(i=0;i<mA;i++)
	{
		for(j=0;j<nB;j++)
		{
			for(k=0,C[i][j]=0.0;k<nA;k++)
			{
				C[i][j] += A[i][k]*B[k][j];
			}
		}
	}
	return 1;
}

void transpose(double A[N][N], double C[N][N], int nA, int mA)
{
}

void prn_matr(double A[N][N], int nA, int mA)
{
	int i,j;

	for(i=0;i<mA;i++)
	{
		printf("[");
		for(j=0;j<nA;j++)
		{
			printf("%6.2lf",A[i][j]);
			if(j<nA-1)
				printf(",");
		}
		printf("]\n");
	}
}