/******************************************************************************
 *
 * Programmer:  David Harbaugh
 * Program:     Serial Version of Matrix Multiply
 * File:        serial.c
 * Course:      CS 575: Parallel Processing (Distance Learning Section)
 * Date:        3 February 2005
 * Description:
 *
 *      Test program that mulitplies two matrices together.
 * 
 *
 ******************************************************************************/

#include <stdio.h>
#include <stdlib.h>

#define  SIZE  16

void print_matrix(int A[][SIZE], int s, char l);
void mult_matrix(int A[][SIZE], int B[][SIZE], int C[][SIZE], int s);
void add_matrix(int A[][SIZE], int B[][SIZE], int C[][SIZE], int s);

int main(void) {

  int A[SIZE][SIZE];                   /*  First matrix                      */
  int B[SIZE][SIZE];                   /*  Second matrix                     */
  int C[SIZE][SIZE];                   /*  Result of multiplication          */
  int i, j;

  /* Intialize A & B matrices */
  for (i=0;i<SIZE;i++)
   for (j=0;j<SIZE;j++) {
     A[i][j] = i+j;
     B[i][j] = i-j;
   }


  print_matrix(A, SIZE, 'A');
  printf("\n\n");
  print_matrix(B, SIZE, 'B');
  printf("\n\n");

  mult_matrix(A, B, C, SIZE);
  print_matrix(C, SIZE, 'C');
  printf("\n\n");

  add_matrix(A, B, C, SIZE);
  print_matrix(C, SIZE, 'C');
  printf("\n\n");

  exit(0);
}



/******************************************************************************
 *
 *  Function:  print_matrix
 *  Purpose:   Display to the screen a s x s integer matrix 
 *
 ******************************************************************************/

void print_matrix(int matrix[][SIZE], int s, char l) {
  int i,j;

  printf("%c: ", l);

  for (i=0;i<s;i++) {
    for (j=0;j<s;j++) {
      printf("%6d", matrix[i][j]);
    }
    printf("\n   ");
  }
}


/******************************************************************************
 *
 *  Function:  mult_matrix
 *  Purpose:   Multiply two SIZE x SIZE integer matrices, A and B, and store
 *             the result in C. 
 *
 ******************************************************************************/

void mult_matrix(int A[][SIZE], int B[][SIZE], int C[][SIZE], int s) {

  int i, j, k;
  int sum;

  for (i=0;i<s;i++) {
    for (j=0;j<s;j++) {
      sum = 0;

      for (k=0;k<s;k++) {
        sum += A[i][k] * B[k][j];
      }

      C[i][j] = sum;
    }
  }
}


/******************************************************************************
 *
 *  Function:  add_matrix
 *  Purpose:   Add two SIZE x SIZE integer matrices, A and B, and store the
 *             result in C. 
 *
 ******************************************************************************/

void add_matrix(int A[][SIZE], int B[][SIZE], int C[][SIZE], int s) {

  int i, j, k;

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