Re: Complex square root of -1 : zsqrt(-1)

But I wouldn't bother really, use some assembly à la
double inline __declspec (naked) __fastcall sqrt(double n)
{
    _asm fld qword ptr [esp+4]
    _asm fsqrt
    _asm ret 8
}