#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 */
double sqrt(double);
double calc_gz_norm(double *phi, unsigned long nvox) {
  unsigned long m = 0, t = 0;
  double varsum = 0.0;
  for (m = 0; m < nvox; m++) {
    if (*(phi + m) > 0) {
      varsum += (*(phi + m)) * (*(phi + m));
      ++t;
    }
  }
  if (t > 0)
    return sqrt((varsum) / ((double)t));
  else
    return 0.0;
}


x86 assembly (clang, -O0):
<calc_gz_norm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movq   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x18(%rbp),%rax
jae    11b6 <calc_gz_norm+0xa6>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    11a0 <calc_gz_norm+0x90>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    11a5 <calc_gz_norm+0x95>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1140 <calc_gz_norm+0x30>
cmpq   $0x0,-0x28(%rbp)
jbe    1200 <calc_gz_norm+0xf0>
movsd  -0x30(%rbp),%xmm0
movaps 0xe33(%rip),%xmm2
movq   -0x28(%rbp),%xmm1
punpckldq %xmm2,%xmm1
movapd 0xe32(%rip),%xmm2
subpd  %xmm2,%xmm1
movaps %xmm1,%xmm2
unpckhpd %xmm1,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    1208 <calc_gz_norm+0xf8>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_gz_norm>:
test   %rsi,%rsi
je     1181 <calc_gz_norm+0x71>
cmp    $0x1,%rsi
jne    1186 <calc_gz_norm+0x76>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
xor    %ecx,%ecx
test   $0x1,%sil
je     1148 <calc_gz_norm+0x38>
movsd  (%rdi,%rcx,8),%xmm1
xorpd  %xmm2,%xmm2
xor    %ecx,%ecx
ucomisd %xmm2,%xmm1
seta   %cl
jbe    1145 <calc_gz_norm+0x35>
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    %rcx,%rax
test   %rax,%rax
je     1181 <calc_gz_norm+0x71>
movq   %rax,%xmm1
punpckldq 0xea6(%rip),%xmm1
subpd  0xeae(%rip),%xmm1
movapd %xmm1,%xmm2
unpckhpd %xmm1,%xmm2
addsd  %xmm1,%xmm2
divsd  %xmm2,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     11ee <calc_gz_norm+0xde>
sqrtsd %xmm0,%xmm0
ret
xorpd  %xmm0,%xmm0
ret
mov    %rsi,%r8
and    $0xfffffffffffffffe,%r8
xorpd  %xmm1,%xmm1
xor    %eax,%eax
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    11b3 <calc_gz_norm+0xa3>
nopl   0x0(%rax,%rax,1)
mov    %r9b,%dl
add    %rdx,%rax
add    $0x2,%rcx
cmp    %rcx,%r8
je     1123 <calc_gz_norm+0x13>
movsd  (%rdi,%rcx,8),%xmm2
xor    %edx,%edx
ucomisd %xmm1,%xmm2
seta   %r9b
jbe    11cc <calc_gz_norm+0xbc>
mulsd  %xmm2,%xmm2
addsd  %xmm2,%xmm0
movsd  0x8(%rdi,%rcx,8),%xmm2
mov    %r9b,%dl
add    %rdx,%rax
xor    %edx,%edx
ucomisd %xmm1,%xmm2
seta   %r9b
jbe    11a0 <calc_gz_norm+0x90>
mulsd  %xmm2,%xmm2
addsd  %xmm2,%xmm0
jmp    11a0 <calc_gz_norm+0x90>
jmp    1030 <sqrt@plt>

