function A = MatrixProdFast(E,F,G,H)
A = (E.*F)'*G'*(H.*F)*G;