#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 *);
size_t get_grid_point(int *, int *);
int /*<<< orphan*/ get_third_q_of_triplets_at_q(int **, int const *,
                                                int const *, 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 weights[], const int mesh[3]) {
  int i, j, k, num_ir;
  int address[3][3], address_double[3], bzmesh[3], bzmesh_double[3];
  int *ir_grid_points;
  for (i = 0; i < 3; i++) {
    bzmesh[i] = mesh[i] * 2;
    bzmesh_double[i] = bzmesh[i] * 2;
  }
  num_ir = 0;
  ir_grid_points = (int *)malloc(sizeof(int) * mesh[0] * mesh[1] * mesh[2]);
  for (i = 0; i < mesh[0] * mesh[1] * mesh[2]; i++) {
    if (weights[i] > 0) {
      ir_grid_points[num_ir] = i;
      num_ir++;
    }
  }
#pragma omp parallel for private(j, k, address, address_double)
  for (i = 0; i < num_ir; i++) {
    for (j = 0; j < 3; j++) {
      address[0][j] = bz_grid_address[grid_point][j];
      address[1][j] = bz_grid_address[ir_grid_points[i]][j];
      address[2][j] = -address[0][j] - address[1][j];
    }
    get_third_q_of_triplets_at_q(address, bz_map, mesh, bzmesh, bzmesh_double);
    for (j = 0; j < 3; j++) {
      for (k = 0; k < 3; k++) {
        address_double[k] = address[j][k] * 2;
        if (address_double[k] < 0) {
          address_double[k] += bzmesh_double[k];
        }
      }
      triplets[i][j] = bz_map[get_grid_point(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    $0xb0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    11b3 <get_BZ_triplets_at_q+0x73>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
shl    $0x1,%ecx
movslq -0x34(%rbp),%rax
mov    %ecx,-0x88(%rbp,%rax,4)
movslq -0x34(%rbp),%rax
mov    -0x88(%rbp,%rax,4),%ecx
shl    $0x1,%ecx
movslq -0x34(%rbp),%rax
mov    %ecx,-0x94(%rbp,%rax,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1169 <get_BZ_triplets_at_q+0x29>
movl   $0x0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
movslq (%rax),%rax
shl    $0x2,%rax
mov    -0x30(%rbp),%rcx
movslq 0x4(%rcx),%rcx
imul   %rcx,%rax
mov    -0x30(%rbp),%rcx
movslq 0x8(%rcx),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0xa0(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x30(%rbp),%rcx
mov    (%rcx),%ecx
mov    -0x30(%rbp),%rdx
imul   0x4(%rdx),%ecx
mov    -0x30(%rbp),%rdx
imul   0x8(%rdx),%ecx
cmp    %ecx,%eax
jge    1252 <get_BZ_triplets_at_q+0x112>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    123f <get_BZ_triplets_at_q+0xff>
mov    -0x34(%rbp),%edx
mov    -0xa0(%rbp),%rax
movslq -0x40(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1244 <get_BZ_triplets_at_q+0x104>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11f2 <get_BZ_triplets_at_q+0xb2>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    13df <get_BZ_triplets_at_q+0x29f>
movl   $0x0,-0x38(%rbp)
cmpl   $0x3,-0x38(%rbp)
jge    12e5 <get_BZ_triplets_at_q+0x1a5>
mov    -0x18(%rbp),%rax
movslq -0xc(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
movslq -0x38(%rbp),%rax
mov    %ecx,-0x70(%rbp,%rax,4)
mov    -0x18(%rbp),%rax
mov    -0xa0(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
movslq -0x38(%rbp),%rax
mov    %ecx,-0x64(%rbp,%rax,4)
movslq -0x38(%rbp),%rax
xor    %ecx,%ecx
sub    -0x70(%rbp,%rax,4),%ecx
movslq -0x38(%rbp),%rax
sub    -0x64(%rbp,%rax,4),%ecx
movslq -0x38(%rbp),%rax
mov    %ecx,-0x58(%rbp,%rax,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    126c <get_BZ_triplets_at_q+0x12c>
lea    -0x70(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x30(%rbp),%rdx
lea    -0x88(%rbp),%rcx
lea    -0x94(%rbp),%r8
call   1030 <get_third_q_of_triplets_at_q@plt>
movl   $0x0,-0x38(%rbp)
cmpl   $0x3,-0x38(%rbp)
jge    13cc <get_BZ_triplets_at_q+0x28c>
movl   $0x0,-0x3c(%rbp)
cmpl   $0x3,-0x3c(%rbp)
jge    1380 <get_BZ_triplets_at_q+0x240>
movslq -0x38(%rbp),%rcx
lea    -0x70(%rbp),%rax
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
shl    $0x1,%ecx
movslq -0x3c(%rbp),%rax
mov    %ecx,-0x7c(%rbp,%rax,4)
movslq -0x3c(%rbp),%rax
cmpl   $0x0,-0x7c(%rbp,%rax,4)
jge    136d <get_BZ_triplets_at_q+0x22d>
movslq -0x3c(%rbp),%rax
mov    -0x94(%rbp,%rax,4),%ecx
movslq -0x3c(%rbp),%rax
add    -0x7c(%rbp,%rax,4),%ecx
mov    %ecx,-0x7c(%rbp,%rax,4)
jmp    1372 <get_BZ_triplets_at_q+0x232>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    131c <get_BZ_triplets_at_q+0x1dc>
mov    -0x20(%rbp),%rax
mov    %rax,-0xa8(%rbp)
lea    -0x7c(%rbp),%rdi
lea    -0x88(%rbp),%rsi
call   1050 <get_grid_point@plt>
mov    %rax,%rcx
mov    -0xa8(%rbp),%rax
mov    (%rax,%rcx,4),%edx
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0xc,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    130b <get_BZ_triplets_at_q+0x1cb>
jmp    13d1 <get_BZ_triplets_at_q+0x291>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1259 <get_BZ_triplets_at_q+0x119>
mov    -0xa0(%rbp),%rdi
call   1040 <free@plt>
mov    -0x40(%rbp),%eax
add    $0xb0,%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    %r9,%r15
mov    %r8,%r14
mov    %rcx,0x18(%rsp)
mov    %rdx,0x50(%rsp)
mov    %esi,%r13d
mov    %rdi,%r12
mov    (%r9),%eax
lea    (%rax,%rax,1),%ecx
mov    %ecx,0x5c(%rsp)
shl    $0x2,%eax
mov    %eax,0xc(%rsp)
mov    0x4(%r9),%ecx
lea    (%rcx,%rcx,1),%edx
mov    %edx,0x60(%rsp)
lea    0x0(,%rcx,4),%edx
mov    %edx,0x10(%rsp)
mov    0x8(%r9),%edi
lea    (%rdi,%rdi,1),%edx
mov    %edx,0x64(%rsp)
lea    0x0(,%rdi,4),%edx
mov    %edx,0x14(%rsp)
imul   %ecx,%edi
imul   %eax,%edi
call   1060 <malloc@plt>
mov    (%r15),%ebp
mov    0x4(%r15),%ecx
mov    %ecx,%edx
imul   %ebp,%edx
mov    0x8(%r15),%esi
imul   %esi,%edx
test   %edx,%edx
jle    1400 <get_BZ_triplets_at_q+0x2c0>
xor    %edx,%edx
xor    %ebx,%ebx
jmp    11e4 <get_BZ_triplets_at_q+0xa4>
add    $0x1,%rdx
mov    %ecx,%edi
imul   %ebp,%edi
imul   %esi,%edi
movslq %edi,%rdi
cmp    %rdi,%rdx
jge    1201 <get_BZ_triplets_at_q+0xc1>
cmpl   $0x0,(%r14,%rdx,4)
jle    11d0 <get_BZ_triplets_at_q+0x90>
movslq %ebx,%rbx
mov    %edx,(%rax,%rbx,4)
add    $0x1,%ebx
mov    (%r15),%ebp
mov    0x4(%r15),%ecx
mov    0x8(%r15),%esi
jmp    11d0 <get_BZ_triplets_at_q+0x90>
test   %ebx,%ebx
jle    1409 <get_BZ_triplets_at_q+0x2c9>
movslq %r13d,%rdx
mov    %rbx,0x68(%rsp)
mov    %ebx,%ecx
mov    %rcx,0x80(%rsp)
lea    (%rdx,%rdx,2),%rdx
mov    0x50(%rsp),%rcx
lea    (%rcx,%rdx,4),%rcx
mov    %rcx,0x70(%rsp)
add    $0x8,%r12
xor    %r14d,%r14d
lea    0x5c(%rsp),%rcx
mov    0x18(%rsp),%rsi
mov    %rax,0x78(%rsp)
mov    0x70(%rsp),%rbp
jmp    1282 <get_BZ_triplets_at_q+0x142>
nopw   0x0(%rax,%rax,1)
mov    %rbx,%rdi
call   1050 <get_grid_point@plt>
mov    %r13,%rcx
mov    0x18(%rsp),%rsi
mov    (%rsi,%rax,4),%eax
mov    %eax,(%r12)
add    $0x1,%r14
add    $0xc,%r12
cmp    %r14,0x80(%rsp)
mov    0x78(%rsp),%rax
je     1404 <get_BZ_triplets_at_q+0x2c4>
mov    0x0(%rbp),%ebx
mov    %ebx,0x20(%rsp)
movslq (%rax,%r14,4),%rdx
lea    (%rdx,%rdx,2),%rdi
mov    0x50(%rsp),%rdx
mov    (%rdx,%rdi,4),%edi
mov    %edi,0x2c(%rsp)
add    %edi,%ebx
neg    %ebx
mov    %ebx,0x38(%rsp)
mov    0x4(%rbp),%ebx
mov    %ebx,0x24(%rsp)
movslq (%rax,%r14,4),%rdi
lea    (%rdi,%rdi,2),%rdi
mov    0x4(%rdx,%rdi,4),%edi
mov    %edi,0x30(%rsp)
add    %edi,%ebx
neg    %ebx
mov    %ebx,0x3c(%rsp)
mov    0x8(%rbp),%ebx
mov    %ebx,0x28(%rsp)
movslq (%rax,%r14,4),%rdi
lea    (%rdi,%rdi,2),%rdi
mov    0x8(%rdx,%rdi,4),%edx
mov    %edx,0x34(%rsp)
add    %edx,%ebx
neg    %ebx
mov    %ebx,0x40(%rsp)
lea    0x20(%rsp),%rdi
mov    %r15,%rdx
mov    %rcx,%r13
lea    0xc(%rsp),%r8
call   1030 <get_third_q_of_triplets_at_q@plt>
mov    0x20(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,(%rsp)
test   %ecx,%ecx
jns    130d <get_BZ_triplets_at_q+0x1cd>
add    0xc(%rsp),%eax
mov    %eax,(%rsp)
mov    0x24(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,0x4(%rsp)
test   %ecx,%ecx
jns    1324 <get_BZ_triplets_at_q+0x1e4>
add    0x10(%rsp),%eax
mov    %eax,0x4(%rsp)
mov    0x28(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,0x8(%rsp)
test   %ecx,%ecx
mov    %r13,%rsi
mov    %rsp,%rbx
jns    1341 <get_BZ_triplets_at_q+0x201>
add    0x14(%rsp),%eax
mov    %eax,0x8(%rsp)
mov    %rbx,%rdi
call   1050 <get_grid_point@plt>
mov    0x18(%rsp),%rcx
mov    (%rcx,%rax,4),%eax
mov    %eax,-0x8(%r12)
mov    0x2c(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,(%rsp)
test   %ecx,%ecx
jns    136b <get_BZ_triplets_at_q+0x22b>
add    0xc(%rsp),%eax
mov    %eax,(%rsp)
mov    0x30(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,0x4(%rsp)
test   %ecx,%ecx
jns    1382 <get_BZ_triplets_at_q+0x242>
add    0x10(%rsp),%eax
mov    %eax,0x4(%rsp)
mov    0x34(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,0x8(%rsp)
test   %ecx,%ecx
mov    %r13,%rsi
jns    139c <get_BZ_triplets_at_q+0x25c>
add    0x14(%rsp),%eax
mov    %eax,0x8(%rsp)
mov    %rbx,%rdi
call   1050 <get_grid_point@plt>
mov    0x18(%rsp),%rcx
mov    (%rcx,%rax,4),%eax
mov    %eax,-0x4(%r12)
mov    0x38(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,(%rsp)
test   %ecx,%ecx
jns    13c6 <get_BZ_triplets_at_q+0x286>
add    0xc(%rsp),%eax
mov    %eax,(%rsp)
mov    0x3c(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,0x4(%rsp)
test   %ecx,%ecx
jns    13dd <get_BZ_triplets_at_q+0x29d>
add    0x10(%rsp),%eax
mov    %eax,0x4(%rsp)
mov    0x40(%rsp),%ecx
lea    (%rcx,%rcx,1),%eax
mov    %eax,0x8(%rsp)
test   %ecx,%ecx
mov    %r13,%rsi
jns    1250 <get_BZ_triplets_at_q+0x110>
add    0x14(%rsp),%eax
mov    %eax,0x8(%rsp)
jmp    1250 <get_BZ_triplets_at_q+0x110>
xor    %ebx,%ebx
jmp    1409 <get_BZ_triplets_at_q+0x2c9>
mov    0x68(%rsp),%rbx
mov    %rax,%rdi
call   1040 <free@plt>
mov    %ebx,%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) {
%7 = alloca [3 x i32]*
%8 = alloca i32
%9 = alloca [3 x i32]*
%10 = alloca i32*
%11 = alloca i32*
%12 = alloca i32*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca [3 x [3 x i32]]
%18 = alloca [3 x i32]
%19 = alloca [3 x i32]
%20 = alloca [3 x i32]
%21 = alloca i32*
store [3 x i32]* %0, [3 x i32]** %7
store i32 %1, i32* %8
store [3 x i32]* %2, [3 x i32]** %9
store i32* %3, i32** %10
store i32* %4, i32** %11
store i32* %5, i32** %12
store i32 0, i32* %13
br label %22
22: ; preds = %43, %6
%23 = load i32, i32* %13
%24 = icmp slt i32 %23, 3
br i1 %24, label %25, label %46
25: ; preds = %22
%26 = load i32*, i32** %12
%27 = load i32, i32* %13
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds i32, i32* %26, i64 %28
%30 = load i32, i32* %29
%31 = mul nsw i32 %30, 2
%32 = load i32, i32* %13
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds [3 x i32], [3 x i32]* %19, i64 0, i64 %33
store i32 %31, i32* %34
%35 = load i32, i32* %13
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [3 x i32], [3 x i32]* %19, i64 0, i64 %36
%38 = load i32, i32* %37
%39 = mul nsw i32 %38, 2
%40 = load i32, i32* %13
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds [3 x i32], [3 x i32]* %20, i64 0, i64 %41
store i32 %39, i32* %42
br label %43
43: ; preds = %25
%44 = load i32, i32* %13
%45 = add nsw i32 %44, 1
store i32 %45, i32* %13
br label %22
46: ; preds = %22
store i32 0, i32* %16
%47 = load i32*, i32** %12
%48 = getelementptr inbounds i32, i32* %47, i64 0
%49 = load i32, i32* %48
%50 = sext i32 %49 to i64
%51 = mul i64 4, %50
%52 = load i32*, i32** %12
%53 = getelementptr inbounds i32, i32* %52, i64 1
%54 = load i32, i32* %53
%55 = sext i32 %54 to i64
%56 = mul i64 %51, %55
%57 = load i32*, i32** %12
%58 = getelementptr inbounds i32, i32* %57, i64 2
%59 = load i32, i32* %58
%60 = sext i32 %59 to i64
%61 = mul i64 %56, %60
%62 = trunc i64 %61 to i32
%63 = call i64 @malloc(i32 %62)
%64 = inttoptr i64 %63 to i32*
store i32* %64, i32** %21
store i32 0, i32* %13
br label %65
65: ; preds = %95, %46
%66 = load i32, i32* %13
%67 = load i32*, i32** %12
%68 = getelementptr inbounds i32, i32* %67, i64 0
%69 = load i32, i32* %68
%70 = load i32*, i32** %12
%71 = getelementptr inbounds i32, i32* %70, i64 1
%72 = load i32, i32* %71
%73 = mul nsw i32 %69, %72
%74 = load i32*, i32** %12
%75 = getelementptr inbounds i32, i32* %74, i64 2
%76 = load i32, i32* %75
%77 = mul nsw i32 %73, %76
%78 = icmp slt i32 %66, %77
br i1 %78, label %79, label %98
79: ; preds = %65
%80 = load i32*, i32** %11
%81 = load i32, i32* %13
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds i32, i32* %80, i64 %82
%84 = load i32, i32* %83
%85 = icmp sgt i32 %84, 0
br i1 %85, label %86, label %94
86: ; preds = %79
%87 = load i32, i32* %13
%88 = load i32*, i32** %21
%89 = load i32, i32* %16
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds i32, i32* %88, i64 %90
store i32 %87, i32* %91
%92 = load i32, i32* %16
%93 = add nsw i32 %92, 1
store i32 %93, i32* %16
br label %94
94: ; preds = %86, %79
br label %95
95: ; preds = %94
%96 = load i32, i32* %13
%97 = add nsw i32 %96, 1
store i32 %97, i32* %13
br label %65
98: ; preds = %65
store i32 0, i32* %13
br label %99
99: ; preds = %219, %98
%100 = load i32, i32* %13
%101 = load i32, i32* %16
%102 = icmp slt i32 %100, %101
br i1 %102, label %103, label %222
103: ; preds = %99
store i32 0, i32* %14
br label %104
104: ; preds = %152, %103
%105 = load i32, i32* %14
%106 = icmp slt i32 %105, 3
br i1 %106, label %107, label %155
107: ; preds = %104
%108 = load [3 x i32]*, [3 x i32]** %9
%109 = load i32, i32* %8
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds [3 x i32], [3 x i32]* %108, i64 %110
%112 = load i32, i32* %14
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds [3 x i32], [3 x i32]* %111, i64 0, i64 %113
%115 = load i32, i32* %114
%116 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %17, i64 0, i64 0
%117 = load i32, i32* %14
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds [3 x i32], [3 x i32]* %116, i64 0, i64 %118
store i32 %115, i32* %119
%120 = load [3 x i32]*, [3 x i32]** %9
%121 = load i32*, i32** %21
%122 = load i32, i32* %13
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds i32, i32* %121, i64 %123
%125 = load i32, i32* %124
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds [3 x i32], [3 x i32]* %120, i64 %126
%128 = load i32, i32* %14
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds [3 x i32], [3 x i32]* %127, i64 0, i64 %129
%131 = load i32, i32* %130
%132 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %17, i64 0, i64 1
%133 = load i32, i32* %14
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds [3 x i32], [3 x i32]* %132, i64 0, i64 %134
store i32 %131, i32* %135
%136 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %17, i64 0, i64 0
%137 = load i32, i32* %14
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds [3 x i32], [3 x i32]* %136, i64 0, i64 %138
%140 = load i32, i32* %139
%141 = sub nsw i32 0, %140
%142 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %17, i64 0, i64 1
%143 = load i32, i32* %14
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds [3 x i32], [3 x i32]* %142, i64 0, i64 %144
%146 = load i32, i32* %145
%147 = sub nsw i32 %141, %146
%148 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %17, i64 0, i64 2
%149 = load i32, i32* %14
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds [3 x i32], [3 x i32]* %148, i64 0, i64 %150
store i32 %147, i32* %151
br label %152
152: ; preds = %107
%153 = load i32, i32* %14
%154 = add nsw i32 %153, 1
store i32 %154, i32* %14
br label %104
155: ; preds = %104
%156 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %17, i64 0, i64 0
%157 = bitcast [3 x i32]* %156 to i32**
%158 = load i32*, i32** %10
%159 = load i32*, i32** %12
%160 = getelementptr inbounds [3 x i32], [3 x i32]* %19, i64 0, i64 0
%161 = getelementptr inbounds [3 x i32], [3 x i32]* %20, i64 0, i64 0
%162 = call i32 @get_third_q_of_triplets_at_q(i32** %157, i32* %158, i32* %159, i32* %160, i32* %161)
store i32 0, i32* %14
br label %163
163: ; preds = %215, %155
%164 = load i32, i32* %14
%165 = icmp slt i32 %164, 3
br i1 %165, label %166, label %218
166: ; preds = %163
store i32 0, i32* %15
br label %167
167: ; preds = %198, %166
%168 = load i32, i32* %15
%169 = icmp slt i32 %168, 3
br i1 %169, label %170, label %201
170: ; preds = %167
%171 = load i32, i32* %14
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds [3 x [3 x i32]], [3 x [3 x i32]]* %17, i64 0, i64 %172
%174 = load i32, i32* %15
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds [3 x i32], [3 x i32]* %173, i64 0, i64 %175
%177 = load i32, i32* %176
%178 = mul nsw i32 %177, 2
%179 = load i32, i32* %15
%180 = sext i32 %179 to i64
%181 = getelementptr inbounds [3 x i32], [3 x i32]* %18, i64 0, i64 %180
store i32 %178, i32* %181
%182 = load i32, i32* %15
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds [3 x i32], [3 x i32]* %18, i64 0, i64 %183
%185 = load i32, i32* %184
%186 = icmp slt i32 %185, 0
br i1 %186, label %187, label %197
187: ; preds = %170
%188 = load i32, i32* %15
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds [3 x i32], [3 x i32]* %20, i64 0, i64 %189
%191 = load i32, i32* %190
%192 = load i32, i32* %15
%193 = sext i32 %192 to i64
%194 = getelementptr inbounds [3 x i32], [3 x i32]* %18, i64 0, i64 %193
%195 = load i32, i32* %194
%196 = add nsw i32 %195, %191
store i32 %196, i32* %194
br label %197
197: ; preds = %187, %170
br label %198
198: ; preds = %197
%199 = load i32, i32* %15
%200 = add nsw i32 %199, 1
store i32 %200, i32* %15
br label %167
201: ; preds = %167
%202 = load i32*, i32** %10
%203 = getelementptr inbounds [3 x i32], [3 x i32]* %18, i64 0, i64 0
%204 = getelementptr inbounds [3 x i32], [3 x i32]* %19, i64 0, i64 0
%205 = call i64 @get_grid_point(i32* %203, i32* %204)
%206 = getelementptr inbounds i32, i32* %202, i64 %205
%207 = load i32, i32* %206
%208 = load [3 x i32]*, [3 x i32]** %7
%209 = load i32, i32* %13
%210 = sext i32 %209 to i64
%211 = getelementptr inbounds [3 x i32], [3 x i32]* %208, i64 %210
%212 = load i32, i32* %14
%213 = sext i32 %212 to i64
%214 = getelementptr inbounds [3 x i32], [3 x i32]* %211, i64 0, i64 %213
store i32 %207, i32* %214
br label %215
215: ; preds = %201
%216 = load i32, i32* %14
%217 = add nsw i32 %216, 1
store i32 %217, i32* %14
br label %163
218: ; preds = %163
br label %219
219: ; preds = %218
%220 = load i32, i32* %13
%221 = add nsw i32 %220, 1
store i32 %221, i32* %13
br label %99
222: ; preds = %99
%223 = load i32*, i32** %21
%224 = call i32 @free(i32* %223)
%225 = load i32, i32* %16
ret i32 %225
}

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