code: def matrix_product(mat1, mat2):
      
      return np.dot(mat2.T, mat1.T).T
docs: Compute the product of two Fortran contiguous matrices.
  
      This is to avoid the overhead of NumPy converting to C-contiguous
      before computing a matrix product.
  
      Does so via ``A B = (B^T A^T)^T`` since ``B^T`` and ``A^T`` will be
      C-contiguous without a copy, then the product ``P = B^T A^T`` will
      be C-contiguous and we can return the view ``P^T`` without a copy.
  
      Args:
          mat1 (numpy.ndarray): The left-hand side matrix.
          mat2 (numpy.ndarray): The right-hand side matrix.
  
      Returns:
          numpy.ndarray: The product of the two matrices.