void matrixmul(Mydata A[N][M],
	       MyData B[M][P],
	       MyData C[N][P])
{
 row: for (int i=0 ; i<N ; i++)
    {
    col: for (int j=0 ; j<P ; j++)
	{
	  MyData ABij = (MyData)0.0;

	product: for (int k=0 ; k<M ; k++)
	    {
	      ABij += (A[i][k] * B[k][j]);
	    }

	  C[i][j] = ABij;
	} /* cols */
    } /* rows */
  
  return;
}
