DiagProd.m
diagonal of the product without creating all the product matrix. Very fast and memory saving. Use:instead of diag(A*B) use diagprod(A,B)