dd-complex library
混合精度演算やベクトル演算など全ての種類の演算はまだ実装していませんが,普通の四則演算であれば大体対応しています.
CやC++からインクルードして使用してください.
構造体で書いているので,dd_compの型を宣言して.reと.imで実部と虚部,
さらにそれぞれに対して.hiと.loで上位ビット下位ビットにアクセスできます.
.hiと.loの型は普通のdoubleです.
aa.re.hi = hogehoge;
b[i].im.lo = 0.0;