#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_ivector(int *, int /*<<< orphan*/, int);
int *ivector(int /*<<< orphan*/, int);
int /*<<< orphan*/ resample_shuffle(int *, int, int, long *);
void permut_locations_of_groups(int n, double *x, double *y, double *z,
                                double **Mdij, int *groupi, double *xmix,
                                double *ymix, double *zmix, double **Mdijmix,
                                long *seed) {
  int i, j, *loci;
  int *locg, *group, newgroup, ng, g;
  loci = ivector(0, n);
  locg = ivector(0, n);
  group = ivector(0, n);
  ng = 0;
  for (i = 1; i <= n; i++) {
    newgroup = 1;
    for (g = ng; g >= 1; g--)
      if (groupi[i] == group[g]) {
        newgroup = 0;
        break;
      }
    if (newgroup) {
      ng++;
      group[ng] = groupi[i];
      locg[ng] = i;
    }
  }
  resample_shuffle(locg, 1, ng, seed);
  for (i = 1; i <= n; i++)
    loci[i] = locg[groupi[i]];
  for (i = 1; i <= n; i++) {
    if (Mdij[0][0] == 1.)
      for (j = 1; j <= n; j++)
        Mdijmix[i][j] = Mdij[loci[i]][loci[j]];
    else {
      xmix[i] = x[loci[i]];
      ymix[i] = y[loci[i]];
      zmix[i] = z[loci[i]];
    }
  }
  free_ivector(loci, 0, n);
  free_ivector(locg, 0, n);
  free_ivector(group, 0, n);
}


