#include "include/define.h"
#include <iostream.h>

void matprint(pel_t *mat[],int size){
  int m,n;
  
  for(m=0;m<size;m++){
    for(n=0;n<size;n++)
      cout<<mat[m][n]<<" ";
    cout<<endl;
  }
  cout<<endl;

}

pel_t transpose(pel_t *A[],int m,int n)
{/* returns mth row, nth column element of transpose(A) */

  return A[n][m];

}

void mat_add(pel_t *sum[],pel_t *A[],pel_t *B[],int width,int height){
  //This function adds two matrices together (sum = A + B)

  int m,n; //m is height, n is width
  
  for(m = 0; m < height; m++)
    for(n = 0; n < width; n++) 	
      sum[m][m]=A[m][n] + B[m][n];

}

void mat_sub(pel_t *diff[],pel_t *A[],pel_t *B[],int width, int height){
  //This function subtracts two matrices together (S = A - B)
  
  int m,n; //m is height, n is width
  
  for(m = 0; m < height; m++)
    for(n = 0; n < width; n++) 	
      diff[m][m]=A[m][n] - B[m][n];

}