#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);
void clr(const double *const mis, const int *const np, double *const res) {
  const int n = *np;
  double ms[n], vs[n];
  {
    int i, j;
    for (i = 0; i < n; i++) {
      double acc;
      for (j = 0, acc = 0; j < n; j++)
        acc += mis[n * i + j];
      ms[i] = acc / n;
      for (j = 0, acc = 0; j < n; j++) {
        const double v = mis[n * i + j] - ms[i];
        acc += v * v;
      }
      vs[i] = acc / n;
    }
    for (i = 1; i < n; i++) {
      res[n * i + i] = 0;
      for (j = 0; j < i; j++) {
        const double m = mis[n * i + j];
        double zi = m - ms[i];
        zi = (zi < 0 || vs[i] == 0) ? 0 : zi * zi / vs[i];
        double zj = m - ms[j];
        zj = (zj < 0 || vs[j] == 0) ? 0 : zj * zj / vs[j];
        res[n * i + j] = res[n * j + i] = sqrt(zi * zi + zj * zj);
      }
    }
  }
}


x86 assembly (clang, -O0):
<clr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %rsp,%rcx
mov    %rcx,-0x28(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x78(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x70(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1284 <clr+0x174>
movl   $0x0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11d3 <clr+0xc3>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x3c(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1199 <clr+0x89>
mov    -0x78(%rbp),%rax
movsd  -0x48(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
divsd  %xmm1,%xmm0
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    125b <clr+0x14b>
mov    -0x78(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x1c(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    -0x40(%rbp),%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11fd <clr+0xed>
mov    -0x70(%rbp),%rax
movsd  -0x48(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
divsd  %xmm1,%xmm0
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    117e <clr+0x6e>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1445 <clr+0x335>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x3c(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1432 <clr+0x322>
mov    -0x78(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    -0x1c(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    -0x40(%rbp),%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x60(%rbp),%xmm0
ja     1323 <clr+0x213>
mov    -0x70(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1330 <clr+0x220>
jp     1330 <clr+0x220>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    134c <clr+0x23c>
mov    -0x70(%rbp),%rax
movsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movslq -0x3c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x78(%rbp),%rax
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x58(%rbp),%xmm0
movslq -0x40(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x68(%rbp),%xmm0
ja     139b <clr+0x28b>
mov    -0x70(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    13ab <clr+0x29b>
jp     13ab <clr+0x29b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
jmp    13ca <clr+0x2ba>
mov    -0x70(%rbp),%rax
movsd  -0x68(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x40(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x3c(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    12b7 <clr+0x1a7>
jmp    1437 <clr+0x327>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    128b <clr+0x17b>
mov    -0x28(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<clr>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %rdx,-0x48(%rbp)
mov    (%rsi),%r11d
mov    %rsp,%rcx
lea    0xf(,%r11,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rcx
mov    %rcx,-0x40(%rbp)
mov    %rcx,%rsp
mov    %rsp,%rcx
sub    %rax,%rcx
mov    %rcx,-0x30(%rbp)
mov    %rcx,%rsp
test   %r11d,%r11d
jle    13ff <clr+0x2ef>
cvtsi2sd %r11d,%xmm0
lea    -0x1(%r11),%r8
mov    %r11d,%ecx
and    $0x3,%ecx
mov    %r11d,%edx
and    $0xfffffffc,%edx
lea    0x18(%rdi),%rsi
lea    0x0(,%r11,8),%rax
mov    %rax,-0x38(%rbp)
xor    %r10d,%r10d
mov    %rdi,%r9
jmp    11b5 <clr+0xa5>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
divsd  %xmm0,%xmm2
mov    -0x30(%rbp),%rax
movsd  %xmm2,(%rax,%r10,8)
add    $0x1,%r10
mov    -0x38(%rbp),%rax
add    %rax,%rsi
add    %rax,%r9
cmp    %r11,%r10
je     12af <clr+0x19f>
xorpd  %xmm1,%xmm1
xor    %ebx,%ebx
cmp    $0x3,%r8
jb     11f0 <clr+0xe0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  -0x18(%rsi,%rbx,8),%xmm1
addsd  -0x10(%rsi,%rbx,8),%xmm1
addsd  -0x8(%rsi,%rbx,8),%xmm1
addsd  (%rsi,%rbx,8),%xmm1
add    $0x4,%rbx
cmp    %rbx,%rdx
jne    11d0 <clr+0xc0>
test   %rcx,%rcx
je     120e <clr+0xfe>
lea    (%r9,%rbx,8),%rax
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
addsd  (%rax,%rbx,8),%xmm1
add    $0x1,%rbx
cmp    %rbx,%rcx
jne    1200 <clr+0xf0>
divsd  %xmm0,%xmm1
mov    -0x40(%rbp),%rax
movsd  %xmm1,(%rax,%r10,8)
xorpd  %xmm2,%xmm2
xor    %ebx,%ebx
cmp    $0x3,%r8
jb     1280 <clr+0x170>
nopl   0x0(%rax,%rax,1)
movsd  -0x18(%rsi,%rbx,8),%xmm3
movsd  -0x10(%rsi,%rbx,8),%xmm4
subsd  %xmm1,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm2,%xmm3
subsd  %xmm1,%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm3,%xmm4
movsd  -0x8(%rsi,%rbx,8),%xmm3
subsd  %xmm1,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm4,%xmm3
movsd  (%rsi,%rbx,8),%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm3,%xmm2
add    $0x4,%rbx
cmp    %rbx,%rdx
jne    1230 <clr+0x120>
test   $0x3,%r11b
je     1190 <clr+0x80>
lea    (%r9,%rbx,8),%rbx
xor    %eax,%eax
movsd  (%rbx,%rax,8),%xmm3
subsd  %xmm1,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm3,%xmm2
add    $0x1,%rax
cmp    %rax,%rcx
jne    1290 <clr+0x180>
jmp    1190 <clr+0x80>
cmp    $0x2,%r11d
jl     13ff <clr+0x2ef>
mov    -0x48(%rbp),%rax
lea    0x8(%rax),%r12
lea    (%rdi,%r11,8),%r14
lea    (%rax,%r11,8),%rdi
mov    $0x1,%r15d
mov    $0x8,%r13d
xorpd  %xmm5,%xmm5
mov    %r11,-0x50(%rbp)
jmp    1307 <clr+0x1f7>
nop
add    $0x1,%r15
mov    -0x58(%rbp),%r12
add    $0x8,%r12
mov    -0x38(%rbp),%rax
add    %rax,%r14
add    %rax,%rdi
add    $0x8,%r13
mov    -0x50(%rbp),%r11
cmp    %r11,%r15
je     13ff <clr+0x2ef>
mov    %r15d,%eax
imul   %r11d,%eax
cltq
add    %r15,%rax
mov    -0x48(%rbp),%rcx
movq   $0x0,(%rcx,%rax,8)
mov    -0x40(%rbp),%rax
movsd  (%rax,%r15,8),%xmm6
xor    %ebx,%ebx
mov    %r12,-0x58(%rbp)
mov    %rdi,-0x68(%rbp)
movsd  %xmm6,-0x60(%rbp)
jmp    136e <clr+0x25e>
nopw   0x0(%rax,%rax,1)
call   1030 <sqrt@plt>
movsd  -0x60(%rbp),%xmm6
xorpd  %xmm5,%xmm5
mov    -0x68(%rbp),%rdi
movsd  %xmm0,(%r12)
movsd  %xmm0,(%rdi,%rbx,1)
add    -0x38(%rbp),%r12
add    $0x8,%rbx
cmp    %rbx,%r13
je     12e0 <clr+0x1d0>
movsd  (%r14,%rbx,1),%xmm1
movapd %xmm1,%xmm3
subsd  %xmm6,%xmm3
xorpd  %xmm2,%xmm2
ucomisd %xmm3,%xmm2
xorpd  %xmm0,%xmm0
ja     13b0 <clr+0x2a0>
mov    -0x30(%rbp),%rax
movsd  (%rax,%r15,8),%xmm4
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm4
jne    13a0 <clr+0x290>
jnp    13b0 <clr+0x2a0>
mulsd  %xmm3,%xmm3
divsd  %xmm4,%xmm3
movapd %xmm3,%xmm0
nopl   0x0(%rax)
mov    -0x40(%rbp),%rax
subsd  (%rax,%rbx,1),%xmm1
ucomisd %xmm1,%xmm2
ja     13e0 <clr+0x2d0>
mov    -0x30(%rbp),%rax
movsd  (%rax,%rbx,1),%xmm3
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm3
jne    13d4 <clr+0x2c4>
jnp    13e0 <clr+0x2d0>
mulsd  %xmm1,%xmm1
divsd  %xmm3,%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm2,%xmm2
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
ucomisd %xmm5,%xmm0
jb     1340 <clr+0x230>
sqrtsd %xmm0,%xmm0
jmp    1352 <clr+0x242>
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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