dgemkmm function

Evaluate kron(A,B) * C without storing kron(A,B)