#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*/ GaussianFiltering3D_float(float *, float *, int *, double,
                                             int);
int /*<<< orphan*/ free(float *);
float *mallocf(int);
void StructureTensor3D(float *ux, float *uy, float *uz, float **J, int *dimsu,
                       double rho) {
  int npixelsu;
  int i;
  float *Jxx, *Jxy, *Jxz, *Jyy, *Jyz, *Jzz;
  npixelsu = dimsu[0] * dimsu[1] * dimsu[2];
  Jxx = mallocf(npixelsu);
  Jyy = mallocf(npixelsu);
  Jzz = mallocf(npixelsu);
  Jxy = mallocf(npixelsu);
  Jxz = mallocf(npixelsu);
  Jyz = mallocf(npixelsu);
#pragma omp parallel for
  for (i = 0; i < npixelsu; i++) {
    Jxx[i] = ux[i] * ux[i];
    Jxy[i] = ux[i] * uy[i];
    Jxz[i] = ux[i] * uz[i];
    Jyy[i] = uy[i] * uy[i];
    Jyz[i] = uy[i] * uz[i];
    Jzz[i] = uz[i] * uz[i];
  }
  J[0] = mallocf(npixelsu);
  GaussianFiltering3D_float(Jxx, J[0], dimsu, rho, 4 * rho);
  free(Jxx);
  J[1] = mallocf(npixelsu);
  GaussianFiltering3D_float(Jyy, J[1], dimsu, rho, 4 * rho);
  free(Jyy);
  J[2] = mallocf(npixelsu);
  GaussianFiltering3D_float(Jzz, J[2], dimsu, rho, 4 * rho);
  free(Jzz);
  J[3] = mallocf(npixelsu);
  GaussianFiltering3D_float(Jxy, J[3], dimsu, rho, 4 * rho);
  free(Jxy);
  J[4] = mallocf(npixelsu);
  GaussianFiltering3D_float(Jxz, J[4], dimsu, rho, 4 * rho);
  free(Jxz);
  J[5] = mallocf(npixelsu);
  GaussianFiltering3D_float(Jyz, J[5], dimsu, rho, 4 * rho);
  free(Jyz);
}