LLVM IR (-O0):
define double @calc_gz_norm(double* %0, i64 %1) {
%3 = alloca double
%4 = alloca double*
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca double
store double* %0, double** %4
store i64 %1, i64* %5
store i64 0, i64* %6
store i64 0, i64* %7
store double 0.000000e+00, double* %8
store i64 0, i64* %6
br label %9
9: ; preds = %33, %2
%10 = load i64, i64* %6
%11 = load i64, i64* %5
%12 = icmp ult i64 %10, %11
br i1 %12, label %13, label %36
13: ; preds = %9
%14 = load double*, double** %4
%15 = load i64, i64* %6
%16 = getelementptr inbounds double, double* %14, i64 %15
%17 = load double, double* %16
%18 = fcmp ogt double %17, 0.000000e+00
br i1 %18, label %19, label %32
19: ; preds = %13
%20 = load double*, double** %4
%21 = load i64, i64* %6
%22 = getelementptr inbounds double, double* %20, i64 %21
%23 = load double, double* %22
%24 = load double*, double** %4
%25 = load i64, i64* %6
%26 = getelementptr inbounds double, double* %24, i64 %25
%27 = load double, double* %26
%28 = load double, double* %8
%29 = call double @llvm.fmuladd.f64(double %23, double %27, double %28)
store double %29, double* %8
%30 = load i64, i64* %7
%31 = add i64 %30, 1
store i64 %31, i64* %7
br label %32
32: ; preds = %19, %13
br label %33
33: ; preds = %32
%34 = load i64, i64* %6
%35 = add i64 %34, 1
store i64 %35, i64* %6
br label %9
36: ; preds = %9
%37 = load i64, i64* %7
%38 = icmp ugt i64 %37, 0
br i1 %38, label %39, label %45
39: ; preds = %36
%40 = load double, double* %8
%41 = load i64, i64* %7
%42 = uitofp i64 %41 to double
%43 = fdiv double %40, %42
%44 = call double @sqrt(double %43)
store double %44, double* %3
br label %46
45: ; preds = %36
store double 0.000000e+00, double* %3
br label %46
46: ; preds = %45, %39
%47 = load double, double* %3
ret double %47
}

LLVM IR (-O2):
define double @calc_gz_norm(double* nocapture %0, i64 %1) {
%3 = icmp eq i64 %1, 0
br i1 %3, label %55, label %4
4: ; preds = %2
%5 = and i64 %1, 1
%6 = icmp eq i64 %1, 1
br i1 %6, label %32, label %7
7: ; preds = %4
%8 = and i64 %1, -2
br label %9
9: ; preds = %9, %7
%10 = phi double [ 0.000000e+00, %7 ], [ %28, %9 ]
%11 = phi i64 [ 0, %7 ], [ %27, %9 ]
%12 = phi i64 [ 0, %7 ], [ %29, %9 ]
%13 = phi i64 [ 0, %7 ], [ %30, %9 ]
%14 = getelementptr inbounds double, double* %0, i64 %12
%15 = load double, double* %14
%16 = fcmp ogt double %15, 0.000000e+00
%17 = tail call double @llvm.fmuladd.f64(double %15, double %15, double %10)
%18 = zext i1 %16 to i64
%19 = add i64 %11, %18
%20 = select i1 %16, double %17, double %10
%21 = or i64 %12, 1
%22 = getelementptr inbounds double, double* %0, i64 %21
%23 = load double, double* %22
%24 = fcmp ogt double %23, 0.000000e+00
%25 = tail call double @llvm.fmuladd.f64(double %23, double %23, double %20)
%26 = zext i1 %24 to i64
%27 = add i64 %19, %26
%28 = select i1 %24, double %25, double %20
%29 = add nuw i64 %12, 2
%30 = add i64 %13, 2
%31 = icmp eq i64 %30, %8
br i1 %31, label %32, label %9
32: ; preds = %9, %4
%33 = phi i64 [ undef, %4 ], [ %27, %9 ]
%34 = phi double [ undef, %4 ], [ %28, %9 ]
%35 = phi double [ 0.000000e+00, %4 ], [ %28, %9 ]
%36 = phi i64 [ 0, %4 ], [ %27, %9 ]
%37 = phi i64 [ 0, %4 ], [ %29, %9 ]
%38 = icmp eq i64 %5, 0
br i1 %38, label %47, label %39
39: ; preds = %32
%40 = getelementptr inbounds double, double* %0, i64 %37
%41 = load double, double* %40
%42 = fcmp ogt double %41, 0.000000e+00
%43 = tail call double @llvm.fmuladd.f64(double %41, double %41, double %35)
%44 = zext i1 %42 to i64
%45 = add i64 %36, %44
%46 = select i1 %42, double %43, double %35
br label %47
47: ; preds = %32, %39
%48 = phi i64 [ %33, %32 ], [ %45, %39 ]
%49 = phi double [ %34, %32 ], [ %46, %39 ]
%50 = icmp eq i64 %48, 0
br i1 %50, label %55, label %51
51: ; preds = %47
%52 = uitofp i64 %48 to double
%53 = fdiv double %49, %52
%54 = tail call double @sqrt(double %53)
br label %55
55: ; preds = %2, %47, %51
%56 = phi double [ %54, %51 ], [ 0.000000e+00, %47 ], [ 0.000000e+00, %2 ]
ret double %56
}