x86 assembly (clang, -O0):
<permut_locations_of_groups>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x4(%rbp),%esi
xor    %edi,%edi
call   1040 <ivector@plt>
mov    %rax,-0x40(%rbp)
mov    -0x4(%rbp),%esi
xor    %edi,%edi
call   1040 <ivector@plt>
mov    %rax,-0x48(%rbp)
mov    -0x4(%rbp),%esi
xor    %edi,%edi
call   1040 <ivector@plt>
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x58(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1243 <permut_locations_of_groups+0x113>
movl   $0x1,-0x54(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x5c(%rbp)
cmpl   $0x1,-0x5c(%rbp)
jl     11f9 <permut_locations_of_groups+0xc9>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x50(%rbp),%rcx
movslq -0x5c(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jne    11e6 <permut_locations_of_groups+0xb6>
movl   $0x0,-0x54(%rbp)
jmp    11f9 <permut_locations_of_groups+0xc9>
jmp    11eb <permut_locations_of_groups+0xbb>
mov    -0x5c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x5c(%rbp)
jmp    11b4 <permut_locations_of_groups+0x84>
cmpl   $0x0,-0x54(%rbp)
je     1230 <permut_locations_of_groups+0x100>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x50(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%edx
mov    -0x48(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    1235 <permut_locations_of_groups+0x105>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    119b <permut_locations_of_groups+0x6b>
mov    -0x48(%rbp),%rdi
mov    -0x58(%rbp),%edx
mov    0x30(%rbp),%rcx
mov    $0x1,%esi
call   1050 <resample_shuffle@plt>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1297 <permut_locations_of_groups+0x167>
mov    -0x48(%rbp),%rax
mov    -0x30(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    125f <permut_locations_of_groups+0x12f>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     13a6 <permut_locations_of_groups+0x276>
mov    -0x28(%rbp),%rax
mov    (%rax),%rax
movsd  (%rax),%xmm0
movsd  0xd43(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    132d <permut_locations_of_groups+0x1fd>
jp     132d <permut_locations_of_groups+0x1fd>
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1328 <permut_locations_of_groups+0x1f8>
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x40(%rbp),%rcx
movslq -0x38(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12d4 <permut_locations_of_groups+0x1a4>
jmp    1393 <permut_locations_of_groups+0x263>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x40(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1398 <permut_locations_of_groups+0x268>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    129e <permut_locations_of_groups+0x16e>
mov    -0x40(%rbp),%rdi
mov    -0x4(%rbp),%edx
xor    %esi,%esi
call   1030 <free_ivector@plt>
mov    -0x48(%rbp),%rdi
mov    -0x4(%rbp),%edx
xor    %esi,%esi
call   1030 <free_ivector@plt>
mov    -0x50(%rbp),%rdi
mov    -0x4(%rbp),%edx
xor    %esi,%esi
call   1030 <free_ivector@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<permut_locations_of_groups>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r9,%r15
mov    %r8,0x10(%rsp)
mov    %rcx,0x48(%rsp)
mov    %rdx,0x40(%rsp)
mov    %rsi,0x38(%rsp)
mov    %edi,%ebp
mov    0xb0(%rsp),%rax
mov    %rax,0x8(%rsp)
mov    0xa8(%rsp),%rax
mov    %rax,0x50(%rsp)
mov    0xa0(%rsp),%rax
mov    %rax,0x30(%rsp)
mov    0x98(%rsp),%rax
mov    %rax,0x28(%rsp)
mov    0x90(%rsp),%rax
mov    %rax,0x20(%rsp)
xor    %r14d,%r14d
xor    %edi,%edi
mov    %ebp,%esi
call   1040 <ivector@plt>
mov    %rax,%rbx
xor    %edi,%edi
mov    %ebp,%esi
call   1040 <ivector@plt>
mov    %rax,%r12
xor    %edi,%edi
mov    %ebp,%esi
call   1040 <ivector@plt>
mov    %rax,%r13
test   %ebp,%ebp
jle    140d <permut_locations_of_groups+0x2dd>
mov    %rbp,0x18(%rsp)
add    $0x1,%ebp
mov    $0x1,%eax
mov    0x8(%rsp),%rcx
jmp    11fa <permut_locations_of_groups+0xca>
nopl   0x0(%rax,%rax,1)
movslq %r14d,%rdx
add    $0x1,%r14d
mov    %edi,0x4(%r13,%rdx,4)
mov    %eax,0x4(%r12,%rdx,4)
add    $0x1,%rax
cmp    %rbp,%rax
je     1225 <permut_locations_of_groups+0xf5>
mov    (%r15,%rax,4),%edi
mov    %r14d,%edx
test   %r14d,%r14d
jle    11e0 <permut_locations_of_groups+0xb0>
cs nopw 0x0(%rax,%rax,1)
mov    %edx,%esi
cmp    0x0(%r13,%rsi,4),%edi
je     11f1 <permut_locations_of_groups+0xc1>
lea    -0x1(%rdx),%esi
cmp    $0x2,%edx
mov    %esi,%edx
jge    1210 <permut_locations_of_groups+0xe0>
jmp    11e0 <permut_locations_of_groups+0xb0>
mov    %r12,%rdi
mov    $0x1,%esi
mov    %r14d,%edx
call   1050 <resample_shuffle@plt>
mov    0x18(%rsp),%r14
test   %r14d,%r14d
jle    1424 <permut_locations_of_groups+0x2f4>
lea    -0x1(%rbp),%r10
lea    -0x2(%rbp),%r8
mov    %r10d,%edx
and    $0x3,%edx
mov    $0x1,%ecx
cmp    $0x3,%r8
jb     12b1 <permut_locations_of_groups+0x181>
mov    %r10,%rsi
and    $0xfffffffffffffffc,%rsi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movslq 0x4(%r15,%rcx,4),%rax
mov    (%r12,%rax,4),%eax
mov    %eax,0x4(%rbx,%rcx,4)
movslq 0x8(%r15,%rcx,4),%rax
mov    (%r12,%rax,4),%eax
mov    %eax,0x8(%rbx,%rcx,4)
movslq 0xc(%r15,%rcx,4),%rax
mov    (%r12,%rax,4),%eax
mov    %eax,0xc(%rbx,%rcx,4)
movslq 0x10(%r15,%rcx,4),%rax
mov    (%r12,%rax,4),%eax
mov    %eax,0x10(%rbx,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%rsi
jne    1270 <permut_locations_of_groups+0x140>
add    $0x1,%rcx
test   %rdx,%rdx
je     12d4 <permut_locations_of_groups+0x1a4>
lea    (%rbx,%rcx,4),%rsi
lea    (%r15,%rcx,4),%rcx
xor    %edi,%edi
movslq (%rcx,%rdi,4),%rax
mov    (%r12,%rax,4),%eax
mov    %eax,(%rsi,%rdi,4)
add    $0x1,%rdi
cmp    %rdi,%rdx
jne    12c0 <permut_locations_of_groups+0x190>
test   %r14d,%r14d
jle    1424 <permut_locations_of_groups+0x2f4>
mov    0x10(%rsp),%rax
mov    (%rax),%r9
mov    %r10d,%esi
and    $0x3,%esi
and    $0xfffffffffffffffc,%r10
mov    $0x1,%r11d
movsd  0xd03(%rip),%xmm0
mov    %rbp,%r15
jmp    135f <permut_locations_of_groups+0x22f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x38(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
mov    0x20(%rsp),%rcx
movsd  %xmm1,(%rcx,%r11,8)
mov    0x40(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
mov    0x28(%rsp),%rcx
movsd  %xmm1,(%rcx,%r11,8)
mov    0x48(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
mov    0x30(%rsp),%rax
movsd  %xmm1,(%rax,%r11,8)
add    $0x1,%r11
mov    %r15,%rax
cmp    %r15,%r11
je     1424 <permut_locations_of_groups+0x2f4>
movsd  (%r9),%xmm1
ucomisd %xmm0,%xmm1
movslq (%rbx,%r11,4),%rax
jne    1310 <permut_locations_of_groups+0x1e0>
jp     1310 <permut_locations_of_groups+0x1e0>
mov    0x10(%rsp),%rcx
mov    (%rcx,%rax,8),%rcx
mov    0x50(%rsp),%rax
mov    (%rax,%r11,8),%rdx
mov    $0x1,%ebp
cmp    $0x3,%r8
jb     13dd <permut_locations_of_groups+0x2ad>
xor    %ebp,%ebp
nop
movslq 0x4(%rbx,%rbp,4),%rax
movsd  (%rcx,%rax,8),%xmm1
movsd  %xmm1,0x8(%rdx,%rbp,8)
movslq 0x8(%rbx,%rbp,4),%rax
movsd  (%rcx,%rax,8),%xmm1
movsd  %xmm1,0x10(%rdx,%rbp,8)
movslq 0xc(%rbx,%rbp,4),%rax
movsd  (%rcx,%rax,8),%xmm1
movsd  %xmm1,0x18(%rdx,%rbp,8)
movslq 0x10(%rbx,%rbp,4),%rax
movsd  (%rcx,%rax,8),%xmm1
movsd  %xmm1,0x20(%rdx,%rbp,8)
add    $0x4,%rbp
cmp    %rbp,%r10
jne    1390 <permut_locations_of_groups+0x260>
add    $0x1,%rbp
test   %rsi,%rsi
je     134f <permut_locations_of_groups+0x21f>
lea    (%rdx,%rbp,8),%rdx
lea    (%rbx,%rbp,4),%rbp
xor    %eax,%eax
movslq 0x0(%rbp,%rax,4),%rdi
movsd  (%rcx,%rdi,8),%xmm1
movsd  %xmm1,(%rdx,%rax,8)
add    $0x1,%rax
cmp    %rax,%rsi
jne    13f0 <permut_locations_of_groups+0x2c0>
jmp    134f <permut_locations_of_groups+0x21f>
mov    %r12,%rdi
mov    $0x1,%esi
xor    %edx,%edx
mov    0x8(%rsp),%rcx
call   1050 <resample_shuffle@plt>
mov    %rbp,%r14
mov    %rbx,%rdi
xor    %esi,%esi
mov    %r14d,%edx
call   1030 <free_ivector@plt>
mov    %r12,%rdi
xor    %esi,%esi
mov    %r14d,%edx
call   1030 <free_ivector@plt>
mov    %r13,%rdi
xor    %esi,%esi
mov    %r14d,%edx
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free_ivector@plt>

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

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