void scpcol(data, n, m, symmat) float **data, **symmat;
int n, m;
{
  int i, j1, j2;
  for (j1 = 1; j1 <= m; j1++) {
    for (j2 = j1; j2 <= m; j2++) {
      symmat[j1][j2] = 0.0;
      for (i = 1; i <= n; i++) {
        symmat[j1][j2] += data[i][j1] * data[i][j2];
      }
      symmat[j2][j1] = symmat[j1][j2];
    }
  }
  return;
}


x86 assembly (clang, -O0):
<scpcol>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1207 <scpcol+0x107>
mov    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     11f4 <scpcol+0xf4>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11bc <scpcol+0xbc>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1156 <scpcol+0x56>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112b <scpcol+0x2b>
jmp    11f9 <scpcol+0xf9>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1119 <scpcol+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<scpcol>:
push   %r15
push   %r14
push   %rbx
test   %edx,%edx
jle    11ec <scpcol+0xec>
lea    0x1(%rsi),%r8d
add    $0x1,%edx
lea    -0x1(%r8),%r9
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
neg    %r10
mov    $0x1,%r11d
jmp    113d <scpcol+0x3d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r11
cmp    %rdx,%r11
je     11ec <scpcol+0xec>
mov    (%rcx,%r11,8),%r14
mov    %r11,%r15
jmp    1166 <scpcol+0x66>
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
mov    (%rcx,%r15,8),%rax
movss  %xmm0,(%rax,%r11,4)
add    $0x1,%r15
cmp    %rdx,%r15
je     1130 <scpcol+0x30>
movl   $0x0,(%r14,%r15,4)
test   %esi,%esi
jle    1150 <scpcol+0x50>
mov    $0x1,%eax
xorps  %xmm0,%xmm0
cmp    $0x2,%r8d
je     11c7 <scpcol+0xc7>
mov    (%rdi,%rax,8),%rbx
movss  (%rbx,%r11,4),%xmm1
mulss  (%rbx,%r15,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r14,%r15,4)
mov    0x8(%rdi,%rax,8),%rbx
movss  (%rbx,%r11,4),%xmm0
mulss  (%rbx,%r15,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r14,%r15,4)
lea    (%r10,%rax,1),%rbx
add    $0x2,%rbx
add    $0x2,%rax
cmp    $0x1,%rbx
jne    1180 <scpcol+0x80>
test   $0x1,%r9b
je     1153 <scpcol+0x53>
mov    (%rdi,%rax,8),%rax
movss  (%rax,%r11,4),%xmm1
mulss  (%rax,%r15,4),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,(%r14,%r15,4)
jmp    1153 <scpcol+0x53>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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