#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ qd_add_d(double *, double, double *);
int /*<<< orphan*/ qd_minus(double *);
int /*<<< orphan*/ qd_mul(double *, double *, double *);
int /*<<< orphan*/ qd_mul_d(double *, double, double *);
int /*<<< orphan*/ qd_npwr(double *, int, double *);
int /*<<< orphan*/ qd_real(double, double *);
int sscanf(char const *, char *, int *);
int qd_read(const char *s, double *a) {
  const char *p = s;
  char ch;
  int sign = 0;
  int point = -1;
  int nd = 0;
  int e = 0;
  int done = 0;
  int nread;
  while (*p == ' ')
    p++;
  qd_real(0.0, a);
  while ((done == 0) && (ch = *p) != '\0') {
    if (ch >= '0' && ch <= '9') {
      int d = ch - '0';
      double cd = (double)d;
      qd_mul_d(a, 10.0, a);
      qd_add_d(a, cd, a);
      nd++;
    } else {
      switch (ch) {
      case '.':
        if (point >= 0)
          return -1;
        point = nd;
        break;
      case '-':
      case '+':
        if (sign != 0 || nd > 0)
          return -1;
        sign = (ch == '-') ? -1 : 1;
        break;
      case 'E':
      case 'e':
        nread = sscanf(p + 1, "%d", &e);
        done = 1;
        if (nread != 1)
          return -1;
        break;
      case ' ':
        done = 1;
        break;
      default:
        return -1;
      }
    }
    p++;
  }
  if (point >= 0)
    e -= (nd - point);
  if (e != 0) {
    double acc[4];
    qd_real(10.0, acc);
    qd_npwr(acc, e, acc);
    qd_mul(a, acc, a);
  }
  if (sign < 0)
    qd_minus(a);
  return 0;
}


