int encode_GDS_R4(buf, dchars, size)
char *buf;
unsigned char *dchars;
int size;
{
  double val = 0.0;
  int i;
  int exponent = 0;
  char isneg;
  sscanf(buf, "%lf", &val);
  if (isneg = (val < 0.0)) {
    val = -val;
  }
  exponent += 64;
  if (val == 0.0) {
  } else if (val >= 1.0) {
    while (val >= 1.0) {
      exponent++;
      val /= 16.0;
    }
  } else if (val < (1.0 / 16.0)) {
    while (val < (1.0 / 16.0)) {
      exponent--;
      val *= 16.0;
    }
  }
  dchars[0] = exponent;
  if (isneg)
    dchars[0] |= 0x80;
  for (i = 1; i < 4; i++) {
    val *= 256.0;
    dchars[i] = (int)val;
    val -= dchars[i];
  }
  return 4;
}


x86 assembly (clang, -O0):
<encode_GDS_R4>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x8(%rbp),%rdi
lea    0xee3(%rip),%rsi
lea    -0x20(%rbp),%rdx
mov    $0x0,%al
call   1030 <sscanf@plt>
xorps  %xmm0,%xmm0
ucomisd -0x20(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %al,-0x29(%rbp)
cmp    $0x0,%al
je     1184 <encode_GDS_R4+0x74>
movsd  -0x20(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x28(%rbp),%eax
add    $0x40,%eax
mov    %eax,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11aa <encode_GDS_R4+0x9a>
jp     11aa <encode_GDS_R4+0x9a>
jmp    1260 <encode_GDS_R4+0x150>
movsd  -0x20(%rbp),%xmm0
movsd  0xe49(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1206 <encode_GDS_R4+0xf6>
jmp    11c6 <encode_GDS_R4+0xb6>
movsd  -0x20(%rbp),%xmm0
movsd  0xe2d(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1201 <encode_GDS_R4+0xf1>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  0xe1d(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11c6 <encode_GDS_R4+0xb6>
jmp    125b <encode_GDS_R4+0x14b>
movsd  0xdfa(%rip),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    1256 <encode_GDS_R4+0x146>
jmp    121e <encode_GDS_R4+0x10e>
movsd  0xde2(%rip),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    1251 <encode_GDS_R4+0x141>
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
movsd  0xdce(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    121e <encode_GDS_R4+0x10e>
jmp    1256 <encode_GDS_R4+0x146>
jmp    125b <encode_GDS_R4+0x14b>
jmp    1260 <encode_GDS_R4+0x150>
mov    -0x28(%rbp),%eax
mov    %al,%cl
mov    -0x10(%rbp),%rax
mov    %cl,(%rax)
cmpb   $0x0,-0x29(%rbp)
je     1284 <encode_GDS_R4+0x174>
mov    -0x10(%rbp),%rax
movzbl (%rax),%ecx
or     $0x80,%ecx
mov    %cl,(%rax)
movl   $0x1,-0x24(%rbp)
cmpl   $0x4,-0x24(%rbp)
jge    12e5 <encode_GDS_R4+0x1d5>
movsd  0xd7b(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
cvttsd2si -0x20(%rbp),%eax
mov    %al,%dl
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movzbl (%rax,%rcx,1),%eax
cvtsi2sd %eax,%xmm1
movsd  -0x20(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    128b <encode_GDS_R4+0x17b>
mov    $0x4,%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<encode_GDS_R4>:
push   %rbx
sub    $0x10,%rsp
mov    %rsi,%rbx
movq   $0x0,0x8(%rsp)
lea    0xf08(%rip),%rsi
lea    0x8(%rsp),%rdx
xor    %eax,%eax
call   1030 <sscanf@plt>
movsd  0x8(%rsp),%xmm1
movapd 0xebe(%rip),%xmm3
xorpd  %xmm1,%xmm3
xorpd  %xmm2,%xmm2
movapd %xmm1,%xmm0
cmpltsd %xmm2,%xmm0
andpd  %xmm0,%xmm3
andnpd %xmm1,%xmm0
orpd   %xmm3,%xmm0
mov    $0x40,%eax
ucomisd %xmm2,%xmm0
jne    116c <encode_GDS_R4+0x5c>
jnp    11dd <encode_GDS_R4+0xcd>
ucomisd 0xe9c(%rip),%xmm0
jae    11af <encode_GDS_R4+0x9f>
movsd  0xe9a(%rip),%xmm3
ucomisd %xmm0,%xmm3
jbe    11dd <encode_GDS_R4+0xcd>
mov    $0x40,%eax
movsd  0xe8f(%rip),%xmm4
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0xffffffff,%eax
mulsd  %xmm4,%xmm0
ucomisd %xmm0,%xmm3
ja     11a0 <encode_GDS_R4+0x90>
jmp    11dd <encode_GDS_R4+0xcd>
mov    $0x40,%eax
movsd  0xe5c(%rip),%xmm3
movsd  0xe4c(%rip),%xmm4
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%eax
mulsd  %xmm3,%xmm0
ucomisd %xmm4,%xmm0
jae    11d0 <encode_GDS_R4+0xc0>
mov    %eax,%ecx
or     $0x80,%cl
ucomisd %xmm1,%xmm2
movzbl %cl,%ecx
cmovbe %eax,%ecx
mov    %cl,(%rbx)
movsd  0xe32(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %al,0x1(%rbx)
movzbl %al,%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
subsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %al,0x2(%rbx)
movzbl %al,%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
subsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %al,0x3(%rbx)
mov    $0x4,%eax
add    $0x10,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@.str = [4 x i8] c"%lf\00"
define i32 @encode_GDS_R4(i8* %0, i8* %1, i32 %2) {
%4 = alloca i8*
%5 = alloca i8*
%6 = alloca i32
%7 = alloca double
%8 = alloca i32
%9 = alloca i32
%10 = alloca i8
store i8* %0, i8** %4
store i8* %1, i8** %5
store i32 %2, i32* %6
store double 0.000000e+00, double* %7
store i32 0, i32* %9
%11 = load i8*, i8** %4
%12 = call i32 (i8*, i8*, ...) @sscanf(i8* %11, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str, i64 0, i64 0), double* %7)
%13 = load double, double* %7
%14 = fcmp olt double %13, 0.000000e+00
%15 = zext i1 %14 to i32
%16 = trunc i32 %15 to i8
store i8 %16, i8* %10
%17 = icmp ne i8 %16, 0
br i1 %17, label %18, label %21
18: ; preds = %3
%19 = load double, double* %7
%20 = fneg double %19
store double %20, double* %7
br label %21
21: ; preds = %18, %3
%22 = load i32, i32* %9
%23 = add nsw i32 %22, 64
store i32 %23, i32* %9
%24 = load double, double* %7
%25 = fcmp oeq double %24, 0.000000e+00
br i1 %25, label %26, label %27
26: ; preds = %21
br label %55
27: ; preds = %21
%28 = load double, double* %7
%29 = fcmp oge double %28, 1.000000e+00
br i1 %29, label %30, label %40
30: ; preds = %27
br label %31
31: ; preds = %34, %30
%32 = load double, double* %7
%33 = fcmp oge double %32, 1.000000e+00
br i1 %33, label %34, label %39
34: ; preds = %31
%35 = load i32, i32* %9
%36 = add nsw i32 %35, 1
store i32 %36, i32* %9
%37 = load double, double* %7
%38 = fdiv double %37, 1.600000e+01
store double %38, double* %7
br label %31
39: ; preds = %31
br label %54
40: ; preds = %27
%41 = load double, double* %7
%42 = fcmp olt double %41, 6.250000e-02
br i1 %42, label %43, label %53
43: ; preds = %40
br label %44
44: ; preds = %47, %43
%45 = load double, double* %7
%46 = fcmp olt double %45, 6.250000e-02
br i1 %46, label %47, label %52
47: ; preds = %44
%48 = load i32, i32* %9
%49 = add nsw i32 %48, -1
store i32 %49, i32* %9
%50 = load double, double* %7
%51 = fmul double %50, 1.600000e+01
store double %51, double* %7
br label %44
52: ; preds = %44
br label %53
53: ; preds = %52, %40
br label %54
54: ; preds = %53, %39
br label %55
55: ; preds = %54, %26
%56 = load i32, i32* %9
%57 = trunc i32 %56 to i8
%58 = load i8*, i8** %5
%59 = getelementptr inbounds i8, i8* %58, i64 0
store i8 %57, i8* %59
%60 = load i8, i8* %10
%61 = icmp ne i8 %60, 0
br i1 %61, label %62, label %69
62: ; preds = %55
%63 = load i8*, i8** %5
%64 = getelementptr inbounds i8, i8* %63, i64 0
%65 = load i8, i8* %64
%66 = zext i8 %65 to i32
%67 = or i32 %66, 128
%68 = trunc i32 %67 to i8
store i8 %68, i8* %64
br label %69
69: ; preds = %62, %55
store i32 1, i32* %8
br label %70
70: ; preds = %92, %69
%71 = load i32, i32* %8
%72 = icmp slt i32 %71, 4
br i1 %72, label %73, label %95
73: ; preds = %70
%74 = load double, double* %7
%75 = fmul double %74, 2.560000e+02
store double %75, double* %7
%76 = load double, double* %7
%77 = fptosi double %76 to i32
%78 = trunc i32 %77 to i8
%79 = load i8*, i8** %5
%80 = load i32, i32* %8
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds i8, i8* %79, i64 %81
store i8 %78, i8* %82
%83 = load i8*, i8** %5
%84 = load i32, i32* %8
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds i8, i8* %83, i64 %85
%87 = load i8, i8* %86
%88 = zext i8 %87 to i32
%89 = sitofp i32 %88 to double
%90 = load double, double* %7
%91 = fsub double %90, %89
store double %91, double* %7
br label %92
92: ; preds = %73
%93 = load i32, i32* %8
%94 = add nsw i32 %93, 1
store i32 %94, i32* %8
br label %70
95: ; preds = %70
ret i32 4
}

LLVM IR (-O2):
@.str = [4 x i8] c"%lf\00"
define i32 @encode_GDS_R4(i8* nocapture %0, i8* nocapture writeonly %1, i32 %2) {
%4 = alloca double
%5 = bitcast double* %4 to i8*
store double 0.000000e+00, double* %4
%6 = call i32 (i8*, i8*, ...) @sscanf(i8* %0, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str, i64 0, i64 0), double* %4)
%7 = load double, double* %4
%8 = fcmp olt double %7, 0.000000e+00
%9 = fneg double %7
%10 = select i1 %8, double %9, double %7
%11 = fcmp oeq double %10, 0.000000e+00
br i1 %11, label %28, label %12
12: ; preds = %3
%13 = fcmp ult double %10, 1.000000e+00
br i1 %13, label %20, label %14
14: ; preds = %12, %14
%15 = phi i32 [ %17, %14 ], [ 64, %12 ]
%16 = phi double [ %18, %14 ], [ %10, %12 ]
%17 = add nuw nsw i32 %15, 1
%18 = fmul double %16, 6.250000e-02
%19 = fcmp ult double %18, 1.000000e+00
br i1 %19, label %28, label %14
20: ; preds = %12
%21 = fcmp olt double %10, 6.250000e-02
br i1 %21, label %22, label %28
22: ; preds = %20, %22
%23 = phi double [ %26, %22 ], [ %10, %20 ]
%24 = phi i32 [ %25, %22 ], [ 64, %20 ]
%25 = add nsw i32 %24, -1
%26 = fmul double %23, 1.600000e+01
%27 = fcmp olt double %26, 6.250000e-02
br i1 %27, label %22, label %28
28: ; preds = %14, %22, %20, %3
%29 = phi double [ %10, %3 ], [ %10, %20 ], [ %26, %22 ], [ %18, %14 ]
%30 = phi i32 [ 64, %3 ], [ 64, %20 ], [ %25, %22 ], [ %17, %14 ]
%31 = trunc i32 %30 to i8
%32 = or i8 %31, -128
%33 = select i1 %8, i8 %32, i8 %31
store i8 %33, i8* %1
%34 = fmul double %29, 2.560000e+02
%35 = fptosi double %34 to i32
%36 = trunc i32 %35 to i8
%37 = getelementptr inbounds i8, i8* %1, i64 1
store i8 %36, i8* %37
%38 = and i32 %35, 255
%39 = sitofp i32 %38 to double
%40 = fsub double %34, %39
%41 = fmul double %40, 2.560000e+02
%42 = fptosi double %41 to i32
%43 = trunc i32 %42 to i8
%44 = getelementptr inbounds i8, i8* %1, i64 2
store i8 %43, i8* %44
%45 = and i32 %42, 255
%46 = sitofp i32 %45 to double
%47 = fsub double %41, %46
%48 = fmul double %47, 2.560000e+02
%49 = fptosi double %48 to i32
%50 = trunc i32 %49 to i8
%51 = getelementptr inbounds i8, i8* %1, i64 3
store i8 %50, i8* %51
ret i32 4
}