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

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ free(unsigned int *);
scalar_t__ malloc(unsigned int);
int /*<<< orphan*/ memcpy(double *, double *, int);
unsigned int rand();
void dual_sample(double *uA, unsigned lenA, unsigned sA, double *uB,
                 unsigned lenB, unsigned sB, double *usam) {
  unsigned idx, pos, i;
  unsigned *ptrA = (unsigned *)malloc(lenA * sizeof(unsigned));
  unsigned *ptrB = (unsigned *)malloc(lenB * sizeof(unsigned));
  for (i = 0; i < lenA; ++i) {
    ptrA[i] = i;
  }
  for (i = 0; i < lenB; ++i) {
    ptrB[i] = i;
  }
  for (pos = 0; pos < sA; ++pos) {
    idx = rand() % lenA;
    i = ptrA[pos];
    ptrA[pos] = ptrA[idx];
    ptrA[idx] = i;
  }
  for (pos = 0; pos < sB; ++pos) {
    idx = rand() % lenB;
    i = ptrB[pos];
    ptrB[pos] = ptrB[idx];
    ptrB[idx] = i;
  }
  for (i = 0; i < sA; ++i) {
    memcpy(usam + 6 * i, uA + 6 * ptrA[i], 6 * sizeof(double));
  }
  for (i = 0; i < sB; ++i) {
    memcpy(usam + 6 * (i + sA), uB + 6 * ptrB[i], 6 * sizeof(double));
  }
  free(ptrA);
  free(ptrB);
}


