Dear Julia users,
Is there a way to improve performance of this kind of code without devectorizing (use loops). All variables are 2d arrays. I solve a PDE problem with finite differences method. Thanks for your help
996988748 divV .= diff(Vx,dims=1)/dx
992096000 divV .+= diff(Vy,dims=2)/dy

992096000 Exxc .= diff(Vx,dims=1)/dx
495600000 Exxc .= 1/2*divV
992096000 Eyyc .= diff(Vy,dims=2)/dy
495600000 Eyyc .= 1/2*divV

507248000 Exyv .= diff(Vx_exp,dims=2)
0 Exyv ./= dy
507248000 Exyv .+= diff(Vy_exp,dims=1)
0 Exyv ./= dx
0 Exyv .*= 0.5

1984192000 @views Exyc .= 0.25 * ( Exyv[1:end1,1:end1]
 + Exyv[2:end ,1:end1]
 + Exyv[1:end1,2:end ]
 + Exyv[2:end ,2:end ])