void rdfmaketable() {
  double slo = exp(-3.0), shi = exp(3.0), sinc = exp(0.2);
  const double blor = exp(0.0), bhir = exp(3.0), bincr = exp(0.2);
  const double blob = 0, bhib = 1.0, bincb = 0.1;
  double *r, *rdfa, *rdfd, dr, s, b, flux, eps;
  int i, n, done;
  char mode, dir, dirs[256];
  printf("Function for calculating radial diffusion functions (rdf) and "
         "reactive\n");
  printf("fluxes for alternating reaction and diffusion steps.  This module\n");
  printf(
      "can operate in several modes.  Enter (i) for irreversible reactions\n");
  printf(
      "(r) for reversible reactions with the unbinding radius larger than\n");
  printf("the binding radius, or (b) for other reversible reactions.  Enter "
         "this\n");
  printf("mode in upper case for machine readable output.\n");
  printf("Operation mode: ");
  scanf("%c", &mode);
  printf("Enter the number of radial points in the rdf (e.g. 200): ");
  scanf("%i", &n);
  if (n < 10) {
    printf("Value is too low.  Function stopped.\n");
    return;
  }
  printf("Enter level of precision (e.g. 1e-4): ");
  scanf("%lf", &eps);
  if (eps <= 0) {
    printf("Impossible precision.  Function stopped.\n");
    return;
  }
  printf("Enter u for increasing step lengths, d for decreasing: ");
  scanf("%s", dirs);
  dir = dirs[0];
  if (dir == 'd') {
    s = slo;
    slo = shi;
    shi = s;
    sinc = 1.0 / sinc;
  }
  r = (double *)calloc(n, sizeof(double));
  rdfa = (double *)calloc(n, sizeof(double));
  rdfd = (double *)calloc(n, sizeof(double));
  if (!r || !rdfa || !rdfd) {
    printf("Out of memory.  Function stopped.\n");
    return;
  }
  if (mode == 'i' || mode == 'I')
    b = -1;
  else if (mode == 'r' || mode == 'R')
    b = blor;
  else
    b = blob;
  done = 0;
  if (mode == 'i')
    printf("step     flux\n");
  else if (mode == 'r' || mode == 'b')
    printf("b      step       flux\n");
  else
    printf("\n");
  while (!done) {
    if (mode == 'i' || mode == 'I')
      dr = 10.0 / n;
    else if (mode == 'r' || mode == 'R')
      dr = (b + 3.0) / n;
    else
      dr = 5.0 / n;
    r[0] = 0;
    for (i = 1; r[i - 1] < 1 && i - 1 < n; i++)
      r[i] = r[i - 1] + dr;
    r[i - 1] = 0.9999;
    r[i++] = 1.0001;
    for (; i < n; i++)
      r[i] = r[i - 1] + dr;
    for (i = 0; i < n && r[i] < 1; i++)
      rdfa[i] = 0;
    if (dir == 'u' && (mode == 'i' || mode == 'I'))
      for (; i < n; i++)
        rdfa[i] = 1.0 - 1.0 / r[i];
    else if (dir == 'u' && (mode == 'r' || mode == 'R'))
      for (; i < n && r[i] < b; i++)
        rdfa[i] = 1.0 - (b - r[i]) / r[i] / (b - 1.0);
    for (; i < n; i++)
      rdfa[i] = 1.0;
    flux = 0;
    for (s = slo; slo < shi ? s < shi * sqrt(sinc) : s > shi * sqrt(sinc);
         s *= sinc) {
      if (flux >= 0)
        flux = rdfsteadystate(r, rdfa, rdfd, n, s, b, eps);
      if (mode == 'i')
        printf("%lf %lf\n", s, flux);
      else if (mode == 'r' || mode == 'b')
        printf("%lf %lf %lf\n", b, s, flux);
      else
        printf("%lf,", fabs(flux));
    }
    printf("\n");
    if (mode == 'i' || mode == 'I')
      done = 1;
    else if (mode == 'r' || mode == 'R') {
      b *= bincr;
      if (b > bhir * sqrt(bincr))
        done = 1;
    } else {
      b += bincb;
      if (b > bhib + bincb / 2.0)
        done = 1;
    }
  }
  free(r);
  free(rdfa);
  free(rdfd);
  return;
}


