Fix a bug in the offset computation code in decompose_domain.

Fix a bug in the offset computation code in decompose_domain.

Left in my debug code (turned off with a #ifdef) to make it easier to exclude the potential computation as the reason zeros showed up in the output (which they no longer do).