Allow different delta values (dx, dy, dz) for "Adaptive Weighted" algorithms?

Currently, all version of "AW" algorithms have the same "delta" parameter controls the adaptive weight exp( -(xi-xj)^2/ delta) for all (x, y, z) directions. I can see the need to specify different "deltas" when the x, y, z axis have different resolutions (for example, some tomosynthesis situations). Could we make changes to allow this? We could allow delta to be either a scalar (likely current situation) or a array of (d1, d2, d3). Need some code changes in POCS_TV.cu and other places as well.