int main(int argc, char *argv[]) {
  int av, a, vmax, N, n, num, den, k, kq, kq2, t, v, s, i;
  double sum;
  if (argc < 2 || (n = atoi(argv[1])) <= 0) {
    printf("This program computes the n'th decimal digit of \\pi\n"
           "usage: pi n , where n is the digit you want\n");
    exit(1);
  }
  N = (int)((n + 20) * log(10) / log(2));
  sum = 0;
  for (a = 3; a <= (2 * N); a = next_prime(a)) {
    vmax = (int)(log(2 * N) / log(a));
    av = 1;
    for (i = 0; i < vmax; i++)
      av = av * a;
    s = 0;
    num = 1;
    den = 1;
    v = 0;
    kq = 1;
    kq2 = 1;
    for (k = 1; k <= N; k++) {
      t = k;
      if (kq >= a) {
        do {
          t = t / a;
          v--;
        } while ((t % a) == 0);
        kq = 0;
      }
      kq++;
      num = fmod((double)num * (double)t, av);
      t = (2 * k - 1);
      if (kq2 >= a) {
        if (kq2 == a) {
          do {
            t = t / a;
            v++;
          } while ((t % a) == 0);
        }
        kq2 -= a;
      }
      den = fmod((double)den * (double)t, av);
      kq2 += 2;
      if (v > 0) {
        t = inv_mod(den, av);
        t = fmod((double)t * (double)num, av);
        t = fmod((double)t * (double)k, av);
        for (i = v; i < vmax; i++)
          t = fmod((double)t * (double)a, av);
        s += t;
        if (s >= av)
          s -= av;
      }
    }
    t = pow_mod(10, n - 1, av);
    s = fmod((double)s * (double)t, av);
    sum = fmod(sum + (double)s / (double)av, 1.0);
  }
  printf("Decimal digits of pi at position %d: %09d\n", n, (int)(sum * 1e9));
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
movl   $0x0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
cmpl   $0x2,-0x8(%rbp)
jl     11bb <main+0x3b>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rdi
mov    $0x0,%al
call   1090 <atoi@plt>
mov    %eax,-0x24(%rbp)
cmp    $0x0,%eax
jg     11d3 <main+0x53>
lea    0xe5e(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    $0x1,%edi
call   10a0 <exit@plt>
mov    -0x24(%rbp),%eax
add    $0x14,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xe1e(%rip),%xmm0
call   1080 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xdf8(%rip),%xmm0
call   1080 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x3,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x20(%rbp),%ecx
shl    $0x1,%ecx
cmp    %ecx,%eax
jg     150e <main+0x38e>
mov    -0x20(%rbp),%eax
shl    $0x1,%eax
cvtsi2sd %eax,%xmm0
call   1080 <log@plt>
movsd  %xmm0,-0x68(%rbp)
cvtsi2sdl -0x18(%rbp),%xmm0
call   1080 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x1c(%rbp)
movl   $0x1,-0x14(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12a3 <main+0x123>
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    127f <main+0xff>
movl   $0x0,-0x44(%rbp)
movl   $0x1,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
movl   $0x0,-0x40(%rbp)
movl   $0x1,-0x34(%rbp)
movl   $0x1,-0x38(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     149c <main+0x31c>
mov    -0x30(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x18(%rbp),%eax
jl     1321 <main+0x1a1>
jmp    12f7 <main+0x177>
mov    -0x3c(%rbp),%eax
cltd
idivl  -0x18(%rbp)
mov    %eax,-0x3c(%rbp)
mov    -0x40(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
cltd
idivl  -0x18(%rbp)
cmp    $0x0,%edx
je     12f7 <main+0x177>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
cvtsi2sdl -0x28(%rbp),%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
call   1040 <fmod@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x28(%rbp)
mov    -0x30(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x18(%rbp),%eax
jl     13a5 <main+0x225>
mov    -0x38(%rbp),%eax
cmp    -0x18(%rbp),%eax
jne    139a <main+0x21a>
jmp    1372 <main+0x1f2>
mov    -0x3c(%rbp),%eax
cltd
idivl  -0x18(%rbp)
mov    %eax,-0x3c(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
cltd
idivl  -0x18(%rbp)
cmp    $0x0,%edx
je     1372 <main+0x1f2>
jmp    139a <main+0x21a>
mov    -0x18(%rbp),%ecx
mov    -0x38(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x38(%rbp)
cvtsi2sdl -0x2c(%rbp),%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
call   1040 <fmod@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x40(%rbp)
jle    1489 <main+0x309>
mov    -0x2c(%rbp),%edi
mov    -0x14(%rbp),%esi
mov    $0x0,%al
call   1070 <inv_mod@plt>
mov    %eax,-0x3c(%rbp)
cvtsi2sdl -0x3c(%rbp),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
call   1040 <fmod@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x3c(%rbp)
cvtsi2sdl -0x3c(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
call   1040 <fmod@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1464 <main+0x2e4>
cvtsi2sdl -0x3c(%rbp),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
call   1040 <fmod@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    142b <main+0x2ab>
mov    -0x3c(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     1484 <main+0x304>
mov    -0x14(%rbp),%ecx
mov    -0x44(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x44(%rbp)
jmp    1489 <main+0x309>
jmp    148e <main+0x30e>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12d4 <main+0x154>
mov    -0x24(%rbp),%esi
sub    $0x1,%esi
mov    -0x14(%rbp),%edx
mov    $0xa,%edi
mov    $0x0,%al
call   1060 <pow_mod@plt>
mov    %eax,-0x3c(%rbp)
cvtsi2sdl -0x44(%rbp),%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
call   1040 <fmod@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x44(%rbp)
movsd  -0x50(%rbp),%xmm0
cvtsi2sdl -0x44(%rbp),%xmm1
cvtsi2sdl -0x14(%rbp),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  0xb26(%rip),%xmm1
call   1040 <fmod@plt>
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%edi
mov    $0x0,%al
call   1050 <next_prime@plt>
mov    %eax,-0x18(%rbp)
jmp    122f <main+0xaf>
mov    -0x24(%rbp),%esi
movsd  0xaf7(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
cvttsd2si %xmm0,%edx
lea    0xb58(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
xor    %eax,%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
cmp    $0x2,%edi
jl     16b1 <main+0x521>
mov    0x8(%rsi),%rdi
xor    %ebx,%ebx
xor    %eax,%eax
call   10a0 <atoi@plt>
test   %eax,%eax
jle    16b1 <main+0x521>
mov    %eax,%esi
lea    0x14(%rsi),%eax
cvtsi2sd %eax,%xmm0
mulsd  0xe33(%rip),%xmm0
divsd  0xe33(%rip),%xmm0
cvttsd2si %xmm0,%r13d
cmp    $0x2,%r13d
jl     1690 <main+0x500>
lea    0x0(%r13,%r13,1),%eax
mov    %eax,0x38(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x50(%rsp)
mov    %rsi,0x48(%rsp)
mov    %esi,%eax
add    $0xffffffff,%eax
mov    %eax,0x34(%rsp)
mov    $0x3,%ebx
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x28(%rsp)
mov    %r13,0x58(%rsp)
jmp    1293 <main+0x103>
xchg   %ax,%ax
mov    $0xa,%edi
mov    0x34(%rsp),%esi
mov    0x1c(%rsp),%edx
xor    %eax,%eax
call   1070 <pow_mod@plt>
xorps  %xmm1,%xmm1
cvtsi2sdl 0x18(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0x8(%rsp),%xmm1
call   1050 <fmod@plt>
cvttpd2dq %xmm0,%xmm0
cvtdq2pd %xmm0,%xmm1
divsd  0x8(%rsp),%xmm1
movsd  0x28(%rsp),%xmm0
addsd  %xmm1,%xmm0
movq   0xd9d(%rip),%xmm1
call   1050 <fmod@plt>
movsd  %xmm0,0x28(%rsp)
mov    %ebx,%edi
xor    %eax,%eax
call   1060 <next_prime@plt>
mov    %eax,%ebx
cmp    0x38(%rsp),%eax
jg     1679 <main+0x4e9>
movsd  0x50(%rsp),%xmm0
call   1090 <log@plt>
movsd  %xmm0,0x8(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm0,0x60(%rsp)
call   1090 <log@plt>
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%eax
mov    $0x1,%r14d
mov    $0x1,%edx
mov    %eax,0x14(%rsp)
test   %eax,%eax
jle    1488 <main+0x2f8>
mov    0x14(%rsp),%eax
cmp    $0x8,%eax
jae    12f0 <main+0x160>
mov    $0x1,%edx
xor    %eax,%eax
jmp    146d <main+0x2dd>
and    $0xfffffff8,%eax
movd   %ebx,%xmm0
pshufd $0x0,%xmm0,%xmm0
lea    -0x8(%rax),%esi
mov    %esi,%edx
shr    $0x3,%edx
add    $0x1,%edx
mov    %edx,%ecx
and    $0x7,%ecx
pshufd $0xf5,%xmm0,%xmm3
movdqa 0xd07(%rip),%xmm1
movdqa %xmm1,%xmm2
cmp    $0x38,%esi
jb     13ef <main+0x25f>
and    $0xfffffff8,%edx
movdqa 0xcef(%rip),%xmm1
movdqa %xmm1,%xmm2
cs nopw 0x0(%rax,%rax,1)
nop
pshufd $0xf5,%xmm2,%xmm5
pmuludq %xmm3,%xmm5
pmuludq %xmm0,%xmm2
pshufd $0xf5,%xmm1,%xmm4
pmuludq %xmm3,%xmm4
pmuludq %xmm0,%xmm1
pmuludq %xmm0,%xmm2
pmuludq %xmm3,%xmm5
pmuludq %xmm0,%xmm1
pmuludq %xmm3,%xmm4
pmuludq %xmm3,%xmm5
pmuludq %xmm0,%xmm2
pmuludq %xmm3,%xmm4
pmuludq %xmm0,%xmm1
pmuludq %xmm0,%xmm2
pmuludq %xmm3,%xmm5
pmuludq %xmm0,%xmm1
pmuludq %xmm3,%xmm4
pmuludq %xmm3,%xmm5
pmuludq %xmm0,%xmm2
pmuludq %xmm3,%xmm4
pmuludq %xmm0,%xmm1
pmuludq %xmm0,%xmm2
pmuludq %xmm3,%xmm5
pmuludq %xmm0,%xmm1
pmuludq %xmm3,%xmm4
pmuludq %xmm3,%xmm5
pmuludq %xmm0,%xmm2
pmuludq %xmm3,%xmm4
pmuludq %xmm0,%xmm1
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pmuludq %xmm3,%xmm5
pshufd $0xe8,%xmm5,%xmm5
punpckldq %xmm5,%xmm2
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm1
pmuludq %xmm3,%xmm4
pshufd $0xe8,%xmm4,%xmm4
punpckldq %xmm4,%xmm1
add    $0xfffffff8,%edx
jne    1340 <main+0x1b0>
test   %ecx,%ecx
je     143b <main+0x2ab>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
pshufd $0xf5,%xmm2,%xmm4
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pmuludq %xmm3,%xmm4
pshufd $0xe8,%xmm4,%xmm4
punpckldq %xmm4,%xmm2
pshufd $0xf5,%xmm1,%xmm4
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm1
pmuludq %xmm3,%xmm4
pshufd $0xe8,%xmm4,%xmm4
punpckldq %xmm4,%xmm1
add    $0xffffffff,%ecx
jne    1400 <main+0x270>
pshufd $0xf5,%xmm2,%xmm0
pshufd $0xf5,%xmm1,%xmm3
pmuludq %xmm0,%xmm3
pmuludq %xmm2,%xmm1
pshufd $0xee,%xmm1,%xmm0
pmuludq %xmm1,%xmm0
pshufd $0xaa,%xmm3,%xmm1
pmuludq %xmm3,%xmm1
pmuludq %xmm0,%xmm1
movd   %xmm1,%edx
cmp    0x14(%rsp),%eax
je     1488 <main+0x2f8>
mov    0x14(%rsp),%ecx
sub    %eax,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
imul   %ebx,%edx
add    $0xffffffff,%ecx
jne    1480 <main+0x2f0>
mov    %edx,0x1c(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
movl   $0x0,0x18(%rsp)
xor    %r15d,%r15d
mov    $0x1,%ecx
mov    $0x1,%ebp
movl   $0x1,0x24(%rsp)
mov    $0x1,%esi
movsd  %xmm1,0x8(%rsp)
jmp    14ea <main+0x35a>
nopl   (%rax)
movq   0x8(%rsp),%xmm1
mov    0x40(%rsp),%r12
mov    0x20(%rsp),%ecx
add    $0x1,%ecx
add    $0x2,%r14d
lea    0x1(%r12),%eax
cmp    %r13d,%r12d
mov    %eax,%ebp
mov    0x3c(%rsp),%esi
je     1220 <main+0x90>
cmp    %ebx,%ecx
jge    1500 <main+0x370>
mov    %ecx,0x20(%rsp)
mov    %ebp,%ecx
jmp    152a <main+0x39a>
cs nopw 0x0(%rax,%rax,1)
mov    %ebp,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %ecx,%eax
cltd
idiv   %ebx
add    $0xffffffff,%r15d
mov    %eax,%ecx
cltd
idiv   %ebx
test   %edx,%edx
je     1510 <main+0x380>
movl   $0x0,0x20(%rsp)
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
mulsd  %xmm2,%xmm0
call   1050 <fmod@plt>
lea    -0x1(%rbp,%rbp,1),%ecx
mov    %r14d,%eax
sub    %ebx,%r14d
mov    %rbp,0x40(%rsp)
jge    1560 <main+0x3d0>
mov    %eax,%r14d
jmp    1590 <main+0x400>
nopw   0x0(%rax,%rax,1)
jne    1590 <main+0x400>
movq   0x8(%rsp),%xmm1
nopl   0x0(%rax,%rax,1)
mov    %ecx,%eax
cltd
idiv   %ebx
add    $0x1,%r15d
mov    %eax,%ecx
cltd
idiv   %ebx
test   %edx,%edx
je     1570 <main+0x3e0>
jmp    1596 <main+0x406>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movq   0x8(%rsp),%xmm1
cvttsd2si %xmm0,%ebp
xorps  %xmm2,%xmm2
cvtsi2sdl 0x24(%rsp),%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
mulsd  %xmm2,%xmm0
call   1050 <fmod@plt>
cvttsd2si %xmm0,%edi
test   %r15d,%r15d
mov    %edi,0x24(%rsp)
mov    %ebp,0x3c(%rsp)
jle    14c0 <main+0x330>
mov    0x1c(%rsp),%esi
xor    %eax,%eax
call   1080 <inv_mod@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
mulsd  %xmm1,%xmm0
movsd  0x8(%rsp),%xmm1
call   1050 <fmod@plt>
cvttpd2dq %xmm0,%xmm0
cvtdq2pd %xmm0,%xmm1
mov    0x40(%rsp),%r12
xorps  %xmm0,%xmm0
cvtsi2sd %r12d,%xmm0
mulsd  %xmm1,%xmm0
movq   0x8(%rsp),%xmm1
call   1050 <fmod@plt>
cvttsd2si %xmm0,%eax
mov    0x14(%rsp),%ebp
cmp    %ebp,%r15d
jge    1651 <main+0x4c1>
sub    %r15d,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0x60(%rsp),%xmm0
movq   0x8(%rsp),%xmm1
call   1050 <fmod@plt>
cvttsd2si %xmm0,%eax
add    $0xffffffff,%ebp
jne    1630 <main+0x4a0>
add    0x18(%rsp),%eax
mov    0x1c(%rsp),%ecx
cmp    %ecx,%eax
mov    $0x0,%edx
cmovl  %edx,%ecx
sub    %ecx,%eax
mov    %eax,0x18(%rsp)
mov    0x58(%rsp),%r13
movq   0x8(%rsp),%xmm1
jmp    14cb <main+0x33b>
movsd  0x28(%rsp),%xmm0
mulsd  0x991(%rip),%xmm0
cvttsd2si %xmm0,%ebx
mov    0x48(%rsp),%rsi
lea    0x999(%rip),%rdi
mov    %ebx,%edx
xor    %eax,%eax
call   1040 <printf@plt>
xor    %eax,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0x9a3(%rip),%rdi
call   1030 <puts@plt>
mov    $0x1,%edi
call   10b0 <exit@plt>

LLVM IR (-O0):
@.str = [97 x i8] c"This program computes the n'th decimal digit of \\pi\0Ausage: pi n , where n is the digit you want\0A\00"
@.str.1 = [43 x i8] c"Decimal digits of pi at position %d: %09d\0A\00"
define i32 @main(i32 %0, i8** %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i8**
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca double
store i32 0, i32* %3
store i32 %0, i32* %4
store i8** %1, i8*** %5
%21 = load i32, i32* %4
%22 = icmp slt i32 %21, 2
br i1 %22, label %29, label %23
23: ; preds = %2
%24 = load i8**, i8*** %5
%25 = getelementptr inbounds i8*, i8** %24, i64 1
%26 = load i8*, i8** %25
%27 = call i32 (i8*, ...) bitcast (i32 (...)* @atoi to i32 (i8*, ...)*)(i8* %26)
store i32 %27, i32* %10
%28 = icmp sle i32 %27, 0
br i1 %28, label %29, label %31
29: ; preds = %23, %2
%30 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([97 x i8], [97 x i8]* @.str, i64 0, i64 0))
unreachable
31: ; preds = %23
%32 = load i32, i32* %10
%33 = add nsw i32 %32, 20
%34 = sitofp i32 %33 to double
%35 = call double @log(double 1.000000e+01)
%36 = fmul double %34, %35
%37 = call double @log(double 2.000000e+00)
%38 = fdiv double %36, %37
%39 = fptosi double %38 to i32
store i32 %39, i32* %9
store double 0.000000e+00, double* %20
store i32 3, i32* %7
br label %40
40: ; preds = %220, %31
%41 = load i32, i32* %7
%42 = load i32, i32* %9
%43 = mul nsw i32 2, %42
%44 = icmp sle i32 %41, %43
br i1 %44, label %45, label %223
45: ; preds = %40
%46 = load i32, i32* %9
%47 = mul nsw i32 2, %46
%48 = sitofp i32 %47 to double
%49 = call double @log(double %48)
%50 = load i32, i32* %7
%51 = sitofp i32 %50 to double
%52 = call double @log(double %51)
%53 = fdiv double %49, %52
%54 = fptosi double %53 to i32
store i32 %54, i32* %8
store i32 1, i32* %6
store i32 0, i32* %19
br label %55
55: ; preds = %63, %45
%56 = load i32, i32* %19
%57 = load i32, i32* %8
%58 = icmp slt i32 %56, %57
br i1 %58, label %59, label %66
59: ; preds = %55
%60 = load i32, i32* %6
%61 = load i32, i32* %7
%62 = mul nsw i32 %60, %61
store i32 %62, i32* %6
br label %63
63: ; preds = %59
%64 = load i32, i32* %19
%65 = add nsw i32 %64, 1
store i32 %65, i32* %19
br label %55
66: ; preds = %55
store i32 0, i32* %18
store i32 1, i32* %11
store i32 1, i32* %12
store i32 0, i32* %17
store i32 1, i32* %14
store i32 1, i32* %15
store i32 1, i32* %13
br label %67
67: ; preds = %195, %66
%68 = load i32, i32* %13
%69 = load i32, i32* %9
%70 = icmp sle i32 %68, %69
br i1 %70, label %71, label %198
71: ; preds = %67
%72 = load i32, i32* %13
store i32 %72, i32* %16
%73 = load i32, i32* %14
%74 = load i32, i32* %7
%75 = icmp sge i32 %73, %74
br i1 %75, label %76, label %89
76: ; preds = %71
br label %77
77: ; preds = %83, %76
%78 = load i32, i32* %16
%79 = load i32, i32* %7
%80 = sdiv i32 %78, %79
store i32 %80, i32* %16
%81 = load i32, i32* %17
%82 = add nsw i32 %81, -1
store i32 %82, i32* %17
br label %83
83: ; preds = %77
%84 = load i32, i32* %16
%85 = load i32, i32* %7
%86 = srem i32 %84, %85
%87 = icmp eq i32 %86, 0
br i1 %87, label %77, label %88
88: ; preds = %83
store i32 0, i32* %14
br label %89
89: ; preds = %88, %71
%90 = load i32, i32* %14
%91 = add nsw i32 %90, 1
store i32 %91, i32* %14
%92 = load i32, i32* %11
%93 = sitofp i32 %92 to double
%94 = load i32, i32* %16
%95 = sitofp i32 %94 to double
%96 = fmul double %93, %95
%97 = load i32, i32* %6
%98 = sitofp i32 %97 to double
%99 = call double @fmod(double %96, double %98)
%100 = fptosi double %99 to i32
store i32 %100, i32* %11
%101 = load i32, i32* %13
%102 = mul nsw i32 2, %101
%103 = sub nsw i32 %102, 1
store i32 %103, i32* %16
%104 = load i32, i32* %15
%105 = load i32, i32* %7
%106 = icmp sge i32 %104, %105
br i1 %106, label %107, label %128
107: ; preds = %89
%108 = load i32, i32* %15
%109 = load i32, i32* %7
%110 = icmp eq i32 %108, %109
br i1 %110, label %111, label %124
111: ; preds = %107
br label %112
112: ; preds = %118, %111
%113 = load i32, i32* %16
%114 = load i32, i32* %7
%115 = sdiv i32 %113, %114
store i32 %115, i32* %16
%116 = load i32, i32* %17
%117 = add nsw i32 %116, 1
store i32 %117, i32* %17
br label %118
118: ; preds = %112
%119 = load i32, i32* %16
%120 = load i32, i32* %7
%121 = srem i32 %119, %120
%122 = icmp eq i32 %121, 0
br i1 %122, label %112, label %123
123: ; preds = %118
br label %124
124: ; preds = %123, %107
%125 = load i32, i32* %7
%126 = load i32, i32* %15
%127 = sub nsw i32 %126, %125
store i32 %127, i32* %15
br label %128
128: ; preds = %124, %89
%129 = load i32, i32* %12
%130 = sitofp i32 %129 to double
%131 = load i32, i32* %16
%132 = sitofp i32 %131 to double
%133 = fmul double %130, %132
%134 = load i32, i32* %6
%135 = sitofp i32 %134 to double
%136 = call double @fmod(double %133, double %135)
%137 = fptosi double %136 to i32
store i32 %137, i32* %12
%138 = load i32, i32* %15
%139 = add nsw i32 %138, 2
store i32 %139, i32* %15
%140 = load i32, i32* %17
%141 = icmp sgt i32 %140, 0
br i1 %141, label %142, label %194
142: ; preds = %128
%143 = load i32, i32* %12
%144 = load i32, i32* %6
%145 = call i32 (i32, i32, ...) bitcast (i32 (...)* @inv_mod to i32 (i32, i32, ...)*)(i32 %143, i32 %144)
store i32 %145, i32* %16
%146 = load i32, i32* %16
%147 = sitofp i32 %146 to double
%148 = load i32, i32* %11
%149 = sitofp i32 %148 to double
%150 = fmul double %147, %149
%151 = load i32, i32* %6
%152 = sitofp i32 %151 to double
%153 = call double @fmod(double %150, double %152)
%154 = fptosi double %153 to i32
store i32 %154, i32* %16
%155 = load i32, i32* %16
%156 = sitofp i32 %155 to double
%157 = load i32, i32* %13
%158 = sitofp i32 %157 to double
%159 = fmul double %156, %158
%160 = load i32, i32* %6
%161 = sitofp i32 %160 to double
%162 = call double @fmod(double %159, double %161)
%163 = fptosi double %162 to i32
store i32 %163, i32* %16
%164 = load i32, i32* %17
store i32 %164, i32* %19
br label %165
165: ; preds = %179, %142
%166 = load i32, i32* %19
%167 = load i32, i32* %8
%168 = icmp slt i32 %166, %167
br i1 %168, label %169, label %182
169: ; preds = %165
%170 = load i32, i32* %16
%171 = sitofp i32 %170 to double
%172 = load i32, i32* %7
%173 = sitofp i32 %172 to double
%174 = fmul double %171, %173
%175 = load i32, i32* %6
%176 = sitofp i32 %175 to double
%177 = call double @fmod(double %174, double %176)
%178 = fptosi double %177 to i32
store i32 %178, i32* %16
br label %179
179: ; preds = %169
%180 = load i32, i32* %19
%181 = add nsw i32 %180, 1
store i32 %181, i32* %19
br label %165
182: ; preds = %165
%183 = load i32, i32* %16
%184 = load i32, i32* %18
%185 = add nsw i32 %184, %183
store i32 %185, i32* %18
%186 = load i32, i32* %18
%187 = load i32, i32* %6
%188 = icmp sge i32 %186, %187
br i1 %188, label %189, label %193
189: ; preds = %182
%190 = load i32, i32* %6
%191 = load i32, i32* %18
%192 = sub nsw i32 %191, %190
store i32 %192, i32* %18
br label %193
193: ; preds = %189, %182
br label %194
194: ; preds = %193, %128
br label %195
195: ; preds = %194
%196 = load i32, i32* %13
%197 = add nsw i32 %196, 1
store i32 %197, i32* %13
br label %67
198: ; preds = %67
%199 = load i32, i32* %10
%200 = sub nsw i32 %199, 1
%201 = load i32, i32* %6
%202 = call i32 (i32, i32, i32, ...) bitcast (i32 (...)* @pow_mod to i32 (i32, i32, i32, ...)*)(i32 10, i32 %200, i32 %201)
store i32 %202, i32* %16
%203 = load i32, i32* %18
%204 = sitofp i32 %203 to double
%205 = load i32, i32* %16
%206 = sitofp i32 %205 to double
%207 = fmul double %204, %206
%208 = load i32, i32* %6
%209 = sitofp i32 %208 to double
%210 = call double @fmod(double %207, double %209)
%211 = fptosi double %210 to i32
store i32 %211, i32* %18
%212 = load double, double* %20
%213 = load i32, i32* %18
%214 = sitofp i32 %213 to double
%215 = load i32, i32* %6
%216 = sitofp i32 %215 to double
%217 = fdiv double %214, %216
%218 = fadd double %212, %217
%219 = call double @fmod(double %218, double 1.000000e+00)
store double %219, double* %20
br label %220
220: ; preds = %198
%221 = load i32, i32* %7
%222 = call i32 (i32, ...) bitcast (i32 (...)* @next_prime to i32 (i32, ...)*)(i32 %221)
store i32 %222, i32* %7
br label %40
223: ; preds = %40
%224 = load i32, i32* %10
%225 = load double, double* %20
%226 = fmul double %225, 1.000000e+09
%227 = fptosi double %226 to i32
%228 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([43 x i8], [43 x i8]* @.str.1, i64 0, i64 0), i32 %224, i32 %227)
ret i32 0
}

LLVM IR (-O2):
@.str.1 = [43 x i8] c"Decimal digits of pi at position %d: %09d\0A\00"
@str = [96 x i8] c"This program computes the n'th decimal digit of \\pi\0Ausage: pi n , where n is the digit you want\00"
define i32 @main(i32 %0, i8** nocapture %1) {
%3 = icmp slt i32 %0, 2
br i1 %3, label %9, label %4
4: ; preds = %2
%5 = getelementptr inbounds i8*, i8** %1, i64 1
%6 = load i8*, i8** %5
%7 = tail call i32 (i8*, ...) bitcast (i32 (...)* @atoi to i32 (i8*, ...)*)(i8* %6)
%8 = icmp slt i32 %7, 1
br i1 %8, label %9, label %11
9: ; preds = %4, %2
%10 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([96 x i8], [96 x i8]* @str, i64 0, i64 0))
tail call void @exit(i32 1)
unreachable
11: ; preds = %4
%12 = add nuw nsw i32 %7, 20
%13 = sitofp i32 %12 to double
%14 = fmul double %13, 0x40026BB1BBB55516
%15 = fdiv double %14, 0x3FE62E42FEFA39EF
%16 = fptosi double %15 to i32
%17 = shl nsw i32 %16, 1
%18 = icmp slt i32 %16, 2
br i1 %18, label %201, label %19
19: ; preds = %11
%20 = sitofp i32 %17 to double
%21 = add nsw i32 %7, -1
br label %22
22: ; preds = %19, %185
%23 = phi double [ 0.000000e+00, %19 ], [ %195, %185 ]
%24 = phi i32 [ 3, %19 ], [ %196, %185 ]
%25 = tail call double @log(double %20)
%26 = sitofp i32 %24 to double
%27 = tail call double @log(double %26)
%28 = fdiv double %25, %27
%29 = fptosi double %28 to i32
%30 = icmp sgt i32 %29, 0
br i1 %30, label %31, label %91
31: ; preds = %22
%32 = icmp ult i32 %29, 8
br i1 %32, label %88, label %33
33: ; preds = %31
%34 = and i32 %29, -8
%35 = insertelement <4 x i32> poison, i32 %24, i64 0
%36 = shufflevector <4 x i32> %35, <4 x i32> poison, <4 x i32> zeroinitializer
%37 = insertelement <4 x i32> poison, i32 %24, i64 0
%38 = shufflevector <4 x i32> %37, <4 x i32> poison, <4 x i32> zeroinitializer
%39 = add i32 %34, -8
%40 = lshr exact i32 %39, 3
%41 = add nuw nsw i32 %40, 1
%42 = and i32 %41, 7
%43 = icmp ult i32 %39, 56
br i1 %43, label %68, label %44
44: ; preds = %33
%45 = and i32 %41, 1073741816
br label %46
46: ; preds = %46, %44
%47 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %44 ], [ %64, %46 ]
%48 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %44 ], [ %65, %46 ]
%49 = phi i32 [ 0, %44 ], [ %66, %46 ]
%50 = mul <4 x i32> %47, %36
%51 = mul <4 x i32> %48, %38
%52 = mul <4 x i32> %50, %36
%53 = mul <4 x i32> %51, %38
%54 = mul <4 x i32> %52, %36
%55 = mul <4 x i32> %53, %38
%56 = mul <4 x i32> %54, %36
%57 = mul <4 x i32> %55, %38
%58 = mul <4 x i32> %56, %36
%59 = mul <4 x i32> %57, %38
%60 = mul <4 x i32> %58, %36
%61 = mul <4 x i32> %59, %38
%62 = mul <4 x i32> %60, %36
%63 = mul <4 x i32> %61, %38
%64 = mul <4 x i32> %62, %36
%65 = mul <4 x i32> %63, %38
%66 = add nuw i32 %49, 8
%67 = icmp eq i32 %66, %45
br i1 %67, label %68, label %46
68: ; preds = %46, %33
%69 = phi <4 x i32> [ undef, %33 ], [ %64, %46 ]
%70 = phi <4 x i32> [ undef, %33 ], [ %65, %46 ]
%71 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %33 ], [ %64, %46 ]
%72 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %33 ], [ %65, %46 ]
%73 = icmp eq i32 %42, 0
br i1 %73, label %82, label %74
74: ; preds = %68, %74
%75 = phi <4 x i32> [ %78, %74 ], [ %71, %68 ]
%76 = phi <4 x i32> [ %79, %74 ], [ %72, %68 ]
%77 = phi i32 [ %80, %74 ], [ 0, %68 ]
%78 = mul <4 x i32> %75, %36
%79 = mul <4 x i32> %76, %38
%80 = add i32 %77, 1
%81 = icmp eq i32 %80, %42
br i1 %81, label %82, label %74
82: ; preds = %74, %68
%83 = phi <4 x i32> [ %69, %68 ], [ %78, %74 ]
%84 = phi <4 x i32> [ %70, %68 ], [ %79, %74 ]
%85 = mul <4 x i32> %84, %83
%86 = call i32 @llvm.vector.reduce.mul.v4i32(<4 x i32> %85)
%87 = icmp eq i32 %34, %29
br i1 %87, label %91, label %88
88: ; preds = %31, %82
%89 = phi i32 [ 0, %31 ], [ %34, %82 ]
%90 = phi i32 [ 1, %31 ], [ %86, %82 ]
br label %94
91: ; preds = %94, %82, %22
%92 = phi i32 [ 1, %22 ], [ %86, %82 ], [ %97, %94 ]
%93 = sitofp i32 %92 to double
br label %100
94: ; preds = %88, %94
%95 = phi i32 [ %98, %94 ], [ %89, %88 ]
%96 = phi i32 [ %97, %94 ], [ %90, %88 ]
%97 = mul nsw i32 %96, %24
%98 = add nuw nsw i32 %95, 1
%99 = icmp eq i32 %98, %29
br i1 %99, label %91, label %94
100: ; preds = %91, %181
%101 = phi i32 [ 0, %91 ], [ %182, %181 ]
%102 = phi i32 [ 0, %91 ], [ %145, %181 ]
%103 = phi i32 [ 1, %91 ], [ %151, %181 ]
%104 = phi i32 [ 1, %91 ], [ %120, %181 ]
%105 = phi i32 [ 1, %91 ], [ %183, %181 ]
%106 = phi i32 [ 1, %91 ], [ %150, %181 ]
%107 = phi i32 [ 1, %91 ], [ %125, %181 ]
%108 = icmp slt i32 %104, %24
br i1 %108, label %116, label %109
109: ; preds = %100, %109
%110 = phi i32 [ %112, %109 ], [ %105, %100 ]
%111 = phi i32 [ %113, %109 ], [ %102, %100 ]
%112 = sdiv i32 %110, %24
%113 = add nsw i32 %111, -1
%114 = srem i32 %112, %24
%115 = icmp eq i32 %114, 0
br i1 %115, label %109, label %116
116: ; preds = %109, %100
%117 = phi i32 [ %104, %100 ], [ 0, %109 ]
%118 = phi i32 [ %105, %100 ], [ %112, %109 ]
%119 = phi i32 [ %102, %100 ], [ %113, %109 ]
%120 = add nsw i32 %117, 1
%121 = sitofp i32 %107 to double
%122 = sitofp i32 %118 to double
%123 = fmul double %121, %122
%124 = tail call double @fmod(double %123, double %93)
%125 = fptosi double %124 to i32
%126 = shl nuw nsw i32 %105, 1
%127 = add nsw i32 %126, -1
%128 = icmp slt i32 %103, %24
br i1 %128, label %142, label %129
129: ; preds = %116
%130 = icmp eq i32 %103, %24
br i1 %130, label %131, label %138
131: ; preds = %129, %131
%132 = phi i32 [ %134, %131 ], [ %127, %129 ]
%133 = phi i32 [ %135, %131 ], [ %119, %129 ]
%134 = sdiv i32 %132, %24
%135 = add nsw i32 %133, 1
%136 = srem i32 %134, %24
%137 = icmp eq i32 %136, 0
br i1 %137, label %131, label %138
138: ; preds = %131, %129
%139 = phi i32 [ %127, %129 ], [ %134, %131 ]
%140 = phi i32 [ %119, %129 ], [ %135, %131 ]
%141 = sub nsw i32 %103, %24
br label %142
142: ; preds = %138, %116
%143 = phi i32 [ %141, %138 ], [ %103, %116 ]
%144 = phi i32 [ %139, %138 ], [ %127, %116 ]
%145 = phi i32 [ %140, %138 ], [ %119, %116 ]
%146 = sitofp i32 %106 to double
%147 = sitofp i32 %144 to double
%148 = fmul double %146, %147
%149 = tail call double @fmod(double %148, double %93)
%150 = fptosi double %149 to i32
%151 = add nsw i32 %143, 2
%152 = icmp sgt i32 %145, 0
br i1 %152, label %153, label %181
153: ; preds = %142
%154 = tail call i32 (i32, i32, ...) bitcast (i32 (...)* @inv_mod to i32 (i32, i32, ...)*)(i32 %150, i32 %92)
%155 = sitofp i32 %154 to double
%156 = sitofp i32 %125 to double
%157 = fmul double %156, %155
%158 = tail call double @fmod(double %157, double %93)
%159 = fptosi double %158 to i32
%160 = sitofp i32 %159 to double
%161 = sitofp i32 %105 to double
%162 = fmul double %161, %160
%163 = tail call double @fmod(double %162, double %93)
%164 = fptosi double %163 to i32
%165 = icmp slt i32 %145, %29
br i1 %165, label %166, label %175
166: ; preds = %153, %166
%167 = phi i32 [ %173, %166 ], [ %164, %153 ]
%168 = phi i32 [ %172, %166 ], [ %145, %153 ]
%169 = sitofp i32 %167 to double
%170 = fmul double %26, %169
%171 = tail call double @fmod(double %170, double %93)
%172 = add nsw i32 %168, 1
%173 = fptosi double %171 to i32
%174 = icmp eq i32 %172, %29
br i1 %174, label %175, label %166
175: ; preds = %166, %153
%176 = phi i32 [ %164, %153 ], [ %173, %166 ]
%177 = add nsw i32 %176, %101
%178 = icmp slt i32 %177, %92
%179 = select i1 %178, i32 0, i32 %92
%180 = sub nsw i32 %177, %179
br label %181
181: ; preds = %175, %142
%182 = phi i32 [ %101, %142 ], [ %180, %175 ]
%183 = add nuw i32 %105, 1
%184 = icmp eq i32 %105, %16
br i1 %184, label %185, label %100
185: ; preds = %181
%186 = tail call i32 (i32, i32, i32, ...) bitcast (i32 (...)* @pow_mod to i32 (i32, i32, i32, ...)*)(i32 10, i32 %21, i32 %92)
%187 = sitofp i32 %182 to double
%188 = sitofp i32 %186 to double
%189 = fmul double %187, %188
%190 = tail call double @fmod(double %189, double %93)
%191 = fptosi double %190 to i32
%192 = sitofp i32 %191 to double
%193 = fdiv double %192, %93
%194 = fadd double %23, %193
%195 = tail call double @fmod(double %194, double 1.000000e+00)
%196 = tail call i32 (i32, ...) bitcast (i32 (...)* @next_prime to i32 (i32, ...)*)(i32 %24)
%197 = icmp sgt i32 %196, %17
br i1 %197, label %198, label %22
198: ; preds = %185
%199 = fmul double %195, 1.000000e+09
%200 = fptosi double %199 to i32
br label %201
201: ; preds = %198, %11
%202 = phi i32 [ 0, %11 ], [ %200, %198 ]
%203 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([43 x i8], [43 x i8]* @.str.1, i64 0, i64 0), i32 %7, i32 %202)
ret i32 0
}