void cfwpminmx2(float qe[], float qbme, float *wpmax, float *wpmin, int nx,
                int ny, int nxe, int nye) {
  int j, k;
  float tpmax, tpmin, at1;
  tpmax = qbme * qe[0];
  tpmin = tpmax;
  for (k = 0; k < ny; k++) {
    for (j = 0; j < nx; j++) {
      at1 = qbme * qe[j + nxe * k];
      tpmax = at1 > tpmax ? at1 : tpmax;
      tpmin = at1 < tpmin ? at1 : tpmin;
    }
  }
  *wpmax = tpmax;
  *wpmin = tpmin;
  return;
}


x86 assembly (clang, -O0):
<cfwpminmx2>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8d,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movss  -0xc(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    1210 <cfwpminmx2+0x110>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11fd <cfwpminmx2+0xfd>
movss  -0xc(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
ucomiss -0x38(%rbp),%xmm0
jbe    11a5 <cfwpminmx2+0xa5>
movss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
jmp    11af <cfwpminmx2+0xaf>
movss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x40(%rbp),%xmm1
movss  -0x3c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    11db <cfwpminmx2+0xdb>
movss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
jmp    11e5 <cfwpminmx2+0xe5>
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1159 <cfwpminmx2+0x59>
jmp    1202 <cfwpminmx2+0x102>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1146 <cfwpminmx2+0x46>
movss  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x3c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<cfwpminmx2>:
push   %r15
push   %r14
push   %r12
push   %rbx
movss  (%rdi),%xmm1
mulss  %xmm0,%xmm1
test   %r8d,%r8d
jle    11c0 <cfwpminmx2+0xc0>
movslq %r9d,%r9
mov    %r8d,%r10d
mov    %ecx,%r8d
mov    %r8d,%r11d
and    $0xfffffffe,%r11d
lea    0x4(%rdi),%r12
lea    0x0(,%r9,4),%r14
xor    %r15d,%r15d
movaps %xmm1,%xmm2
jmp    114c <cfwpminmx2+0x4c>
nopl   0x0(%rax)
add    $0x1,%r15
add    %r14,%r12
cmp    %r10,%r15
je     11c3 <cfwpminmx2+0xc3>
test   %ecx,%ecx
jle    1140 <cfwpminmx2+0x40>
xor    %ebx,%ebx
cmp    $0x1,%ecx
je     1194 <cfwpminmx2+0x94>
nopw   0x0(%rax,%rax,1)
movss  -0x4(%r12,%rbx,4),%xmm3
mulss  %xmm0,%xmm3
movaps %xmm3,%xmm4
maxss  %xmm2,%xmm4
minss  %xmm1,%xmm3
movss  (%r12,%rbx,4),%xmm1
mulss  %xmm0,%xmm1
movaps %xmm1,%xmm2
maxss  %xmm4,%xmm2
minss  %xmm3,%xmm1
add    $0x2,%rbx
cmp    %rbx,%r11
jne    1160 <cfwpminmx2+0x60>
test   $0x1,%r8b
je     1140 <cfwpminmx2+0x40>
mov    %r15,%rax
imul   %r9,%rax
add    %rax,%rbx
movss  (%rdi,%rbx,4),%xmm3
mulss  %xmm0,%xmm3
movaps %xmm3,%xmm4
maxss  %xmm2,%xmm4
minss  %xmm1,%xmm3
movaps %xmm4,%xmm2
movaps %xmm3,%xmm1
jmp    1140 <cfwpminmx2+0x40>
movaps %xmm1,%xmm2
movss  %xmm2,(%rsi)
movss  %xmm1,(%rdx)
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @cfwpminmx2(float* %0, float %1, float* %2, float* %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = alloca float*
%10 = alloca float
%11 = alloca float*
%12 = alloca float*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca float
%20 = alloca float
%21 = alloca float
store float* %0, float** %9
store float %1, float* %10
store float* %2, float** %11
store float* %3, float** %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store i32 %7, i32* %16
%22 = load float, float* %10
%23 = load float*, float** %9
%24 = getelementptr inbounds float, float* %23, i64 0
%25 = load float, float* %24
%26 = fmul float %22, %25
store float %26, float* %19
%27 = load float, float* %19
store float %27, float* %20
store i32 0, i32* %18
br label %28
28: ; preds = %71, %8
%29 = load i32, i32* %18
%30 = load i32, i32* %14
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %74
32: ; preds = %28
store i32 0, i32* %17
br label %33
33: ; preds = %67, %32
%34 = load i32, i32* %17
%35 = load i32, i32* %13
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %70
37: ; preds = %33
%38 = load float, float* %10
%39 = load float*, float** %9
%40 = load i32, i32* %17
%41 = load i32, i32* %15
%42 = load i32, i32* %18
%43 = mul nsw i32 %41, %42
%44 = add nsw i32 %40, %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds float, float* %39, i64 %45
%47 = load float, float* %46
%48 = fmul float %38, %47
store float %48, float* %21
%49 = load float, float* %21
%50 = load float, float* %19
%51 = fcmp ogt float %49, %50
br i1 %51, label %52, label %54
52: ; preds = %37
%53 = load float, float* %21
br label %56
54: ; preds = %37
%55 = load float, float* %19
br label %56
56: ; preds = %54, %52
%57 = phi float [ %53, %52 ], [ %55, %54 ]
store float %57, float* %19
%58 = load float, float* %21
%59 = load float, float* %20
%60 = fcmp olt float %58, %59
br i1 %60, label %61, label %63
61: ; preds = %56
%62 = load float, float* %21
br label %65
63: ; preds = %56
%64 = load float, float* %20
br label %65
65: ; preds = %63, %61
%66 = phi float [ %62, %61 ], [ %64, %63 ]
store float %66, float* %20
br label %67
67: ; preds = %65
%68 = load i32, i32* %17
%69 = add nsw i32 %68, 1
store i32 %69, i32* %17
br label %33
70: ; preds = %33
br label %71
71: ; preds = %70
%72 = load i32, i32* %18
%73 = add nsw i32 %72, 1
store i32 %73, i32* %18
br label %28
74: ; preds = %28
%75 = load float, float* %19
%76 = load float*, float** %11
store float %75, float* %76
%77 = load float, float* %20
%78 = load float*, float** %12
store float %77, float* %78
ret void
}

LLVM IR (-O2):
define void @cfwpminmx2(float* nocapture %0, float %1, float* nocapture writeonly %2, float* nocapture writeonly %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = load float, float* %0
%10 = fmul float %9, %1
%11 = icmp sgt i32 %5, 0
br i1 %11, label %12, label %72
12: ; preds = %8
%13 = icmp sgt i32 %4, 0
%14 = sext i32 %6 to i64
%15 = zext i32 %5 to i64
%16 = zext i32 %4 to i64
%17 = and i64 %16, 1
%18 = icmp eq i32 %4, 1
%19 = and i64 %16, 4294967294
%20 = icmp eq i64 %17, 0
br label %21
21: ; preds = %12, %67
%22 = phi i64 [ 0, %12 ], [ %70, %67 ]
%23 = phi float [ %10, %12 ], [ %69, %67 ]
%24 = phi float [ %10, %12 ], [ %68, %67 ]
br i1 %13, label %25, label %67
25: ; preds = %21
%26 = mul nsw i64 %22, %14
br i1 %18, label %52, label %27
27: ; preds = %25, %27
%28 = phi i64 [ %49, %27 ], [ 0, %25 ]
%29 = phi float [ %48, %27 ], [ %23, %25 ]
%30 = phi float [ %46, %27 ], [ %24, %25 ]
%31 = phi i64 [ %50, %27 ], [ 0, %25 ]
%32 = add nsw i64 %28, %26
%33 = getelementptr inbounds float, float* %0, i64 %32
%34 = load float, float* %33
%35 = fmul float %34, %1
%36 = fcmp ogt float %35, %30
%37 = select i1 %36, float %35, float %30
%38 = fcmp olt float %35, %29
%39 = select i1 %38, float %35, float %29
%40 = or i64 %28, 1
%41 = add nsw i64 %40, %26
%42 = getelementptr inbounds float, float* %0, i64 %41
%43 = load float, float* %42
%44 = fmul float %43, %1
%45 = fcmp ogt float %44, %37
%46 = select i1 %45, float %44, float %37
%47 = fcmp olt float %44, %39
%48 = select i1 %47, float %44, float %39
%49 = add nuw nsw i64 %28, 2
%50 = add i64 %31, 2
%51 = icmp eq i64 %50, %19
br i1 %51, label %52, label %27
52: ; preds = %27, %25
%53 = phi float [ undef, %25 ], [ %46, %27 ]
%54 = phi float [ undef, %25 ], [ %48, %27 ]
%55 = phi i64 [ 0, %25 ], [ %49, %27 ]
%56 = phi float [ %23, %25 ], [ %48, %27 ]
%57 = phi float [ %24, %25 ], [ %46, %27 ]
br i1 %20, label %67, label %58
58: ; preds = %52
%59 = add nsw i64 %55, %26
%60 = getelementptr inbounds float, float* %0, i64 %59
%61 = load float, float* %60
%62 = fmul float %61, %1
%63 = fcmp ogt float %62, %57
%64 = select i1 %63, float %62, float %57
%65 = fcmp olt float %62, %56
%66 = select i1 %65, float %62, float %56
br label %67
67: ; preds = %58, %52, %21
%68 = phi float [ %24, %21 ], [ %53, %52 ], [ %64, %58 ]
%69 = phi float [ %23, %21 ], [ %54, %52 ], [ %66, %58 ]
%70 = add nuw nsw i64 %22, 1
%71 = icmp eq i64 %70, %15
br i1 %71, label %72, label %21
72: ; preds = %67, %8
%73 = phi float [ %10, %8 ], [ %68, %67 ]
%74 = phi float [ %10, %8 ], [ %69, %67 ]
store float %73, float* %2
store float %74, float* %3
ret void
}