void svdvar(double *v, double *w, int n, double *cov) {
  static double wti[(3 * 11)];
  double sum;
  int i, j, k;
  for (i = 0; i < n; i++)
    wti[i] = w[i] ? 1.0 / (w[i] * w[i]) : 0.0;
  for (i = 0; i < n; i++)
    for (j = 0; j <= i; j++) {
      for (sum = 0.0, k = 0; k < n; k++)
        sum += v[k * n + i] * v[k * n + j] * wti[k];
      cov[j * n + i] = cov[i * n + j] = sum;
    }
  return;
}


x86 assembly (clang, -O0):
<svdvar>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11ab <svdvar+0xab>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    114b <svdvar+0x4b>
jp     114b <svdvar+0x4b>
jmp    117b <svdvar+0x7b>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movsd  0xe93(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1188 <svdvar+0x88>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1188 <svdvar+0x88>
movsd  -0x40(%rbp),%xmm0
movslq -0x2c(%rbp),%rcx
lea    0x2e98(%rip),%rax
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    111a <svdvar+0x1a>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    129a <svdvar+0x19a>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1287 <svdvar+0x187>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1248 <svdvar+0x148>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movslq -0x34(%rbp),%rcx
lea    0x2e0d(%rip),%rax
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11e0 <svdvar+0xe0>
movsd  -0x28(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11c5 <svdvar+0xc5>
jmp    128c <svdvar+0x18c>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b2 <svdvar+0xb2>
pop    %rbp
ret

x86 assembly (clang, -O2):
<svdvar>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    1307 <svdvar+0x207>
mov    %edx,%r9d
mov    %r9,%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x8(%rsp)
je     1192 <svdvar+0x92>
mov    %r9d,%ebx
and    $0xfffffffe,%ebx
xor    %eax,%eax
movsd  0xecd(%rip),%xmm0
lea    0x2ef6(%rip),%rbp
jmp    114f <svdvar+0x4f>
nopl   0x0(%rax)
movsd  %xmm1,0x8(%rbp,%rax,8)
add    $0x2,%rax
cmp    %rax,%rbx
je     1194 <svdvar+0x94>
movsd  (%rsi,%rax,8),%xmm2
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm2
xorpd  %xmm3,%xmm3
jne    1164 <svdvar+0x64>
jnp    1170 <svdvar+0x70>
mulsd  %xmm2,%xmm2
movapd %xmm0,%xmm3
divsd  %xmm2,%xmm3
movsd  %xmm3,0x0(%rbp,%rax,8)
movsd  0x8(%rsi,%rax,8),%xmm2
ucomisd %xmm1,%xmm2
jne    1184 <svdvar+0x84>
jnp    1140 <svdvar+0x40>
mulsd  %xmm2,%xmm2
movapd %xmm0,%xmm1
divsd  %xmm2,%xmm1
jmp    1140 <svdvar+0x40>
xor    %eax,%eax
test   $0x1,%r9b
je     11c7 <svdvar+0xc7>
movsd  (%rsi,%rax,8),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    11ab <svdvar+0xab>
jnp    11bb <svdvar+0xbb>
mulsd  %xmm1,%xmm1
movsd  0xe49(%rip),%xmm0
divsd  %xmm1,%xmm0
lea    0x2e6e(%rip),%rsi
movsd  %xmm0,(%rsi,%rax,8)
test   %edx,%edx
jle    1307 <svdvar+0x207>
mov    %r9d,%edx
and    $0xfffffffe,%edx
lea    (%r9,%r9,1),%rsi
mov    $0x1,%eax
mov    %rax,-0x10(%rsp)
xor    %r11d,%r11d
lea    0x2e43(%rip),%rax
mov    %r9,-0x18(%rsp)
jmp    1219 <svdvar+0x119>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r11
addq   $0x1,-0x10(%rsp)
addq   $0x1,-0x18(%rsp)
cmp    %r9,%r11
je     1307 <svdvar+0x207>
mov    %r11,%r15
imul   %r9,%r15
mov    %r9,%r10
xor    %r8d,%r8d
jmp    1257 <svdvar+0x157>
nopl   0x0(%rax,%rax,1)
lea    (%r8,%r15,1),%rbx
movsd  %xmm0,(%rcx,%rbx,8)
mov    %r8,%rbx
imul   %r9,%rbx
add    %r11,%rbx
movsd  %xmm0,(%rcx,%rbx,8)
add    $0x1,%r8
add    $0x1,%r10
cmp    -0x10(%rsp),%r8
je     1200 <svdvar+0x100>
cmpq   $0x0,-0x8(%rsp)
je     12d0 <svdvar+0x1d0>
xorpd  %xmm0,%xmm0
mov    %r11,%rbp
mov    %r8,%r13
mov    -0x18(%rsp),%r14
mov    %r10,%r12
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rdi,%rbp,8),%xmm1
mulsd  (%rdi,%r13,8),%xmm1
mulsd  (%rax,%rbx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  (%rdi,%r14,8),%xmm0
mulsd  (%rdi,%r12,8),%xmm0
mulsd  0x8(%rax,%rbx,8),%xmm0
addsd  %xmm1,%xmm0
add    $0x2,%rbx
add    %rsi,%r12
add    %rsi,%r14
add    %rsi,%r13
add    %rsi,%rbp
cmp    %rbx,%rdx
jne    1280 <svdvar+0x180>
test   $0x1,%r9b
je     1230 <svdvar+0x130>
jmp    12e0 <svdvar+0x1e0>
nopl   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
test   $0x1,%r9b
je     1230 <svdvar+0x130>
mov    %rbx,%r14
imul   %r9,%r14
lea    (%r14,%r11,1),%rbp
movsd  (%rdi,%rbp,8),%xmm1
add    %r8,%r14
mulsd  (%rdi,%r14,8),%xmm1
mulsd  (%rax,%rbx,8),%xmm1
addsd  %xmm1,%xmm0
jmp    1230 <svdvar+0x130>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@svdvar.wti = internal global [33 x double] zeroinitializer
define void @svdvar(double* %0, double* %1, i32 %2, double* %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca double*
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store double* %3, double** %8
store i32 0, i32* %10
br label %13
13: ; preds = %43, %4
%14 = load i32, i32* %10
%15 = load i32, i32* %7
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %46
17: ; preds = %13
%18 = load double*, double** %6
%19 = load i32, i32* %10
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double, double* %18, i64 %20
%22 = load double, double* %21
%23 = fcmp une double %22, 0.000000e+00
br i1 %23, label %24, label %37
24: ; preds = %17
%25 = load double*, double** %6
%26 = load i32, i32* %10
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = load double*, double** %6
%31 = load i32, i32* %10
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
%35 = fmul double %29, %34
%36 = fdiv double 1.000000e+00, %35
br label %38
37: ; preds = %17
br label %38
38: ; preds = %37, %24
%39 = phi double [ %36, %24 ], [ 0.000000e+00, %37 ]
%40 = load i32, i32* %10
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %41
store double %39, double* %42
br label %43
43: ; preds = %38
%44 = load i32, i32* %10
%45 = add nsw i32 %44, 1
store i32 %45, i32* %10
br label %13
46: ; preds = %13
store i32 0, i32* %10
br label %47
47: ; preds = %112, %46
%48 = load i32, i32* %10
%49 = load i32, i32* %7
%50 = icmp slt i32 %48, %49
br i1 %50, label %51, label %115
51: ; preds = %47
store i32 0, i32* %11
br label %52
52: ; preds = %108, %51
%53 = load i32, i32* %11
%54 = load i32, i32* %10
%55 = icmp sle i32 %53, %54
br i1 %55, label %56, label %111
56: ; preds = %52
store double 0.000000e+00, double* %9
store i32 0, i32* %12
br label %57
57: ; preds = %87, %56
%58 = load i32, i32* %12
%59 = load i32, i32* %7
%60 = icmp slt i32 %58, %59
br i1 %60, label %61, label %90
61: ; preds = %57
%62 = load double*, double** %5
%63 = load i32, i32* %12
%64 = load i32, i32* %7
%65 = mul nsw i32 %63, %64
%66 = load i32, i32* %10
%67 = add nsw i32 %65, %66
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %62, i64 %68
%70 = load double, double* %69
%71 = load double*, double** %5
%72 = load i32, i32* %12
%73 = load i32, i32* %7
%74 = mul nsw i32 %72, %73
%75 = load i32, i32* %11
%76 = add nsw i32 %74, %75
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %71, i64 %77
%79 = load double, double* %78
%80 = fmul double %70, %79
%81 = load i32, i32* %12
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %82
%84 = load double, double* %83
%85 = load double, double* %9
%86 = call double @llvm.fmuladd.f64(double %80, double %84, double %85)
store double %86, double* %9
br label %87
87: ; preds = %61
%88 = load i32, i32* %12
%89 = add nsw i32 %88, 1
store i32 %89, i32* %12
br label %57
90: ; preds = %57
%91 = load double, double* %9
%92 = load double*, double** %8
%93 = load i32, i32* %10
%94 = load i32, i32* %7
%95 = mul nsw i32 %93, %94
%96 = load i32, i32* %11
%97 = add nsw i32 %95, %96
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %92, i64 %98
store double %91, double* %99
%100 = load double*, double** %8
%101 = load i32, i32* %11
%102 = load i32, i32* %7
%103 = mul nsw i32 %101, %102
%104 = load i32, i32* %10
%105 = add nsw i32 %103, %104
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds double, double* %100, i64 %106
store double %91, double* %107
br label %108
108: ; preds = %90
%109 = load i32, i32* %11
%110 = add nsw i32 %109, 1
store i32 %110, i32* %11
br label %52
111: ; preds = %52
br label %112
112: ; preds = %111
%113 = load i32, i32* %10
%114 = add nsw i32 %113, 1
store i32 %114, i32* %10
br label %47
115: ; preds = %47
ret void
}

LLVM IR (-O2):
@svdvar.wti = internal unnamed_addr global [33 x double] zeroinitializer
define void @svdvar(double* nocapture %0, double* nocapture %1, i32 %2, double* nocapture writeonly %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %124
6: ; preds = %4
%7 = zext i32 %2 to i64
%8 = add nsw i64 %7, -1
%9 = and i64 %7, 1
%10 = icmp eq i64 %8, 0
br i1 %10, label %13, label %11
11: ; preds = %6
%12 = and i64 %7, 4294967294
br label %34
13: ; preds = %53, %6
%14 = phi i64 [ 0, %6 ], [ %56, %53 ]
%15 = icmp eq i64 %9, 0
br i1 %15, label %26, label %16
16: ; preds = %13
%17 = getelementptr inbounds double, double* %1, i64 %14
%18 = load double, double* %17
%19 = fcmp une double %18, 0.000000e+00
br i1 %19, label %20, label %23
20: ; preds = %16
%21 = fmul double %18, %18
%22 = fdiv double 1.000000e+00, %21
br label %23
23: ; preds = %20, %16
%24 = phi double [ %22, %20 ], [ 0.000000e+00, %16 ]
%25 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %14
store double %24, double* %25
br label %26
26: ; preds = %13, %23
br i1 %5, label %27, label %124
27: ; preds = %26
%28 = zext i32 %2 to i64
%29 = zext i32 %2 to i64
%30 = and i64 %28, 1
%31 = icmp eq i64 %8, 0
%32 = and i64 %28, 4294967294
%33 = icmp eq i64 %30, 0
br label %59
34: ; preds = %53, %11
%35 = phi i64 [ 0, %11 ], [ %56, %53 ]
%36 = phi i64 [ 0, %11 ], [ %57, %53 ]
%37 = getelementptr inbounds double, double* %1, i64 %35
%38 = load double, double* %37
%39 = fcmp une double %38, 0.000000e+00
br i1 %39, label %40, label %43
40: ; preds = %34
%41 = fmul double %38, %38
%42 = fdiv double 1.000000e+00, %41
br label %43
43: ; preds = %34, %40
%44 = phi double [ %42, %40 ], [ 0.000000e+00, %34 ]
%45 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %35
store double %44, double* %45
%46 = or i64 %35, 1
%47 = getelementptr inbounds double, double* %1, i64 %46
%48 = load double, double* %47
%49 = fcmp une double %48, 0.000000e+00
br i1 %49, label %50, label %53
50: ; preds = %43
%51 = fmul double %48, %48
%52 = fdiv double 1.000000e+00, %51
br label %53
53: ; preds = %50, %43
%54 = phi double [ %52, %50 ], [ 0.000000e+00, %43 ]
%55 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %46
store double %54, double* %55
%56 = add nuw nsw i64 %35, 2
%57 = add i64 %36, 2
%58 = icmp eq i64 %57, %12
br i1 %58, label %13, label %34
59: ; preds = %27, %120
%60 = phi i64 [ 0, %27 ], [ %121, %120 ]
%61 = phi i64 [ 1, %27 ], [ %122, %120 ]
%62 = mul nsw i64 %60, %28
br label %63
63: ; preds = %111, %59
%64 = phi i64 [ 0, %59 ], [ %118, %111 ]
br i1 %31, label %95, label %65
65: ; preds = %63, %65
%66 = phi i64 [ %92, %65 ], [ 0, %63 ]
%67 = phi double [ %91, %65 ], [ 0.000000e+00, %63 ]
%68 = phi i64 [ %93, %65 ], [ 0, %63 ]
%69 = mul nsw i64 %66, %28
%70 = add nuw nsw i64 %69, %60
%71 = getelementptr inbounds double, double* %0, i64 %70
%72 = load double, double* %71
%73 = add nuw nsw i64 %69, %64
%74 = getelementptr inbounds double, double* %0, i64 %73
%75 = load double, double* %74
%76 = fmul double %72, %75
%77 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %66
%78 = load double, double* %77
%79 = tail call double @llvm.fmuladd.f64(double %76, double %78, double %67)
%80 = or i64 %66, 1
%81 = mul nsw i64 %80, %28
%82 = add nuw nsw i64 %81, %60
%83 = getelementptr inbounds double, double* %0, i64 %82
%84 = load double, double* %83
%85 = add nuw nsw i64 %81, %64
%86 = getelementptr inbounds double, double* %0, i64 %85
%87 = load double, double* %86
%88 = fmul double %84, %87
%89 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %80
%90 = load double, double* %89
%91 = tail call double @llvm.fmuladd.f64(double %88, double %90, double %79)
%92 = add nuw nsw i64 %66, 2
%93 = add i64 %68, 2
%94 = icmp eq i64 %93, %32
br i1 %94, label %95, label %65
95: ; preds = %65, %63
%96 = phi double [ undef, %63 ], [ %91, %65 ]
%97 = phi i64 [ 0, %63 ], [ %92, %65 ]
%98 = phi double [ 0.000000e+00, %63 ], [ %91, %65 ]
br i1 %33, label %111, label %99
99: ; preds = %95
%100 = mul nsw i64 %97, %28
%101 = add nuw nsw i64 %100, %60
%102 = getelementptr inbounds double, double* %0, i64 %101
%103 = load double, double* %102
%104 = add nuw nsw i64 %100, %64
%105 = getelementptr inbounds double, double* %0, i64 %104
%106 = load double, double* %105
%107 = fmul double %103, %106
%108 = getelementptr inbounds [33 x double], [33 x double]* @svdvar.wti, i64 0, i64 %97
%109 = load double, double* %108
%110 = tail call double @llvm.fmuladd.f64(double %107, double %109, double %98)
br label %111
111: ; preds = %95, %99
%112 = phi double [ %96, %95 ], [ %110, %99 ]
%113 = add nuw nsw i64 %64, %62
%114 = getelementptr inbounds double, double* %3, i64 %113
store double %112, double* %114
%115 = mul nsw i64 %64, %28
%116 = add nuw nsw i64 %115, %60
%117 = getelementptr inbounds double, double* %3, i64 %116
store double %112, double* %117
%118 = add nuw nsw i64 %64, 1
%119 = icmp eq i64 %118, %61
br i1 %119, label %120, label %63
120: ; preds = %111
%121 = add nuw nsw i64 %60, 1
%122 = add nuw nsw i64 %61, 1
%123 = icmp eq i64 %121, %29
br i1 %123, label %124, label %59
124: ; preds = %120, %4, %26
ret void
}