void poly_sub2(double u1[5][5], double u2[5][5], int p1[2], int p2[2],
               double u3[5][5], int p3[2]) {
  int i, j, p11, p12, p21, p22;
  int i1_ok, i2_ok;
  p11 = p1[0];
  p12 = p1[1];
  p21 = p2[0];
  p22 = p2[1];
  p3[0] = ((p11) > (p21)) ? (p11) : (p21);
  p3[1] = ((p12) > (p22)) ? (p12) : (p22);
  for (i = 0; i <= p3[1]; i++) {
    i1_ok = (i > p12);
    i2_ok = (i > p22);
    for (j = 0; j <= p3[0]; j++) {
      if (i2_ok || (j > p21))
        u3[i][j] = u1[i][j];
      else if (i1_ok || (j > p11))
        u3[i][j] = -u2[i][j];
      else
        u3[i][j] = u1[i][j] - u2[i][j];
    }
  }
  return;
}


x86 assembly (clang, -O0):
<poly_sub2>:
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)
mov    %r9,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x44(%rbp)
mov    -0x20(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x44(%rbp),%eax
jle    1159 <poly_sub2+0x59>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x54(%rbp)
jmp    115f <poly_sub2+0x5f>
mov    -0x44(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%ecx
mov    -0x30(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x40(%rbp),%eax
cmp    -0x48(%rbp),%eax
jle    117f <poly_sub2+0x7f>
mov    -0x40(%rbp),%eax
mov    %eax,-0x58(%rbp)
jmp    1185 <poly_sub2+0x85>
mov    -0x48(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x58(%rbp),%ecx
mov    -0x30(%rbp),%rax
mov    %ecx,0x4(%rax)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x30(%rbp),%rcx
cmp    0x4(%rcx),%eax
jg     12fe <poly_sub2+0x1fe>
mov    -0x34(%rbp),%eax
cmp    -0x40(%rbp),%eax
setg   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x48(%rbp),%eax
setg   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x50(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x30(%rbp),%rcx
cmp    (%rcx),%eax
jg     12eb <poly_sub2+0x1eb>
cmpl   $0x0,-0x50(%rbp)
jne    11f4 <poly_sub2+0xf4>
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
jle    1229 <poly_sub2+0x129>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    12d8 <poly_sub2+0x1d8>
cmpl   $0x0,-0x4c(%rbp)
jne    123f <poly_sub2+0x13f>
mov    -0x38(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jle    128b <poly_sub2+0x18b>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    12d3 <poly_sub2+0x1d3>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x28,%rcx,%rcx
add    %rcx,%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    12d8 <poly_sub2+0x1d8>
jmp    12dd <poly_sub2+0x1dd>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11cf <poly_sub2+0xcf>
jmp    12f0 <poly_sub2+0x1f0>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1196 <poly_sub2+0x96>
pop    %rbp
ret

x86 assembly (clang, -O2):
<poly_sub2>:
push   %r15
push   %r14
push   %r12
push   %rbx
movslq (%rdx),%r10
movslq 0x4(%rdx),%r15
movslq (%rcx),%r12
movslq 0x4(%rcx),%rcx
cmp    %r12d,%r10d
mov    %r12d,%r11d
cmovg  %r10d,%r11d
mov    %r11d,(%r9)
cmp    %ecx,%r15d
mov    %ecx,%r14d
cmovg  %r15d,%r14d
mov    %r14d,0x4(%r9)
test   %r14d,%r14d
js     11bb <poly_sub2+0xbb>
lea    0x1(%r11),%ebx
add    $0x1,%r14d
xor    %edx,%edx
movapd 0xeb5(%rip),%xmm0
jmp    1165 <poly_sub2+0x65>
nopl   (%rax)
add    $0x1,%rdx
add    $0x28,%r8
add    $0x28,%rsi
add    $0x28,%rdi
cmp    %r14,%rdx
je     11bb <poly_sub2+0xbb>
test   %r11d,%r11d
js     1150 <poly_sub2+0x50>
xor    %eax,%eax
jmp    1184 <poly_sub2+0x84>
xchg   %ax,%ax
movsd  (%rdi,%rax,8),%xmm1
movsd  %xmm1,(%r8,%rax,8)
add    $0x1,%rax
cmp    %rax,%rbx
je     1150 <poly_sub2+0x50>
cmp    %rcx,%rdx
jg     1170 <poly_sub2+0x70>
cmp    %r12,%rax
jg     1170 <poly_sub2+0x70>
cmp    %r15,%rdx
jg     11b0 <poly_sub2+0xb0>
cmp    %r10,%rax
jg     11b0 <poly_sub2+0xb0>
movsd  (%rdi,%rax,8),%xmm1
subsd  (%rsi,%rax,8),%xmm1
jmp    1175 <poly_sub2+0x75>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rax,8),%xmm1
xorpd  %xmm0,%xmm1
jmp    1175 <poly_sub2+0x75>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @poly_sub2([5 x double]* %0, [5 x double]* %1, i32* %2, i32* %3, [5 x double]* %4, i32* %5) {
%7 = alloca [5 x double]*
%8 = alloca [5 x double]*
%9 = alloca i32*
%10 = alloca i32*
%11 = alloca [5 x double]*
%12 = alloca i32*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
store [5 x double]* %0, [5 x double]** %7
store [5 x double]* %1, [5 x double]** %8
store i32* %2, i32** %9
store i32* %3, i32** %10
store [5 x double]* %4, [5 x double]** %11
store i32* %5, i32** %12
%21 = load i32*, i32** %9
%22 = getelementptr inbounds i32, i32* %21, i64 0
%23 = load i32, i32* %22
store i32 %23, i32* %15
%24 = load i32*, i32** %9
%25 = getelementptr inbounds i32, i32* %24, i64 1
%26 = load i32, i32* %25
store i32 %26, i32* %16
%27 = load i32*, i32** %10
%28 = getelementptr inbounds i32, i32* %27, i64 0
%29 = load i32, i32* %28
store i32 %29, i32* %17
%30 = load i32*, i32** %10
%31 = getelementptr inbounds i32, i32* %30, i64 1
%32 = load i32, i32* %31
store i32 %32, i32* %18
%33 = load i32, i32* %15
%34 = load i32, i32* %17
%35 = icmp sgt i32 %33, %34
br i1 %35, label %36, label %38
36: ; preds = %6
%37 = load i32, i32* %15
br label %40
38: ; preds = %6
%39 = load i32, i32* %17
br label %40
40: ; preds = %38, %36
%41 = phi i32 [ %37, %36 ], [ %39, %38 ]
%42 = load i32*, i32** %12
%43 = getelementptr inbounds i32, i32* %42, i64 0
store i32 %41, i32* %43
%44 = load i32, i32* %16
%45 = load i32, i32* %18
%46 = icmp sgt i32 %44, %45
br i1 %46, label %47, label %49
47: ; preds = %40
%48 = load i32, i32* %16
br label %51
49: ; preds = %40
%50 = load i32, i32* %18
br label %51
51: ; preds = %49, %47
%52 = phi i32 [ %48, %47 ], [ %50, %49 ]
%53 = load i32*, i32** %12
%54 = getelementptr inbounds i32, i32* %53, i64 1
store i32 %52, i32* %54
store i32 0, i32* %13
br label %55
55: ; preds = %154, %51
%56 = load i32, i32* %13
%57 = load i32*, i32** %12
%58 = getelementptr inbounds i32, i32* %57, i64 1
%59 = load i32, i32* %58
%60 = icmp sle i32 %56, %59
br i1 %60, label %61, label %157
61: ; preds = %55
%62 = load i32, i32* %13
%63 = load i32, i32* %16
%64 = icmp sgt i32 %62, %63
%65 = zext i1 %64 to i32
store i32 %65, i32* %19
%66 = load i32, i32* %13
%67 = load i32, i32* %18
%68 = icmp sgt i32 %66, %67
%69 = zext i1 %68 to i32
store i32 %69, i32* %20
store i32 0, i32* %14
br label %70
70: ; preds = %150, %61
%71 = load i32, i32* %14
%72 = load i32*, i32** %12
%73 = getelementptr inbounds i32, i32* %72, i64 0
%74 = load i32, i32* %73
%75 = icmp sle i32 %71, %74
br i1 %75, label %76, label %153
76: ; preds = %70
%77 = load i32, i32* %20
%78 = icmp ne i32 %77, 0
br i1 %78, label %83, label %79
79: ; preds = %76
%80 = load i32, i32* %14
%81 = load i32, i32* %17
%82 = icmp sgt i32 %80, %81
br i1 %82, label %83, label %99
83: ; preds = %79, %76
%84 = load [5 x double]*, [5 x double]** %7
%85 = load i32, i32* %13
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds [5 x double], [5 x double]* %84, i64 %86
%88 = load i32, i32* %14
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds [5 x double], [5 x double]* %87, i64 0, i64 %89
%91 = load double, double* %90
%92 = load [5 x double]*, [5 x double]** %11
%93 = load i32, i32* %13
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds [5 x double], [5 x double]* %92, i64 %94
%96 = load i32, i32* %14
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds [5 x double], [5 x double]* %95, i64 0, i64 %97
store double %91, double* %98
br label %149
99: ; preds = %79
%100 = load i32, i32* %19
%101 = icmp ne i32 %100, 0
br i1 %101, label %106, label %102
102: ; preds = %99
%103 = load i32, i32* %14
%104 = load i32, i32* %15
%105 = icmp sgt i32 %103, %104
br i1 %105, label %106, label %123
106: ; preds = %102, %99
%107 = load [5 x double]*, [5 x double]** %8
%108 = load i32, i32* %13
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds [5 x double], [5 x double]* %107, i64 %109
%111 = load i32, i32* %14
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds [5 x double], [5 x double]* %110, i64 0, i64 %112
%114 = load double, double* %113
%115 = fneg double %114
%116 = load [5 x double]*, [5 x double]** %11
%117 = load i32, i32* %13
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds [5 x double], [5 x double]* %116, i64 %118
%120 = load i32, i32* %14
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds [5 x double], [5 x double]* %119, i64 0, i64 %121
store double %115, double* %122
br label %148
123: ; preds = %102
%124 = load [5 x double]*, [5 x double]** %7
%125 = load i32, i32* %13
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds [5 x double], [5 x double]* %124, i64 %126
%128 = load i32, i32* %14
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds [5 x double], [5 x double]* %127, i64 0, i64 %129
%131 = load double, double* %130
%132 = load [5 x double]*, [5 x double]** %8
%133 = load i32, i32* %13
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds [5 x double], [5 x double]* %132, i64 %134
%136 = load i32, i32* %14
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds [5 x double], [5 x double]* %135, i64 0, i64 %137
%139 = load double, double* %138
%140 = fsub double %131, %139
%141 = load [5 x double]*, [5 x double]** %11
%142 = load i32, i32* %13
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds [5 x double], [5 x double]* %141, i64 %143
%145 = load i32, i32* %14
%146 = sext i32 %145 to i64
%147 = getelementptr inbounds [5 x double], [5 x double]* %144, i64 0, i64 %146
store double %140, double* %147
br label %148
148: ; preds = %123, %106
br label %149
149: ; preds = %148, %83
br label %150
150: ; preds = %149
%151 = load i32, i32* %14
%152 = add nsw i32 %151, 1
store i32 %152, i32* %14
br label %70
153: ; preds = %70
br label %154
154: ; preds = %153
%155 = load i32, i32* %13
%156 = add nsw i32 %155, 1
store i32 %156, i32* %13
br label %55
157: ; preds = %55
ret void
}

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