#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(int *);
scalar_t__ get_third_q_of_triplets_at_q(int **, int, int const *, int const *,
                                        int *);
size_t kpt_get_grid_point_double_mesh(int *, int *);
scalar_t__ malloc(int);
int get_BZ_triplets_at_q(int triplets[][3], const int grid_point,
                         int bz_grid_address[][3], const int bz_map[],
                         const int map_triplets[], const int num_map_triplets,
                         const int mesh[3]) {
  int i, j, k, num_ir;
  int bz_address[3][3], bz_address_double[3], bzmesh[3];
  int *ir_grid_points;
  for (i = 0; i < 3; i++) {
    bzmesh[i] = mesh[i] * 2;
  }
  num_ir = 0;
  ir_grid_points = (int *)malloc(sizeof(int) * num_map_triplets);
  for (i = 0; i < num_map_triplets; i++) {
    if (map_triplets[i] == i) {
      ir_grid_points[num_ir] = i;
      num_ir++;
    }
  }
#pragma omp parallel for private(j, k, bz_address, bz_address_double)
  for (i = 0; i < num_ir; i++) {
    for (j = 0; j < 3; j++) {
      bz_address[0][j] = bz_grid_address[grid_point][j];
      bz_address[1][j] = bz_grid_address[ir_grid_points[i]][j];
      bz_address[2][j] = -bz_address[0][j] - bz_address[1][j];
    }
    for (j = 2; j > -1; j--) {
      if (get_third_q_of_triplets_at_q(bz_address, j, bz_map, mesh, bzmesh) ==
          0) {
        break;
      }
    }
    for (j = 0; j < 3; j++) {
      for (k = 0; k < 3; k++) {
        bz_address_double[k] = bz_address[j][k] * 2;
      }
      triplets[i][j] =
          bz_map[kpt_get_grid_point_double_mesh(bz_address_double, bzmesh)];
    }
  }
  free(ir_grid_points);
  return num_ir;
}


