2. loop: (x[i,...,n] + y[i,...,n]) * z[i,...,n] - w[i,...,n]
although I don't think the loops are going to be the main overhead in this case.