x86 assembly (clang, -O0):
<rdfmaketable>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1e0,%rsp
movsd  0xe9d(%rip),%xmm0
call   1060 <exp@plt>
movsd  %xmm0,-0x8(%rbp)
movsd  0xe83(%rip),%xmm0
call   1060 <exp@plt>
movsd  %xmm0,-0x10(%rbp)
movsd  0xe69(%rip),%xmm0
call   1060 <exp@plt>
movsd  %xmm0,-0x18(%rbp)
xorps  %xmm0,%xmm0
call   1060 <exp@plt>
movsd  %xmm0,-0x20(%rbp)
movsd  0xe52(%rip),%xmm0
call   1060 <exp@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  0xe38(%rip),%xmm0
call   1060 <exp@plt>
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe16(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe01(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
lea    0xe65(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xe9f(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xed4(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xf0b(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xf41(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xf7b(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xf9e(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xfa1(%rip),%rdi
lea    -0x95(%rbp),%rsi
mov    $0x0,%al
call   1070 <scanf@plt>
lea    0xf8f(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xfbb(%rip),%rdi
lea    -0x90(%rbp),%rsi
mov    $0x0,%al
call   1070 <scanf@plt>
cmpl   $0xa,-0x90(%rbp)
jge    12be <rdfmaketable+0x14e>
lea    0xf9c(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
jmp    1bf0 <rdfmaketable+0xa80>
lea    0xfaf(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xfc8(%rip),%rdi
lea    -0x88(%rbp),%rsi
mov    $0x0,%al
call   1070 <scanf@plt>
xorps  %xmm0,%xmm0
ucomisd -0x88(%rbp),%xmm0
ff
jb     1305 <rdfmaketable+0x195>
lea    0xfa6(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
jmp    1bf0 <rdfmaketable+0xa80>
lea    0xfbd(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    -0x1a0(%rbp),%rsi
lea    0xfe0(%rip),%rdi
mov    $0x0,%al
call   1070 <scanf@plt>
mov    -0x1a0(%rbp),%al
mov    %al,-0x96(%rbp)
movsbl -0x96(%rbp),%eax
cmp    $0x64,%eax
jne    1374 <rdfmaketable+0x204>
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  0xc9e(%rip),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movslq -0x90(%rbp),%rdi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x50(%rbp)
movslq -0x90(%rbp),%rdi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x58(%rbp)
movslq -0x90(%rbp),%rdi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x60(%rbp)
cmpq   $0x0,-0x50(%rbp)
je     13d4 <rdfmaketable+0x264>
cmpq   $0x0,-0x58(%rbp)
je     13d4 <rdfmaketable+0x264>
cmpq   $0x0,-0x60(%rbp)
jne    13e7 <rdfmaketable+0x277>
lea    0xf29(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
jmp    1bf0 <rdfmaketable+0xa80>
movsbl -0x95(%rbp),%eax
cmp    $0x69,%eax
je     1407 <rdfmaketable+0x297>
movsbl -0x95(%rbp),%eax
cmp    $0x49,%eax
jne    1419 <rdfmaketable+0x2a9>
movsd  0xc19(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    1455 <rdfmaketable+0x2e5>
movsbl -0x95(%rbp),%eax
cmp    $0x72,%eax
je     1439 <rdfmaketable+0x2c9>
movsbl -0x95(%rbp),%eax
cmp    $0x52,%eax
jne    1448 <rdfmaketable+0x2d8>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    1450 <rdfmaketable+0x2e0>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    1455 <rdfmaketable+0x2e5>
movl   $0x0,-0x94(%rbp)
movsbl -0x95(%rbp),%eax
cmp    $0x69,%eax
jne    1482 <rdfmaketable+0x312>
lea    0xeb1(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
jmp    14c8 <rdfmaketable+0x358>
movsbl -0x95(%rbp),%eax
cmp    $0x72,%eax
je     14a2 <rdfmaketable+0x332>
movsbl -0x95(%rbp),%eax
cmp    $0x62,%eax
jne    14b5 <rdfmaketable+0x345>
lea    0xe8d(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
jmp    14c3 <rdfmaketable+0x353>
lea    0xe0b(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
jmp    14c8 <rdfmaketable+0x358>
jmp    14cd <rdfmaketable+0x35d>
cmpl   $0x0,-0x94(%rbp)
setne  %al
xor    $0xff,%al
test   $0x1,%al
jne    14e6 <rdfmaketable+0x376>
jmp    1bd5 <rdfmaketable+0xa65>
movsbl -0x95(%rbp),%eax
cmp    $0x69,%eax
je     1506 <rdfmaketable+0x396>
movsbl -0x95(%rbp),%eax
cmp    $0x49,%eax
jne    1524 <rdfmaketable+0x3b4>
cvtsi2sdl -0x90(%rbp),%xmm1
ff
movsd  0xb22(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    1585 <rdfmaketable+0x415>
movsbl -0x95(%rbp),%eax
cmp    $0x72,%eax
je     1544 <rdfmaketable+0x3d4>
movsbl -0x95(%rbp),%eax
cmp    $0x52,%eax
jne    1567 <rdfmaketable+0x3f7>
movsd  0xacc(%rip),%xmm0
addsd  -0x78(%rbp),%xmm0
cvtsi2sdl -0x90(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    1580 <rdfmaketable+0x410>
cvtsi2sdl -0x90(%rbp),%xmm1
ff
movsd  0xab9(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    1585 <rdfmaketable+0x415>
mov    -0x50(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x8c(%rbp)
mov    -0x50(%rbp),%rcx
mov    -0x8c(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rdx
xor    %eax,%eax
movsd  0xa54(%rip),%xmm0
ucomisd (%rcx,%rdx,8),%xmm0
mov    %al,-0x1a1(%rbp)
jbe    15dd <rdfmaketable+0x46d>
mov    -0x8c(%rbp),%eax
sub    $0x1,%eax
cmp    -0x90(%rbp),%eax
setl   %al
mov    %al,-0x1a1(%rbp)
mov    -0x1a1(%rbp),%al
test   $0x1,%al
jne    15f0 <rdfmaketable+0x480>
jmp    162e <rdfmaketable+0x4be>
mov    -0x50(%rbp),%rax
mov    -0x8c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    159a <rdfmaketable+0x42a>
mov    -0x50(%rbp),%rax
mov    -0x8c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  0xa02(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
mov    -0x8c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x8c(%rbp)
movslq %ecx,%rcx
movsd  0x9d5(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
cmp    -0x90(%rbp),%eax
jge    16c0 <rdfmaketable+0x550>
mov    -0x50(%rbp),%rax
mov    -0x8c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    1670 <rdfmaketable+0x500>
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%ecx
xor    %eax,%eax
cmp    -0x90(%rbp),%ecx
mov    %al,-0x1a2(%rbp)
jge    1705 <rdfmaketable+0x595>
mov    -0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  0x911(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
seta   %al
mov    %al,-0x1a2(%rbp)
mov    -0x1a2(%rbp),%al
test   $0x1,%al
jne    1718 <rdfmaketable+0x5a8>
jmp    173f <rdfmaketable+0x5cf>
mov    -0x58(%rbp),%rax
movslq -0x8c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    16ca <rdfmaketable+0x55a>
movsbl -0x96(%rbp),%eax
cmp    $0x75,%eax
jne    17d3 <rdfmaketable+0x663>
movsbl -0x95(%rbp),%eax
cmp    $0x69,%eax
je     176f <rdfmaketable+0x5ff>
movsbl -0x95(%rbp),%eax
cmp    $0x49,%eax
jne    17d3 <rdfmaketable+0x663>
jmp    1774 <rdfmaketable+0x604>
mov    -0x8c(%rbp),%eax
cmp    -0x90(%rbp),%eax
jge    17ce <rdfmaketable+0x65e>
mov    -0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  0x86f(%rip),%xmm1
divsd  (%rax,%rcx,8),%xmm1
movsd  0x862(%rip),%xmm0
subsd  %xmm1,%xmm0
mov    -0x58(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    1774 <rdfmaketable+0x604>
jmp    18cb <rdfmaketable+0x75b>
movsbl -0x96(%rbp),%eax
cmp    $0x75,%eax
jne    18c6 <rdfmaketable+0x756>
movsbl -0x95(%rbp),%eax
cmp    $0x72,%eax
je     1803 <rdfmaketable+0x693>
movsbl -0x95(%rbp),%eax
cmp    $0x52,%eax
jne    18c6 <rdfmaketable+0x756>
jmp    1808 <rdfmaketable+0x698>
mov    -0x8c(%rbp),%ecx
xor    %eax,%eax
cmp    -0x90(%rbp),%ecx
mov    %al,-0x1a3(%rbp)
jge    1844 <rdfmaketable+0x6d4>
mov    -0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x78(%rbp),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x1a3(%rbp)
mov    -0x1a3(%rbp),%al
test   $0x1,%al
jne    1857 <rdfmaketable+0x6e7>
jmp    18c1 <rdfmaketable+0x751>
movsd  -0x78(%rbp),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
movsd  -0x78(%rbp),%xmm0
movsd  0x77f(%rip),%xmm2
subsd  %xmm2,%xmm0
divsd  %xmm0,%xmm1
movsd  0x76f(%rip),%xmm0
subsd  %xmm1,%xmm0
mov    -0x58(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    1808 <rdfmaketable+0x698>
jmp    18c6 <rdfmaketable+0x756>
jmp    18cb <rdfmaketable+0x75b>
jmp    18d0 <rdfmaketable+0x760>
mov    -0x8c(%rbp),%eax
cmp    -0x90(%rbp),%eax
jge    190e <rdfmaketable+0x79e>
mov    -0x58(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  0x713(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    18d0 <rdfmaketable+0x760>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x8(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1986 <rdfmaketable+0x816>
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x1b0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
call   1080 <sqrt@plt>
movsd  -0x1b8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x1b0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x1a8(%rbp)
jmp    19d3 <rdfmaketable+0x863>
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
call   1080 <sqrt@plt>
movsd  -0x1c8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x1c0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x1a8(%rbp)
mov    -0x1a8(%rbp),%eax
cmp    $0x0,%eax
je     1ad3 <rdfmaketable+0x963>
movsd  -0x80(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1a28 <rdfmaketable+0x8b8>
mov    -0x50(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0x60(%rbp),%rdx
mov    -0x90(%rbp),%ecx
movsd  -0x70(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm1
movsd  -0x88(%rbp),%xmm2
ff
mov    $0x3,%al
call   1090 <rdfsteadystate@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsbl -0x95(%rbp),%eax
cmp    $0x69,%eax
jne    1a55 <rdfmaketable+0x8e5>
movsd  -0x70(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm1
lea    0x909(%rip),%rdi
mov    $0x2,%al
call   1040 <printf@plt>
jmp    1aba <rdfmaketable+0x94a>
movsbl -0x95(%rbp),%eax
cmp    $0x72,%eax
je     1a75 <rdfmaketable+0x905>
movsbl -0x95(%rbp),%eax
cmp    $0x62,%eax
jne    1a97 <rdfmaketable+0x927>
movsd  -0x78(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm2
lea    0x8c3(%rip),%rdi
mov    $0x3,%al
call   1040 <printf@plt>
jmp    1ab5 <rdfmaketable+0x945>
movsd  -0x80(%rbp),%xmm0
movaps 0x5bd(%rip),%xmm1
pand   %xmm1,%xmm0
lea    0x8ad(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
jmp    1aba <rdfmaketable+0x94a>
jmp    1abf <rdfmaketable+0x94f>
movsd  -0x18(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    1920 <rdfmaketable+0x7b0>
lea    0x7ed(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movsbl -0x95(%rbp),%eax
cmp    $0x69,%eax
je     1b01 <rdfmaketable+0x991>
movsbl -0x95(%rbp),%eax
cmp    $0x49,%eax
jne    1b10 <rdfmaketable+0x9a0>
movl   $0x1,-0x94(%rbp)
jmp    1bd0 <rdfmaketable+0xa60>
movsbl -0x95(%rbp),%eax
cmp    $0x72,%eax
je     1b30 <rdfmaketable+0x9c0>
movsbl -0x95(%rbp),%eax
cmp    $0x52,%eax
jne    1b93 <rdfmaketable+0xa23>
movsd  -0x30(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
call   1080 <sqrt@plt>
movsd  -0x1d8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x1d0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    1b8e <rdfmaketable+0xa1e>
movl   $0x1,-0x94(%rbp)
jmp    1bcb <rdfmaketable+0xa5b>
movsd  0x465(%rip),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  0x49e(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1bc6 <rdfmaketable+0xa56>
movl   $0x1,-0x94(%rbp)
jmp    1bcb <rdfmaketable+0xa5b>
jmp    1bd0 <rdfmaketable+0xa60>
jmp    14cd <rdfmaketable+0x35d>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
add    $0x1e0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rdfmaketable>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x178,%rsp
lea    0xfc9(%rip),%rdi
call   1050 <puts@plt>
lea    0x1004(%rip),%rdi
call   1050 <puts@plt>
lea    0x103a(%rip),%rdi
call   1050 <puts@plt>
lea    0x1072(%rip),%rdi
call   1050 <puts@plt>
lea    0x10a9(%rip),%rdi
call   1050 <puts@plt>
lea    0x10e4(%rip),%rdi
call   1050 <puts@plt>
lea    0xeb8(%rip),%rdi
xor    %eax,%eax
call   1060 <printf@plt>
lea    0xebb(%rip),%rdi
lea    0x6(%rsp),%rsi
xor    %eax,%eax
call   1080 <scanf@plt>
lea    0xeab(%rip),%rdi
xor    %eax,%eax
call   1060 <printf@plt>
lea    0xed7(%rip),%rdi
lea    0x14(%rsp),%rsi
xor    %eax,%eax
call   1080 <scanf@plt>
cmpl   $0x9,0x14(%rsp)
jg     1233 <rdfmaketable+0xb3>
lea    0x112f(%rip),%rdi
call   1050 <puts@plt>
jmp    19bd <rdfmaketable+0x83d>
lea    0xeaf(%rip),%rdi
xor    %eax,%eax
call   1060 <printf@plt>
lea    0xec8(%rip),%rdi
lea    0x40(%rsp),%rsi
xor    %eax,%eax
call   1080 <scanf@plt>
xorpd  %xmm0,%xmm0
ucomisd 0x40(%rsp),%xmm0
jae    12b2 <rdfmaketable+0x132>
lea    0xead(%rip),%rdi
xor    %eax,%eax
call   1060 <printf@plt>
lea    0xed7(%rip),%rdi
lea    0x70(%rsp),%rsi
xor    %eax,%eax
call   1080 <scanf@plt>
mov    0x70(%rsp),%al
xor    %r12d,%r12d
movsd  0xda0(%rip),%xmm1
movsd  0xda0(%rip),%xmm0
mov    %al,0x7(%rsp)
cmp    $0x64,%al
sete   %r13b
je     12c3 <rdfmaketable+0x143>
movsd  %xmm1,0x28(%rsp)
movsd  %xmm0,0x38(%rsp)
jmp    12cf <rdfmaketable+0x14f>
lea    0x1076(%rip),%rdi
call   1050 <puts@plt>
jmp    19bd <rdfmaketable+0x83d>
movsd  %xmm1,0x38(%rsp)
movsd  %xmm0,0x28(%rsp)
movslq 0x14(%rsp),%rbp
mov    $0x8,%esi
mov    %rbp,%rdi
call   1070 <calloc@plt>
mov    %rax,%rbx
mov    $0x8,%esi
mov    %rbp,%rdi
call   1070 <calloc@plt>
mov    %rax,%r14
mov    $0x8,%esi
mov    %rbp,%rdi
call   1070 <calloc@plt>
test   %rbx,%rbx
je     132a <rdfmaketable+0x1aa>
test   %r14,%r14
je     132a <rdfmaketable+0x1aa>
mov    %rax,%r15
test   %rax,%rax
je     132a <rdfmaketable+0x1aa>
mov    %r13b,%r12b
mov    0x6(%rsp),%al
mov    %eax,%ecx
and    $0xdf,%cl
cmp    $0x52,%cl
je     133b <rdfmaketable+0x1bb>
xorpd  %xmm0,%xmm0
jmp    1343 <rdfmaketable+0x1c3>
lea    0xfb7(%rip),%rdi
call   1050 <puts@plt>
jmp    19bd <rdfmaketable+0x83d>
movsd  0xcfd(%rip),%xmm0
lea    0xcb6(%rip),%rdx
cmp    $0x49,%cl
jne    1357 <rdfmaketable+0x1d7>
movsd  0xcf1(%rip),%xmm0
movsd  (%rdx,%r12,8),%xmm1
movsd  %xmm1,0x20(%rsp)
cmp    $0x69,%al
movsd  %xmm0,0x18(%rsp)
jne    137b <rdfmaketable+0x1fb>
lea    0xfad(%rip),%rdi
call   1050 <puts@plt>
jmp    1399 <rdfmaketable+0x219>
and    $0xef,%al
cmp    $0x62,%al
jne    138f <rdfmaketable+0x20f>
lea    0xf82(%rip),%rdi
call   1050 <puts@plt>
jmp    1399 <rdfmaketable+0x219>
mov    $0xa,%edi
call   1040 <putchar@plt>
mov    0x6(%rsp),%dl
lea    0x10(%r14),%rax
mov    %rax,0x30(%rsp)
lea    0x10(%rbx),%rax
mov    %rax,0x48(%rsp)
lea    0x70(%r14),%rax
mov    %rax,0x50(%rsp)
movsd  0xc80(%rip),%xmm2
movabs $0x3ff0000000000000,%r13
movsd  0x20(%rsp),%xmm0
sqrtsd %xmm0,%xmm0
movsd  %xmm0,0x58(%rsp)
lea    0xd72(%rip),%r12
movsd  0x18(%rsp),%xmm3
jmp    1414 <rdfmaketable+0x294>
nopl   0x0(%rax)
movsd  0x18(%rsp),%xmm3
mulsd  0xc72(%rip),%xmm3
ucomisd 0xc72(%rip),%xmm3
movsd  0xc32(%rip),%xmm2
ja     19a5 <rdfmaketable+0x825>
and    $0xdf,%dl
movapd %xmm3,%xmm1
addsd  0xc2d(%rip),%xmm1
cmp    $0x52,%dl
je     1430 <rdfmaketable+0x2b0>
movsd  0xc28(%rip),%xmm1
movsd  0xc28(%rip),%xmm0
cmp    $0x49,%dl
je     1441 <rdfmaketable+0x2c1>
movapd %xmm1,%xmm0
movslq 0x14(%rsp),%r10
xorps  %xmm1,%xmm1
cvtsi2sd %r10d,%xmm1
mov    %r10d,%ecx
divsd  %xmm1,%xmm0
movq   $0x0,(%rbx)
test   %r10,%r10
jle    14a0 <rdfmaketable+0x320>
xorpd  %xmm1,%xmm1
mov    $0x1,%eax
nopw   0x0(%rax,%rax,1)
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
lea    0x1(%rax),%rsi
ucomisd %xmm1,%xmm2
jbe    148b <rdfmaketable+0x30b>
cmp    %rcx,%rax
mov    %rsi,%rax
jb     1470 <rdfmaketable+0x2f0>
lea    -0x1(%rsi),%rax
jmp    14a7 <rdfmaketable+0x327>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    $0x1,%esi
xor    %eax,%eax
movabs $0x3fefff2e48e8a71e,%rdi
ff ef 3f
mov    %rdi,(%rbx,%rax,8)
mov    %esi,%eax
movabs $0x3ff00068db8bac71,%rdi
mov    %rdi,(%rbx,%rax,8)
add    $0x1,%esi
cmp    %r10d,%esi
jge    1560 <rdfmaketable+0x3e0>
movapd %xmm0,%xmm1
addsd  0xb8b(%rip),%xmm1
movsd  %xmm1,0x8(%rbx,%rax,8)
lea    0x2(%rax),%rsi
cmp    %esi,%ecx
jle    1560 <rdfmaketable+0x3e0>
movsd  0x8(%rbx,%rax,8),%xmm1
mov    %r10d,%ebp
sub    %eax,%ebp
add    $0x2,%ebp
mov    %ecx,%edi
sub    %eax,%edi
add    $0xfffffffd,%edi
and    $0x3,%ebp
je     1522 <rdfmaketable+0x3a2>
cs nopw 0x0(%rax,%rax,1)
nop
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
add    $0x1,%rsi
add    $0xffffffff,%ebp
jne    1510 <rdfmaketable+0x390>
cmp    $0x3,%edi
jb     1560 <rdfmaketable+0x3e0>
nopw   0x0(%rax,%rax,1)
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rbx,%rsi,8)
addsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rbx,%rsi,8)
addsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rbx,%rsi,8)
add    $0x4,%rsi
cmp    %esi,%ecx
jg     1530 <rdfmaketable+0x3b0>
nop
xor    %esi,%esi
test   %ecx,%ecx
jle    158a <rdfmaketable+0x40a>
cs nopw 0x0(%rax,%rax,1)
ucomisd (%rbx,%rsi,8),%xmm2
jbe    158a <rdfmaketable+0x40a>
movq   $0x0,(%r14,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    1570 <rdfmaketable+0x3f0>
mov    %ecx,%esi
cmpb   $0x75,0x7(%rsp)
jne    1720 <rdfmaketable+0x5a0>
cmp    $0x49,%dl
jne    16b0 <rdfmaketable+0x530>
cmp    %ecx,%esi
jge    1720 <rdfmaketable+0x5a0>
mov    %esi,%edx
not    %esi
add    %r10d,%esi
je     1680 <rdfmaketable+0x500>
add    $0x1,%rsi
mov    %rsi,%r8
and    $0xfffffffffffffffe,%r8
lea    -0x2(%r8),%rax
mov    %rax,%r9
shr    %r9
add    $0x1,%r9
test   %rax,%rax
je     1994 <rdfmaketable+0x814>
mov    %r9,%rdi
and    $0xfffffffffffffffe,%rdi
mov    0x30(%rsp),%rax
lea    (%rax,%rdx,8),%r11
mov    0x48(%rsp),%rax
lea    (%rax,%rdx,8),%rax
xor    %ebp,%ebp
movapd 0xa18(%rip),%xmm4
nopl   0x0(%rax,%rax,1)
movupd -0x10(%rax,%rbp,8),%xmm0
movapd %xmm4,%xmm1
divpd  %xmm0,%xmm1
movapd %xmm4,%xmm0
subpd  %xmm1,%xmm0
movupd %xmm0,-0x10(%r11,%rbp,8)
movupd (%rax,%rbp,8),%xmm0
movapd %xmm4,%xmm1
divpd  %xmm0,%xmm1
movapd %xmm4,%xmm0
subpd  %xmm1,%xmm0
movupd %xmm0,(%r11,%rbp,8)
add    $0x4,%rbp
add    $0xfffffffffffffffe,%rdi
jne    1600 <rdfmaketable+0x480>
test   $0x1,%r9b
je     166e <rdfmaketable+0x4ee>
add    %rdx,%rbp
movupd (%rbx,%rbp,8),%xmm0
movapd 0x9b8(%rip),%xmm4
movapd %xmm4,%xmm1
divpd  %xmm0,%xmm1
movapd %xmm4,%xmm0
subpd  %xmm1,%xmm0
movupd %xmm0,(%r14,%rbp,8)
add    %r8,%rdx
cmp    %r8,%rsi
je     169f <rdfmaketable+0x51f>
cs nopw 0x0(%rax,%rax,1)
movapd %xmm2,%xmm0
divsd  (%rbx,%rdx,8),%xmm0
movapd %xmm2,%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rdx,8)
add    $0x1,%rdx
cmp    %edx,%ecx
jg     1680 <rdfmaketable+0x500>
mov    %edx,%esi
jmp    1720 <rdfmaketable+0x5a0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmp    $0x52,%dl
jne    1720 <rdfmaketable+0x5a0>
cmp    %ecx,%esi
jge    1720 <rdfmaketable+0x5a0>
movapd %xmm3,%xmm0
addsd  0x983(%rip),%xmm0
mov    %esi,%esi
lea    (%r14,%rsi,8),%rdi
lea    (%rbx,%rsi,8),%rbp
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rdx,8),%xmm1
ucomisd %xmm1,%xmm3
jbe    171d <rdfmaketable+0x59d>
movapd %xmm3,%xmm2
subsd  %xmm1,%xmm2
divsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
movsd  0x93c(%rip),%xmm1
subsd  %xmm2,%xmm1
movsd  %xmm1,(%rdi,%rdx,8)
add    $0x1,%rdx
lea    (%rdx,%rsi,1),%eax
cmp    %eax,%ecx
jg     16e0 <rdfmaketable+0x560>
jmp    17fd <rdfmaketable+0x67d>
add    %edx,%esi
nop
cmp    %ecx,%esi
jge    17fd <rdfmaketable+0x67d>
movslq %esi,%rcx
mov    %r10,%rdx
sub    %rcx,%rdx
cmp    $0x4,%rdx
jb     17f0 <rdfmaketable+0x670>
mov    %rdx,%r8
and    $0xfffffffffffffffc,%r8
lea    -0x4(%r8),%rax
mov    %rax,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%ebp
and    $0x3,%ebp
cmp    $0xc,%rax
jae    1768 <rdfmaketable+0x5e8>
xor    %esi,%esi
movapd 0x8aa(%rip),%xmm0
jmp    17b9 <rdfmaketable+0x639>
and    $0xfffffffffffffffc,%rdi
mov    0x50(%rsp),%rax
lea    (%rax,%rcx,8),%rax
xor    %esi,%esi
movapd 0x891(%rip),%xmm0
nop
movupd %xmm0,-0x70(%rax,%rsi,8)
movupd %xmm0,-0x60(%rax,%rsi,8)
movupd %xmm0,-0x50(%rax,%rsi,8)
movupd %xmm0,-0x40(%rax,%rsi,8)
movupd %xmm0,-0x30(%rax,%rsi,8)
movupd %xmm0,-0x20(%rax,%rsi,8)
movupd %xmm0,-0x10(%rax,%rsi,8)
movupd %xmm0,(%rax,%rsi,8)
add    $0x10,%rsi
add    $0xfffffffffffffffc,%rdi
jne    1780 <rdfmaketable+0x600>
test   %rbp,%rbp
je     17e4 <rdfmaketable+0x664>
add    %rcx,%rsi
mov    0x30(%rsp),%rax
lea    (%rax,%rsi,8),%rax
shl    $0x5,%rbp
xor    %esi,%esi
movupd %xmm0,-0x10(%rax,%rsi,1)
movupd %xmm0,(%rax,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%rbp
jne    17d0 <rdfmaketable+0x650>
cmp    %r8,%rdx
je     17fd <rdfmaketable+0x67d>
add    %r8,%rcx
nopl   0x0(%rax)
mov    %r13,(%r14,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r10
jne    17f0 <rdfmaketable+0x670>
movsd  %xmm3,0x18(%rsp)
xorpd  %xmm0,%xmm0
movapd %xmm0,0x60(%rsp)
movsd  0x28(%rsp),%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  0x20(%rsp),%xmm0
lea    0x929(%rip),%rbp
jmp    1860 <rdfmaketable+0x6e0>
nopl   0x0(%rax,%rax,1)
mov    %r12,%rdi
movsd  0x8(%rsp),%xmm0
movapd %xmm2,%xmm1
mov    $0x2,%al
call   1060 <printf@plt>
movsd  0x28(%rsp),%xmm1
movsd  0x20(%rsp),%xmm0
movsd  0x8(%rsp),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rsp)
ucomisd 0x828(%rip),%xmm0
movsd  0x58(%rsp),%xmm0
jae    1881 <rdfmaketable+0x701>
movsd  0x20(%rsp),%xmm0
call   1090 <sqrt@plt>
movsd  0x28(%rsp),%xmm1
movsd  0x38(%rsp),%xmm2
ucomisd %xmm1,%xmm2
mulsd  %xmm2,%xmm0
jbe    18a0 <rdfmaketable+0x720>
ucomisd 0x8(%rsp),%xmm0
ja     18b0 <rdfmaketable+0x730>
jmp    1950 <rdfmaketable+0x7d0>
xchg   %ax,%ax
movsd  0x8(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jbe    1950 <rdfmaketable+0x7d0>
movapd 0x60(%rsp),%xmm2
ucomisd 0x7d2(%rip),%xmm2
jb     18ed <rdfmaketable+0x76d>
mov    0x14(%rsp),%ecx
movsd  0x40(%rsp),%xmm2
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r15,%rdx
movsd  0x8(%rsp),%xmm0
movsd  0x18(%rsp),%xmm1
mov    $0x3,%al
call   10a0 <rdfsteadystate@plt>
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movzbl 0x6(%rsp),%eax
cmp    $0x69,%al
movapd %xmm2,0x60(%rsp)
je     1830 <rdfmaketable+0x6b0>
and    $0xef,%al
cmp    $0x62,%al
jne    1930 <rdfmaketable+0x7b0>
mov    %rbp,%rdi
movsd  0x18(%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
mov    $0x3,%al
call   1060 <printf@plt>
jmp    1844 <rdfmaketable+0x6c4>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
andpd  0x6e8(%rip),%xmm2
lea    0x81d(%rip),%rdi
movapd %xmm2,%xmm0
mov    $0x1,%al
call   1060 <printf@plt>
jmp    1844 <rdfmaketable+0x6c4>
nop
mov    $0xa,%edi
call   1040 <putchar@plt>
mov    0x6(%rsp),%dl
mov    %edx,%eax
and    $0xdf,%al
cmp    $0x52,%al
je     13f0 <rdfmaketable+0x270>
cmp    $0x49,%al
movsd  0x18(%rsp),%xmm3
movsd  0x6c6(%rip),%xmm2
je     19a5 <rdfmaketable+0x825>
addsd  0x6fc(%rip),%xmm3
ucomisd 0x6fc(%rip),%xmm3
jbe    1414 <rdfmaketable+0x294>
jmp    19a5 <rdfmaketable+0x825>
xor    %ebp,%ebp
test   $0x1,%r9b
jne    1648 <rdfmaketable+0x4c8>
jmp    166e <rdfmaketable+0x4ee>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
add    $0x178,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [72 x i8] c"Function for calculating radial diffusion functions (rdf) and reactive\0A\00"
@.str.1 = [67 x i8] c"fluxes for alternating reaction and diffusion steps. This module\0A\00"
@.str.2 = [69 x i8] c"can operate in several modes. Enter (i) for irreversible reactions\0A\00"
@.str.3 = [68 x i8] c"(r) for reversible reactions with the unbinding radius larger than\0A\00"
@.str.4 = [72 x i8] c"the binding radius, or (b) for other reversible reactions. Enter this\0A\00"
@.str.5 = [49 x i8] c"mode in upper case for machine readable output.\0A\00"
@.str.6 = [17 x i8] c"Operation mode: \00"
@.str.7 = [3 x i8] c"%c\00"
@.str.8 = [58 x i8] c"Enter the number of radial points in the rdf (e.g. 200): \00"
@.str.9 = [3 x i8] c"%i\00"
@.str.10 = [38 x i8] c"Value is too low. Function stopped.\0A\00"
@.str.11 = [39 x i8] c"Enter level of precision (e.g. 1e-4): \00"
@.str.12 = [4 x i8] c"%lf\00"
@.str.13 = [42 x i8] c"Impossible precision. Function stopped.\0A\00"
@.str.14 = [56 x i8] c"Enter u for increasing step lengths, d for decreasing: \00"
@.str.15 = [3 x i8] c"%s\00"
@.str.16 = [35 x i8] c"Out of memory. Function stopped.\0A\00"
@.str.17 = [15 x i8] c"step flux\0A\00"
@.str.18 = [24 x i8] c"b step flux\0A\00"
@.str.19 = [2 x i8] c"\0A\00"
@.str.20 = [9 x i8] c"%lf %lf\0A\00"
@.str.21 = [13 x i8] c"%lf %lf %lf\0A\00"
@.str.22 = [5 x i8] c"%lf,\00"
define void @rdfmaketable() {
%1 = alloca double
%2 = alloca double
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double*
%11 = alloca double*
%12 = alloca double*
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i8
%22 = alloca i8
%23 = alloca [256 x i8]
%24 = call double @exp(double -3.000000e+00)
store double %24, double* %1
%25 = call double @exp(double 3.000000e+00)
store double %25, double* %2
%26 = call double @exp(double 2.000000e-01)
store double %26, double* %3
%27 = call double @exp(double 0.000000e+00)
store double %27, double* %4
%28 = call double @exp(double 3.000000e+00)
store double %28, double* %5
%29 = call double @exp(double 2.000000e-01)
store double %29, double* %6
store double 0.000000e+00, double* %7
store double 1.000000e+00, double* %8
store double 1.000000e-01, double* %9
%30 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([72 x i8], [72 x i8]* @.str, i64 0, i64 0))
%31 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([67 x i8], [67 x i8]* @.str.1, i64 0, i64 0))
%32 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([69 x i8], [69 x i8]* @.str.2, i64 0, i64 0))
%33 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([68 x i8], [68 x i8]* @.str.3, i64 0, i64 0))
%34 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([72 x i8], [72 x i8]* @.str.4, i64 0, i64 0))
%35 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str.5, i64 0, i64 0))
%36 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str.6, i64 0, i64 0))
%37 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.7, i64 0, i64 0), i8* %21)
%38 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([58 x i8], [58 x i8]* @.str.8, i64 0, i64 0))
%39 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.9, i64 0, i64 0), i32* %19)
%40 = load i32, i32* %19
%41 = icmp slt i32 %40, 10
br i1 %41, label %42, label %44
42: ; preds = %0
%43 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([38 x i8], [38 x i8]* @.str.10, i64 0, i64 0))
br label %467
44: ; preds = %0
%45 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str.11, i64 0, i64 0))
%46 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.12, i64 0, i64 0), double* %17)
%47 = load double, double* %17
%48 = fcmp ole double %47, 0.000000e+00
br i1 %48, label %49, label %51
49: ; preds = %44
%50 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([42 x i8], [42 x i8]* @.str.13, i64 0, i64 0))
br label %467
51: ; preds = %44
%52 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([56 x i8], [56 x i8]* @.str.14, i64 0, i64 0))
%53 = getelementptr inbounds [256 x i8], [256 x i8]* %23, i64 0, i64 0
%54 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.15, i64 0, i64 0), i8* %53)
%55 = getelementptr inbounds [256 x i8], [256 x i8]* %23, i64 0, i64 0
%56 = load i8, i8* %55
store i8 %56, i8* %22
%57 = load i8, i8* %22
%58 = sext i8 %57 to i32
%59 = icmp eq i32 %58, 100
br i1 %59, label %60, label %66
60: ; preds = %51
%61 = load double, double* %1
store double %61, double* %14
%62 = load double, double* %2
store double %62, double* %1
%63 = load double, double* %14
store double %63, double* %2
%64 = load double, double* %3
%65 = fdiv double 1.000000e+00, %64
store double %65, double* %3
br label %66
66: ; preds = %60, %51
%67 = load i32, i32* %19
%68 = sext i32 %67 to i64
%69 = call i8* @calloc(i64 %68, i64 8)
%70 = bitcast i8* %69 to double*
store double* %70, double** %10
%71 = load i32, i32* %19
%72 = sext i32 %71 to i64
%73 = call i8* @calloc(i64 %72, i64 8)
%74 = bitcast i8* %73 to double*
store double* %74, double** %11
%75 = load i32, i32* %19
%76 = sext i32 %75 to i64
%77 = call i8* @calloc(i64 %76, i64 8)
%78 = bitcast i8* %77 to double*
store double* %78, double** %12
%79 = load double*, double** %10
%80 = icmp ne double* %79, null
br i1 %80, label %81, label %87
81: ; preds = %66
%82 = load double*, double** %11
%83 = icmp ne double* %82, null
br i1 %83, label %84, label %87
84: ; preds = %81
%85 = load double*, double** %12
%86 = icmp ne double* %85, null
br i1 %86, label %89, label %87
87: ; preds = %84, %81, %66
%88 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str.16, i64 0, i64 0))
br label %467
89: ; preds = %84
%90 = load i8, i8* %21
%91 = sext i8 %90 to i32
%92 = icmp eq i32 %91, 105
br i1 %92, label %97, label %93
93: ; preds = %89
%94 = load i8, i8* %21
%95 = sext i8 %94 to i32
%96 = icmp eq i32 %95, 73
br i1 %96, label %97, label %98
97: ; preds = %93, %89
store double -1.000000e+00, double* %15
br label %110
98: ; preds = %93
%99 = load i8, i8* %21
%100 = sext i8 %99 to i32
%101 = icmp eq i32 %100, 114
br i1 %101, label %106, label %102
102: ; preds = %98
%103 = load i8, i8* %21
%104 = sext i8 %103 to i32
%105 = icmp eq i32 %104, 82
br i1 %105, label %106, label %108
106: ; preds = %102, %98
%107 = load double, double* %4
store double %107, double* %15
br label %109
108: ; preds = %102
store double 0.000000e+00, double* %15
br label %109
109: ; preds = %108, %106
br label %110
110: ; preds = %109, %97
store i32 0, i32* %20
%111 = load i8, i8* %21
%112 = sext i8 %111 to i32
%113 = icmp eq i32 %112, 105
br i1 %113, label %114, label %116
114: ; preds = %110
%115 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str.17, i64 0, i64 0))
br label %129
116: ; preds = %110
%117 = load i8, i8* %21
%118 = sext i8 %117 to i32
%119 = icmp eq i32 %118, 114
br i1 %119, label %124, label %120
120: ; preds = %116
%121 = load i8, i8* %21
%122 = sext i8 %121 to i32
%123 = icmp eq i32 %122, 98
br i1 %123, label %124, label %126
124: ; preds = %120, %116
%125 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str.18, i64 0, i64 0))
br label %128
126: ; preds = %120
%127 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.19, i64 0, i64 0))
br label %128
128: ; preds = %126, %124
br label %129
129: ; preds = %128, %114
br label %130
130: ; preds = %459, %129
%131 = load i32, i32* %20
%132 = icmp ne i32 %131, 0
%133 = xor i1 %132, true
br i1 %133, label %134, label %460
134: ; preds = %130
%135 = load i8, i8* %21
%136 = sext i8 %135 to i32
%137 = icmp eq i32 %136, 105
br i1 %137, label %142, label %138
138: ; preds = %134
%139 = load i8, i8* %21
%140 = sext i8 %139 to i32
%141 = icmp eq i32 %140, 73
br i1 %141, label %142, label %146
142: ; preds = %138, %134
%143 = load i32, i32* %19
%144 = sitofp i32 %143 to double
%145 = fdiv double 1.000000e+01, %144
store double %145, double* %13
br label %165
146: ; preds = %138
%147 = load i8, i8* %21
%148 = sext i8 %147 to i32
%149 = icmp eq i32 %148, 114
br i1 %149, label %154, label %150
150: ; preds = %146
%151 = load i8, i8* %21
%152 = sext i8 %151 to i32
%153 = icmp eq i32 %152, 82
br i1 %153, label %154, label %160
154: ; preds = %150, %146
%155 = load double, double* %15
%156 = fadd double %155, 3.000000e+00
%157 = load i32, i32* %19
%158 = sitofp i32 %157 to double
%159 = fdiv double %156, %158
store double %159, double* %13
br label %164
160: ; preds = %150
%161 = load i32, i32* %19
%162 = sitofp i32 %161 to double
%163 = fdiv double 5.000000e+00, %162
store double %163, double* %13
br label %164
164: ; preds = %160, %154
br label %165
165: ; preds = %164, %142
%166 = load double*, double** %10
%167 = getelementptr inbounds double, double* %166, i64 0
store double 0.000000e+00, double* %167
store i32 1, i32* %18
br label %168
168: ; preds = %196, %165
%169 = load double*, double** %10
%170 = load i32, i32* %18
%171 = sub nsw i32 %170, 1
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %169, i64 %172
%174 = load double, double* %173
%175 = fcmp olt double %174, 1.000000e+00
br i1 %175, label %176, label %181
176: ; preds = %168
%177 = load i32, i32* %18
%178 = sub nsw i32 %177, 1
%179 = load i32, i32* %19
%180 = icmp slt i32 %178, %179
br label %181
181: ; preds = %176, %168
%182 = phi i1 [ false, %168 ], [ %180, %176 ]
br i1 %182, label %183, label %199
183: ; preds = %181
%184 = load double*, double** %10
%185 = load i32, i32* %18
%186 = sub nsw i32 %185, 1
%187 = sext i32 %186 to i64
%188 = getelementptr inbounds double, double* %184, i64 %187
%189 = load double, double* %188
%190 = load double, double* %13
%191 = fadd double %189, %190
%192 = load double*, double** %10
%193 = load i32, i32* %18
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds double, double* %192, i64 %194
store double %191, double* %195
br label %196
196: ; preds = %183
%197 = load i32, i32* %18
%198 = add nsw i32 %197, 1
store i32 %198, i32* %18
br label %168
199: ; preds = %181
%200 = load double*, double** %10
%201 = load i32, i32* %18
%202 = sub nsw i32 %201, 1
%203 = sext i32 %202 to i64
%204 = getelementptr inbounds double, double* %200, i64 %203
store double 9.999000e-01, double* %204
%205 = load double*, double** %10
%206 = load i32, i32* %18
%207 = add nsw i32 %206, 1
store i32 %207, i32* %18
%208 = sext i32 %206 to i64
%209 = getelementptr inbounds double, double* %205, i64 %208
store double 1.000100e+00, double* %209
br label %210
210: ; preds = %227, %199
%211 = load i32, i32* %18
%212 = load i32, i32* %19
%213 = icmp slt i32 %211, %212
br i1 %213, label %214, label %230
214: ; preds = %210
%215 = load double*, double** %10
%216 = load i32, i32* %18
%217 = sub nsw i32 %216, 1
%218 = sext i32 %217 to i64
%219 = getelementptr inbounds double, double* %215, i64 %218
%220 = load double, double* %219
%221 = load double, double* %13
%222 = fadd double %220, %221
%223 = load double*, double** %10
%224 = load i32, i32* %18
%225 = sext i32 %224 to i64
%226 = getelementptr inbounds double, double* %223, i64 %225
store double %222, double* %226
br label %227
227: ; preds = %214
%228 = load i32, i32* %18
%229 = add nsw i32 %228, 1
store i32 %229, i32* %18
br label %210
230: ; preds = %210
store i32 0, i32* %18
br label %231
231: ; preds = %249, %230
%232 = load i32, i32* %18
%233 = load i32, i32* %19
%234 = icmp slt i32 %232, %233
br i1 %234, label %235, label %242
235: ; preds = %231
%236 = load double*, double** %10
%237 = load i32, i32* %18
%238 = sext i32 %237 to i64
%239 = getelementptr inbounds double, double* %236, i64 %238
%240 = load double, double* %239
%241 = fcmp olt double %240, 1.000000e+00
br label %242
242: ; preds = %235, %231
%243 = phi i1 [ false, %231 ], [ %241, %235 ]
br i1 %243, label %244, label %252
244: ; preds = %242
%245 = load double*, double** %11
%246 = load i32, i32* %18
%247 = sext i32 %246 to i64
%248 = getelementptr inbounds double, double* %245, i64 %247
store double 0.000000e+00, double* %248
br label %249
249: ; preds = %244
%250 = load i32, i32* %18
%251 = add nsw i32 %250, 1
store i32 %251, i32* %18
br label %231
252: ; preds = %242
%253 = load i8, i8* %22
%254 = sext i8 %253 to i32
%255 = icmp eq i32 %254, 117
br i1 %255, label %256, label %285
256: ; preds = %252
%257 = load i8, i8* %21
%258 = sext i8 %257 to i32
%259 = icmp eq i32 %258, 105
br i1 %259, label %264, label %260
260: ; preds = %256
%261 = load i8, i8* %21
%262 = sext i8 %261 to i32
%263 = icmp eq i32 %262, 73
br i1 %263, label %264, label %285
264: ; preds = %260, %256
br label %265
265: ; preds = %281, %264
%266 = load i32, i32* %18
%267 = load i32, i32* %19
%268 = icmp slt i32 %266, %267
br i1 %268, label %269, label %284
269: ; preds = %265
%270 = load double*, double** %10
%271 = load i32, i32* %18
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds double, double* %270, i64 %272
%274 = load double, double* %273
%275 = fdiv double 1.000000e+00, %274
%276 = fsub double 1.000000e+00, %275
%277 = load double*, double** %11
%278 = load i32, i32* %18
%279 = sext i32 %278 to i64
%280 = getelementptr inbounds double, double* %277, i64 %279
store double %276, double* %280
br label %281
281: ; preds = %269
%282 = load i32, i32* %18
%283 = add nsw i32 %282, 1
store i32 %283, i32* %18
br label %265
284: ; preds = %265
br label %339
285: ; preds = %260, %252
%286 = load i8, i8* %22
%287 = sext i8 %286 to i32
%288 = icmp eq i32 %287, 117
br i1 %288, label %289, label %338
289: ; preds = %285
%290 = load i8, i8* %21
%291 = sext i8 %290 to i32
%292 = icmp eq i32 %291, 114
br i1 %292, label %297, label %293
293: ; preds = %289
%294 = load i8, i8* %21
%295 = sext i8 %294 to i32
%296 = icmp eq i32 %295, 82
br i1 %296, label %297, label %338
297: ; preds = %293, %289
br label %298
298: ; preds = %334, %297
%299 = load i32, i32* %18
%300 = load i32, i32* %19
%301 = icmp slt i32 %299, %300
br i1 %301, label %302, label %310
302: ; preds = %298
%303 = load double*, double** %10
%304 = load i32, i32* %18
%305 = sext i32 %304 to i64
%306 = getelementptr inbounds double, double* %303, i64 %305
%307 = load double, double* %306
%308 = load double, double* %15
%309 = fcmp olt double %307, %308
br label %310
310: ; preds = %302, %298
%311 = phi i1 [ false, %298 ], [ %309, %302 ]
br i1 %311, label %312, label %337
312: ; preds = %310
%313 = load double, double* %15
%314 = load double*, double** %10
%315 = load i32, i32* %18
%316 = sext i32 %315 to i64
%317 = getelementptr inbounds double, double* %314, i64 %316
%318 = load double, double* %317
%319 = fsub double %313, %318
%320 = load double*, double** %10
%321 = load i32, i32* %18
%322 = sext i32 %321 to i64
%323 = getelementptr inbounds double, double* %320, i64 %322
%324 = load double, double* %323
%325 = fdiv double %319, %324
%326 = load double, double* %15
%327 = fsub double %326, 1.000000e+00
%328 = fdiv double %325, %327
%329 = fsub double 1.000000e+00, %328
%330 = load double*, double** %11
%331 = load i32, i32* %18
%332 = sext i32 %331 to i64
%333 = getelementptr inbounds double, double* %330, i64 %332
store double %329, double* %333
br label %334
334: ; preds = %312
%335 = load i32, i32* %18
%336 = add nsw i32 %335, 1
store i32 %336, i32* %18
br label %298
337: ; preds = %310
br label %338
338: ; preds = %337, %293, %285
br label %339
339: ; preds = %338, %284
br label %340
340: ; preds = %349, %339
%341 = load i32, i32* %18
%342 = load i32, i32* %19
%343 = icmp slt i32 %341, %342
br i1 %343, label %344, label %352
344: ; preds = %340
%345 = load double*, double** %11
%346 = load i32, i32* %18
%347 = sext i32 %346 to i64
%348 = getelementptr inbounds double, double* %345, i64 %347
store double 1.000000e+00, double* %348
br label %349
349: ; preds = %344
%350 = load i32, i32* %18
%351 = add nsw i32 %350, 1
store i32 %351, i32* %18
br label %340
352: ; preds = %340
store double 0.000000e+00, double* %16
%353 = load double, double* %1
store double %353, double* %14
br label %354
354: ; preds = %417, %352
%355 = load double, double* %1
%356 = load double, double* %2
%357 = fcmp olt double %355, %356
br i1 %357, label %358, label %366
358: ; preds = %354
%359 = load double, double* %14
%360 = load double, double* %2
%361 = load double, double* %3
%362 = call double @sqrt(double %361)
%363 = fmul double %360, %362
%364 = fcmp olt double %359, %363
%365 = zext i1 %364 to i32
br label %374
366: ; preds = %354
%367 = load double, double* %14
%368 = load double, double* %2
%369 = load double, double* %3
%370 = call double @sqrt(double %369)
%371 = fmul double %368, %370
%372 = fcmp ogt double %367, %371
%373 = zext i1 %372 to i32
br label %374
374: ; preds = %366, %358
%375 = phi i32 [ %365, %358 ], [ %373, %366 ]
%376 = icmp ne i32 %375, 0
br i1 %376, label %377, label %421
377: ; preds = %374
%378 = load double, double* %16
%379 = fcmp oge double %378, 0.000000e+00
br i1 %379, label %380, label %390
380: ; preds = %377
%381 = load double*, double** %10
%382 = load double*, double** %11
%383 = load double*, double** %12
%384 = load i32, i32* %19
%385 = load double, double* %14
%386 = load double, double* %15
%387 = load double, double* %17
%388 = call i32 (double*, double*, double*, i32, double, double, double, ...) bitcast (i32 (...)* @rdfsteadystate to i32 (double*, double*, double*, i32, double, double, double, ...)*)(double* %381, double* %382, double* %383, i32 %384, double %385, double %386, double %387)
%389 = sitofp i32 %388 to double
store double %389, double* %16
br label %390
390: ; preds = %380, %377
%391 = load i8, i8* %21
%392 = sext i8 %391 to i32
%393 = icmp eq i32 %392, 105
br i1 %393, label %394, label %398
394: ; preds = %390
%395 = load double, double* %14
%396 = load double, double* %16
%397 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([9 x i8], [9 x i8]* @.str.20, i64 0, i64 0), double %395, double %396)
br label %416
398: ; preds = %390
%399 = load i8, i8* %21
%400 = sext i8 %399 to i32
%401 = icmp eq i32 %400, 114
br i1 %401, label %406, label %402
402: ; preds = %398
%403 = load i8, i8* %21
%404 = sext i8 %403 to i32
%405 = icmp eq i32 %404, 98
br i1 %405, label %406, label %411
406: ; preds = %402, %398
%407 = load double, double* %15
%408 = load double, double* %14
%409 = load double, double* %16
%410 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.21, i64 0, i64 0), double %407, double %408, double %409)
br label %415
411: ; preds = %402
%412 = load double, double* %16
%413 = call double @llvm.fabs.f64(double %412)
%414 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([5 x i8], [5 x i8]* @.str.22, i64 0, i64 0), double %413)
br label %415
415: ; preds = %411, %406
br label %416
416: ; preds = %415, %394
br label %417
417: ; preds = %416
%418 = load double, double* %3
%419 = load double, double* %14
%420 = fmul double %419, %418
store double %420, double* %14
br label %354
421: ; preds = %374
%422 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.19, i64 0, i64 0))
%423 = load i8, i8* %21
%424 = sext i8 %423 to i32
%425 = icmp eq i32 %424, 105
br i1 %425, label %430, label %426
426: ; preds = %421
%427 = load i8, i8* %21
%428 = sext i8 %427 to i32
%429 = icmp eq i32 %428, 73
br i1 %429, label %430, label %431
430: ; preds = %426, %421
store i32 1, i32* %20
br label %459
431: ; preds = %426
%432 = load i8, i8* %21
%433 = sext i8 %432 to i32
%434 = icmp eq i32 %433, 114
br i1 %434, label %439, label %435
435: ; preds = %431
%436 = load i8, i8* %21
%437 = sext i8 %436 to i32
%438 = icmp eq i32 %437, 82
br i1 %438, label %439, label %451
439: ; preds = %435, %431
%440 = load double, double* %6
%441 = load double, double* %15
%442 = fmul double %441, %440
store double %442, double* %15
%443 = load double, double* %15
%444 = load double, double* %5
%445 = load double, double* %6
%446 = call double @sqrt(double %445)
%447 = fmul double %444, %446
%448 = fcmp ogt double %443, %447
br i1 %448, label %449, label %450
449: ; preds = %439
store i32 1, i32* %20
br label %450
450: ; preds = %449, %439
br label %458
451: ; preds = %435
%452 = load double, double* %15
%453 = fadd double %452, 1.000000e-01
store double %453, double* %15
%454 = load double, double* %15
%455 = fcmp ogt double %454, 1.050000e+00
br i1 %455, label %456, label %457
456: ; preds = %451
store i32 1, i32* %20
br label %457
457: ; preds = %456, %451
br label %458
458: ; preds = %457, %450
br label %459
459: ; preds = %458, %430
br label %130
460: ; preds = %130
%461 = load double*, double** %10
%462 = bitcast double* %461 to i8*
%463 = load double*, double** %11
%464 = bitcast double* %463 to i8*
%465 = load double*, double** %12
%466 = bitcast double* %465 to i8*
br label %467
467: ; preds = %460, %87, %49, %42
ret void
}

LLVM IR (-O2):
@.str.6 = [17 x i8] c"Operation mode: \00"
@.str.7 = [3 x i8] c"%c\00"
@.str.8 = [58 x i8] c"Enter the number of radial points in the rdf (e.g. 200): \00"
@.str.9 = [3 x i8] c"%i\00"
@.str.11 = [39 x i8] c"Enter level of precision (e.g. 1e-4): \00"
@.str.12 = [4 x i8] c"%lf\00"
@.str.14 = [56 x i8] c"Enter u for increasing step lengths, d for decreasing: \00"
@.str.15 = [3 x i8] c"%s\00"
@.str.20 = [9 x i8] c"%lf %lf\0A\00"
@.str.21 = [13 x i8] c"%lf %lf %lf\0A\00"
@.str.22 = [5 x i8] c"%lf,\00"
@str = [71 x i8] c"Function for calculating radial diffusion functions (rdf) and reactive\00"
@str.23 = [66 x i8] c"fluxes for alternating reaction and diffusion steps. This module\00"
@str.24 = [68 x i8] c"can operate in several modes. Enter (i) for irreversible reactions\00"
@str.25 = [67 x i8] c"(r) for reversible reactions with the unbinding radius larger than\00"
@str.26 = [71 x i8] c"the binding radius, or (b) for other reversible reactions. Enter this\00"
@str.27 = [48 x i8] c"mode in upper case for machine readable output.\00"
@str.28 = [34 x i8] c"Out of memory. Function stopped.\00"
@str.29 = [23 x i8] c"b step flux\00"
@str.30 = [14 x i8] c"step flux\00"
@str.31 = [41 x i8] c"Impossible precision. Function stopped.\00"
@str.32 = [37 x i8] c"Value is too low. Function stopped.\00"
define void @rdfmaketable() {
%1 = alloca double
%2 = alloca i32
%3 = alloca i8
%4 = alloca [256 x i8]
%5 = bitcast double* %1 to i8*
%6 = bitcast i32* %2 to i8*
%7 = getelementptr inbounds [256 x i8], [256 x i8]* %4, i64 0, i64 0
%8 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([71 x i8], [71 x i8]* @str, i64 0, i64 0))
%9 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([66 x i8], [66 x i8]* @str.23, i64 0, i64 0))
%10 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([68 x i8], [68 x i8]* @str.24, i64 0, i64 0))
%11 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([67 x i8], [67 x i8]* @str.25, i64 0, i64 0))
%12 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([71 x i8], [71 x i8]* @str.26, i64 0, i64 0))
%13 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([48 x i8], [48 x i8]* @str.27, i64 0, i64 0))
%14 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([17 x i8], [17 x i8]* @.str.6, i64 0, i64 0))
%15 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.7, i64 0, i64 0), i8* %3)
%16 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([58 x i8], [58 x i8]* @.str.8, i64 0, i64 0))
%17 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.9, i64 0, i64 0), i32* %2)
%18 = load i32, i32* %2
%19 = icmp slt i32 %18, 10
br i1 %19, label %20, label %22
20: ; preds = %0
%21 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([37 x i8], [37 x i8]* @str.32, i64 0, i64 0))
br label %383
22: ; preds = %0
%23 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([39 x i8], [39 x i8]* @.str.11, i64 0, i64 0))
%24 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.12, i64 0, i64 0), double* %1)
%25 = load double, double* %1
%26 = fcmp ugt double %25, 0.000000e+00
br i1 %26, label %29, label %27
27: ; preds = %22
%28 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([41 x i8], [41 x i8]* @str.31, i64 0, i64 0))
br label %383
29: ; preds = %22
%30 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([56 x i8], [56 x i8]* @.str.14, i64 0, i64 0))
%31 = call i32 (i8*, ...) @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.15, i64 0, i64 0), i8* %7)
%32 = load i8, i8* %7
%33 = icmp eq i8 %32, 100
%34 = select i1 %33, double 0x3FEA330AD6166159, double 0x3FF38ADD9E58AC8D
%35 = select i1 %33, double 0x3FA97DB0CCCEB0AF, double 0x403415E5BF6FB106
%36 = select i1 %33, double 0x403415E5BF6FB106, double 0x3FA97DB0CCCEB0AF
%37 = load i32, i32* %2
%38 = sext i32 %37 to i64
%39 = call i8* @calloc(i64 %38, i64 8)
%40 = bitcast i8* %39 to double*
%41 = call i8* @calloc(i64 %38, i64 8)
%42 = bitcast i8* %41 to double*
%43 = call i8* @calloc(i64 %38, i64 8)
%44 = bitcast i8* %43 to double*
%45 = icmp ne i8* %39, null
%46 = icmp ne i8* %41, null
%47 = and i1 %45, %46
%48 = icmp ne i8* %43, null
%49 = and i1 %47, %48
br i1 %49, label %52, label %50
50: ; preds = %29
%51 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([34 x i8], [34 x i8]* @str.28, i64 0, i64 0))
br label %383
52: ; preds = %29
%53 = load i8, i8* %3
%54 = and i8 %53, -33
%55 = icmp eq i8 %54, 73
%56 = icmp eq i8 %54, 82
%57 = select i1 %56, double 1.000000e+00, double 0.000000e+00
%58 = select i1 %55, double -1.000000e+00, double %57
%59 = icmp eq i8 %53, 105
br i1 %59, label %60, label %62
60: ; preds = %52
%61 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([14 x i8], [14 x i8]* @str.30, i64 0, i64 0))
br label %69
62: ; preds = %52
%63 = and i8 %53, -17
%64 = icmp eq i8 %63, 98
br i1 %64, label %65, label %67
65: ; preds = %62
%66 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([23 x i8], [23 x i8]* @str.29, i64 0, i64 0))
br label %69
67: ; preds = %62
%68 = call i32 @putchar(i32 10)
br label %69
69: ; preds = %65, %67, %60
%70 = icmp eq i8 %32, 117
%71 = fcmp olt double %36, %35
%72 = load i8, i8* %3
%73 = getelementptr i8, i8* %39, i64 8
br label %74
74: ; preds = %377, %69
%75 = phi i8 [ %72, %69 ], [ %372, %377 ]
%76 = phi double [ %58, %69 ], [ %378, %377 ]
%77 = and i8 %75, -33
%78 = icmp eq i8 %77, 73
%79 = icmp eq i8 %77, 82
%80 = fadd double %76, 3.000000e+00
%81 = select i1 %79, double %80, double 5.000000e+00
%82 = select i1 %78, double 1.000000e+01, double %81
%83 = load i32, i32* %2
%84 = sitofp i32 %83 to double
%85 = fdiv double %82, %84
store double 0.000000e+00, double* %40
%86 = icmp sgt i32 %83, 0
br i1 %86, label %87, label %100
87: ; preds = %74
%88 = zext i32 %83 to i64
br label %89
89: ; preds = %87, %89
%90 = phi i64 [ 1, %87 ], [ %94, %89 ]
%91 = phi double [ 0.000000e+00, %87 ], [ %92, %89 ]
%92 = fadd double %85, %91
%93 = getelementptr inbounds double, double* %40, i64 %90
store double %92, double* %93
%94 = add nuw nsw i64 %90, 1
%95 = fcmp olt double %92, 1.000000e+00
%96 = icmp ult i64 %90, %88
%97 = select i1 %95, i1 %96, i1 false
br i1 %97, label %89, label %98
98: ; preds = %89
%99 = trunc i64 %94 to i32
br label %100
100: ; preds = %98, %74
%101 = phi i64 [ 0, %74 ], [ %90, %98 ]
%102 = phi i32 [ 1, %74 ], [ %99, %98 ]
%103 = getelementptr inbounds double, double* %40, i64 %101
store double 9.999000e-01, double* %103
%104 = zext i32 %102 to i64
%105 = getelementptr inbounds double, double* %40, i64 %104
store double 1.000100e+00, double* %105
%106 = add nuw nsw i32 %102, 1
%107 = icmp slt i32 %106, %83
br i1 %107, label %108, label %139
108: ; preds = %100
%109 = add nuw nsw i64 %104, 1
%110 = fadd double %85, 1.000100e+00
%111 = getelementptr inbounds double, double* %40, i64 %109
store double %110, double* %111
%112 = add nuw nsw i64 %104, 2
%113 = trunc i64 %112 to i32
%114 = icmp sgt i32 %83, %113
br i1 %114, label %115, label %139
115: ; preds = %108
%116 = shl nuw nsw i64 %104, 3
%117 = getelementptr i8, i8* %73, i64 %116
%118 = bitcast i8* %117 to double*
%119 = load double, double* %118
%120 = add i32 %83, 2
%121 = sub i32 %120, %102
%122 = add i32 %83, -3
%123 = sub i32 %122, %102
%124 = and i32 %121, 3
%125 = icmp eq i32 %124, 0
br i1 %125, label %135, label %126
126: ; preds = %115, %126
%127 = phi double [ %130, %126 ], [ %119, %115 ]
%128 = phi i64 [ %132, %126 ], [ %112, %115 ]
%129 = phi i32 [ %133, %126 ], [ 0, %115 ]
%130 = fadd double %85, %127
%131 = getelementptr inbounds double, double* %40, i64 %128
store double %130, double* %131
%132 = add nuw nsw i64 %128, 1
%133 = add i32 %129, 1
%134 = icmp eq i32 %133, %124
br i1 %134, label %135, label %126
135: ; preds = %126, %115
%136 = phi double [ %119, %115 ], [ %130, %126 ]
%137 = phi i64 [ %112, %115 ], [ %132, %126 ]
%138 = icmp ult i32 %123, 3
br i1 %138, label %139, label %142
139: ; preds = %135, %142, %108, %100
br i1 %86, label %140, label %170
140: ; preds = %139
%141 = zext i32 %83 to i64
br label %159
142: ; preds = %135, %142
%143 = phi double [ %154, %142 ], [ %136, %135 ]
%144 = phi i64 [ %156, %142 ], [ %137, %135 ]
%145 = fadd double %85, %143
%146 = getelementptr inbounds double, double* %40, i64 %144
store double %145, double* %146
%147 = add nuw nsw i64 %144, 1
%148 = fadd double %85, %145
%149 = getelementptr inbounds double, double* %40, i64 %147
store double %148, double* %149
%150 = add nuw nsw i64 %144, 2
%151 = fadd double %85, %148
%152 = getelementptr inbounds double, double* %40, i64 %150
store double %151, double* %152
%153 = add nuw nsw i64 %144, 3
%154 = fadd double %85, %151
%155 = getelementptr inbounds double, double* %40, i64 %153
store double %154, double* %155
%156 = add nuw nsw i64 %144, 4
%157 = trunc i64 %156 to i32
%158 = icmp sgt i32 %83, %157
br i1 %158, label %142, label %139
159: ; preds = %140, %164
%160 = phi i64 [ 0, %140 ], [ %166, %164 ]
%161 = getelementptr inbounds double, double* %40, i64 %160
%162 = load double, double* %161
%163 = fcmp olt double %162, 1.000000e+00
br i1 %163, label %164, label %168
164: ; preds = %159
%165 = getelementptr inbounds double, double* %42, i64 %160
store double 0.000000e+00, double* %165
%166 = add nuw nsw i64 %160, 1
%167 = icmp eq i64 %166, %141
br i1 %167, label %170, label %159
168: ; preds = %159
%169 = trunc i64 %160 to i32
br label %170
170: ; preds = %164, %168, %139
%171 = phi i32 [ 0, %139 ], [ %169, %168 ], [ %83, %164 ]
br i1 %70, label %172, label %267
172: ; preds = %170
br i1 %78, label %173, label %241
173: ; preds = %172
%174 = icmp slt i32 %171, %83
br i1 %174, label %175, label %267
175: ; preds = %173
%176 = zext i32 %171 to i64
%177 = xor i32 %171, -1
%178 = add i32 %83, %177
%179 = zext i32 %178 to i64
%180 = add nuw nsw i64 %179, 1
%181 = icmp eq i32 %178, 0
br i1 %181, label %229, label %182
182: ; preds = %175
%183 = and i64 %180, 8589934590
%184 = add nuw nsw i64 %183, %176
%185 = add nsw i64 %183, -2
%186 = lshr exact i64 %185, 1
%187 = add nuw i64 %186, 1
%188 = and i64 %187, 1
%189 = icmp eq i64 %185, 0
br i1 %189, label %215, label %190
190: ; preds = %182
%191 = and i64 %187, -2
br label %192
192: ; preds = %192, %190
%193 = phi i64 [ 0, %190 ], [ %212, %192 ]
%194 = phi i64 [ 0, %190 ], [ %213, %192 ]
%195 = add i64 %193, %176
%196 = getelementptr inbounds double, double* %40, i64 %195
%197 = bitcast double* %196 to <2 x double>*
%198 = load <2 x double>, <2 x double>* %197
%199 = fdiv <2 x double> <double 1.000000e+00, double 1.000000e+00>, %198
%200 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %199
%201 = getelementptr inbounds double, double* %42, i64 %195
%202 = bitcast double* %201 to <2 x double>*
store <2 x double> %200, <2 x double>* %202
%203 = or i64 %193, 2
%204 = add i64 %203, %176
%205 = getelementptr inbounds double, double* %40, i64 %204
%206 = bitcast double* %205 to <2 x double>*
%207 = load <2 x double>, <2 x double>* %206
%208 = fdiv <2 x double> <double 1.000000e+00, double 1.000000e+00>, %207
%209 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %208
%210 = getelementptr inbounds double, double* %42, i64 %204
%211 = bitcast double* %210 to <2 x double>*
store <2 x double> %209, <2 x double>* %211
%212 = add nuw i64 %193, 4
%213 = add nuw i64 %194, 2
%214 = icmp eq i64 %213, %191
br i1 %214, label %215, label %192
215: ; preds = %192, %182
%216 = phi i64 [ 0, %182 ], [ %212, %192 ]
%217 = icmp eq i64 %188, 0
br i1 %217, label %227, label %218
218: ; preds = %215
%219 = add i64 %216, %176
%220 = getelementptr inbounds double, double* %40, i64 %219
%221 = bitcast double* %220 to <2 x double>*
%222 = load <2 x double>, <2 x double>* %221
%223 = fdiv <2 x double> <double 1.000000e+00, double 1.000000e+00>, %222
%224 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %223
%225 = getelementptr inbounds double, double* %42, i64 %219
%226 = bitcast double* %225 to <2 x double>*
store <2 x double> %224, <2 x double>* %226
br label %227
227: ; preds = %215, %218
%228 = icmp eq i64 %180, %183
br i1 %228, label %262, label %229
229: ; preds = %175, %227
%230 = phi i64 [ %176, %175 ], [ %184, %227 ]
br label %231
231: ; preds = %229, %231
%232 = phi i64 [ %238, %231 ], [ %230, %229 ]
%233 = getelementptr inbounds double, double* %40, i64 %232
%234 = load double, double* %233
%235 = fdiv double 1.000000e+00, %234
%236 = fsub double 1.000000e+00, %235
%237 = getelementptr inbounds double, double* %42, i64 %232
store double %236, double* %237
%238 = add nuw nsw i64 %232, 1
%239 = trunc i64 %238 to i32
%240 = icmp sgt i32 %83, %239
br i1 %240, label %231, label %262
241: ; preds = %172
%242 = icmp eq i8 %77, 82
%243 = icmp slt i32 %171, %83
%244 = select i1 %242, i1 %243, i1 false
br i1 %244, label %245, label %267
245: ; preds = %241
%246 = fadd double %76, -1.000000e+00
%247 = zext i32 %171 to i64
br label %248
248: ; preds = %245, %253
%249 = phi i64 [ %247, %245 ], [ %259, %253 ]
%250 = getelementptr inbounds double, double* %40, i64 %249
%251 = load double, double* %250
%252 = fcmp olt double %251, %76
br i1 %252, label %253, label %265
253: ; preds = %248
%254 = fsub double %76, %251
%255 = fdiv double %254, %251
%256 = fdiv double %255, %246
%257 = fsub double 1.000000e+00, %256
%258 = getelementptr inbounds double, double* %42, i64 %249
store double %257, double* %258
%259 = add nuw nsw i64 %249, 1
%260 = trunc i64 %259 to i32
%261 = icmp sgt i32 %83, %260
br i1 %261, label %248, label %337
262: ; preds = %231, %227
%263 = phi i64 [ %184, %227 ], [ %238, %231 ]
%264 = trunc i64 %263 to i32
br label %267
265: ; preds = %248
%266 = trunc i64 %249 to i32
br label %267
267: ; preds = %265, %262, %173, %170, %241
%268 = phi i32 [ %171, %241 ], [ %171, %170 ], [ %171, %173 ], [ %264, %262 ], [ %266, %265 ]
%269 = icmp slt i32 %268, %83
br i1 %269, label %270, label %337
270: ; preds = %267
%271 = sext i32 %268 to i64
%272 = sext i32 %83 to i64
%273 = sub nsw i64 %272, %271
%274 = icmp ult i64 %273, 4
br i1 %274, label %330, label %275
275: ; preds = %270
%276 = and i64 %273, -4
%277 = add nsw i64 %276, %271
%278 = add nsw i64 %276, -4
%279 = lshr exact i64 %278, 2
%280 = add nuw nsw i64 %279, 1
%281 = and i64 %280, 3
%282 = icmp ult i64 %278, 12
br i1 %282, label %314, label %283
283: ; preds = %275
%284 = and i64 %280, 9223372036854775804
br label %285
285: ; preds = %285, %283
%286 = phi i64 [ 0, %283 ], [ %311, %285 ]
%287 = phi i64 [ 0, %283 ], [ %312, %285 ]
%288 = add i64 %286, %271
%289 = getelementptr inbounds double, double* %42, i64 %288
%290 = bitcast double* %289 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %290
%291 = getelementptr inbounds double, double* %289, i64 2
%292 = bitcast double* %291 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %292
%293 = or i64 %286, 4
%294 = add i64 %293, %271
%295 = getelementptr inbounds double, double* %42, i64 %294
%296 = bitcast double* %295 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %296
%297 = getelementptr inbounds double, double* %295, i64 2
%298 = bitcast double* %297 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %298
%299 = or i64 %286, 8
%300 = add i64 %299, %271
%301 = getelementptr inbounds double, double* %42, i64 %300
%302 = bitcast double* %301 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %302
%303 = getelementptr inbounds double, double* %301, i64 2
%304 = bitcast double* %303 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %304
%305 = or i64 %286, 12
%306 = add i64 %305, %271
%307 = getelementptr inbounds double, double* %42, i64 %306
%308 = bitcast double* %307 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %308
%309 = getelementptr inbounds double, double* %307, i64 2
%310 = bitcast double* %309 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %310
%311 = add nuw i64 %286, 16
%312 = add nuw i64 %287, 4
%313 = icmp eq i64 %312, %284
br i1 %313, label %314, label %285
314: ; preds = %285, %275
%315 = phi i64 [ 0, %275 ], [ %311, %285 ]
%316 = icmp eq i64 %281, 0
br i1 %316, label %328, label %317
317: ; preds = %314, %317
%318 = phi i64 [ %325, %317 ], [ %315, %314 ]
%319 = phi i64 [ %326, %317 ], [ 0, %314 ]
%320 = add i64 %318, %271
%321 = getelementptr inbounds double, double* %42, i64 %320
%322 = bitcast double* %321 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %322
%323 = getelementptr inbounds double, double* %321, i64 2
%324 = bitcast double* %323 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %324
%325 = add nuw i64 %318, 4
%326 = add i64 %319, 1
%327 = icmp eq i64 %326, %281
br i1 %327, label %328, label %317
328: ; preds = %317, %314
%329 = icmp eq i64 %273, %276
br i1 %329, label %337, label %330
330: ; preds = %270, %328
%331 = phi i64 [ %271, %270 ], [ %277, %328 ]
br label %332
332: ; preds = %330, %332
%333 = phi i64 [ %335, %332 ], [ %331, %330 ]
%334 = getelementptr inbounds double, double* %42, i64 %333
store double 1.000000e+00, double* %334
%335 = add nsw i64 %333, 1
%336 = icmp eq i64 %335, %272
br i1 %336, label %337, label %332
337: ; preds = %253, %332, %328, %267
br label %338
338: ; preds = %337, %368
%339 = phi double [ %369, %368 ], [ %36, %337 ]
%340 = phi double [ %355, %368 ], [ 0.000000e+00, %337 ]
%341 = call double @sqrt(double %34)
%342 = fmul double %35, %341
br i1 %71, label %343, label %345
343: ; preds = %338
%344 = fcmp olt double %339, %342
br i1 %344, label %347, label %370
345: ; preds = %338
%346 = fcmp ogt double %339, %342
br i1 %346, label %347, label %370
347: ; preds = %343, %345
%348 = fcmp ult double %340, 0.000000e+00
br i1 %348, label %354, label %349
349: ; preds = %347
%350 = load i32, i32* %2
%351 = load double, double* %1
%352 = call i32 (double*, double*, double*, i32, double, double, double, ...) bitcast (i32 (...)* @rdfsteadystate to i32 (double*, double*, double*, i32, double, double, double, ...)*)(double* %40, double* %42, double* %44, i32 %350, double %339, double %76, double %351)
%353 = sitofp i32 %352 to double
br label %354
354: ; preds = %349, %347
%355 = phi double [ %353, %349 ], [ %340, %347 ]
%356 = load i8, i8* %3
%357 = icmp eq i8 %356, 105
br i1 %357, label %358, label %360
358: ; preds = %354
%359 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([9 x i8], [9 x i8]* @.str.20, i64 0, i64 0), double %339, double %355)
br label %368
360: ; preds = %354
%361 = and i8 %356, -17
%362 = icmp eq i8 %361, 98
br i1 %362, label %363, label %365
363: ; preds = %360
%364 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([13 x i8], [13 x i8]* @.str.21, i64 0, i64 0), double %76, double %339, double %355)
br label %368
365: ; preds = %360
%366 = call double @llvm.fabs.f64(double %355)
%367 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([5 x i8], [5 x i8]* @.str.22, i64 0, i64 0), double %366)
br label %368
368: ; preds = %358, %365, %363
%369 = fmul double %34, %339
br label %338
370: ; preds = %343, %345
%371 = call i32 @putchar(i32 10)
%372 = load i8, i8* %3
%373 = and i8 %372, -33
switch i8 %373, label %379 [
i8 73, label %382
i8 82, label %374
]
374: ; preds = %370
%375 = fmul double %76, 0x3FF38ADD9E58AC8D
%376 = fcmp ogt double %375, 0x403632ACEF67FE3D
br i1 %376, label %382, label %377
377: ; preds = %374, %379
%378 = phi double [ %375, %374 ], [ %380, %379 ]
br label %74
379: ; preds = %370
%380 = fadd double %76, 1.000000e-01
%381 = fcmp ogt double %380, 1.050000e+00
br i1 %381, label %382, label %377
382: ; preds = %379, %374, %370
br label %383
383: ; preds = %382, %50, %27, %20
ret void
}