[julia-users] Efficient way to compute X' diag(w) Y

[Matthieu]

A lot of statistical algorithms require to compute X' diag(w) Y where X and 
Y are two matrices and w is a vector of weight (Y may or may not equal X). 
Is there a way to compute this product efficiently in Julia?