Add support for (Ax) * y in elementwise mult, by internally forming a matrix mapping the global variable to y.