void fillDPtable(unsigned long *dptable, const char *u, unsigned long n,
                 const char *v, unsigned long m) {
  unsigned long i, j, nw, we;
  for (i = 0; i <= n; i++)
    dptable[i] = i;
  for (j = 1; j <= m; j++) {
    nw = dptable[0];
    dptable[0] = j;
    for (i = 1; i <= n; i++) {
      we = dptable[i];
      dptable[i] = nw + (u[i - 1] == v[j - 1] ? 0 : 1);
      if (dptable[i - 1] + 1 < dptable[i])
        dptable[i] = dptable[i - 1] + 1;
      if (we + 1 < dptable[i])
        dptable[i] = we + 1;
      nw = we;
    }
  }
}


x86 assembly (clang, -O0):
<fillDPtable>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     114f <fillDPtable+0x4f>
mov    -0x30(%rbp),%rdx
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1120 <fillDPtable+0x20>
movq   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     1285 <fillDPtable+0x185>
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rcx,(%rax)
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     126f <fillDPtable+0x16f>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x48(%rbp)
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
movsbl (%rax,%rcx,1),%esi
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
movsbl (%rax,%rcx,1),%edi
mov    $0x1,%eax
xor    %ecx,%ecx
cmp    %edi,%esi
cmove  %ecx,%eax
cltq
add    %rax,%rdx
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rax
add    $0x1,%rax
mov    -0x8(%rbp),%rcx
mov    -0x30(%rbp),%rdx
cmp    (%rcx,%rdx,8),%rax
jae    1228 <fillDPtable+0x128>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rdx
add    $0x1,%rdx
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    -0x8(%rbp),%rcx
mov    -0x30(%rbp),%rdx
cmp    (%rcx,%rdx,8),%rax
jae    1256 <fillDPtable+0x156>
mov    -0x48(%rbp),%rdx
add    $0x1,%rdx
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1183 <fillDPtable+0x83>
jmp    1274 <fillDPtable+0x174>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1157 <fillDPtable+0x57>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fillDPtable>:
push   %r14
push   %rbx
lea    0x1(%rdx),%r9
cmp    $0x4,%r9
jae    1114 <fillDPtable+0x14>
xor    %eax,%eax
jmp    1264 <fillDPtable+0x164>
mov    %r9,%rax
and    $0xfffffffffffffffc,%rax
lea    -0x4(%rax),%rbx
mov    %rbx,%r11
shr    $0x2,%r11
add    $0x1,%r11
mov    %r11d,%r10d
and    $0x3,%r10d
cmp    $0xc,%rbx
jae    1146 <fillDPtable+0x46>
movdqa 0xec1(%rip),%xmm0
xor    %ebx,%ebx
jmp    1218 <fillDPtable+0x118>
and    $0xfffffffffffffffc,%r11
movdqa 0xeae(%rip),%xmm0
xor    %ebx,%ebx
movdqa 0xeb3(%rip),%xmm8
movdqa 0xeba(%rip),%xmm9
movdqa 0xec1(%rip),%xmm10
movdqa 0xec9(%rip),%xmm4
movdqa 0xed1(%rip),%xmm5
movdqa 0xed9(%rip),%xmm6
movdqa 0xee1(%rip),%xmm7
movdqa 0xee9(%rip),%xmm1
nopw   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm2
paddq  %xmm8,%xmm2
movdqu %xmm0,(%rdi,%rbx,8)
movdqu %xmm2,0x10(%rdi,%rbx,8)
movdqa %xmm0,%xmm2
paddq  %xmm9,%xmm2
movdqa %xmm0,%xmm3
paddq  %xmm10,%xmm3
movdqu %xmm2,0x20(%rdi,%rbx,8)
movdqu %xmm3,0x30(%rdi,%rbx,8)
movdqa %xmm0,%xmm2
paddq  %xmm4,%xmm2
movdqa %xmm0,%xmm3
paddq  %xmm5,%xmm3
movdqu %xmm2,0x40(%rdi,%rbx,8)
movdqu %xmm3,0x50(%rdi,%rbx,8)
movdqa %xmm0,%xmm2
paddq  %xmm6,%xmm2
movdqa %xmm0,%xmm3
paddq  %xmm7,%xmm3
movdqu %xmm2,0x60(%rdi,%rbx,8)
movdqu %xmm3,0x70(%rdi,%rbx,8)
add    $0x10,%rbx
paddq  %xmm1,%xmm0
add    $0xfffffffffffffffc,%r11
jne    11a0 <fillDPtable+0xa0>
test   %r10,%r10
je     126c <fillDPtable+0x16c>
lea    (%rdi,%rbx,8),%r11
add    $0x10,%r11
shl    $0x5,%r10
xor    %ebx,%ebx
movdqa 0xddd(%rip),%xmm1
movdqa 0xde5(%rip),%xmm2
nopl   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm3
paddq  %xmm1,%xmm3
movdqu %xmm0,-0x10(%r11,%rbx,1)
movdqu %xmm3,(%r11,%rbx,1)
paddq  %xmm2,%xmm0
add    $0x20,%rbx
cmp    %rbx,%r10
jne    1240 <fillDPtable+0x140>
jmp    126c <fillDPtable+0x16c>
mov    %rax,(%rdi,%rax,8)
add    $0x1,%rax
cmp    %rax,%r9
jne    1264 <fillDPtable+0x164>
test   %r8,%r8
je     12f1 <fillDPtable+0x1f1>
mov    (%rdi),%r10
mov    $0x1,%r9d
jmp    129f <fillDPtable+0x19f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r9,%r10
lea    0x1(%r9),%rax
cmp    %r8,%r9
mov    %rax,%r9
je     12f1 <fillDPtable+0x1f1>
mov    %r9,(%rdi)
test   %rdx,%rdx
je     1290 <fillDPtable+0x190>
xor    %r14d,%r14d
mov    %r9,%rax
jmp    12c1 <fillDPtable+0x1c1>
nop
mov    %r11,0x8(%rdi,%r14,8)
add    $0x1,%r14
mov    %r11,%rax
cmp    %r14,%rdx
je     1290 <fillDPtable+0x190>
movzbl (%rsi,%r14,1),%r11d
xor    %ebx,%ebx
cmp    -0x1(%rcx,%r9,1),%r11b
setne  %bl
add    %r10,%rbx
mov    0x8(%rdi,%r14,8),%r10
add    $0x1,%rax
cmp    %rbx,%rax
cmovae %rbx,%rax
lea    0x1(%r10),%r11
cmp    %rax,%r11
jb     12b0 <fillDPtable+0x1b0>
mov    %rax,%r11
jmp    12b0 <fillDPtable+0x1b0>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @fillDPtable(i64* %0, i8* %1, i64 %2, i8* %3, i64 %4) {
%6 = alloca i64*
%7 = alloca i8*
%8 = alloca i64
%9 = alloca i8*
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
store i64* %0, i64** %6
store i8* %1, i8** %7
store i64 %2, i64* %8
store i8* %3, i8** %9
store i64 %4, i64* %10
store i64 0, i64* %11
br label %15
15: ; preds = %24, %5
%16 = load i64, i64* %11
%17 = load i64, i64* %8
%18 = icmp ule i64 %16, %17
br i1 %18, label %19, label %27
19: ; preds = %15
%20 = load i64, i64* %11
%21 = load i64*, i64** %6
%22 = load i64, i64* %11
%23 = getelementptr inbounds i64, i64* %21, i64 %22
store i64 %20, i64* %23
br label %24
24: ; preds = %19
%25 = load i64, i64* %11
%26 = add i64 %25, 1
store i64 %26, i64* %11
br label %15
27: ; preds = %15
store i64 1, i64* %12
br label %28
28: ; preds = %110, %27
%29 = load i64, i64* %12
%30 = load i64, i64* %10
%31 = icmp ule i64 %29, %30
br i1 %31, label %32, label %113
32: ; preds = %28
%33 = load i64*, i64** %6
%34 = getelementptr inbounds i64, i64* %33, i64 0
%35 = load i64, i64* %34
store i64 %35, i64* %13
%36 = load i64, i64* %12
%37 = load i64*, i64** %6
%38 = getelementptr inbounds i64, i64* %37, i64 0
store i64 %36, i64* %38
store i64 1, i64* %11
br label %39
39: ; preds = %106, %32
%40 = load i64, i64* %11
%41 = load i64, i64* %8
%42 = icmp ule i64 %40, %41
br i1 %42, label %43, label %109
43: ; preds = %39
%44 = load i64*, i64** %6
%45 = load i64, i64* %11
%46 = getelementptr inbounds i64, i64* %44, i64 %45
%47 = load i64, i64* %46
store i64 %47, i64* %14
%48 = load i64, i64* %13
%49 = load i8*, i8** %7
%50 = load i64, i64* %11
%51 = sub i64 %50, 1
%52 = getelementptr inbounds i8, i8* %49, i64 %51
%53 = load i8, i8* %52
%54 = sext i8 %53 to i32
%55 = load i8*, i8** %9
%56 = load i64, i64* %12
%57 = sub i64 %56, 1
%58 = getelementptr inbounds i8, i8* %55, i64 %57
%59 = load i8, i8* %58
%60 = sext i8 %59 to i32
%61 = icmp eq i32 %54, %60
%62 = zext i1 %61 to i64
%63 = select i1 %61, i32 0, i32 1
%64 = sext i32 %63 to i64
%65 = add i64 %48, %64
%66 = load i64*, i64** %6
%67 = load i64, i64* %11
%68 = getelementptr inbounds i64, i64* %66, i64 %67
store i64 %65, i64* %68
%69 = load i64*, i64** %6
%70 = load i64, i64* %11
%71 = sub i64 %70, 1
%72 = getelementptr inbounds i64, i64* %69, i64 %71
%73 = load i64, i64* %72
%74 = add i64 %73, 1
%75 = load i64*, i64** %6
%76 = load i64, i64* %11
%77 = getelementptr inbounds i64, i64* %75, i64 %76
%78 = load i64, i64* %77
%79 = icmp ult i64 %74, %78
br i1 %79, label %80, label %90
80: ; preds = %43
%81 = load i64*, i64** %6
%82 = load i64, i64* %11
%83 = sub i64 %82, 1
%84 = getelementptr inbounds i64, i64* %81, i64 %83
%85 = load i64, i64* %84
%86 = add i64 %85, 1
%87 = load i64*, i64** %6
%88 = load i64, i64* %11
%89 = getelementptr inbounds i64, i64* %87, i64 %88
store i64 %86, i64* %89
br label %90
90: ; preds = %80, %43
%91 = load i64, i64* %14
%92 = add i64 %91, 1
%93 = load i64*, i64** %6
%94 = load i64, i64* %11
%95 = getelementptr inbounds i64, i64* %93, i64 %94
%96 = load i64, i64* %95
%97 = icmp ult i64 %92, %96
br i1 %97, label %98, label %104
98: ; preds = %90
%99 = load i64, i64* %14
%100 = add i64 %99, 1
%101 = load i64*, i64** %6
%102 = load i64, i64* %11
%103 = getelementptr inbounds i64, i64* %101, i64 %102
store i64 %100, i64* %103
br label %104
104: ; preds = %98, %90
%105 = load i64, i64* %14
store i64 %105, i64* %13
br label %106
106: ; preds = %104
%107 = load i64, i64* %11
%108 = add i64 %107, 1
store i64 %108, i64* %11
br label %39
109: ; preds = %39
br label %110
110: ; preds = %109
%111 = load i64, i64* %12
%112 = add i64 %111, 1
store i64 %112, i64* %12
br label %28
113: ; preds = %28
ret void
}

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