void weight_matrix_max(float **weight_matrix, int motif_length,
                       float *max_weight_matrix) {
  int i;
  int j;
  float current_max;
  for (i = 1; i <= motif_length; i++) {
    current_max = -1000000;
    for (j = 1; j <= 4; j++) {
      if (weight_matrix[j][i] > current_max)
        current_max = weight_matrix[j][i];
    }
    max_weight_matrix[i] = current_max;
  }
}


x86 assembly (clang, -O0):
<weight_matrix_max>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11ac <weight_matrix_max+0xac>
movss  0xed6(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x1,-0x20(%rbp)
cmpl   $0x4,-0x20(%rbp)
jg     118c <weight_matrix_max+0x8c>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x24(%rbp),%xmm0
jbe    1179 <weight_matrix_max+0x79>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
jmp    117e <weight_matrix_max+0x7e>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1136 <weight_matrix_max+0x36>
movss  -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1116 <weight_matrix_max+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<weight_matrix_max>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %esi,%esi
jle    1244 <weight_matrix_max+0x144>
mov    %esi,%ebp
lea    0x1(%rbp),%r8
mov    0x8(%rdi),%r9
mov    0x10(%rdi),%r10
mov    0x18(%rdi),%rcx
mov    0x20(%rdi),%rdi
mov    $0x1,%eax
cmp    $0x4,%esi
jb     1205 <weight_matrix_max+0x105>
mov    %rbp,-0x8(%rsp)
lea    0x4(%rdx),%rbp
lea    (%rdx,%r8,4),%rsi
lea    0x4(%r9),%r14
lea    (%r9,%r8,4),%rbx
lea    0x4(%r10),%r15
lea    (%r10,%r8,4),%r12
cmp    %rbx,%rbp
setb   -0x9(%rsp)
cmp    %rsi,%r14
lea    0x4(%rcx),%r11
setb   %r14b
cmp    %r12,%rbp
setb   %r12b
cmp    %rsi,%r15
lea    (%rcx,%r8,4),%rbx
setb   %r15b
cmp    %rbx,%rbp
setb   %r13b
cmp    %rsi,%r11
lea    (%rdi,%r8,4),%rbx
setb   %r11b
cmp    %rbx,%rbp
lea    0x4(%rdi),%rbx
setb   %bpl
cmp    %rsi,%rbx
setb   %sil
test   %r14b,-0x9(%rsp)
jne    1205 <weight_matrix_max+0x105>
and    %r15b,%r12b
jne    1205 <weight_matrix_max+0x105>
and    %r11b,%r13b
jne    1205 <weight_matrix_max+0x105>
and    %sil,%bpl
jne    1205 <weight_matrix_max+0x105>
mov    -0x8(%rsp),%rbp
mov    %ebp,%esi
and    $0xfffffffc,%esi
lea    0x1(%rsi),%rax
xor    %ebx,%ebx
movaps 0xe37(%rip),%xmm0
nopl   0x0(%rax)
movups 0x4(%r9,%rbx,4),%xmm1
maxps  %xmm0,%xmm1
movups 0x4(%r10,%rbx,4),%xmm2
maxps  %xmm1,%xmm2
movups 0x4(%rcx,%rbx,4),%xmm1
maxps  %xmm2,%xmm1
movups 0x4(%rdi,%rbx,4),%xmm2
maxps  %xmm1,%xmm2
movups %xmm2,0x4(%rdx,%rbx,4)
add    $0x4,%rbx
cmp    %rbx,%rsi
jne    11d0 <weight_matrix_max+0xd0>
cmp    %rsi,%rbp
je     1244 <weight_matrix_max+0x144>
movss  0xe03(%rip),%xmm0
nopl   (%rax)
movss  (%r9,%rax,4),%xmm1
maxss  %xmm0,%xmm1
movss  (%r10,%rax,4),%xmm2
maxss  %xmm1,%xmm2
movss  (%rcx,%rax,4),%xmm1
maxss  %xmm2,%xmm1
movss  (%rdi,%rax,4),%xmm2
maxss  %xmm1,%xmm2
movss  %xmm2,(%rdx,%rax,4)
add    $0x1,%rax
cmp    %rax,%r8
jne    1210 <weight_matrix_max+0x110>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @weight_matrix_max(float** %0, i32 %1, float* %2) {
%4 = alloca float**
%5 = alloca i32
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca float
store float** %0, float*** %4
store i32 %1, i32* %5
store float* %2, float** %6
store i32 1, i32* %7
br label %10
10: ; preds = %50, %3
%11 = load i32, i32* %7
%12 = load i32, i32* %5
%13 = icmp sle i32 %11, %12
br i1 %13, label %14, label %53
14: ; preds = %10
store float -1.000000e+06, float* %9
store i32 1, i32* %8
br label %15
15: ; preds = %41, %14
%16 = load i32, i32* %8
%17 = icmp sle i32 %16, 4
br i1 %17, label %18, label %44
18: ; preds = %15
%19 = load float**, float*** %4
%20 = load i32, i32* %8
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds float*, float** %19, i64 %21
%23 = load float*, float** %22
%24 = load i32, i32* %7
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds float, float* %23, i64 %25
%27 = load float, float* %26
%28 = load float, float* %9
%29 = fcmp ogt float %27, %28
br i1 %29, label %30, label %40
30: ; preds = %18
%31 = load float**, float*** %4
%32 = load i32, i32* %8
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float*, float** %31, i64 %33
%35 = load float*, float** %34
%36 = load i32, i32* %7
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds float, float* %35, i64 %37
%39 = load float, float* %38
store float %39, float* %9
br label %40
40: ; preds = %30, %18
br label %41
41: ; preds = %40
%42 = load i32, i32* %8
%43 = add nsw i32 %42, 1
store i32 %43, i32* %8
br label %15
44: ; preds = %15
%45 = load float, float* %9
%46 = load float*, float** %6
%47 = load i32, i32* %7
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %46, i64 %48
store float %45, float* %49
br label %50
50: ; preds = %44
%51 = load i32, i32* %7
%52 = add nsw i32 %51, 1
store i32 %52, i32* %7
br label %10
53: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @weight_matrix_max(float** nocapture %0, i32 %1, float* nocapture writeonly %2) {
%4 = icmp slt i32 %1, 1
br i1 %4, label %99, label %5
5: ; preds = %3
%6 = getelementptr inbounds float*, float** %0, i64 1
%7 = add nuw i32 %1, 1
%8 = zext i32 %7 to i64
%9 = load float*, float** %6
%10 = getelementptr inbounds float*, float** %0, i64 2
%11 = load float*, float** %10
%12 = getelementptr inbounds float*, float** %0, i64 3
%13 = load float*, float** %12
%14 = getelementptr inbounds float*, float** %0, i64 4
%15 = load float*, float** %14
%16 = add nsw i64 %8, -1
%17 = icmp ult i64 %16, 4
br i1 %17, label %76, label %18
18: ; preds = %5
%19 = getelementptr float, float* %2, i64 1
%20 = getelementptr float, float* %2, i64 %8
%21 = getelementptr float, float* %9, i64 1
%22 = getelementptr float, float* %9, i64 %8
%23 = getelementptr float, float* %11, i64 1
%24 = getelementptr float, float* %11, i64 %8
%25 = getelementptr float, float* %13, i64 1
%26 = getelementptr float, float* %13, i64 %8
%27 = getelementptr float, float* %15, i64 1
%28 = getelementptr float, float* %15, i64 %8
%29 = icmp ult float* %19, %22
%30 = icmp ult float* %21, %20
%31 = and i1 %29, %30
%32 = icmp ult float* %19, %24
%33 = icmp ult float* %23, %20
%34 = and i1 %32, %33
%35 = or i1 %31, %34
%36 = icmp ult float* %19, %26
%37 = icmp ult float* %25, %20
%38 = and i1 %36, %37
%39 = or i1 %35, %38
%40 = icmp ult float* %19, %28
%41 = icmp ult float* %27, %20
%42 = and i1 %40, %41
%43 = or i1 %39, %42
br i1 %43, label %76, label %44
44: ; preds = %18
%45 = and i64 %16, -4
%46 = or i64 %45, 1
br label %47
47: ; preds = %47, %44
%48 = phi i64 [ 0, %44 ], [ %72, %47 ]
%49 = or i64 %48, 1
%50 = getelementptr inbounds float, float* %9, i64 %49
%51 = bitcast float* %50 to <4 x float>*
%52 = load <4 x float>, <4 x float>* %51
%53 = fcmp ogt <4 x float> %52, <float -1.000000e+06, float -1.000000e+06, float -1.000000e+06, float -1.000000e+06>
%54 = select <4 x i1> %53, <4 x float> %52, <4 x float> <float -1.000000e+06, float -1.000000e+06, float -1.000000e+06, float -1.000000e+06>
%55 = getelementptr inbounds float, float* %11, i64 %49
%56 = bitcast float* %55 to <4 x float>*
%57 = load <4 x float>, <4 x float>* %56
%58 = fcmp ogt <4 x float> %57, %54
%59 = select <4 x i1> %58, <4 x float> %57, <4 x float> %54
%60 = getelementptr inbounds float, float* %13, i64 %49
%61 = bitcast float* %60 to <4 x float>*
%62 = load <4 x float>, <4 x float>* %61
%63 = fcmp ogt <4 x float> %62, %59
%64 = select <4 x i1> %63, <4 x float> %62, <4 x float> %59
%65 = getelementptr inbounds float, float* %15, i64 %49
%66 = bitcast float* %65 to <4 x float>*
%67 = load <4 x float>, <4 x float>* %66
%68 = fcmp ogt <4 x float> %67, %64
%69 = select <4 x i1> %68, <4 x float> %67, <4 x float> %64
%70 = getelementptr inbounds float, float* %2, i64 %49
%71 = bitcast float* %70 to <4 x float>*
store <4 x float> %69, <4 x float>* %71
%72 = add nuw i64 %48, 4
%73 = icmp eq i64 %72, %45
br i1 %73, label %74, label %47
74: ; preds = %47
%75 = icmp eq i64 %16, %45
br i1 %75, label %99, label %76
76: ; preds = %18, %5, %74
%77 = phi i64 [ 1, %18 ], [ 1, %5 ], [ %46, %74 ]
br label %78
78: ; preds = %76, %78
%79 = phi i64 [ %97, %78 ], [ %77, %76 ]
%80 = getelementptr inbounds float, float* %9, i64 %79
%81 = load float, float* %80
%82 = fcmp ogt float %81, -1.000000e+06
%83 = select i1 %82, float %81, float -1.000000e+06
%84 = getelementptr inbounds float, float* %11, i64 %79
%85 = load float, float* %84
%86 = fcmp ogt float %85, %83
%87 = select i1 %86, float %85, float %83
%88 = getelementptr inbounds float, float* %13, i64 %79
%89 = load float, float* %88
%90 = fcmp ogt float %89, %87
%91 = select i1 %90, float %89, float %87
%92 = getelementptr inbounds float, float* %15, i64 %79
%93 = load float, float* %92
%94 = fcmp ogt float %93, %91
%95 = select i1 %94, float %93, float %91
%96 = getelementptr inbounds float, float* %2, i64 %79
store float %95, float* %96
%97 = add nuw nsw i64 %79, 1
%98 = icmp eq i64 %97, %8
br i1 %98, label %99, label %78
99: ; preds = %78, %74, %3
ret void
}