This answer uses vectors.
Let $\vec{CB}=\vec b,\vec{CD}=\vec d$.
Then, we can write $\vec{CP}=m\vec b+n\vec d$ where $m,n\in\mathbb R$ and have $$|\vec{CB}|^2=|\vec{CP}|^2\iff |\vec b|^2=m^2|\vec b|^2+n^2|\vec d|^2+2mn\vec b\cdot\vec d\tag1$$
Now we have
$\vec{FG}=\vec{CG}-\vec{CF}=\frac 12(\vec{CP}+\vec{CA})-\frac 12\vec d=\frac 12(m\vec b+n\vec d+\vec b+\vec d)-\frac 12\vec d=\frac 12(m\vec b+n\vec d+\vec b)$
and
$\vec{BP}=\vec{CP}-\vec{CB}=m\vec b+n\vec d-\vec b$
So,
$$\begin{align}\vec{FG}\cdot\vec{BP}&=\frac 12(m\vec b+n\vec d+\vec b)\cdot (m\vec b+n\vec d-\vec b)\\\\&=\frac 12(m|\vec b|^2+n^2|\vec d|^2+2mn\vec b\cdot\vec d-|\vec b|^2)\end{align}$$
We see that this equals $0$ from $(1)$.