void
square_dgemm(const unsigned M, 
             const double *A, const double *B, double *C)
{
  extern int sdgemm_(const unsigned* M, const double* A, const double* B, 
                    double* C);

  sdgemm_(&M, A, B, C);
}