>>11066263Using arctan(x) = (i/2)log((1-ix)/(1+ix)) and a bit of algebra makes it work. No cleverness and a bunch of real/imaginary parts.
(-i)sqrt(i)arctan(sqrt(tan(x))/sqrt(i)) + (i)sqrt(-i)arctan(sqrt(tan(x))/sqrt(-i)) + C
= Re[2(-i)sqrt(i)arctan(sqrt(tan(x))/sqrt(i))] + C (LOOKS REAL)
= Re[sqrt(i)log((1-sqrt(i)sqrt(tan(x)))/(1+sqrt(i)sqrt(tan(x))))] + C
= (sqrt(2)/2)Re[log((1-sqrt(i)sqrt(tan(x)))/(1+sqrt(i)sqrt(tan(x))))] - (sqrt(2)/2)Im[log((1-sqrt(i)sqrt(tan(x)))/(1+sqrt(i)sqrt(tan(x))))] + C
= (sqrt(2)/2)log(|(1-sqrt(i)sqrt(tan(x)))/(1+sqrt(i)sqrt(tan(x)))|) - (sqrt(2)/2)Arg((1-sqrt(i)sqrt(tan(x)))/(1+sqrt(i)sqrt(tan(x)))) + C
log(|(1-sqrt(i)sqrt(tan(x)))/(1+sqrt(i)sqrt(tan(x)))|)
= log(|(1-(sqrt(2)/2)sqrt(tan(x))-i(sqrt(2)/2)sqrt(tan(x)))/(1+(sqrt(2)/2)sqrt(tan(x))+i(sqrt(2)/2)sqrt(tan(x)))|)
= (1/2)log(|(1-sqrt(2tan(x))+ tan(x))/(1+sqrt(2tan(x))+tan(x))|)
Arg((1-sqrt(i)sqrt(tan(x)))/(1+sqrt(i)sqrt(tan(x))))
= Arg((1-sqrt(i)sqrt(tan(x)))(1+ sqrt(-i)sqrt(tan(x)))/(1+sqrt(2)sqrt(tan(x))+tan(x)))
= Arg((1-sqrt(i)sqrt(tan(x)))(1+ sqrt(-i)sqrt(tan(x))))
= Arg(1-tan(x)-(i)sqrt(2)sqrt(tan(x)))
= arctan(sqrt(2tan(x))/(tan(x)-1))
The final result is
-(sqrt(2)/2)arctan(sqrt(2tan(x))/(tan(x)-1)) + (sqrt(2)/4)log(|(1-sqrt(2tan(x))+ tan(x))/(1+sqrt(2tan(x))+tan(x))|) + C
which agrees with OP.
It looks nasty but it's not too bad on paper.