x86 assembly (clang, -O0):
<qd_read>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0xffffffff,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
movsbl (%rax),%eax
cmp    $0x20,%eax
jne    11cf <qd_read+0x5f>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    11ae <qd_read+0x3e>
mov    -0x18(%rbp),%rdi
xorps  %xmm0,%xmm0
call   1080 <qd_real@plt>
xor    %eax,%eax
cmpl   $0x0,-0x38(%rbp)
mov    %al,-0x71(%rbp)
jne    11ff <qd_read+0x8f>
mov    -0x20(%rbp),%rax
mov    (%rax),%al
mov    %al,-0x21(%rbp)
movsbl %al,%eax
cmp    $0x0,%eax
setne  %al
mov    %al,-0x71(%rbp)
mov    -0x71(%rbp),%al
test   $0x1,%al
jne    120f <qd_read+0x9f>
jmp    13a4 <qd_read+0x234>
movsbl -0x21(%rbp),%eax
cmp    $0x30,%eax
jl     1272 <qd_read+0x102>
movsbl -0x21(%rbp),%eax
cmp    $0x39,%eax
jg     1272 <qd_read+0x102>
movsbl -0x21(%rbp),%eax
sub    $0x30,%eax
mov    %eax,-0x40(%rbp)
cvtsi2sdl -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x18(%rbp),%rsi
movsd  0xdb3(%rip),%xmm0
call   1060 <qd_mul_d@plt>
mov    -0x18(%rbp),%rdi
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rsi
call   1090 <qd_add_d@plt>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1393 <qd_read+0x223>
movsbl -0x21(%rbp),%eax
mov    %eax,-0x78(%rbp)
sub    $0x20,%eax
je     1376 <qd_read+0x206>
jmp    1287 <qd_read+0x117>
mov    -0x78(%rbp),%eax
sub    $0x2b,%eax
je     12fd <qd_read+0x18d>
jmp    1298 <qd_read+0x128>
mov    -0x78(%rbp),%eax
sub    $0x2d,%eax
je     12fd <qd_read+0x18d>
jmp    12a9 <qd_read+0x139>
mov    -0x78(%rbp),%eax
sub    $0x2e,%eax
je     12dc <qd_read+0x16c>
jmp    12ba <qd_read+0x14a>
mov    -0x78(%rbp),%eax
sub    $0x45,%eax
je     1339 <qd_read+0x1c9>
jmp    12cb <qd_read+0x15b>
mov    -0x78(%rbp),%eax
sub    $0x65,%eax
je     1339 <qd_read+0x1c9>
jmp    1382 <qd_read+0x212>
cmpl   $0x0,-0x2c(%rbp)
jl     12f2 <qd_read+0x182>
movl   $0xffffffff,-0x4(%rbp)
jmp    1412 <qd_read+0x2a2>
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    138e <qd_read+0x21e>
cmpl   $0x0,-0x28(%rbp)
jne    1311 <qd_read+0x1a1>
cmpl   $0x0,-0x30(%rbp)
jle    131d <qd_read+0x1ad>
movl   $0xffffffff,-0x4(%rbp)
jmp    1412 <qd_read+0x2a2>
movsbl -0x21(%rbp),%edx
mov    $0x1,%eax
mov    $0xffffffff,%ecx
cmp    $0x2d,%edx
cmove  %ecx,%eax
mov    %eax,-0x28(%rbp)
jmp    138e <qd_read+0x21e>
mov    -0x20(%rbp),%rdi
add    $0x1,%rdi
lea    0xcc0(%rip),%rsi
lea    -0x34(%rbp),%rdx
call   1050 <sscanf@plt>
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x38(%rbp)
cmpl   $0x1,-0x3c(%rbp)
je     1371 <qd_read+0x201>
movl   $0xffffffff,-0x4(%rbp)
jmp    1412 <qd_read+0x2a2>
jmp    138e <qd_read+0x21e>
movl   $0x1,-0x38(%rbp)
jmp    138e <qd_read+0x21e>
movl   $0xffffffff,-0x4(%rbp)
jmp    1412 <qd_read+0x2a2>
jmp    1393 <qd_read+0x223>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    11db <qd_read+0x6b>
cmpl   $0x0,-0x2c(%rbp)
jl     13bc <qd_read+0x24c>
mov    -0x30(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
je     13f8 <qd_read+0x288>
lea    -0x70(%rbp),%rdi
movsd  0xc2e(%rip),%xmm0
call   1080 <qd_real@plt>
lea    -0x70(%rbp),%rdi
mov    -0x34(%rbp),%esi
lea    -0x70(%rbp),%rdx
call   1040 <qd_npwr@plt>
mov    -0x18(%rbp),%rdi
lea    -0x70(%rbp),%rsi
mov    -0x18(%rbp),%rdx
call   1070 <qd_mul@plt>
cmpl   $0x0,-0x28(%rbp)
jge    140b <qd_read+0x29b>
mov    -0x18(%rbp),%rdi
call   1030 <qd_minus@plt>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qd_read>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rsi,%r14
mov    %rdi,%rbx
movl   $0x0,0x4(%rsp)
nopl   0x0(%rax)
cmpb   $0x20,(%rbx)
lea    0x1(%rbx),%rbx
je     1190 <qd_read+0x20>
xorps  %xmm0,%xmm0
mov    %r14,%rdi
call   1080 <qd_real@plt>
mov    -0x1(%rbx),%al
test   %al,%al
je     1263 <qd_read+0xf3>
xor    %r13d,%r13d
mov    $0xffffffff,%r12d
xor    %r15d,%r15d
mov    $0xffffffff,%ebp
jmp    1213 <qd_read+0xa3>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0xffffffd0,%ecx
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
movsd  %xmm0,0x8(%rsp)
mov    %r14,%rdi
movsd  0xe15(%rip),%xmm0
mov    %r14,%rsi
call   1060 <qd_mul_d@plt>
mov    %r14,%rdi
movsd  0x8(%rsp),%xmm0
mov    %r14,%rsi
call   1090 <qd_add_d@plt>
add    $0x1,%r13d
movzbl (%rbx),%eax
add    $0x1,%rbx
test   %al,%al
je     126f <qd_read+0xff>
movsbl %al,%ecx
lea    -0x30(%rcx),%edx
cmp    $0x9,%dl
jbe    11d0 <qd_read+0x60>
lea    -0x20(%rcx),%edx
cmp    $0x25,%edx
ja     12e4 <qd_read+0x174>
lea    0xdd7(%rip),%rsi
movslq (%rsi,%rdx,4),%rcx
add    %rsi,%rcx
jmp    *%rcx
test   %r15d,%r15d
jne    12c7 <qd_read+0x157>
test   %r13d,%r13d
jg     12c7 <qd_read+0x157>
xor    %ecx,%ecx
cmp    $0x2d,%al
setne  %cl
lea    (%rcx,%rcx,1),%r15d
add    $0xffffffff,%r15d
jmp    1208 <qd_read+0x98>
test   %r12d,%r12d
mov    %r13d,%r12d
js     1208 <qd_read+0x98>
jmp    12c7 <qd_read+0x157>
mov    $0xffffffff,%r12d
xor    %r15d,%r15d
xor    %r13d,%r13d
test   %r12d,%r12d
js     12d8 <qd_read+0x168>
sub    %r13d,%r12d
add    0x4(%rsp),%r12d
mov    %r12d,0x4(%rsp)
test   %r12d,%r12d
je     12b8 <qd_read+0x148>
movsd  0xd72(%rip),%xmm0
lea    0x10(%rsp),%rbx
mov    %rbx,%rdi
call   1080 <qd_real@plt>
mov    0x4(%rsp),%esi
mov    %rbx,%rdi
mov    %rbx,%rdx
call   1040 <qd_npwr@plt>
mov    %r14,%rdi
mov    %rbx,%rsi
mov    %r14,%rdx
call   1070 <qd_mul@plt>
xor    %ebp,%ebp
test   %r15d,%r15d
jns    12c7 <qd_read+0x157>
mov    %r14,%rdi
call   1030 <qd_minus@plt>
mov    %ebp,%eax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x4(%rsp),%r12d
test   %r12d,%r12d
jne    1286 <qd_read+0x116>
jmp    12b8 <qd_read+0x148>
cmp    $0x65,%ecx
jne    12c7 <qd_read+0x157>
lea    0xdb0(%rip),%rsi
lea    0x4(%rsp),%rdx
mov    %rbx,%rdi
call   1050 <sscanf@plt>
cmp    $0x1,%eax
je     126f <qd_read+0xff>
jmp    12c7 <qd_read+0x157>

LLVM IR (-O0):
@.str = [3 x i8] c"%d\00"
define i32 @qd_read(i8* %0, double* %1) {
%3 = alloca i32
%4 = alloca i8*
%5 = alloca double*
%6 = alloca i8*
%7 = alloca i8
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double
%16 = alloca [4 x double]
store i8* %0, i8** %4
store double* %1, double** %5
%17 = load i8*, i8** %4
store i8* %17, i8** %6
store i32 0, i32* %8
store i32 -1, i32* %9
store i32 0, i32* %10
store i32 0, i32* %11
store i32 0, i32* %12
br label %18
18: ; preds = %23, %2
%19 = load i8*, i8** %6
%20 = load i8, i8* %19
%21 = sext i8 %20 to i32
%22 = icmp eq i32 %21, 32
br i1 %22, label %23, label %26
23: ; preds = %18
%24 = load i8*, i8** %6
%25 = getelementptr inbounds i8, i8* %24, i32 1
store i8* %25, i8** %6
br label %18
26: ; preds = %18
%27 = load double*, double** %5
%28 = call i32 @qd_real(double 0.000000e+00, double* %27)
br label %29
29: ; preds = %95, %26
%30 = load i32, i32* %12
%31 = icmp eq i32 %30, 0
br i1 %31, label %32, label %37
32: ; preds = %29
%33 = load i8*, i8** %6
%34 = load i8, i8* %33
store i8 %34, i8* %7
%35 = sext i8 %34 to i32
%36 = icmp ne i32 %35, 0
br label %37
37: ; preds = %32, %29
%38 = phi i1 [ false, %29 ], [ %36, %32 ]
br i1 %38, label %39, label %98
39: ; preds = %37
%40 = load i8, i8* %7
%41 = sext i8 %40 to i32
%42 = icmp sge i32 %41, 48
br i1 %42, label %43, label %62
43: ; preds = %39
%44 = load i8, i8* %7
%45 = sext i8 %44 to i32
%46 = icmp sle i32 %45, 57
br i1 %46, label %47, label %62
47: ; preds = %43
%48 = load i8, i8* %7
%49 = sext i8 %48 to i32
%50 = sub nsw i32 %49, 48
store i32 %50, i32* %14
%51 = load i32, i32* %14
%52 = sitofp i32 %51 to double
store double %52, double* %15
%53 = load double*, double** %5
%54 = load double*, double** %5
%55 = call i32 @qd_mul_d(double* %53, double 1.000000e+01, double* %54)
%56 = load double*, double** %5
%57 = load double, double* %15
%58 = load double*, double** %5
%59 = call i32 @qd_add_d(double* %56, double %57, double* %58)
%60 = load i32, i32* %10
%61 = add nsw i32 %60, 1
store i32 %61, i32* %10
br label %95
62: ; preds = %43, %39
%63 = load i8, i8* %7
%64 = sext i8 %63 to i32
switch i32 %64, label %93 [
i32 46, label %65
i32 45, label %71
i32 43, label %71
i32 69, label %84
i32 101, label %84
i32 32, label %92
]
65: ; preds = %62
%66 = load i32, i32* %9
%67 = icmp sge i32 %66, 0
br i1 %67, label %68, label %69
68: ; preds = %65
store i32 -1, i32* %3
br label %128
69: ; preds = %65
%70 = load i32, i32* %10
store i32 %70, i32* %9
br label %94
71: ; preds = %62, %62
%72 = load i32, i32* %8
%73 = icmp ne i32 %72, 0
br i1 %73, label %77, label %74
74: ; preds = %71
%75 = load i32, i32* %10
%76 = icmp sgt i32 %75, 0
br i1 %76, label %77, label %78
77: ; preds = %74, %71
store i32 -1, i32* %3
br label %128
78: ; preds = %74
%79 = load i8, i8* %7
%80 = sext i8 %79 to i32
%81 = icmp eq i32 %80, 45
%82 = zext i1 %81 to i64
%83 = select i1 %81, i32 -1, i32 1
store i32 %83, i32* %8
br label %94
84: ; preds = %62, %62
%85 = load i8*, i8** %6
%86 = getelementptr inbounds i8, i8* %85, i64 1
%87 = call i32 @sscanf(i8* %86, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %11)
store i32 %87, i32* %13
store i32 1, i32* %12
%88 = load i32, i32* %13
%89 = icmp ne i32 %88, 1
br i1 %89, label %90, label %91
90: ; preds = %84
store i32 -1, i32* %3
br label %128
91: ; preds = %84
br label %94
92: ; preds = %62
store i32 1, i32* %12
br label %94
93: ; preds = %62
store i32 -1, i32* %3
br label %128
94: ; preds = %92, %91, %78, %69
br label %95
95: ; preds = %94, %47
%96 = load i8*, i8** %6
%97 = getelementptr inbounds i8, i8* %96, i32 1
store i8* %97, i8** %6
br label %29
98: ; preds = %37
%99 = load i32, i32* %9
%100 = icmp sge i32 %99, 0
br i1 %100, label %101, label %107
101: ; preds = %98
%102 = load i32, i32* %10
%103 = load i32, i32* %9
%104 = sub nsw i32 %102, %103
%105 = load i32, i32* %11
%106 = sub nsw i32 %105, %104
store i32 %106, i32* %11
br label %107
107: ; preds = %101, %98
%108 = load i32, i32* %11
%109 = icmp ne i32 %108, 0
br i1 %109, label %110, label %121
110: ; preds = %107
%111 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%112 = call i32 @qd_real(double 1.000000e+01, double* %111)
%113 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%114 = load i32, i32* %11
%115 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%116 = call i32 @qd_npwr(double* %113, i32 %114, double* %115)
%117 = load double*, double** %5
%118 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%119 = load double*, double** %5
%120 = call i32 @qd_mul(double* %117, double* %118, double* %119)
br label %121
121: ; preds = %110, %107
%122 = load i32, i32* %8
%123 = icmp slt i32 %122, 0
br i1 %123, label %124, label %127
124: ; preds = %121
%125 = load double*, double** %5
%126 = call i32 @qd_minus(double* %125)
br label %127
127: ; preds = %124, %121
store i32 0, i32* %3
br label %128
128: ; preds = %127, %93, %90, %77, %68
%129 = load i32, i32* %3
ret i32 %129
}

LLVM IR (-O2):
@.str = [3 x i8] c"%d\00"
define i32 @qd_read(i8* nocapture %0, double* %1) {
%3 = alloca i32
%4 = alloca [4 x double]
%5 = bitcast i32* %3 to i8*
store i32 0, i32* %3
br label %6
6: ; preds = %6, %2
%7 = phi i8* [ %0, %2 ], [ %10, %6 ]
%8 = load i8, i8* %7
%9 = icmp eq i8 %8, 32
%10 = getelementptr inbounds i8, i8* %7, i64 1
br i1 %9, label %6, label %11
11: ; preds = %6
%12 = tail call i32 @qd_real(double 0.000000e+00, double* %1)
%13 = load i8, i8* %7
%14 = icmp eq i8 %13, 0
br i1 %14, label %51, label %15
15: ; preds = %11, %44
%16 = phi i8 [ %49, %44 ], [ %13, %11 ]
%17 = phi i32 [ %45, %44 ], [ -1, %11 ]
%18 = phi i32 [ %46, %44 ], [ 0, %11 ]
%19 = phi i32 [ %47, %44 ], [ 0, %11 ]
%20 = phi i8* [ %48, %44 ], [ %7, %11 ]
%21 = sext i8 %16 to i32
%22 = add i8 %16, -48
%23 = icmp ult i8 %22, 10
br i1 %23, label %24, label %30
24: ; preds = %15
%25 = add nsw i32 %21, -48
%26 = sitofp i32 %25 to double
%27 = call i32 @qd_mul_d(double* %1, double 1.000000e+01, double* %1)
%28 = call i32 @qd_add_d(double* %1, double %26, double* %1)
%29 = add nsw i32 %18, 1
br label %44
30: ; preds = %15
switch i32 %21, label %76 [
i32 46, label %31
i32 45, label %33
i32 43, label %33
i32 69, label %40
i32 101, label %40
i32 32, label %51
]
31: ; preds = %30
%32 = icmp sgt i32 %17, -1
br i1 %32, label %76, label %44
33: ; preds = %30, %30
%34 = icmp ne i32 %19, 0
%35 = icmp sgt i32 %18, 0
%36 = select i1 %34, i1 true, i1 %35
br i1 %36, label %76, label %37
37: ; preds = %33
%38 = icmp eq i8 %16, 45
%39 = select i1 %38, i32 -1, i32 1
br label %44
40: ; preds = %30, %30
%41 = getelementptr inbounds i8, i8* %20, i64 1
%42 = call i32 @sscanf(i8* %41, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %3)
%43 = icmp eq i32 %42, 1
br i1 %43, label %51, label %76
44: ; preds = %31, %37, %24
%45 = phi i32 [ %17, %24 ], [ %17, %37 ], [ %18, %31 ]
%46 = phi i32 [ %29, %24 ], [ %18, %37 ], [ %18, %31 ]
%47 = phi i32 [ %19, %24 ], [ %39, %37 ], [ %19, %31 ]
%48 = getelementptr inbounds i8, i8* %20, i64 1
%49 = load i8, i8* %48
%50 = icmp eq i8 %49, 0
br i1 %50, label %51, label %15
51: ; preds = %44, %30, %11, %40
%52 = phi i32 [ %19, %40 ], [ 0, %11 ], [ %19, %30 ], [ %47, %44 ]
%53 = phi i32 [ %18, %40 ], [ 0, %11 ], [ %18, %30 ], [ %46, %44 ]
%54 = phi i32 [ %17, %40 ], [ -1, %11 ], [ %17, %30 ], [ %45, %44 ]
%55 = icmp sgt i32 %54, -1
br i1 %55, label %56, label %60
56: ; preds = %51
%57 = sub i32 %54, %53
%58 = load i32, i32* %3
%59 = add i32 %58, %57
store i32 %59, i32* %3
br label %62
60: ; preds = %51
%61 = load i32, i32* %3
br label %62
62: ; preds = %60, %56
%63 = phi i32 [ %61, %60 ], [ %59, %56 ]
%64 = icmp eq i32 %63, 0
br i1 %64, label %72, label %65
65: ; preds = %62
%66 = bitcast [4 x double]* %4 to i8*
%67 = getelementptr inbounds [4 x double], [4 x double]* %4, i64 0, i64 0
%68 = call i32 @qd_real(double 1.000000e+01, double* %67)
%69 = load i32, i32* %3
%70 = call i32 @qd_npwr(double* %67, i32 %69, double* %67)
%71 = call i32 @qd_mul(double* %1, double* %67, double* %1)
br label %72
72: ; preds = %65, %62
%73 = icmp slt i32 %52, 0
br i1 %73, label %74, label %76
74: ; preds = %72
%75 = call i32 @qd_minus(double* %1)
br label %76
76: ; preds = %30, %33, %31, %40, %72, %74
%77 = phi i32 [ 0, %74 ], [ 0, %72 ], [ -1, %40 ], [ -1, %31 ], [ -1, %33 ], [ -1, %30 ]
ret i32 %77
}