void cgpost2l(float part[], float q[], float qm, int nop, int idimp, int nxv,
              int nyv) {
  int j, nn, mm, np, mp;
  float dxp, dyp, amx, amy;
  for (j = 0; j < nop; j++) {
    nn = part[idimp * j];
    mm = part[1 + idimp * j];
    dxp = qm * (part[idimp * j] - (float)nn);
    dyp = part[1 + idimp * j] - (float)mm;
    mm = nxv * mm;
    amx = qm - dxp;
    mp = mm + nxv;
    amy = 1.0 - dyp;
    np = nn + 1;
    q[np + mp] += dxp * dyp;
    q[nn + mp] += amx * dyp;
    q[np + mm] += dxp * amy;
    q[nn + mm] += amx * amy;
  }
  return;
}


x86 assembly (clang, -O0):
<cgpost2l>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    12aa <cgpost2l+0x1aa>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
movslq %ecx,%rcx
cvttss2si (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
cvttss2si (%rax,%rcx,4),%eax
mov    %eax,-0x30(%rbp)
movss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtsi2ssl -0x2c(%rbp),%xmm2
subss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtsi2ssl -0x30(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%eax
imul   -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
movss  -0x14(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x38(%rbp)
movss  -0x40(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe1e(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
movss  -0x3c(%rbp),%xmm0
movss  -0x40(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x44(%rbp),%xmm0
movss  -0x40(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x3c(%rbp),%xmm0
movss  -0x48(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x44(%rbp),%xmm0
movss  -0x48(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1126 <cgpost2l+0x26>
pop    %rbp
ret

x86 assembly (clang, -O2):
<cgpost2l>:
test   %edx,%edx
jle    11c4 <cgpost2l+0xc4>
movslq %ecx,%r9
mov    %edx,%r10d
add    $0x4,%rdi
shl    $0x2,%r9
movss  0xee2(%rip),%xmm1
xchg   %ax,%ax
movss  -0x4(%rdi),%xmm2
movss  (%rdi),%xmm3
cvttss2si %xmm3,%edx
cvttps2dq %xmm2,%xmm4
cvtdq2ps %xmm4,%xmm4
cvttss2si %xmm2,%r11d
subss  %xmm4,%xmm2
cvttps2dq %xmm3,%xmm4
cvtdq2ps %xmm4,%xmm4
mulss  %xmm0,%xmm2
subss  %xmm4,%xmm3
imul   %r8d,%edx
movaps %xmm0,%xmm4
lea    (%rdx,%r8,1),%ecx
movaps %xmm1,%xmm5
subss  %xmm2,%xmm4
subss  %xmm3,%xmm5
lea    0x1(%r11,%rcx,1),%eax
cltq
movaps %xmm2,%xmm6
mulss  %xmm3,%xmm6
addss  (%rsi,%rax,4),%xmm6
add    %r11d,%ecx
movss  %xmm6,(%rsi,%rax,4)
movslq %ecx,%rax
mulss  %xmm4,%xmm3
addss  (%rsi,%rax,4),%xmm3
movss  %xmm3,(%rsi,%rax,4)
lea    0x1(%rdx,%r11,1),%eax
cltq
mulss  %xmm5,%xmm2
addss  (%rsi,%rax,4),%xmm2
add    %r11d,%edx
movss  %xmm2,(%rsi,%rax,4)
movslq %edx,%rax
mulss  %xmm4,%xmm5
addss  (%rsi,%rax,4),%xmm5
movss  %xmm5,(%rsi,%rax,4)
add    %r9,%rdi
add    $0xffffffffffffffff,%r10
jne    1120 <cgpost2l+0x20>
ret

LLVM IR (-O0):
define void @cgpost2l(float* %0, float* %1, float %2, i32 %3, i32 %4, i32 %5, i32 %6) {
%8 = alloca float*
%9 = alloca float*
%10 = alloca float
%11 = alloca i32
%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 float
%21 = alloca float
%22 = alloca float
%23 = alloca float
store float* %0, float** %8
store float* %1, float** %9
store float %2, float* %10
store i32 %3, i32* %11
store i32 %4, i32* %12
store i32 %5, i32* %13
store i32 %6, i32* %14
store i32 0, i32* %15
br label %24
24: ; preds = %124, %7
%25 = load i32, i32* %15
%26 = load i32, i32* %11
%27 = icmp slt i32 %25, %26
br i1 %27, label %28, label %127
28: ; preds = %24
%29 = load float*, float** %8
%30 = load i32, i32* %12
%31 = load i32, i32* %15
%32 = mul nsw i32 %30, %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %29, i64 %33
%35 = load float, float* %34
%36 = fptosi float %35 to i32
store i32 %36, i32* %16
%37 = load float*, float** %8
%38 = load i32, i32* %12
%39 = load i32, i32* %15
%40 = mul nsw i32 %38, %39
%41 = add nsw i32 1, %40
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %37, i64 %42
%44 = load float, float* %43
%45 = fptosi float %44 to i32
store i32 %45, i32* %17
%46 = load float, float* %10
%47 = load float*, float** %8
%48 = load i32, i32* %12
%49 = load i32, i32* %15
%50 = mul nsw i32 %48, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %47, i64 %51
%53 = load float, float* %52
%54 = load i32, i32* %16
%55 = sitofp i32 %54 to float
%56 = fsub float %53, %55
%57 = fmul float %46, %56
store float %57, float* %20
%58 = load float*, float** %8
%59 = load i32, i32* %12
%60 = load i32, i32* %15
%61 = mul nsw i32 %59, %60
%62 = add nsw i32 1, %61
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds float, float* %58, i64 %63
%65 = load float, float* %64
%66 = load i32, i32* %17
%67 = sitofp i32 %66 to float
%68 = fsub float %65, %67
store float %68, float* %21
%69 = load i32, i32* %13
%70 = load i32, i32* %17
%71 = mul nsw i32 %69, %70
store i32 %71, i32* %17
%72 = load float, float* %10
%73 = load float, float* %20
%74 = fsub float %72, %73
store float %74, float* %22
%75 = load i32, i32* %17
%76 = load i32, i32* %13
%77 = add nsw i32 %75, %76
store i32 %77, i32* %19
%78 = load float, float* %21
%79 = fpext float %78 to double
%80 = fsub double 1.000000e+00, %79
%81 = fptrunc double %80 to float
store float %81, float* %23
%82 = load i32, i32* %16
%83 = add nsw i32 %82, 1
store i32 %83, i32* %18
%84 = load float, float* %20
%85 = load float, float* %21
%86 = load float*, float** %9
%87 = load i32, i32* %18
%88 = load i32, i32* %19
%89 = add nsw i32 %87, %88
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds float, float* %86, i64 %90
%92 = load float, float* %91
%93 = call float @llvm.fmuladd.f32(float %84, float %85, float %92)
store float %93, float* %91
%94 = load float, float* %22
%95 = load float, float* %21
%96 = load float*, float** %9
%97 = load i32, i32* %16
%98 = load i32, i32* %19
%99 = add nsw i32 %97, %98
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %96, i64 %100
%102 = load float, float* %101
%103 = call float @llvm.fmuladd.f32(float %94, float %95, float %102)
store float %103, float* %101
%104 = load float, float* %20
%105 = load float, float* %23
%106 = load float*, float** %9
%107 = load i32, i32* %18
%108 = load i32, i32* %17
%109 = add nsw i32 %107, %108
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds float, float* %106, i64 %110
%112 = load float, float* %111
%113 = call float @llvm.fmuladd.f32(float %104, float %105, float %112)
store float %113, float* %111
%114 = load float, float* %22
%115 = load float, float* %23
%116 = load float*, float** %9
%117 = load i32, i32* %16
%118 = load i32, i32* %17
%119 = add nsw i32 %117, %118
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %116, i64 %120
%122 = load float, float* %121
%123 = call float @llvm.fmuladd.f32(float %114, float %115, float %122)
store float %123, float* %121
br label %124
124: ; preds = %28
%125 = load i32, i32* %15
%126 = add nsw i32 %125, 1
store i32 %126, i32* %15
br label %24
127: ; preds = %24
ret void
}

LLVM IR (-O2):
define void @cgpost2l(float* nocapture %0, float* nocapture %1, float %2, i32 %3, i32 %4, i32 %5, i32 %6) {
%8 = icmp sgt i32 %3, 0
br i1 %8, label %9, label %54
9: ; preds = %7
%10 = sext i32 %4 to i64
%11 = zext i32 %3 to i64
br label %12
12: ; preds = %9, %12
%13 = phi i64 [ 0, %9 ], [ %52, %12 ]
%14 = mul nsw i64 %13, %10
%15 = getelementptr inbounds float, float* %0, i64 %14
%16 = load float, float* %15
%17 = fptosi float %16 to i32
%18 = add nsw i64 %14, 1
%19 = getelementptr inbounds float, float* %0, i64 %18
%20 = load float, float* %19
%21 = fptosi float %20 to i32
%22 = sitofp i32 %17 to float
%23 = fsub float %16, %22
%24 = fmul float %23, %2
%25 = sitofp i32 %21 to float
%26 = fsub float %20, %25
%27 = mul nsw i32 %21, %5
%28 = fsub float %2, %24
%29 = add nsw i32 %27, %5
%30 = fsub float 1.000000e+00, %26
%31 = add nsw i32 %17, 1
%32 = add nsw i32 %29, %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %1, i64 %33
%35 = load float, float* %34
%36 = tail call float @llvm.fmuladd.f32(float %24, float %26, float %35)
store float %36, float* %34
%37 = add nsw i32 %29, %17
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds float, float* %1, i64 %38
%40 = load float, float* %39
%41 = tail call float @llvm.fmuladd.f32(float %28, float %26, float %40)
store float %41, float* %39
%42 = add nsw i32 %27, %31
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %1, i64 %43
%45 = load float, float* %44
%46 = tail call float @llvm.fmuladd.f32(float %24, float %30, float %45)
store float %46, float* %44
%47 = add nsw i32 %27, %17
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %1, i64 %48
%50 = load float, float* %49
%51 = tail call float @llvm.fmuladd.f32(float %28, float %30, float %50)
store float %51, float* %49
%52 = add nuw nsw i64 %13, 1
%53 = icmp eq i64 %52, %11
br i1 %53, label %54, label %12
54: ; preds = %12, %7
ret void
}