import llvmdc.intrinsics;

extern(C) int printf(char*,...);
extern(C) int scanf(char*,...);

void main()
{
    float f;
    printf("Enter float: ");
    //scanf("%f", &f);
    f = 1.22345;
    float sf = llvm_sqrt_f32(f);
    printf("sqrt(%f) = %f\n", f, sf);

    double d;
    printf("Enter double: ");
    //scanf("%lf", &d);
    d = 2.2311167895435245;
    double sd = llvm_sqrt_f64(d);
    printf("sqrt(%lf) = %lf\n", d, sd);

    real r;
    printf("Enter real: ");
    //scanf("%lf", &d);
    r = 3.2311167891231231234754764576;
    version(X86)
    {
        real sr = llvm_sqrt_f80(r);
        printf("sqrt(%llf) = %llf\n", r, sr);
    }
    else
    {
        real sr = llvm_sqrt_f64(r);
        printf("sqrt(%f) = %lf\n", r, sr);
    }
}