/*file MatrixMult.c
	C program that calls an assembly routine. It asks the
	user for the size and initial values for the matrices
	and then shows the results of their multiplication.
	EMR 990407
      Modified:
	MRR 000327
*/
#include "support.h"

int *mata; /* First Matrix */
int *matb; /* Second Matrix */
int *matc; /* Resultant Matrix */
int size; /* Matrices Size */

void main()
{
      set_up_transfer_base();
	int i, l;
	int temp;

	/* Ask the user for the size */
	printf("Enter size of matrices > ");
	scanf("%d",&size);
	printf("\n");

	/* Separate memory for the matrices */
	mata = malloc(size*size);
	matb = malloc(size*size);
	matc = malloc(size*size);
	
	/* Ask user for initial values */
	for(int j=0; j<size; j++)
	{
		for(int m=0; m<size; m++)
		{
			printf("A%d",j+1);
			printf("%d = ",m+1);
			scanf("%d", &temp);
			printf("\n");
			mata[j*size+m]=temp;
			
			printf("B%d",j+1);
			printf("%d = ",m+1);
			scanf("%d", &temp);
			printf("\n");
			matb[j*size+m]=temp;
			
			/* Clear resultant matrix memory */
			matc[j*size+m]=0;
		}
	}
	
	/* Calling assembly routine */
	matmult(size, matc, mata, matb);
	
	/* Display results */
	for(i=0; i<size; i++)
	{
		printf("| ");
		for(l=0; l<size; l++)
			printf("%d ",mata[i*size+l]);
		printf("| ");
		if((i+1)==(size/2))
		{
			printf("*");
		}
		else
		{
			printf(" ");
		}
		printf(" | ");
		for(l=0; l<size; l++)	
			printf("%d ",matb[i*size+l]);
		printf("| ");
		if((i+1)==(size/2))
		{
			printf("=");
		}
		else
		{
			printf(" ");
		}
		printf(" | ");
		for(l=0; l<size; l++)
			printf("%d ",matc[i*size+l]);
		printf("|");
		printf("\n");
	}
}