x86 assembly (clang, -O0):
<dual_sample>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    11b4 <dual_sample+0x74>
mov    -0x2c(%rbp),%edx
mov    -0x38(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    118d <dual_sample+0x4d>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jae    11e2 <dual_sample+0xa2>
mov    -0x2c(%rbp),%edx
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11bb <dual_sample+0x7b>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x10(%rbp),%eax
jae    1240 <dual_sample+0x100>
mov    $0x0,%al
call   1060 <rand@plt>
xor    %edx,%edx
divl   -0xc(%rbp)
mov    %edx,-0x24(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    (%rax,%rcx,4),%edx
mov    -0x38(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%edx
mov    -0x38(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11e9 <dual_sample+0xa9>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jae    129e <dual_sample+0x15e>
mov    $0x0,%al
call   1060 <rand@plt>
xor    %edx,%edx
divl   -0x1c(%rbp)
mov    %edx,-0x24(%rbp)
mov    -0x40(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x40(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    (%rax,%rcx,4),%edx
mov    -0x40(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%edx
mov    -0x40(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1247 <dual_sample+0x107>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jae    12f2 <dual_sample+0x1b2>
mov    0x10(%rbp),%rdi
imul   $0x6,-0x2c(%rbp),%eax
mov    %eax,%eax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x8(%rbp),%rsi
mov    -0x38(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x6,(%rax,%rcx,4),%eax
mov    %eax,%eax
shl    $0x3,%rax
add    %rax,%rsi
mov    $0x30,%edx
call   1040 <memcpy@plt>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12a5 <dual_sample+0x165>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jae    134b <dual_sample+0x20b>
mov    0x10(%rbp),%rdi
mov    -0x2c(%rbp),%eax
add    -0x10(%rbp),%eax
imul   $0x6,%eax,%eax
mov    %eax,%eax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x18(%rbp),%rsi
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x6,(%rax,%rcx,4),%eax
mov    %eax,%eax
shl    $0x3,%rax
add    %rax,%rsi
mov    $0x30,%edx
call   1040 <memcpy@plt>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12f9 <dual_sample+0x1b9>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dual_sample>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r9d,0x14(%rsp)
mov    %r8d,%ebx
mov    %rcx,0x30(%rsp)
mov    %rdx,0x18(%rsp)
mov    %esi,%r12d
mov    %rdi,%r13
mov    0x70(%rsp),%rbp
lea    0x0(,%r12,4),%edi
call   1050 <malloc@plt>
mov    %rax,%r15
mov    %rbx,0x8(%rsp)
lea    0x0(,%rbx,4),%edi
call   1050 <malloc@plt>
mov    %rax,%r14
test   %r12d,%r12d
je     131d <dual_sample+0x1dd>
mov    %r12d,%eax
cmp    $0x8,%r12d
jae    11a8 <dual_sample+0x68>
xor    %ecx,%ecx
jmp    1310 <dual_sample+0x1d0>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0x18,%rdi
jae    11d6 <dual_sample+0x96>
movdqa 0xe31(%rip),%xmm0
xor    %edi,%edi
jmp    12b0 <dual_sample+0x170>
and    $0xfffffffffffffffc,%rsi
movdqa 0xe1e(%rip),%xmm0
xor    %edi,%edi
movdqa 0xe23(%rip),%xmm8
movdqa 0xe2a(%rip),%xmm9
movdqa 0xe31(%rip),%xmm10
movdqa 0xe39(%rip),%xmm4
movdqa 0xe41(%rip),%xmm5
movdqa 0xe49(%rip),%xmm6
movdqa 0xe51(%rip),%xmm7
movdqa 0xe59(%rip),%xmm1
nopw   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm2
paddd  %xmm8,%xmm2
movdqu %xmm0,(%r15,%rdi,4)
movdqu %xmm2,0x10(%r15,%rdi,4)
movdqa %xmm0,%xmm2
paddd  %xmm9,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm10,%xmm3
movdqu %xmm2,0x20(%r15,%rdi,4)
movdqu %xmm3,0x30(%r15,%rdi,4)
movdqa %xmm0,%xmm2
paddd  %xmm4,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm5,%xmm3
movdqu %xmm2,0x40(%r15,%rdi,4)
movdqu %xmm3,0x50(%r15,%rdi,4)
movdqa %xmm0,%xmm2
paddd  %xmm6,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm7,%xmm3
movdqu %xmm2,0x60(%r15,%rdi,4)
movdqu %xmm3,0x70(%r15,%rdi,4)
add    $0x20,%rdi
paddd  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rsi
jne    1230 <dual_sample+0xf0>
test   %rdx,%rdx
je     1300 <dual_sample+0x1c0>
lea    (%r15,%rdi,4),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %edi,%edi
movdqa 0xd45(%rip),%xmm1
movdqa 0xd4d(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movdqa %xmm0,%xmm3
paddd  %xmm1,%xmm3
movdqu %xmm0,-0x10(%rsi,%rdi,1)
movdqu %xmm3,(%rsi,%rdi,1)
paddd  %xmm2,%xmm0
add    $0x20,%rdi
cmp    %rdi,%rdx
jne    12e0 <dual_sample+0x1a0>
cmp    %rax,%rcx
je     131d <dual_sample+0x1dd>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ecx,(%r15,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1310 <dual_sample+0x1d0>
cmpl   $0x0,0x8(%rsp)
je     14ad <dual_sample+0x36d>
mov    0x8(%rsp),%rcx
mov    %ecx,%eax
cmp    $0x8,%ecx
jae    133b <dual_sample+0x1fb>
xor    %ecx,%ecx
jmp    14a0 <dual_sample+0x360>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0x18,%rdi
jae    1369 <dual_sample+0x229>
movdqa 0xc9e(%rip),%xmm0
xor    %edi,%edi
jmp    1440 <dual_sample+0x300>
and    $0xfffffffffffffffc,%rsi
movdqa 0xc8b(%rip),%xmm0
xor    %edi,%edi
movdqa 0xc90(%rip),%xmm8
movdqa 0xc97(%rip),%xmm9
movdqa 0xc9e(%rip),%xmm10
movdqa 0xca6(%rip),%xmm4
movdqa 0xcae(%rip),%xmm5
movdqa 0xcb6(%rip),%xmm6
movdqa 0xcbe(%rip),%xmm7
movdqa 0xcc6(%rip),%xmm1
nopw   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm2
paddd  %xmm8,%xmm2
movdqu %xmm0,(%r14,%rdi,4)
movdqu %xmm2,0x10(%r14,%rdi,4)
movdqa %xmm0,%xmm2
paddd  %xmm9,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm10,%xmm3
movdqu %xmm2,0x20(%r14,%rdi,4)
movdqu %xmm3,0x30(%r14,%rdi,4)
movdqa %xmm0,%xmm2
paddd  %xmm4,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm5,%xmm3
movdqu %xmm2,0x40(%r14,%rdi,4)
movdqu %xmm3,0x50(%r14,%rdi,4)
movdqa %xmm0,%xmm2
paddd  %xmm6,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm7,%xmm3
movdqu %xmm2,0x60(%r14,%rdi,4)
movdqu %xmm3,0x70(%r14,%rdi,4)
add    $0x20,%rdi
paddd  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rsi
jne    13c0 <dual_sample+0x280>
test   %rdx,%rdx
je     1490 <dual_sample+0x350>
lea    (%r14,%rdi,4),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %edi,%edi
movdqa 0xbb5(%rip),%xmm1
movdqa 0xbbd(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movdqa %xmm0,%xmm3
paddd  %xmm1,%xmm3
movdqu %xmm0,-0x10(%rsi,%rdi,1)
movdqu %xmm3,(%rsi,%rdi,1)
paddd  %xmm2,%xmm0
add    $0x20,%rdi
cmp    %rdi,%rdx
jne    1470 <dual_sample+0x330>
cmp    %rax,%rcx
je     14ad <dual_sample+0x36d>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ecx,(%r14,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    14a0 <dual_sample+0x360>
mov    %r13,0x28(%rsp)
mov    0x18(%rsp),%rax
mov    %eax,%r13d
test   %eax,%eax
je     14e5 <dual_sample+0x3a5>
xor    %ebx,%ebx
xor    %eax,%eax
call   1060 <rand@plt>
xor    %edx,%edx
div    %r12d
mov    (%r15,%rbx,4),%eax
mov    (%r15,%rdx,4),%ecx
mov    %ecx,(%r15,%rbx,4)
mov    %eax,(%r15,%rdx,4)
add    $0x1,%rbx
cmp    %rbx,%r13
jne    14c0 <dual_sample+0x380>
mov    %r15,0x20(%rsp)
mov    0x14(%rsp),%eax
mov    %eax,%r12d
test   %eax,%eax
mov    0x8(%rsp),%r15
je     1525 <dual_sample+0x3e5>
xor    %ebx,%ebx
nopl   0x0(%rax)
xor    %eax,%eax
call   1060 <rand@plt>
xor    %edx,%edx
div    %r15d
mov    (%r14,%rbx,4),%eax
mov    (%r14,%rdx,4),%ecx
mov    %ecx,(%r14,%rbx,4)
mov    %eax,(%r14,%rdx,4)
add    $0x1,%rbx
cmp    %rbx,%r12
jne    1500 <dual_sample+0x3c0>
mov    %r12,0x8(%rsp)
cmpl   $0x0,0x18(%rsp)
mov    %rbp,%r12
mov    0x28(%rsp),%r15
je     157c <dual_sample+0x43c>
add    %r13,%r13
lea    0x0(%r13,%r13,2),%r13
xor    %ebx,%ebx
mov    0x20(%rsp),%rbp
nopl   0x0(%rax,%rax,1)
mov    %ebx,%eax
and    $0xfffffffe,%eax
lea    (%r12,%rax,8),%rdi
mov    0x0(%rbp),%eax
add    %eax,%eax
lea    (%rax,%rax,2),%eax
lea    (%r15,%rax,8),%rsi
mov    $0x30,%edx
call   1040 <memcpy@plt>
add    $0x6,%rbx
add    $0x4,%rbp
cmp    %rbx,%r13
jne    1550 <dual_sample+0x410>
cmpl   $0x0,0x14(%rsp)
mov    0x30(%rsp),%r13
mov    0x8(%rsp),%r15
je     15d0 <dual_sample+0x490>
mov    0x18(%rsp),%rax
add    %eax,%eax
lea    (%rax,%rax,2),%ebx
xor    %ebp,%ebp
nopl   0x0(%rax)
mov    %ebx,%eax
lea    (%r12,%rax,8),%rdi
mov    (%r14,%rbp,4),%eax
add    %eax,%eax
lea    (%rax,%rax,2),%eax
lea    0x0(,%rax,8),%rsi
add    %r13,%rsi
mov    $0x30,%edx
call   1040 <memcpy@plt>
add    $0x1,%rbp
add    $0x6,%ebx
cmp    %rbp,%r15
jne    15a0 <dual_sample+0x460>
mov    0x20(%rsp),%rdi
call   1030 <free@plt>
mov    %r14,%rdi
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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