int main(int argc, char *argv[]) {
  float a, b, c, sd, alpha, beta, z;
  double a2, b2, c2, sd2, alpha2, beta2, z2;
  a = a2 = 1.00000;
  b = b2 = -333.333;
  c = c2 = 1.23456;
  printf("a = %-12f : ", a);
  print_bit(a);
  printf("b = %-12f : ", b);
  print_bit(b);
  printf("c = %-12f : ", c);
  print_bit(c);
  sd = sqrt(b * b - 4 * a * c);
  sd2 = sqrt(b2 * b2 - 4 * a2 * c2);
  alpha = (-b + sd) / (2 * a);
  beta = (-b - sd) / (2 * a);
  alpha2 = (-b2 + sd2) / (2 * a2);
  beta2 = (-b2 - sd2) / (2 * a2);
  printf("alpha = %15e    beta = %15e\n", alpha, beta);
  printf("alpha2= %15e    beta2= %15e\n", alpha2, beta2);
  printf("( sd= %f )\n", sd);
  printf("( sd= %lf )\n", sd2);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
movl   $0x0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  0xed2(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movss  0xedd(%rip),%xmm0
movss  %xmm0,-0x14(%rbp)
movsd  0xeb0(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movss  0xebf(%rip),%xmm0
movss  %xmm0,-0x18(%rbp)
movsd  0xe8e(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movss  0xea1(%rip),%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xe98(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1040 <print_bit@plt>
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xe7e(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1040 <print_bit@plt>
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xe64(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1040 <print_bit@plt>
movss  -0x18(%rbp),%xmm0
movss  -0x18(%rbp),%xmm2
movss  0xe11(%rip),%xmm1
mulss  -0x14(%rbp),%xmm1
mulss  -0x1c(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  0xdaa(%rip),%xmm1
mulsd  -0x38(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  %xmm0,-0x50(%rbp)
movss  -0x18(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
addss  -0x20(%rbp),%xmm0
movss  0xd78(%rip),%xmm1
mulss  -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x18(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
subss  -0x20(%rbp),%xmm0
movss  0xd4b(%rip),%xmm1
mulss  -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  0xcec(%rip),%xmm1
mulsd  -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  0xcb5(%rip),%xmm1
mulsd  -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x28(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
lea    0xcf1(%rip),%rdi
mov    $0x2,%al
call   1030 <printf@plt>
movsd  -0x58(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm1
lea    0xcf6(%rip),%rdi
mov    $0x2,%al
call   1030 <printf@plt>
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xcfc(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movsd  -0x50(%rbp),%xmm0
lea    0xcf5(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
xor    %eax,%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rax
lea    0xf20(%rip),%rdi
movsd  0xed0(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
movsd  0xec1(%rip),%xmm0
mov    $0x1,%al
call   1040 <print_bit@plt>
lea    0xf08(%rip),%rdi
movsd  0xeb3(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
movsd  0xea4(%rip),%xmm0
mov    $0x1,%al
call   1040 <print_bit@plt>
lea    0xef0(%rip),%rdi
movsd  0xe96(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
movsd  0xe87(%rip),%xmm0
mov    $0x1,%al
call   1040 <print_bit@plt>
lea    0xed8(%rip),%rdi
movsd  0xe79(%rip),%xmm0
movsd  0xe79(%rip),%xmm1
mov    $0x2,%al
call   1030 <printf@plt>
lea    0xed7(%rip),%rdi
movsd  0xe6b(%rip),%xmm0
movsd  0xe6b(%rip),%xmm1
mov    $0x2,%al
call   1030 <printf@plt>
lea    0xed6(%rip),%rdi
movsd  0xe5d(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
lea    0xecc(%rip),%rdi
movsd  0xe4f(%rip),%xmm0
mov    $0x1,%al
call   1030 <printf@plt>
xor    %eax,%eax
pop    %rcx
ret

LLVM IR (-O0):
@.str = [13 x i8] c"a = %-12f : \00"
@.str.1 = [13 x i8] c"b = %-12f : \00"
@.str.2 = [13 x i8] c"c = %-12f : \00"
@.str.3 = [29 x i8] c"alpha = %15e beta = %15e\0A\00"
@.str.4 = [29 x i8] c"alpha2= %15e beta2= %15e\0A\00"
@.str.5 = [12 x i8] c"( sd= %f )\0A\00"
@.str.6 = [13 x i8] c"( sd= %lf )\0A\00"
define i32 @main(i32 %0, i8** %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i8**
%6 = alloca float
%7 = alloca float
%8 = alloca float
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
store i32 0, i32* %3
store i32 %0, i32* %4
store i8** %1, i8*** %5
store double 1.000000e+00, double* %13
store float 1.000000e+00, float* %6
store double -3.333330e+02, double* %14
store float 0xC074D55400000000, float* %7
store double 1.234560e+00, double* %15
store float 0x3FF3C0C200000000, float* %8
%20 = load float, float* %6
%21 = fpext float %20 to double
%22 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str, i64 0, i64 0), double %21)
%23 = load float, float* %6
%24 = fpext float %23 to double
%25 = call i32 (double, ...) bitcast (i32 (...)* @print_bit to i32 (double, ...)*)(double %24)
%26 = load float, float* %7
%27 = fpext float %26 to double
%28 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.1, i64 0, i64 0), double %27)
%29 = load float, float* %7
%30 = fpext float %29 to double
%31 = call i32 (double, ...) bitcast (i32 (...)* @print_bit to i32 (double, ...)*)(double %30)
%32 = load float, float* %8
%33 = fpext float %32 to double
%34 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.2, i64 0, i64 0), double %33)
%35 = load float, float* %8
%36 = fpext float %35 to double
%37 = call i32 (double, ...) bitcast (i32 (...)* @print_bit to i32 (double, ...)*)(double %36)
%38 = load float, float* %7
%39 = load float, float* %7
%40 = load float, float* %6
%41 = fmul float 4.000000e+00, %40
%42 = load float, float* %8
%43 = fmul float %41, %42
%44 = fneg float %43
%45 = call float @llvm.fmuladd.f32(float %38, float %39, float %44)
%46 = fpext float %45 to double
%47 = call double @sqrt(double %46)
%48 = fptrunc double %47 to float
store float %48, float* %9
%49 = load double, double* %14
%50 = load double, double* %14
%51 = load double, double* %13
%52 = fmul double 4.000000e+00, %51
%53 = load double, double* %15
%54 = fmul double %52, %53
%55 = fneg double %54
%56 = call double @llvm.fmuladd.f64(double %49, double %50, double %55)
%57 = call double @sqrt(double %56)
store double %57, double* %16
%58 = load float, float* %7
%59 = fneg float %58
%60 = load float, float* %9
%61 = fadd float %59, %60
%62 = load float, float* %6
%63 = fmul float 2.000000e+00, %62
%64 = fdiv float %61, %63
store float %64, float* %10
%65 = load float, float* %7
%66 = fneg float %65
%67 = load float, float* %9
%68 = fsub float %66, %67
%69 = load float, float* %6
%70 = fmul float 2.000000e+00, %69
%71 = fdiv float %68, %70
store float %71, float* %11
%72 = load double, double* %14
%73 = fneg double %72
%74 = load double, double* %16
%75 = fadd double %73, %74
%76 = load double, double* %13
%77 = fmul double 2.000000e+00, %76
%78 = fdiv double %75, %77
store double %78, double* %17
%79 = load double, double* %14
%80 = fneg double %79
%81 = load double, double* %16
%82 = fsub double %80, %81
%83 = load double, double* %13
%84 = fmul double 2.000000e+00, %83
%85 = fdiv double %82, %84
store double %85, double* %18
%86 = load float, float* %10
%87 = fpext float %86 to double
%88 = load float, float* %11
%89 = fpext float %88 to double
%90 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([29 x i8], [29 x i8]* @.str.3, i64 0, i64 0), double %87, double %89)
%91 = load double, double* %17
%92 = load double, double* %18
%93 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([29 x i8], [29 x i8]* @.str.4, i64 0, i64 0), double %91, double %92)
%94 = load float, float* %9
%95 = fpext float %94 to double
%96 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([12 x i8], [12 x i8]* @.str.5, i64 0, i64 0), double %95)
%97 = load double, double* %16
%98 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.6, i64 0, i64 0), double %97)
ret i32 0
}

LLVM IR (-O2):
@.str = [13 x i8] c"a = %-12f : \00"
@.str.1 = [13 x i8] c"b = %-12f : \00"
@.str.2 = [13 x i8] c"c = %-12f : \00"
@.str.3 = [29 x i8] c"alpha = %15e beta = %15e\0A\00"
@.str.4 = [29 x i8] c"alpha2= %15e beta2= %15e\0A\00"
@.str.5 = [12 x i8] c"( sd= %f )\0A\00"
@.str.6 = [13 x i8] c"( sd= %lf )\0A\00"
define i32 @main(i32 %0, i8** nocapture readnone %1) {
%3 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([13 x i8], [13 x i8]* @.str, i64 0, i64 0), double 1.000000e+00)
%4 = tail call i32 (double, ...) bitcast (i32 (...)* @print_bit to i32 (double, ...)*)(double 1.000000e+00)
%5 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([13 x i8], [13 x i8]* @.str.1, i64 0, i64 0), double 0xC074D55400000000)
%6 = tail call i32 (double, ...) bitcast (i32 (...)* @print_bit to i32 (double, ...)*)(double 0xC074D55400000000)
%7 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([13 x i8], [13 x i8]* @.str.2, i64 0, i64 0), double 0x3FF3C0C200000000)
%8 = tail call i32 (double, ...) bitcast (i32 (...)* @print_bit to i32 (double, ...)*)(double 0x3FF3C0C200000000)
%9 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([29 x i8], [29 x i8]* @.str.3, i64 0, i64 0), double 0x4074D544C0000000, double 0x3F6E600000000000)
%10 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([29 x i8], [29 x i8]* @.str.4, i64 0, i64 0), double 0x4074D544CC2BC6FB, double 0x3F6E574624380000)
%11 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([12 x i8], [12 x i8]* @.str.5, i64 0, i64 0), double 0x4074D535A0000000)
%12 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([13 x i8], [13 x i8]* @.str.6, i64 0, i64 0), double 0x4074D535A088B4DF)
ret i32 0
}