x86 assembly (clang, -O0):
<get_BZ_triplets_at_q>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
cmpl   $0x3,-0x30(%rbp)
jge    119b <get_BZ_triplets_at_q+0x5b>
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
shl    $0x1,%ecx
movslq -0x30(%rbp),%rax
mov    %ecx,-0x78(%rbp,%rax,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    116d <get_BZ_triplets_at_q+0x2d>
movl   $0x0,-0x3c(%rbp)
movslq -0x2c(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x80(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1206 <get_BZ_triplets_at_q+0xc6>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x30(%rbp),%eax
jne    11f3 <get_BZ_triplets_at_q+0xb3>
mov    -0x30(%rbp),%edx
mov    -0x80(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11f8 <get_BZ_triplets_at_q+0xb8>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11bc <get_BZ_triplets_at_q+0x7c>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    138e <get_BZ_triplets_at_q+0x24e>
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    1296 <get_BZ_triplets_at_q+0x156>
mov    -0x18(%rbp),%rax
movslq -0xc(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
movslq -0x34(%rbp),%rax
mov    %ecx,-0x60(%rbp,%rax,4)
mov    -0x18(%rbp),%rax
mov    -0x80(%rbp),%rcx
movslq -0x30(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
movslq -0x34(%rbp),%rax
mov    %ecx,-0x54(%rbp,%rax,4)
movslq -0x34(%rbp),%rax
xor    %ecx,%ecx
sub    -0x60(%rbp,%rax,4),%ecx
movslq -0x34(%rbp),%rax
sub    -0x54(%rbp,%rax,4),%ecx
movslq -0x34(%rbp),%rax
mov    %ecx,-0x48(%rbp,%rax,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1220 <get_BZ_triplets_at_q+0xe0>
movl   $0x2,-0x34(%rbp)
cmpl   $0xffffffff,-0x34(%rbp)
jle    12e1 <get_BZ_triplets_at_q+0x1a1>
lea    -0x60(%rbp),%rdi
mov    -0x34(%rbp),%esi
mov    -0x20(%rbp),%rdx
mov    0x10(%rbp),%rcx
lea    -0x78(%rbp),%r8
call   1030 <get_third_q_of_triplets_at_q@plt>
cmp    $0x0,%rax
jne    12ce <get_BZ_triplets_at_q+0x18e>
jmp    12e1 <get_BZ_triplets_at_q+0x1a1>
jmp    12d3 <get_BZ_triplets_at_q+0x193>
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    129d <get_BZ_triplets_at_q+0x15d>
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    137b <get_BZ_triplets_at_q+0x23b>
movl   $0x0,-0x38(%rbp)
cmpl   $0x3,-0x38(%rbp)
jge    1332 <get_BZ_triplets_at_q+0x1f2>
movslq -0x34(%rbp),%rcx
lea    -0x60(%rbp),%rax
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
shl    $0x1,%ecx
movslq -0x38(%rbp),%rax
mov    %ecx,-0x6c(%rbp,%rax,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12f9 <get_BZ_triplets_at_q+0x1b9>
mov    -0x20(%rbp),%rax
mov    %rax,-0x88(%rbp)
lea    -0x6c(%rbp),%rdi
lea    -0x78(%rbp),%rsi
call   1050 <kpt_get_grid_point_double_mesh@plt>
mov    %rax,%rcx
mov    -0x88(%rbp),%rax
mov    (%rax,%rcx,4),%edx
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12e8 <get_BZ_triplets_at_q+0x1a8>
jmp    1380 <get_BZ_triplets_at_q+0x240>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    120d <get_BZ_triplets_at_q+0xcd>
mov    -0x80(%rbp),%rdi
call   1040 <free@plt>
mov    -0x3c(%rbp),%eax
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_BZ_triplets_at_q>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9d,%r12d
mov    %r8,%r14
mov    %rcx,%rbx
mov    %rdx,0x50(%rsp)
mov    %esi,%r15d
mov    %rdi,0x80(%rsp)
mov    0xc0(%rsp),%rbp
mov    0x0(%rbp),%eax
add    %eax,%eax
mov    %eax,0x5c(%rsp)
mov    0x4(%rbp),%eax
add    %eax,%eax
mov    %eax,0x60(%rsp)
mov    0x8(%rbp),%eax
add    %eax,%eax
mov    %eax,0x64(%rsp)
lea    0x0(,%r12,4),%edi
call   1060 <malloc@plt>
mov    %rax,0x18(%rsp)
test   %r12d,%r12d
jle    11bd <get_BZ_triplets_at_q+0x7d>
mov    %r12d,%r8d
cmp    $0x1,%r12d
jne    11c5 <get_BZ_triplets_at_q+0x85>
xor    %ecx,%ecx
xor    %r12d,%r12d
test   $0x1,%r8b
jne    122d <get_BZ_triplets_at_q+0xed>
jmp    1246 <get_BZ_triplets_at_q+0x106>
xor    %r12d,%r12d
jmp    13f2 <get_BZ_triplets_at_q+0x2b2>
mov    %r8d,%edx
and    $0xfffffffe,%edx
xor    %esi,%esi
xor    %r12d,%r12d
mov    0x18(%rsp),%rax
jmp    11e9 <get_BZ_triplets_at_q+0xa9>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rsi
cmp    %rdx,%rsi
je     1223 <get_BZ_triplets_at_q+0xe3>
mov    %rsi,%rcx
mov    (%r14,%rsi,4),%esi
cmp    %rsi,%rcx
jne    1205 <get_BZ_triplets_at_q+0xc5>
movslq %r12d,%r12
mov    0x18(%rsp),%rsi
mov    %ecx,(%rsi,%r12,4)
add    $0x1,%r12d
mov    0x4(%r14,%rcx,4),%edi
lea    0x1(%rcx),%rsi
cmp    %rdi,%rsi
jne    11e0 <get_BZ_triplets_at_q+0xa0>
movslq %r12d,%r12
lea    0x1(%rcx),%edi
mov    %edi,(%rax,%r12,4)
add    $0x1,%r12d
jmp    11e0 <get_BZ_triplets_at_q+0xa0>
add    $0x2,%rcx
test   $0x1,%r8b
je     1246 <get_BZ_triplets_at_q+0x106>
mov    (%r14,%rcx,4),%eax
cmp    %rax,%rcx
jne    1246 <get_BZ_triplets_at_q+0x106>
movslq %r12d,%r12
mov    0x18(%rsp),%rax
mov    %ecx,(%rax,%r12,4)
add    $0x1,%r12d
test   %r12d,%r12d
jle    13f2 <get_BZ_triplets_at_q+0x2b2>
movslq %r15d,%rax
mov    %r12,0x68(%rsp)
mov    %r12d,%ecx
mov    %rcx,0x78(%rsp)
lea    (%rax,%rax,2),%rax
mov    0x50(%rsp),%rcx
lea    (%rcx,%rax,4),%rax
mov    %rax,0x70(%rsp)
xor    %r12d,%r12d
lea    0x20(%rsp),%r15
lea    0x5c(%rsp),%r14
jmp    134c <get_BZ_triplets_at_q+0x20c>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0x20(%rsp),%eax
mov    0x24(%rsp),%ecx
add    %eax,%eax
mov    %eax,0xc(%rsp)
add    %ecx,%ecx
mov    %ecx,0x10(%rsp)
mov    0x28(%rsp),%eax
add    %eax,%eax
mov    %eax,0x14(%rsp)
lea    0xc(%rsp),%rdi
mov    %r14,%rsi
call   1050 <kpt_get_grid_point_double_mesh@plt>
mov    (%rbx,%rax,4),%eax
lea    (%r12,%r12,2),%rbp
mov    0x80(%rsp),%r13
mov    %eax,0x0(%r13,%rbp,4)
mov    0x2c(%rsp),%eax
mov    0x30(%rsp),%ecx
add    %eax,%eax
mov    %eax,0xc(%rsp)
add    %ecx,%ecx
mov    %ecx,0x10(%rsp)
mov    0x34(%rsp),%eax
add    %eax,%eax
mov    %eax,0x14(%rsp)
lea    0xc(%rsp),%rdi
mov    %r14,%rsi
call   1050 <kpt_get_grid_point_double_mesh@plt>
mov    (%rbx,%rax,4),%eax
mov    %eax,0x4(%r13,%rbp,4)
mov    0x38(%rsp),%eax
mov    0x3c(%rsp),%ecx
add    %eax,%eax
mov    %eax,0xc(%rsp)
add    %ecx,%ecx
mov    %ecx,0x10(%rsp)
mov    0x40(%rsp),%eax
add    %eax,%eax
mov    %eax,0x14(%rsp)
lea    0xc(%rsp),%rdi
mov    %r14,%rsi
call   1050 <kpt_get_grid_point_double_mesh@plt>
mov    (%rbx,%rax,4),%eax
mov    %eax,0x8(%r13,%rbp,4)
mov    0xc0(%rsp),%rbp
add    $0x1,%r12
cmp    0x78(%rsp),%r12
je     13ed <get_BZ_triplets_at_q+0x2ad>
mov    0x70(%rsp),%rdi
mov    (%rdi),%eax
mov    %eax,0x20(%rsp)
mov    0x18(%rsp),%rsi
movslq (%rsi,%r12,4),%rcx
lea    (%rcx,%rcx,2),%rcx
mov    0x50(%rsp),%rdx
mov    (%rdx,%rcx,4),%ecx
mov    %ecx,0x2c(%rsp)
add    %ecx,%eax
neg    %eax
mov    %eax,0x38(%rsp)
mov    0x4(%rdi),%eax
mov    %eax,0x24(%rsp)
movslq (%rsi,%r12,4),%rcx
lea    (%rcx,%rcx,2),%rcx
mov    0x4(%rdx,%rcx,4),%ecx
mov    %ecx,0x30(%rsp)
add    %ecx,%eax
neg    %eax
mov    %eax,0x3c(%rsp)
mov    0x8(%rdi),%eax
mov    %eax,0x28(%rsp)
movslq (%rsi,%r12,4),%rcx
lea    (%rcx,%rcx,2),%rcx
mov    0x8(%rdx,%rcx,4),%ecx
mov    %ecx,0x34(%rsp)
add    %ecx,%eax
neg    %eax
mov    %eax,0x40(%rsp)
mov    $0x2,%r13d
nopl   0x0(%rax)
mov    %r15,%rdi
mov    %r13d,%esi
mov    %rbx,%rdx
mov    %rbp,%rcx
mov    %r14,%r8
call   1030 <get_third_q_of_triplets_at_q@plt>
add    $0xffffffff,%r13d
setb   %cl
test   %rax,%rax
je     1290 <get_BZ_triplets_at_q+0x150>
test   %cl,%cl
jne    13c0 <get_BZ_triplets_at_q+0x280>
jmp    1290 <get_BZ_triplets_at_q+0x150>
mov    0x68(%rsp),%r12
mov    0x18(%rsp),%rdi
call   1040 <free@plt>
mov    %r12d,%eax
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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