void cradau(
    int n, void fcn(int *, double *, double *, double *, double *, int *),
    double x, double *y, double xend, double h, double rtol, double atol,
    void jac(int *, double *, double *, double *, int *, double *, double *),
    int ijac, int mljac, int mujac,
    void mas(int *n, double *am, int *lmas, double *rpar, int *ipar), int imas,
    int mlmas, int mumas,
    void solout(int *, double *, double *, double *, double *, int *, int *,
                double *, int *, int *),
    int iout, double *work, int *iwork, double *rpar, int *ipar, int *idid) {
  int N = n;
  double X = x;
  double XEND = xend;
  double H = h;
  int ITOL = 0;
  double RTOL = rtol;
  double ATOL = atol;
  int IJAC = ijac;
  int MLJAC = mljac;
  int MUJAC = mujac;
  int IMAS = imas;
  int MLMAS = mlmas;
  int MUMAS = mumas;
  int IOUT = iout;
  double *WORK;
  int LWORK;
  int *IWORK;
  int LIWORK = 1;
  int nsmax = iwork[11];
  int ljac, lmas, le;
  if (mljac == n) {
    ljac = n;
    le = n;
  } else {
    ljac = mljac + mujac + 1;
    le = 2 * mljac + mujac + 1;
  }
  if (imas == 0)
    lmas = 0;
  else if (mlmas = n)
    lmas = n;
  else
    lmas = mlmas + mumas + 1;
  if (nsmax == 0)
    nsmax = 7;
  LWORK = n * (ljac + lmas + nsmax * le + 3 * nsmax + 3) + 20;
  LIWORK = (2 + (nsmax - 1) / 2) * n + 20;
  WORK = malloc((LWORK) * sizeof(double));
  IWORK = malloc(LIWORK * sizeof(int));
  for (le = 0; le < 20; ++le) {
    WORK[le] = work[le];
    IWORK[le] = iwork[le];
  }
  radau_(&N, fcn, &X, y, &XEND, &H, &RTOL, &ATOL, &ITOL, jac, &IJAC, &MLJAC,
         &MUJAC, mas, &IMAS, &MLMAS, &MUMAS, solout, &IOUT, WORK, &LWORK, IWORK,
         &LIWORK, rpar, ipar, idid);
  for (le = 0; le < 20; ++le) {
    work[le] = WORK[le];
    iwork[le] = IWORK[le];
  }
  free(WORK);
  free(IWORK);
}


