Let us define $K(x)$ as complete elliptic integral of the first kind, where $x$ is elliptic modulus. A possible closed-form is ($G$ denotes Catalan's constant.) $$ \int_{0}^{\pi/2} \ln\left[ \tan\left ( \frac{\theta}{2}\right) \right ]^2 K\left ( \sin\theta \right )\text{d}\theta =\frac{\Gamma\left ( \frac14 \right )^4G }{8\pi}. $$ It looks like a "product" of two solvable integrals (both are elementary): $$\int_{0}^{\pi/2} \ln\left[ \tan\left ( \frac{\theta}{2}\right) \right ]^2\text{d}\theta =\frac{\pi^3 }{8}$$ and $$ \int_{0}^{\pi/2} K\left ( \sin\theta \right )\text{d}\theta =\frac{\Gamma\left ( \frac14 \right )^4 }{16\pi}. $$
Question: How can we evaluate the integral? I try to utilize the Fourier series $K(\sin\theta) =\pi\sum_{n\ge0} \frac{\left ( \frac12 \right )_n^2 }{(n!)^2} \sin\left ( \left ( 4n+1 \right )\theta \right )$ to prove, but seems not to go well. I appreciate for your help.
An Interesting Observation: We find $$ \int_{0}^{\pi/2} \ln\left[ \tan\left ( \frac{\theta}{2}\right) \right ]^4 K\left ( \sin\theta \right )\text{d}\theta =\frac{3\,\Gamma\left ( \frac14 \right )^4}{4\pi}(G^2+\beta(4)) $$ where $\beta(.)$ is Dirichlet's $\beta$ function.