x86 assembly (clang, -O0):
<StructureTensor3D>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    -0x28(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x28(%rbp),%rcx
imul   0x8(%rcx),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,-0x40(%rbp)
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,-0x58(%rbp)
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,-0x68(%rbp)
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,-0x48(%rbp)
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,-0x50(%rbp)
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12bd <StructureTensor3D+0x18d>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x68(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11b9 <StructureTensor3D+0x89>
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x40(%rbp),%rdi
mov    -0x20(%rbp),%rax
mov    (%rax),%rsi
mov    -0x28(%rbp),%rdx
movsd  -0x30(%rbp),%xmm0
movsd  0xd15(%rip),%xmm1
mulsd  -0x30(%rbp),%xmm1
cvttsd2si %xmm1,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,0x8(%rax)
mov    -0x58(%rbp),%rdi
mov    -0x20(%rbp),%rax
mov    0x8(%rax),%rsi
mov    -0x28(%rbp),%rdx
movsd  -0x30(%rbp),%xmm0
movsd  0xcce(%rip),%xmm1
mulsd  -0x30(%rbp),%xmm1
cvttsd2si %xmm1,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,0x10(%rax)
mov    -0x68(%rbp),%rdi
mov    -0x20(%rbp),%rax
mov    0x10(%rax),%rsi
mov    -0x28(%rbp),%rdx
movsd  -0x30(%rbp),%xmm0
movsd  0xc87(%rip),%xmm1
mulsd  -0x30(%rbp),%xmm1
cvttsd2si %xmm1,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    -0x68(%rbp),%rdi
call   1030 <free@plt>
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,0x18(%rax)
mov    -0x48(%rbp),%rdi
mov    -0x20(%rbp),%rax
mov    0x18(%rax),%rsi
mov    -0x28(%rbp),%rdx
movsd  -0x30(%rbp),%xmm0
movsd  0xc40(%rip),%xmm1
mulsd  -0x30(%rbp),%xmm1
cvttsd2si %xmm1,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,0x20(%rax)
mov    -0x50(%rbp),%rdi
mov    -0x20(%rbp),%rax
mov    0x20(%rax),%rsi
mov    -0x28(%rbp),%rdx
movsd  -0x30(%rbp),%xmm0
movsd  0xbf9(%rip),%xmm1
mulsd  -0x30(%rbp),%xmm1
cvttsd2si %xmm1,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
mov    -0x34(%rbp),%edi
call   1040 <mallocf@plt>
mov    %rax,%rcx
mov    -0x20(%rbp),%rax
mov    %rcx,0x28(%rax)
mov    -0x60(%rbp),%rdi
mov    -0x20(%rbp),%rax
mov    0x28(%rax),%rsi
mov    -0x28(%rbp),%rdx
movsd  -0x30(%rbp),%xmm0
movsd  0xbb2(%rip),%xmm1
mulsd  -0x30(%rbp),%xmm1
cvttsd2si %xmm1,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<StructureTensor3D>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
movsd  %xmm0,(%rsp)
mov    %rcx,0x28(%rsp)
mov    %rdx,%r13
mov    %rsi,%rbp
mov    %rdi,%rbx
mov    0x4(%r8),%r14d
imul   (%r8),%r14d
mov    %r8,0x30(%rsp)
imul   0x8(%r8),%r14d
mov    %r14d,%edi
call   1040 <mallocf@plt>
mov    %rax,%r12
mov    %r14d,%edi
call   1040 <mallocf@plt>
mov    %rax,%r15
mov    %r14d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x8(%rsp)
mov    %r14d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x10(%rsp)
mov    %r14d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x18(%rsp)
mov    %r14d,%edi
call   1040 <mallocf@plt>
mov    0x8(%rsp),%r8
mov    0x10(%rsp),%rdi
mov    0x18(%rsp),%rsi
test   %r14d,%r14d
jle    123c <StructureTensor3D+0x10c>
mov    %r14d,%edx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rbx,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
movss  %xmm0,(%r12,%rcx,4)
movss  (%rbx,%rcx,4),%xmm0
mulss  0x0(%rbp,%rcx,4),%xmm0
movss  %xmm0,(%rdi,%rcx,4)
movss  (%rbx,%rcx,4),%xmm0
mulss  0x0(%r13,%rcx,4),%xmm0
movss  %xmm0,(%rsi,%rcx,4)
movss  0x0(%rbp,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
movss  %xmm0,(%r15,%rcx,4)
movss  0x0(%rbp,%rcx,4),%xmm0
mulss  0x0(%r13,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  0x0(%r13,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
movss  %xmm0,(%r8,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rdx
jne    11d0 <StructureTensor3D+0xa0>
mov    %r14d,%r13d
mov    %rax,0x20(%rsp)
mov    %r14d,%edi
call   1040 <mallocf@plt>
mov    0x28(%rsp),%rcx
mov    %rax,(%rcx)
mov    %rcx,%rbx
movsd  0xda1(%rip),%xmm1
movsd  (%rsp),%xmm0
mulsd  %xmm0,%xmm1
cvttsd2si %xmm1,%ecx
mov    %r12,%rdi
mov    %rax,%rsi
mov    0x30(%rsp),%r14
mov    %r14,%rdx
mov    %ecx,%ebp
call   1050 <GaussianFiltering3D_float@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %r13d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x8(%rbx)
mov    %r15,%rdi
mov    %rax,%rsi
mov    %r14,%rdx
movsd  (%rsp),%xmm0
mov    %ebp,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %r13d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x10(%rbx)
mov    %rbx,%r12
mov    0x8(%rsp),%rbx
mov    %rbx,%rdi
mov    %rax,%rsi
mov    %r14,%rdx
movsd  (%rsp),%xmm0
mov    %ebp,%ecx
mov    %ebp,%r15d
call   1050 <GaussianFiltering3D_float@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r13d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x18(%r12)
mov    0x10(%rsp),%rbp
mov    %rbp,%rdi
mov    %rax,%rsi
mov    %r14,%rdx
movsd  (%rsp),%xmm0
mov    %r15d,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %r13d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x20(%r12)
mov    0x18(%rsp),%rbp
mov    %rbp,%rdi
mov    %rax,%rsi
mov    %r14,%rdx
movsd  (%rsp),%xmm0
mov    %r15d,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %r13d,%edi
call   1040 <mallocf@plt>
mov    %rax,0x28(%r12)
mov    0x20(%rsp),%rbp
mov    %rbp,%rdi
mov    %rax,%rsi
mov    %r14,%rdx
movsd  (%rsp),%xmm0
mov    %r15d,%ecx
call   1050 <GaussianFiltering3D_float@plt>
mov    %rbp,%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 @StructureTensor3D(float* %0, float* %1, float* %2, float** %3, i32* %4, double %5) {
%7 = alloca float*
%8 = alloca float*
%9 = alloca float*
%10 = alloca float**
%11 = alloca i32*
%12 = alloca double
%13 = alloca i32
%14 = alloca i32
%15 = alloca float*
%16 = alloca float*
%17 = alloca float*
%18 = alloca float*
%19 = alloca float*
%20 = alloca float*
store float* %0, float** %7
store float* %1, float** %8
store float* %2, float** %9
store float** %3, float*** %10
store i32* %4, i32** %11
store double %5, double* %12
%21 = load i32*, i32** %11
%22 = getelementptr inbounds i32, i32* %21, i64 0
%23 = load i32, i32* %22
%24 = load i32*, i32** %11
%25 = getelementptr inbounds i32, i32* %24, i64 1
%26 = load i32, i32* %25
%27 = mul nsw i32 %23, %26
%28 = load i32*, i32** %11
%29 = getelementptr inbounds i32, i32* %28, i64 2
%30 = load i32, i32* %29
%31 = mul nsw i32 %27, %30
store i32 %31, i32* %13
%32 = load i32, i32* %13
%33 = call float* @mallocf(i32 %32)
store float* %33, float** %15
%34 = load i32, i32* %13
%35 = call float* @mallocf(i32 %34)
store float* %35, float** %18
%36 = load i32, i32* %13
%37 = call float* @mallocf(i32 %36)
store float* %37, float** %20
%38 = load i32, i32* %13
%39 = call float* @mallocf(i32 %38)
store float* %39, float** %16
%40 = load i32, i32* %13
%41 = call float* @mallocf(i32 %40)
store float* %41, float** %17
%42 = load i32, i32* %13
%43 = call float* @mallocf(i32 %42)
store float* %43, float** %19
store i32 0, i32* %14
br label %44
44: ; preds = %139, %6
%45 = load i32, i32* %14
%46 = load i32, i32* %13
%47 = icmp slt i32 %45, %46
br i1 %47, label %48, label %142
48: ; preds = %44
%49 = load float*, float** %7
%50 = load i32, i32* %14
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %49, i64 %51
%53 = load float, float* %52
%54 = load float*, float** %7
%55 = load i32, i32* %14
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds float, float* %54, i64 %56
%58 = load float, float* %57
%59 = fmul float %53, %58
%60 = load float*, float** %15
%61 = load i32, i32* %14
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %60, i64 %62
store float %59, float* %63
%64 = load float*, float** %7
%65 = load i32, i32* %14
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds float, float* %64, i64 %66
%68 = load float, float* %67
%69 = load float*, float** %8
%70 = load i32, i32* %14
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds float, float* %69, i64 %71
%73 = load float, float* %72
%74 = fmul float %68, %73
%75 = load float*, float** %16
%76 = load i32, i32* %14
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float, float* %75, i64 %77
store float %74, float* %78
%79 = load float*, float** %7
%80 = load i32, i32* %14
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float, float* %79, i64 %81
%83 = load float, float* %82
%84 = load float*, float** %9
%85 = load i32, i32* %14
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds float, float* %84, i64 %86
%88 = load float, float* %87
%89 = fmul float %83, %88
%90 = load float*, float** %17
%91 = load i32, i32* %14
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds float, float* %90, i64 %92
store float %89, float* %93
%94 = load float*, float** %8
%95 = load i32, i32* %14
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds float, float* %94, i64 %96
%98 = load float, float* %97
%99 = load float*, float** %8
%100 = load i32, i32* %14
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds float, float* %99, i64 %101
%103 = load float, float* %102
%104 = fmul float %98, %103
%105 = load float*, float** %18
%106 = load i32, i32* %14
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds float, float* %105, i64 %107
store float %104, float* %108
%109 = load float*, float** %8
%110 = load i32, i32* %14
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds float, float* %109, i64 %111
%113 = load float, float* %112
%114 = load float*, float** %9
%115 = load i32, i32* %14
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds float, float* %114, i64 %116
%118 = load float, float* %117
%119 = fmul float %113, %118
%120 = load float*, float** %19
%121 = load i32, i32* %14
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds float, float* %120, i64 %122
store float %119, float* %123
%124 = load float*, float** %9
%125 = load i32, i32* %14
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds float, float* %124, i64 %126
%128 = load float, float* %127
%129 = load float*, float** %9
%130 = load i32, i32* %14
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds float, float* %129, i64 %131
%133 = load float, float* %132
%134 = fmul float %128, %133
%135 = load float*, float** %20
%136 = load i32, i32* %14
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds float, float* %135, i64 %137
store float %134, float* %138
br label %139
139: ; preds = %48
%140 = load i32, i32* %14
%141 = add nsw i32 %140, 1
store i32 %141, i32* %14
br label %44
142: ; preds = %44
%143 = load i32, i32* %13
%144 = call float* @mallocf(i32 %143)
%145 = load float**, float*** %10
%146 = getelementptr inbounds float*, float** %145, i64 0
store float* %144, float** %146
%147 = load float*, float** %15
%148 = load float**, float*** %10
%149 = getelementptr inbounds float*, float** %148, i64 0
%150 = load float*, float** %149
%151 = load i32*, i32** %11
%152 = load double, double* %12
%153 = load double, double* %12
%154 = fmul double 4.000000e+00, %153
%155 = fptosi double %154 to i32
%156 = call i32 @GaussianFiltering3D_float(float* %147, float* %150, i32* %151, double %152, i32 %155)
%157 = load float*, float** %15
%158 = call i32 @free(float* %157)
%159 = load i32, i32* %13
%160 = call float* @mallocf(i32 %159)
%161 = load float**, float*** %10
%162 = getelementptr inbounds float*, float** %161, i64 1
store float* %160, float** %162
%163 = load float*, float** %18
%164 = load float**, float*** %10
%165 = getelementptr inbounds float*, float** %164, i64 1
%166 = load float*, float** %165
%167 = load i32*, i32** %11
%168 = load double, double* %12
%169 = load double, double* %12
%170 = fmul double 4.000000e+00, %169
%171 = fptosi double %170 to i32
%172 = call i32 @GaussianFiltering3D_float(float* %163, float* %166, i32* %167, double %168, i32 %171)
%173 = load float*, float** %18
%174 = call i32 @free(float* %173)
%175 = load i32, i32* %13
%176 = call float* @mallocf(i32 %175)
%177 = load float**, float*** %10
%178 = getelementptr inbounds float*, float** %177, i64 2
store float* %176, float** %178
%179 = load float*, float** %20
%180 = load float**, float*** %10
%181 = getelementptr inbounds float*, float** %180, i64 2
%182 = load float*, float** %181
%183 = load i32*, i32** %11
%184 = load double, double* %12
%185 = load double, double* %12
%186 = fmul double 4.000000e+00, %185
%187 = fptosi double %186 to i32
%188 = call i32 @GaussianFiltering3D_float(float* %179, float* %182, i32* %183, double %184, i32 %187)
%189 = load float*, float** %20
%190 = call i32 @free(float* %189)
%191 = load i32, i32* %13
%192 = call float* @mallocf(i32 %191)
%193 = load float**, float*** %10
%194 = getelementptr inbounds float*, float** %193, i64 3
store float* %192, float** %194
%195 = load float*, float** %16
%196 = load float**, float*** %10
%197 = getelementptr inbounds float*, float** %196, i64 3
%198 = load float*, float** %197
%199 = load i32*, i32** %11
%200 = load double, double* %12
%201 = load double, double* %12
%202 = fmul double 4.000000e+00, %201
%203 = fptosi double %202 to i32
%204 = call i32 @GaussianFiltering3D_float(float* %195, float* %198, i32* %199, double %200, i32 %203)
%205 = load float*, float** %16
%206 = call i32 @free(float* %205)
%207 = load i32, i32* %13
%208 = call float* @mallocf(i32 %207)
%209 = load float**, float*** %10
%210 = getelementptr inbounds float*, float** %209, i64 4
store float* %208, float** %210
%211 = load float*, float** %17
%212 = load float**, float*** %10
%213 = getelementptr inbounds float*, float** %212, i64 4
%214 = load float*, float** %213
%215 = load i32*, i32** %11
%216 = load double, double* %12
%217 = load double, double* %12
%218 = fmul double 4.000000e+00, %217
%219 = fptosi double %218 to i32
%220 = call i32 @GaussianFiltering3D_float(float* %211, float* %214, i32* %215, double %216, i32 %219)
%221 = load float*, float** %17
%222 = call i32 @free(float* %221)
%223 = load i32, i32* %13
%224 = call float* @mallocf(i32 %223)
%225 = load float**, float*** %10
%226 = getelementptr inbounds float*, float** %225, i64 5
store float* %224, float** %226
%227 = load float*, float** %19
%228 = load float**, float*** %10
%229 = getelementptr inbounds float*, float** %228, i64 5
%230 = load float*, float** %229
%231 = load i32*, i32** %11
%232 = load double, double* %12
%233 = load double, double* %12
%234 = fmul double 4.000000e+00, %233
%235 = fptosi double %234 to i32
%236 = call i32 @GaussianFiltering3D_float(float* %227, float* %230, i32* %231, double %232, i32 %235)
%237 = load float*, float** %19
%238 = call i32 @free(float* %237)
ret void
}

LLVM IR (-O2):
define void @StructureTensor3D(float* nocapture %0, float* nocapture %1, float* nocapture %2, float** nocapture writeonly %3, i32* %4, double %5) {
%7 = load i32, i32* %4
%8 = getelementptr inbounds i32, i32* %4, i64 1
%9 = load i32, i32* %8
%10 = mul nsw i32 %9, %7
%11 = getelementptr inbounds i32, i32* %4, i64 2
%12 = load i32, i32* %11
%13 = mul nsw i32 %10, %12
%14 = tail call float* @mallocf(i32 %13)
%15 = tail call float* @mallocf(i32 %13)
%16 = tail call float* @mallocf(i32 %13)
%17 = tail call float* @mallocf(i32 %13)
%18 = tail call float* @mallocf(i32 %13)
%19 = tail call float* @mallocf(i32 %13)
%20 = icmp sgt i32 %13, 0
br i1 %20, label %21, label %51
21: ; preds = %6
%22 = zext i32 %13 to i64
br label %23
23: ; preds = %21, %23
%24 = phi i64 [ 0, %21 ], [ %49, %23 ]
%25 = getelementptr inbounds float, float* %0, i64 %24
%26 = load float, float* %25
%27 = fmul float %26, %26
%28 = getelementptr inbounds float, float* %14, i64 %24
store float %27, float* %28
%29 = load float, float* %25
%30 = getelementptr inbounds float, float* %1, i64 %24
%31 = load float, float* %30
%32 = fmul float %29, %31
%33 = getelementptr inbounds float, float* %17, i64 %24
store float %32, float* %33
%34 = load float, float* %25
%35 = getelementptr inbounds float, float* %2, i64 %24
%36 = load float, float* %35
%37 = fmul float %34, %36
%38 = getelementptr inbounds float, float* %18, i64 %24
store float %37, float* %38
%39 = load float, float* %30
%40 = fmul float %39, %39
%41 = getelementptr inbounds float, float* %15, i64 %24
store float %40, float* %41
%42 = load float, float* %30
%43 = load float, float* %35
%44 = fmul float %42, %43
%45 = getelementptr inbounds float, float* %19, i64 %24
store float %44, float* %45
%46 = load float, float* %35
%47 = fmul float %46, %46
%48 = getelementptr inbounds float, float* %16, i64 %24
store float %47, float* %48
%49 = add nuw nsw i64 %24, 1
%50 = icmp eq i64 %49, %22
br i1 %50, label %51, label %23
51: ; preds = %23, %6
%52 = tail call float* @mallocf(i32 %13)
store float* %52, float** %3
%53 = fmul double %5, 4.000000e+00
%54 = fptosi double %53 to i32
%55 = tail call i32 @GaussianFiltering3D_float(float* %14, float* %52, i32* %4, double %5, i32 %54)
%56 = tail call i32 @free(float* %14)
%57 = tail call float* @mallocf(i32 %13)
%58 = getelementptr inbounds float*, float** %3, i64 1
store float* %57, float** %58
%59 = tail call i32 @GaussianFiltering3D_float(float* %15, float* %57, i32* %4, double %5, i32 %54)
%60 = tail call i32 @free(float* %15)
%61 = tail call float* @mallocf(i32 %13)
%62 = getelementptr inbounds float*, float** %3, i64 2
store float* %61, float** %62
%63 = tail call i32 @GaussianFiltering3D_float(float* %16, float* %61, i32* %4, double %5, i32 %54)
%64 = tail call i32 @free(float* %16)
%65 = tail call float* @mallocf(i32 %13)
%66 = getelementptr inbounds float*, float** %3, i64 3
store float* %65, float** %66
%67 = tail call i32 @GaussianFiltering3D_float(float* %17, float* %65, i32* %4, double %5, i32 %54)
%68 = tail call i32 @free(float* %17)
%69 = tail call float* @mallocf(i32 %13)
%70 = getelementptr inbounds float*, float** %3, i64 4
store float* %69, float** %70
%71 = tail call i32 @GaussianFiltering3D_float(float* %18, float* %69, i32* %4, double %5, i32 %54)
%72 = tail call i32 @free(float* %18)
%73 = tail call float* @mallocf(i32 %13)
%74 = getelementptr inbounds float*, float** %3, i64 5
store float* %73, float** %74
%75 = tail call i32 @GaussianFiltering3D_float(float* %19, float* %73, i32* %4, double %5, i32 %54)
%76 = tail call i32 @free(float* %19)
ret void
}