Consider the Fuchsian group $\Gamma:=\Big\langle\begin{pmatrix}1&1\\1&2\end{pmatrix}, \begin{pmatrix}1&-1\\-1&2\end{pmatrix}\Big\rangle$. A commonly studied fact is that when this acts upon the upper half-plane $\mathbb{H}^2$,
the quotient $T:=\mathbb{H}^2/\Gamma$ is a hyperbolic once-punctured torus. For instance taking the Ford domain of $\Gamma$, we get as a fundamental region the hyperbolic quadrilateral in the center of the following tessellation.
I am interested in getting similar pictures using Dirichlet domains instead. For $p\in\mathbb{H}^3$ let $D(p)$ denote the Dirichlet domain for $\Gamma$ centered at $p$. Note that since $\Gamma$ is torsion-free, there are no restrictions on which $p$ can be used.
What is the effect on $D(p)$ as we vary $p$? Is $D(p)$ always a four-sided polygon with ideal vertices? Also, for a given $p$ what elements of $\Gamma$ contribute sides to $D(p)$, and which ones occur in gluing the sides?
UPDATE: Even an investigation of the simplest scenario here lead me to complications, see https://mathoverflow.net/q/219139/14835 for more on that.
Here is a partial answer, which may point the way to a reasonable conjecture. Along the way I do answer one of your several questions.
First a correction: the Ford domain is independent of $p$ (up to covering transformation). It depends on the Fuchsian group alone. The Dirichlet domain $D(p)$ depends, of course, on $p$.
Next some setup, designed to learn as much about $D(p)$ as possible without doing any actual work.
The boundary $\partial D(p)$, when projected to the quotient punctured torus $T$, becomes an embedded 1-complex $S(p)\subset T$. $S(p)$ is a closed subset of $T$, and it is not compact, indeed $S(p)$ accumulates on the puncture of $T$.
Let $\infty$ denote the "puncture" of $T$, aka the end of $T$. Denote $\overline T = T \cup \infty$ to be the one-point compactification, aka the end compactification in this case, so $\overline T$ is an ordinary, unpunctured, compact torus. Denote the closure of $S(p)$ in $\overline T$ as $\overline S(p) = S(p) \cup \infty$. Then $\overline S(p)$ is itself a compact 1-complex. The sides of $\partial D(p)$ glue up in pairs to form the 1-cells of $\overline S(p)$, and so the number of $1$-cells of $\overline S(p)$ is always half the number of sides of $D(p)$.
The topology of $\overline S(p)$ has the following restrictions. First, $\overline S(p)$ is a spine of the once-punctured torus $\overline T - \{p\}$, and so $\overline S(p)$ is a graph of rank $2$ and Euler characteristic $-1$. Its vertex $\infty$ has valence $\ge 1$. All finite vertices have valence $\ge 3$. So by an Euler characteristic argument on the finite graph $\overline S(p)$ the possibilities for $\overline S(p)$ can be classified into the following seven types, in order of "decreasing genericity":
With that, you have asked:
which is equivalent to asking
For example, in your picture if you take $p=i$ (not $p=2i$) then $D(i)$ is the ideal quadrilateral with endpoints $-1,0,1,\infty$ as your picture shows, and this is an an "ideal square" meaning that its isometric symmetry group is the order 8 dihedral group. The "center" of this ideal square is $p=i$. So, with the choice $p=i$, type 7 indeed occurs.
However, the answer to your question is "no". For example, if $p$ is near the puncture, e.g. if $p=\lambda i$ for very large $\lambda$, then it's not hard to see that the valence of $\infty$ equals $1$ and so either type 1, 2, or 3 occurs.
Perhaps a more natural conjecture is that on an open dense set of values of $p$, what occurs is type 1, which is the most generic type. An alternative conjecture is that, since the shape of $T$ itself is rather nongeneric in the Teichmuller space, on an open dense set of values of $p$ what occurs is type 2 or 3.