unsigned work_size(const unsigned nr, const unsigned ns, const unsigned nt,
                   const unsigned npt_max) {
  unsigned n1, n2, wsize;
  if (nr > ns) {
    if (nr > nt)
      n1 = nr, n2 = (ns > nt ? ns : nt);
    else
      n1 = nt, n2 = nr;
  } else {
    if (ns > nt)
      n1 = ns, n2 = (nr > nt ? nr : nt);
    else
      n1 = nt, n2 = ns;
  }
  wsize = (12 + 2 * (nr + ns + nt + nr * ns)) * npt_max +
          (2 * nr > ns ? 2 * nr : ns);
  do {
    const unsigned temp = (2 * (nt + nr * ns));
    wsize = temp > wsize ? temp : wsize;
  } while (0);
  do {
    const unsigned temp = (4 * (nr + ns * nt));
    wsize = temp > wsize ? temp : wsize;
  } while (0);
  do {
    const unsigned temp = (4 * (n1 + n2));
    wsize = temp > wsize ? temp : wsize;
  } while (0);
  do {
    const unsigned temp = (npt_max * (15 + 3 * (2 * n1 + n2)));
    wsize = temp > wsize ? temp : wsize;
  } while (0);
  return wsize;
}


x86 assembly (clang, -O0):
<work_size>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    -0x4(%rbp),%eax
cmp    -0x8(%rbp),%eax
jbe    1167 <work_size+0x67>
mov    -0x4(%rbp),%eax
cmp    -0xc(%rbp),%eax
jbe    1156 <work_size+0x56>
mov    -0x4(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jbe    1145 <work_size+0x45>
mov    -0x8(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    114b <work_size+0x4b>
mov    -0xc(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    1162 <work_size+0x62>
mov    -0xc(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    11b2 <work_size+0xb2>
mov    -0x8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jbe    11a1 <work_size+0xa1>
mov    -0x8(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x4(%rbp),%eax
cmp    -0xc(%rbp),%eax
jbe    1190 <work_size+0x90>
mov    -0x4(%rbp),%eax
mov    %eax,-0x34(%rbp)
jmp    1196 <work_size+0x96>
mov    -0xc(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    11ad <work_size+0xad>
mov    -0xc(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    11b2 <work_size+0xb2>
mov    -0x4(%rbp),%eax
add    -0x8(%rbp),%eax
add    -0xc(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    %ecx,%eax
shl    $0x1,%eax
add    $0xc,%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
cmp    -0x8(%rbp),%eax
jbe    11ee <work_size+0xee>
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11f4 <work_size+0xf4>
mov    -0x8(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x3c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x1c(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    %ecx,%eax
shl    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jbe    1228 <work_size+0x128>
mov    -0x20(%rbp),%eax
mov    %eax,-0x40(%rbp)
jmp    122e <work_size+0x12e>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    1239 <work_size+0x139>
mov    -0x4(%rbp),%eax
mov    -0x8(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    %ecx,%eax
shl    $0x2,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jbe    1262 <work_size+0x162>
mov    -0x24(%rbp),%eax
mov    %eax,-0x44(%rbp)
jmp    1268 <work_size+0x168>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    1273 <work_size+0x173>
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jbe    1296 <work_size+0x196>
mov    -0x28(%rbp),%eax
mov    %eax,-0x48(%rbp)
jmp    129c <work_size+0x19c>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    12a7 <work_size+0x1a7>
mov    -0x10(%rbp),%eax
mov    -0x14(%rbp),%ecx
shl    $0x1,%ecx
add    -0x18(%rbp),%ecx
imul   $0x3,%ecx,%ecx
add    $0xf,%ecx
imul   %ecx,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jbe    12d6 <work_size+0x1d6>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
jmp    12dc <work_size+0x1dc>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<work_size>:
cmp    %esi,%edi
jbe    111c <work_size+0x1c>
mov    %edi,%r9d
mov    %edx,%r8d
cmp    %edx,%edi
jbe    1132 <work_size+0x32>
cmp    %edx,%esi
mov    %edx,%r9d
cmova  %esi,%r9d
mov    %edi,%r8d
jmp    1132 <work_size+0x32>
mov    %esi,%r9d
mov    %edx,%r8d
cmp    %edx,%esi
jbe    1132 <work_size+0x32>
cmp    %edx,%edi
mov    %edx,%r9d
cmova  %edi,%r9d
mov    %esi,%r8d
mov    %esi,%r10d
imul   %edi,%r10d
add    %edx,%r10d
lea    (%r10,%rdi,1),%eax
add    %esi,%eax
lea    (%rax,%rax,1),%r11d
add    $0xc,%r11d
imul   %ecx,%r11d
lea    (%rdi,%rdi,1),%eax
cmp    %esi,%eax
cmovbe %esi,%eax
add    %r11d,%eax
add    %r10d,%r10d
cmp    %eax,%r10d
cmovbe %eax,%r10d
imul   %esi,%edx
add    %edi,%edx
shl    $0x2,%edx
cmp    %r10d,%edx
cmovbe %r10d,%edx
lea    (%r8,%r9,1),%esi
shl    $0x2,%esi
cmp    %edx,%esi
cmovbe %edx,%esi
lea    (%r9,%r8,2),%eax
lea    (%rax,%rax,2),%eax
add    $0xf,%eax
imul   %ecx,%eax
cmp    %esi,%eax
cmovbe %esi,%eax
ret

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

LLVM IR (-O2):
define i32 @work_size(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = icmp ugt i32 %0, %1
br i1 %5, label %6, label %11
6: ; preds = %4
%7 = icmp ugt i32 %0, %2
br i1 %7, label %8, label %16
8: ; preds = %6
%9 = icmp ugt i32 %1, %2
%10 = select i1 %9, i32 %1, i32 %2
br label %16
11: ; preds = %4
%12 = icmp ugt i32 %1, %2
br i1 %12, label %13, label %16
13: ; preds = %11
%14 = icmp ugt i32 %0, %2
%15 = select i1 %14, i32 %0, i32 %2
br label %16
16: ; preds = %11, %6, %13, %8
%17 = phi i32 [ %10, %8 ], [ %15, %13 ], [ %0, %6 ], [ %1, %11 ]
%18 = phi i32 [ %0, %8 ], [ %1, %13 ], [ %2, %6 ], [ %2, %11 ]
%19 = mul i32 %1, %0
%20 = add i32 %19, %2
%21 = add i32 %20, %0
%22 = add i32 %21, %1
%23 = shl i32 %22, 1
%24 = add i32 %23, 12
%25 = mul i32 %24, %3
%26 = shl i32 %0, 1
%27 = icmp ugt i32 %26, %1
%28 = select i1 %27, i32 %26, i32 %1
%29 = add i32 %25, %28
%30 = shl i32 %20, 1
%31 = icmp ugt i32 %30, %29
%32 = select i1 %31, i32 %30, i32 %29
%33 = mul i32 %2, %1
%34 = add i32 %33, %0
%35 = shl i32 %34, 2
%36 = icmp ugt i32 %35, %32
%37 = select i1 %36, i32 %35, i32 %32
%38 = add i32 %18, %17
%39 = shl i32 %38, 2
%40 = icmp ugt i32 %39, %37
%41 = select i1 %40, i32 %39, i32 %37
%42 = shl i32 %18, 1
%43 = add i32 %42, %17
%44 = mul i32 %43, 3
%45 = add i32 %44, 15
%46 = mul i32 %45, %3
%47 = icmp ugt i32 %46, %41
%48 = select i1 %47, i32 %46, i32 %41
ret i32 %48
}