#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 fabs(double);
double sqrt(double);
double dnrm2_(int *n, double *x, int *incx) {
  long int ix, nn, iincx;
  double norm, scale, absxi, ssq, temp;
  nn = *n;
  iincx = *incx;
  if (nn > 0 && iincx > 0) {
    if (nn == 1) {
      norm = fabs(x[0]);
    } else {
      scale = 0.0;
      ssq = 1.0;
      for (ix = (nn - 1) * iincx; ix >= 0; ix -= iincx) {
        if (x[ix] != 0.0) {
          absxi = fabs(x[ix]);
          if (scale < absxi) {
            temp = scale / absxi;
            ssq = ssq * (temp * temp) + 1.0;
            scale = absxi;
          } else {
            temp = absxi / scale;
            ssq += temp * temp;
          }
        }
      }
      norm = scale * sqrt(ssq);
    }
  } else
    norm = 0.0;
  return norm;
}


x86 assembly (clang, -O0):
<dnrm2_>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movslq (%rax),%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movslq (%rax),%rax
mov    %rax,-0x30(%rbp)
cmpq   $0x0,-0x28(%rbp)
jle    12b4 <dnrm2_+0x1a4>
cmpq   $0x0,-0x30(%rbp)
jle    12b4 <dnrm2_+0x1a4>
cmpq   $0x1,-0x28(%rbp)
jne    1178 <dnrm2_+0x68>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movaps 0xea6(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    12af <dnrm2_+0x19f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe78(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x28(%rbp),%rax
sub    $0x1,%rax
imul   -0x30(%rbp),%rax
mov    %rax,-0x20(%rbp)
cmpq   $0x0,-0x20(%rbp)
jl     128a <dnrm2_+0x17a>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11ce <dnrm2_+0xbe>
jp     11ce <dnrm2_+0xbe>
jmp    1271 <dnrm2_+0x161>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe2e(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1241 <dnrm2_+0x131>
movsd  -0x40(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xdd7(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    126c <dnrm2_+0x15c>
movsd  -0x48(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1271 <dnrm2_+0x161>
jmp    1276 <dnrm2_+0x166>
mov    -0x30(%rbp),%rcx
mov    -0x20(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0x20(%rbp)
jmp    119e <dnrm2_+0x8e>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    12bc <dnrm2_+0x1ac>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dnrm2_>:
push   %rax
movslq (%rdi),%rax
xorpd  %xmm6,%xmm6
test   %rax,%rax
jle    11d9 <dnrm2_+0xc9>
mov    (%rdx),%ecx
test   %ecx,%ecx
jle    11d9 <dnrm2_+0xc9>
cmp    $0x1,%eax
jne    113d <dnrm2_+0x2d>
movsd  (%rsi),%xmm0
andps  0xed5(%rip),%xmm0
pop    %rax
ret
movslq %ecx,%rcx
add    $0xffffffffffffffff,%rax
imul   %rcx,%rax
xorpd  %xmm1,%xmm1
test   %rax,%rax
js     11ae <dnrm2_+0x9e>
xorpd  %xmm2,%xmm2
movsd  0xea3(%rip),%xmm3
movapd 0xeab(%rip),%xmm4
movapd %xmm3,%xmm0
xorpd  %xmm6,%xmm6
jmp    1181 <dnrm2_+0x71>
nop
divsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm5
addsd  %xmm5,%xmm0
sub    %rcx,%rax
js     11ba <dnrm2_+0xaa>
movsd  (%rsi,%rax,8),%xmm5
ucomisd %xmm2,%xmm5
jne    118e <dnrm2_+0x7e>
jnp    117c <dnrm2_+0x6c>
andpd  %xmm4,%xmm5
ucomisd %xmm6,%xmm5
jbe    1170 <dnrm2_+0x60>
divsd  %xmm5,%xmm6
mulsd  %xmm6,%xmm6
mulsd  %xmm6,%xmm0
addsd  %xmm3,%xmm0
movapd %xmm5,%xmm6
jmp    117c <dnrm2_+0x6c>
movsd  0xe4a(%rip),%xmm0
xorpd  %xmm6,%xmm6
ucomisd %xmm1,%xmm0
jb     11c6 <dnrm2_+0xb6>
sqrtsd %xmm0,%xmm0
jmp    11d5 <dnrm2_+0xc5>
movsd  %xmm6,(%rsp)
call   1030 <sqrt@plt>
movsd  (%rsp),%xmm6
mulsd  %xmm0,%xmm6
movapd %xmm6,%xmm0
pop    %rax
ret

LLVM IR (-O0):
define double @dnrm2_(i32* %0, double* %1, i32* %2) {
%4 = alloca i32*
%5 = alloca double*
%6 = alloca i32*
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
store i32* %0, i32** %4
store double* %1, double** %5
store i32* %2, i32** %6
%15 = load i32*, i32** %4
%16 = load i32, i32* %15
%17 = sext i32 %16 to i64
store i64 %17, i64* %8
%18 = load i32*, i32** %6
%19 = load i32, i32* %18
%20 = sext i32 %19 to i64
store i64 %20, i64* %9
%21 = load i64, i64* %8
%22 = icmp sgt i64 %21, 0
br i1 %22, label %23, label %87
23: ; preds = %3
%24 = load i64, i64* %9
%25 = icmp sgt i64 %24, 0
br i1 %25, label %26, label %87
26: ; preds = %23
%27 = load i64, i64* %8
%28 = icmp eq i64 %27, 1
br i1 %28, label %29, label %34
29: ; preds = %26
%30 = load double*, double** %5
%31 = getelementptr inbounds double, double* %30, i64 0
%32 = load double, double* %31
%33 = call double @llvm.fabs.f64(double %32)
store double %33, double* %10
br label %86
34: ; preds = %26
store double 0.000000e+00, double* %11
store double 1.000000e+00, double* %13
%35 = load i64, i64* %8
%36 = sub nsw i64 %35, 1
%37 = load i64, i64* %9
%38 = mul nsw i64 %36, %37
store i64 %38, i64* %7
br label %39
39: ; preds = %77, %34
%40 = load i64, i64* %7
%41 = icmp sge i64 %40, 0
br i1 %41, label %42, label %81
42: ; preds = %39
%43 = load double*, double** %5
%44 = load i64, i64* %7
%45 = getelementptr inbounds double, double* %43, i64 %44
%46 = load double, double* %45
%47 = fcmp une double %46, 0.000000e+00
br i1 %47, label %48, label %76
48: ; preds = %42
%49 = load double*, double** %5
%50 = load i64, i64* %7
%51 = getelementptr inbounds double, double* %49, i64 %50
%52 = load double, double* %51
%53 = call double @llvm.fabs.f64(double %52)
store double %53, double* %12
%54 = load double, double* %11
%55 = load double, double* %12
%56 = fcmp olt double %54, %55
br i1 %56, label %57, label %67
57: ; preds = %48
%58 = load double, double* %11
%59 = load double, double* %12
%60 = fdiv double %58, %59
store double %60, double* %14
%61 = load double, double* %13
%62 = load double, double* %14
%63 = load double, double* %14
%64 = fmul double %62, %63
%65 = call double @llvm.fmuladd.f64(double %61, double %64, double 1.000000e+00)
store double %65, double* %13
%66 = load double, double* %12
store double %66, double* %11
br label %75
67: ; preds = %48
%68 = load double, double* %12
%69 = load double, double* %11
%70 = fdiv double %68, %69
store double %70, double* %14
%71 = load double, double* %14
%72 = load double, double* %14
%73 = load double, double* %13
%74 = call double @llvm.fmuladd.f64(double %71, double %72, double %73)
store double %74, double* %13
br label %75
75: ; preds = %67, %57
br label %76
76: ; preds = %75, %42
br label %77
77: ; preds = %76
%78 = load i64, i64* %9
%79 = load i64, i64* %7
%80 = sub nsw i64 %79, %78
store i64 %80, i64* %7
br label %39
81: ; preds = %39
%82 = load double, double* %11
%83 = load double, double* %13
%84 = call double @sqrt(double %83)
%85 = fmul double %82, %84
store double %85, double* %10
br label %86
86: ; preds = %81, %29
br label %88
87: ; preds = %23, %3
store double 0.000000e+00, double* %10
br label %88
88: ; preds = %87, %86
%89 = load double, double* %10
ret double %89
}

LLVM IR (-O2):
define double @dnrm2_(i32* nocapture %0, double* nocapture %1, i32* nocapture %2) {
%4 = load i32, i32* %0
%5 = sext i32 %4 to i64
%6 = load i32, i32* %2
%7 = sext i32 %6 to i64
%8 = icmp sgt i32 %4, 0
%9 = icmp sgt i32 %6, 0
%10 = select i1 %8, i1 %9, i1 false
br i1 %10, label %11, label %47
11: ; preds = %3
%12 = icmp eq i32 %4, 1
br i1 %12, label %13, label %16
13: ; preds = %11
%14 = load double, double* %1
%15 = tail call double @llvm.fabs.f64(double %14)
br label %47
16: ; preds = %11
%17 = add nsw i64 %5, -1
%18 = mul nsw i64 %17, %7
%19 = icmp sgt i64 %18, -1
br i1 %19, label %20, label %42
20: ; preds = %16, %37
%21 = phi double [ %39, %37 ], [ 1.000000e+00, %16 ]
%22 = phi i64 [ %40, %37 ], [ %18, %16 ]
%23 = phi double [ %38, %37 ], [ 0.000000e+00, %16 ]
%24 = getelementptr inbounds double, double* %1, i64 %22
%25 = load double, double* %24
%26 = fcmp une double %25, 0.000000e+00
br i1 %26, label %27, label %37
27: ; preds = %20
%28 = tail call double @llvm.fabs.f64(double %25)
%29 = fcmp olt double %23, %28
br i1 %29, label %30, label %34
30: ; preds = %27
%31 = fdiv double %23, %28
%32 = fmul double %31, %31
%33 = tail call double @llvm.fmuladd.f64(double %21, double %32, double 1.000000e+00)
br label %37
34: ; preds = %27
%35 = fdiv double %28, %23
%36 = tail call double @llvm.fmuladd.f64(double %35, double %35, double %21)
br label %37
37: ; preds = %20, %34, %30
%38 = phi double [ %28, %30 ], [ %23, %34 ], [ %23, %20 ]
%39 = phi double [ %33, %30 ], [ %36, %34 ], [ %21, %20 ]
%40 = sub nsw i64 %22, %7
%41 = icmp sgt i64 %40, -1
br i1 %41, label %20, label %42
42: ; preds = %37, %16
%43 = phi double [ 0.000000e+00, %16 ], [ %38, %37 ]
%44 = phi double [ 1.000000e+00, %16 ], [ %39, %37 ]
%45 = tail call double @sqrt(double %44)
%46 = fmul double %43, %45
br label %47
47: ; preds = %3, %13, %42
%48 = phi double [ %15, %13 ], [ %46, %42 ], [ 0.000000e+00, %3 ]
ret double %48
}