x86 assembly (clang, -O0):
<cradau>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1d8,%rsp
mov    0x68(%rbp),%rax
mov    0x60(%rbp),%rax
mov    0x58(%rbp),%rax
mov    0x50(%rbp),%rax
mov    0x48(%rbp),%rax
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0x2c(%rbp)
mov    %rsi,-0x38(%rbp)
movsd  %xmm0,-0x40(%rbp)
mov    %rdx,-0x48(%rbp)
movsd  %xmm1,-0x50(%rbp)
movsd  %xmm2,-0x58(%rbp)
movsd  %xmm3,-0x60(%rbp)
movsd  %xmm4,-0x68(%rbp)
mov    %rcx,-0x70(%rbp)
mov    %r8d,-0x74(%rbp)
mov    %r9d,-0x78(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x7c(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x0,-0x9c(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x74(%rbp),%eax
mov    %eax,-0xb4(%rbp)
mov    -0x78(%rbp),%eax
mov    %eax,-0xb8(%rbp)
mov    0x10(%rbp),%eax
mov    %eax,-0xbc(%rbp)
mov    0x20(%rbp),%eax
mov    %eax,-0xc0(%rbp)
mov    0x28(%rbp),%eax
mov    %eax,-0xc4(%rbp)
mov    0x30(%rbp),%eax
mov    %eax,-0xc8(%rbp)
mov    0x40(%rbp),%eax
mov    %eax,-0xcc(%rbp)
movl   $0x1,-0xec(%rbp)
mov    0x50(%rbp),%rax
mov    0x2c(%rax),%eax
mov    %eax,-0xf0(%rbp)
mov    -0x78(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jne    1269 <cradau+0x139>
mov    -0x2c(%rbp),%eax
mov    %eax,-0xf4(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0xfc(%rbp)
jmp    128a <cradau+0x15a>
mov    -0x78(%rbp),%eax
add    0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xf4(%rbp)
mov    -0x78(%rbp),%eax
shl    $0x1,%eax
add    0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xfc(%rbp)
cmpl   $0x0,0x20(%rbp)
jne    12a3 <cradau+0x173>
movl   $0x0,-0xf8(%rbp)
jmp    12d4 <cradau+0x1a4>
mov    -0x2c(%rbp),%eax
mov    %eax,0x28(%rbp)
cmp    $0x0,%eax
je     12c0 <cradau+0x190>
mov    -0x2c(%rbp),%eax
mov    %eax,-0xf8(%rbp)
jmp    12cf <cradau+0x19f>
mov    0x28(%rbp),%eax
add    0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xf8(%rbp)
jmp    12d4 <cradau+0x1a4>
cmpl   $0x0,-0xf0(%rbp)
jne    12eb <cradau+0x1bb>
movl   $0x7,-0xf0(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0xf4(%rbp),%ecx
add    -0xf8(%rbp),%ecx
mov    -0xf0(%rbp),%edx
imul   -0xfc(%rbp),%edx
add    %edx,%ecx
imul   $0x3,-0xf0(%rbp),%edx
add    %edx,%ecx
add    $0x3,%ecx
imul   %ecx,%eax
add    $0x14,%eax
mov    %eax,-0xdc(%rbp)
mov    -0xf0(%rbp),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    $0x2,%eax
imul   -0x2c(%rbp),%eax
add    $0x14,%eax
mov    %eax,-0xec(%rbp)
movslq -0xdc(%rbp),%rdi
shl    $0x3,%rdi
call   1040 <malloc@plt>
mov    %rax,-0xd8(%rbp)
movslq -0xec(%rbp),%rdi
shl    $0x2,%rdi
call   1040 <malloc@plt>
mov    %rax,-0xe8(%rbp)
movl   $0x0,-0xfc(%rbp)
cmpl   $0x14,-0xfc(%rbp)
jge    13dd <cradau+0x2ad>
mov    0x48(%rbp),%rax
movslq -0xfc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xd8(%rbp),%rax
movslq -0xfc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x50(%rbp),%rax
movslq -0xfc(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0xe8(%rbp),%rax
movslq -0xfc(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0xfc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xfc(%rbp)
jmp    137a <cradau+0x24a>
mov    -0x38(%rbp),%rsi
mov    -0x48(%rbp),%rcx
mov    -0x70(%rbp),%rax
mov    %rax,-0x150(%rbp)
mov    0x18(%rbp),%rax
mov    %rax,-0x130(%rbp)
mov    0x38(%rbp),%rax
mov    %rax,-0x110(%rbp)
mov    -0xd8(%rbp),%r12
mov    -0xe8(%rbp),%r14
mov    0x58(%rbp),%r11
mov    0x60(%rbp),%r10
mov    0x68(%rbp),%rax
mov    %rax,-0x108(%rbp)
lea    -0x7c(%rbp),%rdi
lea    -0x88(%rbp),%rdx
lea    -0x90(%rbp),%r8
lea    -0x98(%rbp),%r9
lea    -0xa8(%rbp),%rax
lea    -0xb0(%rbp),%rbx
mov    %rbx,-0x160(%rbp)
lea    -0x9c(%rbp),%rbx
mov    %rbx,-0x158(%rbp)
lea    -0xb4(%rbp),%rbx
mov    %rbx,-0x148(%rbp)
lea    -0xb8(%rbp),%rbx
mov    %rbx,-0x140(%rbp)
lea    -0xbc(%rbp),%rbx
mov    %rbx,-0x138(%rbp)
lea    -0xc0(%rbp),%rbx
mov    %rbx,-0x128(%rbp)
lea    -0xc4(%rbp),%rbx
mov    %rbx,-0x120(%rbp)
lea    -0xc8(%rbp),%rbx
mov    %rbx,-0x118(%rbp)
lea    -0xcc(%rbp),%r13
lea    -0xdc(%rbp),%r15
lea    -0xec(%rbp),%rbx
mov    %rax,(%rsp)
mov    -0x160(%rbp),%rax
mov    %rax,0x8(%rsp)
mov    -0x158(%rbp),%rax
mov    %rax,0x10(%rsp)
mov    -0x150(%rbp),%rax
mov    %rax,0x18(%rsp)
mov    -0x148(%rbp),%rax
mov    %rax,0x20(%rsp)
mov    -0x140(%rbp),%rax
mov    %rax,0x28(%rsp)
mov    -0x138(%rbp),%rax
mov    %rax,0x30(%rsp)
mov    -0x130(%rbp),%rax
mov    %rax,0x38(%rsp)
mov    -0x128(%rbp),%rax
mov    %rax,0x40(%rsp)
mov    -0x120(%rbp),%rax
mov    %rax,0x48(%rsp)
mov    -0x118(%rbp),%rax
mov    %rax,0x50(%rsp)
mov    -0x110(%rbp),%rax
mov    %rax,0x58(%rsp)
mov    -0x108(%rbp),%rax
mov    %r13,0x60(%rsp)
mov    %r12,0x68(%rsp)
mov    %r15,0x70(%rsp)
mov    %r14,0x78(%rsp)
mov    %rbx,0x80(%rsp)
mov    %r11,0x88(%rsp)
mov    %r10,0x90(%rsp)
mov    %rax,0x98(%rsp)
mov    $0x0,%al
call   1050 <radau_@plt>
movl   $0x0,-0xfc(%rbp)
cmpl   $0x14,-0xfc(%rbp)
jge    1603 <cradau+0x4d3>
mov    -0xd8(%rbp),%rax
movslq -0xfc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
movslq -0xfc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xe8(%rbp),%rax
movslq -0xfc(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    0x50(%rbp),%rax
movslq -0xfc(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0xfc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xfc(%rbp)
jmp    15a0 <cradau+0x470>
mov    -0xd8(%rbp),%rdi
call   1030 <free@plt>
mov    -0xe8(%rbp),%rdi
call   1030 <free@plt>
add    $0x1d8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

x86 assembly (clang, -O2):
<cradau>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    0x180(%rsp),%r14
mov    0x170(%rsp),%r11d
mov    0x160(%rsp),%eax
mov    0x158(%rsp),%ebp
mov    0x150(%rsp),%r10d
mov    %r8d,0xc4(%rsp)
mov    0x140(%rsp),%r15d
mov    %edi,0xcc(%rsp)
movsd  %xmm0,0x100(%rsp)
movsd  %xmm1,0xf8(%rsp)
movsd  %xmm2,0xf0(%rsp)
movl   $0x0,0xc8(%rsp)
movsd  %xmm3,0xe8(%rsp)
movsd  %xmm4,0xe0(%rsp)
mov    %r9d,0xc0(%rsp)
mov    %r15d,0xbc(%rsp)
mov    %r10d,0xb8(%rsp)
mov    %ebp,0xb4(%rsp)
mov    %eax,0xb0(%rsp)
mov    %r11d,0xac(%rsp)
lea    0x1(%r15),%r11d
xor    %r8d,%r8d
cmp    %edi,%r9d
cmove  %r8d,%r11d
mov    0x2c(%r14),%r8d
lea    0x1(%r15,%r9,2),%ebx
cmove  %r9d,%ebx
add    %r9d,%r11d
add    $0x1,%eax
test   %edi,%edi
cmovne %edi,%eax
test   %r8d,%r8d
mov    $0x7,%ebp
cmovne %r8d,%ebp
mov    %rcx,%r15
mov    %rdx,0xd8(%rsp)
mov    %rsi,0xd0(%rsp)
add    $0x3,%ebx
imul   %ebp,%ebx
add    $0x3,%eax
test   %r10d,%r10d
mov    $0x3,%ecx
cmovne %eax,%ecx
add    %r11d,%ecx
add    %ebx,%ecx
imul   %edi,%ecx
lea    0x14(%rcx),%eax
mov    %eax,0xa8(%rsp)
mov    %ebp,%eax
add    $0xffffffff,%eax
shr    $0x1f,%eax
add    %eax,%ebp
add    $0xffffffff,%ebp
sar    %ebp
add    $0x2,%ebp
imul   %edi,%ebp
lea    0x14(%rbp),%eax
mov    %eax,0xa4(%rsp)
mov    0x178(%rsp),%r12
movslq %ecx,%rax
lea    0xa0(,%rax,8),%rdi
call   1050 <malloc@plt>
mov    %rax,%rbx
movslq %ebp,%rax
lea    0x50(,%rax,4),%rdi
call   1050 <malloc@plt>
mov    %rax,%r13
mov    $0xa0,%edx
mov    %rbx,%rdi
mov    %r12,%rsi
call   1040 <memcpy@plt>
movups (%r14),%xmm0
movups 0x10(%r14),%xmm1
movups 0x20(%r14),%xmm2
movups 0x30(%r14),%xmm3
movups %xmm0,0x0(%r13)
movups %xmm1,0x10(%r13)
movups %xmm2,0x20(%r13)
movups %xmm3,0x30(%r13)
movups 0x40(%r14),%xmm0
movups %xmm0,0x40(%r13)
mov    0x198(%rsp),%rax
mov    %rax,0x98(%rsp)
movups 0x188(%rsp),%xmm0
movups %xmm0,0x88(%rsp)
lea    0xa4(%rsp),%rax
mov    %rax,0x80(%rsp)
mov    %r13,0x78(%rsp)
lea    0xa8(%rsp),%rax
mov    %rax,0x70(%rsp)
mov    %rbx,0x68(%rsp)
lea    0xac(%rsp),%rax
mov    %rax,0x60(%rsp)
mov    0x168(%rsp),%rax
mov    %rax,0x58(%rsp)
lea    0xb0(%rsp),%rax
mov    %rax,0x50(%rsp)
lea    0xb4(%rsp),%rax
mov    %rax,0x48(%rsp)
lea    0xb8(%rsp),%rax
mov    %rax,0x40(%rsp)
mov    0x148(%rsp),%rax
mov    %rax,0x38(%rsp)
lea    0xbc(%rsp),%rax
mov    %rax,0x30(%rsp)
lea    0xc0(%rsp),%rax
mov    %rax,0x28(%rsp)
lea    0xc4(%rsp),%rax
mov    %rax,0x20(%rsp)
mov    %r15,0x18(%rsp)
lea    0xc8(%rsp),%rax
mov    %rax,0x10(%rsp)
lea    0xe0(%rsp),%rax
mov    %rax,0x8(%rsp)
lea    0xe8(%rsp),%rax
mov    %rax,(%rsp)
lea    0xcc(%rsp),%rdi
lea    0x100(%rsp),%rdx
lea    0xf8(%rsp),%r8
lea    0xf0(%rsp),%r9
mov    0xd0(%rsp),%rsi
mov    0xd8(%rsp),%rcx
xor    %eax,%eax
call   1060 <radau_@plt>
movups (%rbx),%xmm0
movups %xmm0,(%r12)
movups 0x10(%rbx),%xmm0
movups %xmm0,0x10(%r12)
movups 0x0(%r13),%xmm0
movups %xmm0,(%r14)
movups 0x20(%rbx),%xmm0
movups %xmm0,0x20(%r12)
movups 0x30(%rbx),%xmm0
movups %xmm0,0x30(%r12)
movups 0x10(%r13),%xmm0
movups %xmm0,0x10(%r14)
movups 0x40(%rbx),%xmm0
movups %xmm0,0x40(%r12)
movups 0x50(%rbx),%xmm0
movups %xmm0,0x50(%r12)
movups 0x20(%r13),%xmm0
movups %xmm0,0x20(%r14)
movups 0x60(%rbx),%xmm0
movups %xmm0,0x60(%r12)
movups 0x70(%rbx),%xmm0
movups %xmm0,0x70(%r12)
movups 0x30(%r13),%xmm0
movups %xmm0,0x30(%r14)
movups 0x80(%rbx),%xmm0
movups %xmm0,0x80(%r12)
movups 0x90(%rbx),%xmm0
movups %xmm0,0x90(%r12)
movups 0x40(%r13),%xmm0
movups %xmm0,0x40(%r14)
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @cradau(i32 %0, void (i32*, double*, double*, double*, double*, i32*)* %1, double %2, double* %3, double %4, double %5, double %6, double %7, void (i32*, double*, double*, double*, i32*, double*, double*)* %8, i32 %9, i32 %10, i32 %11, void (i32*, double*, i32*, double*, i32*)* %12, i32 %13, i32 %14, i32 %15, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)* %16, i32 %17, double* %18, i32* %19, double* %20, i32* %21, i32* %22) {
%24 = alloca i32
%25 = alloca void (i32*, double*, double*, double*, double*, i32*)*
%26 = alloca double
%27 = alloca double*
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca void (i32*, double*, double*, double*, i32*, double*, double*)*
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
%36 = alloca void (i32*, double*, i32*, double*, i32*)*
%37 = alloca i32
%38 = alloca i32
%39 = alloca i32
%40 = alloca void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)*
%41 = alloca i32
%42 = alloca double*
%43 = alloca i32*
%44 = alloca double*
%45 = alloca i32*
%46 = alloca i32*
%47 = alloca i32
%48 = alloca double
%49 = alloca double
%50 = alloca double
%51 = alloca i32
%52 = alloca double
%53 = alloca double
%54 = alloca i32
%55 = alloca i32
%56 = alloca i32
%57 = alloca i32
%58 = alloca i32
%59 = alloca i32
%60 = alloca i32
%61 = alloca double*
%62 = alloca i32
%63 = alloca i32*
%64 = alloca i32
%65 = alloca i32
%66 = alloca i32
%67 = alloca i32
%68 = alloca i32
store i32 %0, i32* %24
store void (i32*, double*, double*, double*, double*, i32*)* %1, void (i32*, double*, double*, double*, double*, i32*)** %25
store double %2, double* %26
store double* %3, double** %27
store double %4, double* %28
store double %5, double* %29
store double %6, double* %30
store double %7, double* %31
store void (i32*, double*, double*, double*, i32*, double*, double*)* %8, void (i32*, double*, double*, double*, i32*, double*, double*)** %32
store i32 %9, i32* %33
store i32 %10, i32* %34
store i32 %11, i32* %35
store void (i32*, double*, i32*, double*, i32*)* %12, void (i32*, double*, i32*, double*, i32*)** %36
store i32 %13, i32* %37
store i32 %14, i32* %38
store i32 %15, i32* %39
store void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)* %16, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)** %40
store i32 %17, i32* %41
store double* %18, double** %42
store i32* %19, i32** %43
store double* %20, double** %44
store i32* %21, i32** %45
store i32* %22, i32** %46
%69 = load i32, i32* %24
store i32 %69, i32* %47
%70 = load double, double* %26
store double %70, double* %48
%71 = load double, double* %28
store double %71, double* %49
%72 = load double, double* %29
store double %72, double* %50
store i32 0, i32* %51
%73 = load double, double* %30
store double %73, double* %52
%74 = load double, double* %31
store double %74, double* %53
%75 = load i32, i32* %33
store i32 %75, i32* %54
%76 = load i32, i32* %34
store i32 %76, i32* %55
%77 = load i32, i32* %35
store i32 %77, i32* %56
%78 = load i32, i32* %37
store i32 %78, i32* %57
%79 = load i32, i32* %38
store i32 %79, i32* %58
%80 = load i32, i32* %39
store i32 %80, i32* %59
%81 = load i32, i32* %41
store i32 %81, i32* %60
store i32 1, i32* %64
%82 = load i32*, i32** %43
%83 = getelementptr inbounds i32, i32* %82, i64 11
%84 = load i32, i32* %83
store i32 %84, i32* %65
%85 = load i32, i32* %34
%86 = load i32, i32* %24
%87 = icmp eq i32 %85, %86
br i1 %87, label %88, label %91
88: ; preds = %23
%89 = load i32, i32* %24
store i32 %89, i32* %66
%90 = load i32, i32* %24
store i32 %90, i32* %68
br label %101
91: ; preds = %23
%92 = load i32, i32* %34
%93 = load i32, i32* %35
%94 = add nsw i32 %92, %93
%95 = add nsw i32 %94, 1
store i32 %95, i32* %66
%96 = load i32, i32* %34
%97 = mul nsw i32 2, %96
%98 = load i32, i32* %35
%99 = add nsw i32 %97, %98
%100 = add nsw i32 %99, 1
store i32 %100, i32* %68
br label %101
101: ; preds = %91, %88
%102 = load i32, i32* %37
%103 = icmp eq i32 %102, 0
br i1 %103, label %104, label %105
104: ; preds = %101
store i32 0, i32* %67
br label %116
105: ; preds = %101
%106 = load i32, i32* %24
store i32 %106, i32* %38
%107 = icmp ne i32 %106, 0
br i1 %107, label %108, label %110
108: ; preds = %105
%109 = load i32, i32* %24
store i32 %109, i32* %67
br label %115
110: ; preds = %105
%111 = load i32, i32* %38
%112 = load i32, i32* %39
%113 = add nsw i32 %111, %112
%114 = add nsw i32 %113, 1
store i32 %114, i32* %67
br label %115
115: ; preds = %110, %108
br label %116
116: ; preds = %115, %104
%117 = load i32, i32* %65
%118 = icmp eq i32 %117, 0
br i1 %118, label %119, label %120
119: ; preds = %116
store i32 7, i32* %65
br label %120
120: ; preds = %119, %116
%121 = load i32, i32* %24
%122 = load i32, i32* %66
%123 = load i32, i32* %67
%124 = add nsw i32 %122, %123
%125 = load i32, i32* %65
%126 = load i32, i32* %68
%127 = mul nsw i32 %125, %126
%128 = add nsw i32 %124, %127
%129 = load i32, i32* %65
%130 = mul nsw i32 3, %129
%131 = add nsw i32 %128, %130
%132 = add nsw i32 %131, 3
%133 = mul nsw i32 %121, %132
%134 = add nsw i32 %133, 20
store i32 %134, i32* %62
%135 = load i32, i32* %65
%136 = sub nsw i32 %135, 1
%137 = sdiv i32 %136, 2
%138 = add nsw i32 2, %137
%139 = load i32, i32* %24
%140 = mul nsw i32 %138, %139
%141 = add nsw i32 %140, 20
store i32 %141, i32* %64
%142 = load i32, i32* %62
%143 = sext i32 %142 to i64
%144 = mul i64 %143, 8
%145 = call i8* @malloc(i64 %144)
%146 = bitcast i8* %145 to double*
store double* %146, double** %61
%147 = load i32, i32* %64
%148 = sext i32 %147 to i64
%149 = mul i64 %148, 4
%150 = call i8* @malloc(i64 %149)
%151 = bitcast i8* %150 to i32*
store i32* %151, i32** %63
store i32 0, i32* %68
br label %152
152: ; preds = %174, %120
%153 = load i32, i32* %68
%154 = icmp slt i32 %153, 20
br i1 %154, label %155, label %177
155: ; preds = %152
%156 = load double*, double** %42
%157 = load i32, i32* %68
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %156, i64 %158
%160 = load double, double* %159
%161 = load double*, double** %61
%162 = load i32, i32* %68
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %161, i64 %163
store double %160, double* %164
%165 = load i32*, i32** %43
%166 = load i32, i32* %68
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds i32, i32* %165, i64 %167
%169 = load i32, i32* %168
%170 = load i32*, i32** %63
%171 = load i32, i32* %68
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds i32, i32* %170, i64 %172
store i32 %169, i32* %173
br label %174
174: ; preds = %155
%175 = load i32, i32* %68
%176 = add nsw i32 %175, 1
store i32 %176, i32* %68
br label %152
177: ; preds = %152
%178 = load void (i32*, double*, double*, double*, double*, i32*)*, void (i32*, double*, double*, double*, double*, i32*)** %25
%179 = load double*, double** %27
%180 = load void (i32*, double*, double*, double*, i32*, double*, double*)*, void (i32*, double*, double*, double*, i32*, double*, double*)** %32
%181 = load void (i32*, double*, i32*, double*, i32*)*, void (i32*, double*, i32*, double*, i32*)** %36
%182 = load void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)*, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)** %40
%183 = load double*, double** %61
%184 = load i32*, i32** %63
%185 = load double*, double** %44
%186 = load i32*, i32** %45
%187 = load i32*, i32** %46
%188 = call i32 (i32*, void (i32*, double*, double*, double*, double*, i32*)*, double*, double*, double*, double*, double*, double*, i32*, void (i32*, double*, double*, double*, i32*, double*, double*)*, i32*, i32*, i32*, void (i32*, double*, i32*, double*, i32*)*, i32*, i32*, i32*, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)*, i32*, double*, i32*, i32*, i32*, double*, i32*, i32*, ...) bitcast (i32 (...)* @radau_ to i32 (i32*, void (i32*, double*, double*, double*, double*, i32*)*, double*, double*, double*, double*, double*, double*, i32*, void (i32*, double*, double*, double*, i32*, double*, double*)*, i32*, i32*, i32*, void (i32*, double*, i32*, double*, i32*)*, i32*, i32*, i32*, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)*, i32*, double*, i32*, i32*, i32*, double*, i32*, i32*, ...)*)(i32* %47, void (i32*, double*, double*, double*, double*, i32*)* %178, double* %48, double* %179, double* %49, double* %50, double* %52, double* %53, i32* %51, void (i32*, double*, double*, double*, i32*, double*, double*)* %180, i32* %54, i32* %55, i32* %56, void (i32*, double*, i32*, double*, i32*)* %181, i32* %57, i32* %58, i32* %59, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)* %182, i32* %60, double* %183, i32* %62, i32* %184, i32* %64, double* %185, i32* %186, i32* %187)
store i32 0, i32* %68
br label %189
189: ; preds = %211, %177
%190 = load i32, i32* %68
%191 = icmp slt i32 %190, 20
br i1 %191, label %192, label %214
192: ; preds = %189
%193 = load double*, double** %61
%194 = load i32, i32* %68
%195 = sext i32 %194 to i64
%196 = getelementptr inbounds double, double* %193, i64 %195
%197 = load double, double* %196
%198 = load double*, double** %42
%199 = load i32, i32* %68
%200 = sext i32 %199 to i64
%201 = getelementptr inbounds double, double* %198, i64 %200
store double %197, double* %201
%202 = load i32*, i32** %63
%203 = load i32, i32* %68
%204 = sext i32 %203 to i64
%205 = getelementptr inbounds i32, i32* %202, i64 %204
%206 = load i32, i32* %205
%207 = load i32*, i32** %43
%208 = load i32, i32* %68
%209 = sext i32 %208 to i64
%210 = getelementptr inbounds i32, i32* %207, i64 %209
store i32 %206, i32* %210
br label %211
211: ; preds = %192
%212 = load i32, i32* %68
%213 = add nsw i32 %212, 1
store i32 %213, i32* %68
br label %189
214: ; preds = %189
%215 = load double*, double** %61
%216 = bitcast double* %215 to i8*
%217 = load i32*, i32** %63
%218 = bitcast i32* %217 to i8*
ret void
}

LLVM IR (-O2):
define void @cradau(i32 %0, void (i32*, double*, double*, double*, double*, i32*)* %1, double %2, double* %3, double %4, double %5, double %6, double %7, void (i32*, double*, double*, double*, i32*, double*, double*)* %8, i32 %9, i32 %10, i32 %11, void (i32*, double*, i32*, double*, i32*)* %12, i32 %13, i32 %14, i32 %15, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)* %16, i32 %17, double* nocapture %18, i32* nocapture %19, double* %20, i32* %21, i32* %22) {
%24 = alloca i32
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca i32
%29 = alloca double
%30 = alloca double
%31 = alloca i32
%32 = alloca i32
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
%36 = alloca i32
%37 = alloca i32
%38 = alloca i32
%39 = alloca i32
%40 = bitcast i32* %24 to i8*
store i32 %0, i32* %24
%41 = bitcast double* %25 to i8*
store double %2, double* %25
%42 = bitcast double* %26 to i8*
store double %4, double* %26
%43 = bitcast double* %27 to i8*
store double %5, double* %27
%44 = bitcast i32* %28 to i8*
store i32 0, i32* %28
%45 = bitcast double* %29 to i8*
store double %6, double* %29
%46 = bitcast double* %30 to i8*
store double %7, double* %30
%47 = bitcast i32* %31 to i8*
store i32 %9, i32* %31
%48 = bitcast i32* %32 to i8*
store i32 %10, i32* %32
%49 = bitcast i32* %33 to i8*
store i32 %11, i32* %33
%50 = bitcast i32* %34 to i8*
store i32 %13, i32* %34
%51 = bitcast i32* %35 to i8*
store i32 %14, i32* %35
%52 = bitcast i32* %36 to i8*
store i32 %15, i32* %36
%53 = bitcast i32* %37 to i8*
store i32 %17, i32* %37
%54 = bitcast i32* %38 to i8*
%55 = bitcast i32* %39 to i8*
%56 = getelementptr inbounds i32, i32* %19, i64 11
%57 = load i32, i32* %56
%58 = icmp eq i32 %10, %0
%59 = add i32 %11, 1
%60 = shl nsw i32 %10, 1
%61 = add i32 %59, %60
%62 = select i1 %58, i32 0, i32 %59
%63 = add i32 %62, %10
%64 = select i1 %58, i32 %10, i32 %61
%65 = bitcast i32* %19 to i8*
%66 = bitcast double* %18 to i8*
%67 = icmp eq i32 %13, 0
%68 = icmp eq i32 %0, 0
%69 = add nsw i32 %15, 1
%70 = select i1 %68, i32 %69, i32 %0
%71 = icmp eq i32 %57, 0
%72 = select i1 %71, i32 7, i32 %57
%73 = add i32 %64, 3
%74 = mul i32 %73, %72
%75 = add i32 %70, 3
%76 = select i1 %67, i32 3, i32 %75
%77 = add i32 %76, %63
%78 = add i32 %77, %74
%79 = mul nsw i32 %78, %0
%80 = add nsw i32 %79, 20
store i32 %80, i32* %38
%81 = add nsw i32 %72, -1
%82 = sdiv i32 %81, 2
%83 = add nsw i32 %82, 2
%84 = mul nsw i32 %83, %0
%85 = add nsw i32 %84, 20
store i32 %85, i32* %39
%86 = sext i32 %80 to i64
%87 = shl nsw i64 %86, 3
%88 = tail call i8* @malloc(i64 %87)
%89 = sext i32 %85 to i64
%90 = shl nsw i64 %89, 2
%91 = tail call i8* @malloc(i64 %90)
%92 = bitcast i8* %88 to double*
%93 = bitcast i8* %91 to i32*
%94 = call i32 (i32*, void (i32*, double*, double*, double*, double*, i32*)*, double*, double*, double*, double*, double*, double*, i32*, void (i32*, double*, double*, double*, i32*, double*, double*)*, i32*, i32*, i32*, void (i32*, double*, i32*, double*, i32*)*, i32*, i32*, i32*, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)*, i32*, double*, i32*, i32*, i32*, double*, i32*, i32*, ...) bitcast (i32 (...)* @radau_ to i32 (i32*, void (i32*, double*, double*, double*, double*, i32*)*, double*, double*, double*, double*, double*, double*, i32*, void (i32*, double*, double*, double*, i32*, double*, double*)*, i32*, i32*, i32*, void (i32*, double*, i32*, double*, i32*)*, i32*, i32*, i32*, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)*, i32*, double*, i32*, i32*, i32*, double*, i32*, i32*, ...)*)(i32* %24, void (i32*, double*, double*, double*, double*, i32*)* %1, double* %25, double* %3, double* %26, double* %27, double* %29, double* %30, i32* %28, void (i32*, double*, double*, double*, i32*, double*, double*)* %8, i32* %31, i32* %32, i32* %33, void (i32*, double*, i32*, double*, i32*)* %12, i32* %34, i32* %35, i32* %36, void (i32*, double*, double*, double*, double*, i32*, i32*, double*, i32*, i32*)* %16, i32* %37, double* %92, i32* %38, i32* %93, i32* %39, double* %20, i32* %21, i32* %22)
%95 = bitcast i8* %88 to <2 x double>*
%96 = load <2 x double>, <2 x double>* %95
%97 = bitcast double* %18 to <2 x double>*
store <2 x double> %96, <2 x double>* %97
%98 = getelementptr inbounds double, double* %92, i64 2
%99 = getelementptr inbounds double, double* %18, i64 2
%100 = bitcast double* %98 to <2 x double>*
%101 = load <2 x double>, <2 x double>* %100
%102 = bitcast double* %99 to <2 x double>*
store <2 x double> %101, <2 x double>* %102
%103 = bitcast i8* %91 to <4 x i32>*
%104 = load <4 x i32>, <4 x i32>* %103
%105 = bitcast i32* %19 to <4 x i32>*
store <4 x i32> %104, <4 x i32>* %105
%106 = getelementptr inbounds double, double* %92, i64 4
%107 = getelementptr inbounds double, double* %18, i64 4
%108 = getelementptr inbounds i32, i32* %93, i64 4
%109 = getelementptr inbounds i32, i32* %19, i64 4
%110 = bitcast double* %106 to <2 x double>*
%111 = load <2 x double>, <2 x double>* %110
%112 = bitcast double* %107 to <2 x double>*
store <2 x double> %111, <2 x double>* %112
%113 = getelementptr inbounds double, double* %92, i64 6
%114 = getelementptr inbounds double, double* %18, i64 6
%115 = bitcast double* %113 to <2 x double>*
%116 = load <2 x double>, <2 x double>* %115
%117 = bitcast double* %114 to <2 x double>*
store <2 x double> %116, <2 x double>* %117
%118 = bitcast i32* %108 to <4 x i32>*
%119 = load <4 x i32>, <4 x i32>* %118
%120 = bitcast i32* %109 to <4 x i32>*
store <4 x i32> %119, <4 x i32>* %120
%121 = getelementptr inbounds double, double* %92, i64 8
%122 = getelementptr inbounds double, double* %18, i64 8
%123 = getelementptr inbounds i32, i32* %93, i64 8
%124 = getelementptr inbounds i32, i32* %19, i64 8
%125 = bitcast double* %121 to <2 x double>*
%126 = load <2 x double>, <2 x double>* %125
%127 = bitcast double* %122 to <2 x double>*
store <2 x double> %126, <2 x double>* %127
%128 = getelementptr inbounds double, double* %92, i64 10
%129 = getelementptr inbounds double, double* %18, i64 10
%130 = bitcast double* %128 to <2 x double>*
%131 = load <2 x double>, <2 x double>* %130
%132 = bitcast double* %129 to <2 x double>*
store <2 x double> %131, <2 x double>* %132
%133 = bitcast i32* %123 to <4 x i32>*
%134 = load <4 x i32>, <4 x i32>* %133
%135 = bitcast i32* %124 to <4 x i32>*
store <4 x i32> %134, <4 x i32>* %135
%136 = getelementptr inbounds double, double* %92, i64 12
%137 = getelementptr inbounds double, double* %18, i64 12
%138 = getelementptr inbounds i32, i32* %93, i64 12
%139 = getelementptr inbounds i32, i32* %19, i64 12
%140 = bitcast double* %136 to <2 x double>*
%141 = load <2 x double>, <2 x double>* %140
%142 = bitcast double* %137 to <2 x double>*
store <2 x double> %141, <2 x double>* %142
%143 = getelementptr inbounds double, double* %92, i64 14
%144 = getelementptr inbounds double, double* %18, i64 14
%145 = bitcast double* %143 to <2 x double>*
%146 = load <2 x double>, <2 x double>* %145
%147 = bitcast double* %144 to <2 x double>*
store <2 x double> %146, <2 x double>* %147
%148 = bitcast i32* %138 to <4 x i32>*
%149 = load <4 x i32>, <4 x i32>* %148
%150 = bitcast i32* %139 to <4 x i32>*
store <4 x i32> %149, <4 x i32>* %150
%151 = getelementptr inbounds double, double* %92, i64 16
%152 = getelementptr inbounds double, double* %18, i64 16
%153 = getelementptr inbounds i32, i32* %93, i64 16
%154 = getelementptr inbounds i32, i32* %19, i64 16
%155 = bitcast double* %151 to <2 x double>*
%156 = load <2 x double>, <2 x double>* %155
%157 = bitcast double* %152 to <2 x double>*
store <2 x double> %156, <2 x double>* %157
%158 = getelementptr inbounds double, double* %92, i64 18
%159 = getelementptr inbounds double, double* %18, i64 18
%160 = bitcast double* %158 to <2 x double>*
%161 = load <2 x double>, <2 x double>* %160
%162 = bitcast double* %159 to <2 x double>*
store <2 x double> %161, <2 x double>* %162
%163 = bitcast i32* %153 to <4 x i32>*
%164 = load <4 x i32>, <4 x i32>* %163
%165 = bitcast i32* %154 to <4 x i32>*
store <4 x i32> %164, <4 x i32>* %165
ret void
}