Given a Gaussian integer $z = a + bi$, where $a, b \in \mathbb{Z}$, $i = \sqrt{-1}$, iterate the function $$f(z) = \frac{z}{1 + i}$$ if $z$ has even Gaussian norm (that is, both $a$ and $b$ are odd, or they're both even), otherwise $f(z) = 3z + i$.
I conjecture that iterating this function eventually leads, if not to $1$, to one of the other Gaussian units ($-1, i, -i$).
For example, starting with $z = 14$, we get $$7 - 7i, -7i, -22i, -11 - 11i, -11, -33 + i, -16 + 17i, \ldots$$
(this is wrong, see edit below)
I have tried a few different values of $z$ with small norm, some purely real, with pencil and paper, haven't gotten far. Also I have tried it in Mathematica, but either I've made some mistakes in my programming that crash the program, or there really are a lot of values that escape to some infinity.
Surely someone else has studied this variant? If so, have they been able to determine anything (like finding a periodic orbit that doesn't include any units)?
EDIT: I made a mistake. Mathematica did save my notebook file at some point before it crashed, and I could have gotten the correct sequence from there instead of having to recalculate it anew. It should go like this: $$7 - 7i, -7i, -20i, 10 + 10i, 10, 5 - 5i, -5i, -14i, -7 - 7i, -7, -20, \ldots$$
Thanks to Mr. Cortek for pointing this out.
The hope here is apparently to find that all Gaussian integers eventually lead to a number of the form $(1 + i)^n u$, where $n$ is a positive real integer and $u$ is a Gaussian unit.
Looking at a list of powers of $1 + i$, it seems to me like they take one of these forms: $\pm 2^n$, $\pm 2^n i$, $\pm 2^n \pm 2^n i$.
Then it's definitely possible for, say, $-3i$ to be followed by $-8i$ and form there it's a spiral drop to 1.
But when the real and imaginary parts are both nonzero, we run into the problem that $3z + i$ affects the real and imaginary parts unevenly.
For example, $5 + 5i$ multiplied by 3 is $15 + 15i$, but then you're only adding $i$ and so $15 + 16i$ falls short of $16 + 16i$ by 1.
Since $\mathbb Z[i]$ is a unique factorization domain and 3 is prime, a number with real and imaginary parts both nonzero is divisible by 3 only if each part is itself also divisible by 3.
Hence $3z$ fails to lead us to a number of the form $2^n + (2^n - 1)i$. A similar problem occurs with $3z + 1$ instead of $3z + i$.