int sub_mat_score(const short int s1i, const short int s2i, int sm[][5][5],
                  const int row, const int len) {
  int score;
  if (row < (15)) {
    score = sm[row][s1i][s2i];
    return score;
  }
  if (len - (row + 1) < (15)) {
    score = sm[((15) * 2) - (len - (row + 1))][s1i][s2i];
    return score;
  }
  score = sm[(15)][s1i][s2i];
  return score;
}


x86 assembly (clang, -O0):
<sub_mat_score>:
push   %rbp
mov    %rsp,%rbp
mov    %si,%ax
mov    %di,%si
mov    %si,-0x6(%rbp)
mov    %ax,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
cmpl   $0xf,-0x14(%rbp)
jge    1158 <sub_mat_score+0x58>
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
imul   $0x64,%rcx,%rcx
add    %rcx,%rax
movswq -0x6(%rbp),%rcx
imul   $0x14,%rcx,%rcx
add    %rcx,%rax
movswq -0x8(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    11d5 <sub_mat_score+0xd5>
mov    -0x18(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
sub    %ecx,%eax
cmp    $0xf,%eax
jge    11ae <sub_mat_score+0xae>
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%edx
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
sub    %ecx,%edx
mov    $0x1e,%ecx
sub    %edx,%ecx
movslq %ecx,%rcx
imul   $0x64,%rcx,%rcx
add    %rcx,%rax
movswq -0x6(%rbp),%rcx
imul   $0x14,%rcx,%rcx
add    %rcx,%rax
movswq -0x8(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    11d5 <sub_mat_score+0xd5>
mov    -0x10(%rbp),%rax
add    $0x5dc,%rax
movswq -0x6(%rbp),%rcx
imul   $0x14,%rcx,%rcx
add    %rcx,%rax
movswq -0x8(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<sub_mat_score>:
cmp    $0xe,%ecx
jg     110a <sub_mat_score+0xa>
movslq %ecx,%rax
jmp    1120 <sub_mat_score+0x20>
not    %ecx
add    %r8d,%ecx
mov    $0xf,%eax
cmp    $0xe,%ecx
jg     1120 <sub_mat_score+0x20>
mov    $0x1e,%eax
sub    %ecx,%eax
movswq %di,%rcx
movswq %si,%rsi
imul   $0x64,%rax,%rax
add    %rax,%rdx
lea    (%rcx,%rcx,4),%rax
lea    (%rdx,%rax,4),%rax
mov    (%rax,%rsi,4),%eax
ret

LLVM IR (-O0):
define i32 @sub_mat_score(i16 signext %0, i16 signext %1, [5 x [5 x i32]]* %2, i32 %3, i32 %4) {
%6 = alloca i32
%7 = alloca i16
%8 = alloca i16
%9 = alloca [5 x [5 x i32]]*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store i16 %0, i16* %7
store i16 %1, i16* %8
store [5 x [5 x i32]]* %2, [5 x [5 x i32]]** %9
store i32 %3, i32* %10
store i32 %4, i32* %11
%13 = load i32, i32* %10
%14 = icmp slt i32 %13, 15
br i1 %14, label %15, label %28
15: ; preds = %5
%16 = load [5 x [5 x i32]]*, [5 x [5 x i32]]** %9
%17 = load i32, i32* %10
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds [5 x [5 x i32]], [5 x [5 x i32]]* %16, i64 %18
%20 = load i16, i16* %7
%21 = sext i16 %20 to i64
%22 = getelementptr inbounds [5 x [5 x i32]], [5 x [5 x i32]]* %19, i64 0, i64 %21
%23 = load i16, i16* %8
%24 = sext i16 %23 to i64
%25 = getelementptr inbounds [5 x i32], [5 x i32]* %22, i64 0, i64 %24
%26 = load i32, i32* %25
store i32 %26, i32* %12
%27 = load i32, i32* %12
store i32 %27, i32* %6
br label %62
28: ; preds = %5
%29 = load i32, i32* %11
%30 = load i32, i32* %10
%31 = add nsw i32 %30, 1
%32 = sub nsw i32 %29, %31
%33 = icmp slt i32 %32, 15
br i1 %33, label %34, label %51
34: ; preds = %28
%35 = load [5 x [5 x i32]]*, [5 x [5 x i32]]** %9
%36 = load i32, i32* %11
%37 = load i32, i32* %10
%38 = add nsw i32 %37, 1
%39 = sub nsw i32 %36, %38
%40 = sub nsw i32 30, %39
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds [5 x [5 x i32]], [5 x [5 x i32]]* %35, i64 %41
%43 = load i16, i16* %7
%44 = sext i16 %43 to i64
%45 = getelementptr inbounds [5 x [5 x i32]], [5 x [5 x i32]]* %42, i64 0, i64 %44
%46 = load i16, i16* %8
%47 = sext i16 %46 to i64
%48 = getelementptr inbounds [5 x i32], [5 x i32]* %45, i64 0, i64 %47
%49 = load i32, i32* %48
store i32 %49, i32* %12
%50 = load i32, i32* %12
store i32 %50, i32* %6
br label %62
51: ; preds = %28
%52 = load [5 x [5 x i32]]*, [5 x [5 x i32]]** %9
%53 = getelementptr inbounds [5 x [5 x i32]], [5 x [5 x i32]]* %52, i64 15
%54 = load i16, i16* %7
%55 = sext i16 %54 to i64
%56 = getelementptr inbounds [5 x [5 x i32]], [5 x [5 x i32]]* %53, i64 0, i64 %55
%57 = load i16, i16* %8
%58 = sext i16 %57 to i64
%59 = getelementptr inbounds [5 x i32], [5 x i32]* %56, i64 0, i64 %58
%60 = load i32, i32* %59
store i32 %60, i32* %12
%61 = load i32, i32* %12
store i32 %61, i32* %6
br label %62
62: ; preds = %51, %34, %15
%63 = load i32, i32* %6
ret i32 %63
}

LLVM IR (-O2):
define i32 @sub_mat_score(i16 signext %0, i16 signext %1, [5 x [5 x i32]]* nocapture %2, i32 %3, i32 %4) {
%6 = icmp slt i32 %3, 15
br i1 %6, label %7, label %9
7: ; preds = %5
%8 = sext i32 %3 to i64
br label %16
9: ; preds = %5
%10 = xor i32 %3, -1
%11 = add i32 %10, %4
%12 = icmp slt i32 %11, 15
br i1 %12, label %13, label %16
13: ; preds = %9
%14 = sub nsw i32 30, %11
%15 = zext i32 %14 to i64
br label %16
16: ; preds = %9, %13, %7
%17 = phi i64 [ %15, %13 ], [ %8, %7 ], [ 15, %9 ]
%18 = sext i16 %0 to i64
%19 = sext i16 %1 to i64
%20 = getelementptr inbounds [5 x [5 x i32]], [5 x [5 x i32]]* %2, i64 %17, i64 %18, i64 %19
%21 = load i32, i32* %20
ret i32 %21
}