order_poly2(poly, ord1, ord2) double poly[25][25];
int *ord1, *ord2;
{
  int i, j;
  *ord1 = *ord2 = 0;
  for (i = 0; i < 25; ++i)
    for (j = 0; j < 25; ++j)
      if (poly[i][j] != 0.0) {
        *ord1 = *ord1 > i ? *ord1 : i;
        ;
        *ord2 = *ord2 > j ? *ord2 : j;
        ;
      }
}


x86 assembly (clang, -O0):
<order_poly2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    -0x20(%rbp),%rax
movl   $0x0,(%rax)
mov    -0x18(%rbp),%rax
movl   $0x0,(%rax)
movl   $0x0,-0x24(%rbp)
cmpl   $0x19,-0x24(%rbp)
jge    11f7 <order_poly2+0xf7>
movl   $0x0,-0x28(%rbp)
cmpl   $0x19,-0x28(%rbp)
jge    11e4 <order_poly2+0xe4>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0xc8,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1179 <order_poly2+0x79>
jp     1179 <order_poly2+0x79>
jmp    11d1 <order_poly2+0xd1>
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
cmp    -0x24(%rbp),%eax
jle    1196 <order_poly2+0x96>
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x2c(%rbp)
jmp    119c <order_poly2+0x9c>
mov    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
cmp    -0x28(%rbp),%eax
jle    11c2 <order_poly2+0xc2>
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x30(%rbp)
jmp    11c8 <order_poly2+0xc8>
mov    -0x28(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
jmp    11d6 <order_poly2+0xd6>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113c <order_poly2+0x3c>
jmp    11e9 <order_poly2+0xe9>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112b <order_poly2+0x2b>
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<order_poly2>:
movl   $0x0,(%rdx)
movl   $0x0,(%rsi)
xor    %r8d,%r8d
xorpd  %xmm0,%xmm0
jmp    1131 <order_poly2+0x31>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r8
add    $0xc8,%rdi
cmp    $0x19,%r8
je     1176 <order_poly2+0x76>
xor    %ecx,%ecx
jmp    114a <order_poly2+0x4a>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rcx
cmp    $0x19,%rcx
je     1120 <order_poly2+0x20>
movsd  (%rdi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1157 <order_poly2+0x57>
jnp    1140 <order_poly2+0x40>
movslq (%rsi),%r9
cmp    %r9,%r8
mov    %r8d,%eax
cmovl  %r9d,%eax
mov    %eax,(%rsi)
movslq (%rdx),%r9
cmp    %r9,%rcx
mov    %ecx,%eax
cmovl  %r9d,%eax
mov    %eax,(%rdx)
jmp    1140 <order_poly2+0x40>
ret

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

LLVM IR (-O2):
define i32 @order_poly2([25 x double]* nocapture %0, i32* nocapture %1, i32* nocapture %2) {
store i32 0, i32* %2
store i32 0, i32* %1
br label %4
4: ; preds = %3, %25
%5 = phi i64 [ 0, %3 ], [ %26, %25 ]
%6 = trunc i64 %5 to i32
br label %7
7: ; preds = %4, %22
%8 = phi i64 [ 0, %4 ], [ %23, %22 ]
%9 = getelementptr inbounds [25 x double], [25 x double]* %0, i64 %5, i64 %8
%10 = load double, double* %9
%11 = fcmp une double %10, 0.000000e+00
br i1 %11, label %12, label %22
12: ; preds = %7
%13 = load i32, i32* %1
%14 = sext i32 %13 to i64
%15 = icmp slt i64 %5, %14
%16 = select i1 %15, i32 %13, i32 %6
store i32 %16, i32* %1
%17 = load i32, i32* %2
%18 = sext i32 %17 to i64
%19 = icmp slt i64 %8, %18
%20 = trunc i64 %8 to i32
%21 = select i1 %19, i32 %17, i32 %20
store i32 %21, i32* %2
br label %22
22: ; preds = %7, %12
%23 = add nuw nsw i64 %8, 1
%24 = icmp eq i64 %23, 25
br i1 %24, label %25, label %7
25: ; preds = %22
%26 = add nuw nsw i64 %5, 1
%27 = icmp eq i64 %26, 25
br i1 %27, label %28, label %4
28: ; preds = %25
ret i32 undef
}