#include "SparseMat.h" 
float * MultiplySparseMats(struct mat_sp A, struct mat_sp B)
{
  float *	C;
  struct mat_sp Bt;
  int		i, j;
	
  /* Initialize product matrix */
  C = malloc(A.m*B.n*sizeof(float));
  
  /* Create Transpose */
  Bt = TransposeSparseMat(B);
  
  /* Fill Matrix */
  for(i=0; i<A.m; i++)
  { 
    for(j=0; j<B.n; j++)
    {
    C[i*A.m + j] = MultiplySparseRows(A.row[i],Bt.row[i]);
    }
  }
  
  FreeSparseMat(Bt);
  
  /* Returning unformated pointer */
  return C;
}