void scalecomplex(float *rdata, float *idata, float rmult, float imult,
                  long npts) {
  long count;
  float rval, ival;
  for (count = 0; count < npts; count++) {
    rval = (*rdata) * rmult - (*idata) * imult;
    ival = (*rdata) * imult + (*idata) * rmult;
    *rdata++ = rval;
    *idata++ = ival;
  }
  return;
}


x86 assembly (clang, -O0):
<scalecomplex>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
movss  %xmm1,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    11cc <scalecomplex+0xcc>
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  -0x14(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm1
mulss  -0x18(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  -0x18(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm1
mulss  -0x14(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  %xmm0,(%rax)
movss  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x10(%rbp)
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1122 <scalecomplex+0x22>
pop    %rbp
ret

x86 assembly (clang, -O2):
<scalecomplex>:
test   %rdx,%rdx
jle    12c1 <scalecomplex+0x1c1>
cmp    $0x4,%rdx
jb     1129 <scalecomplex+0x29>
lea    (%rsi,%rdx,4),%rax
cmp    %rdi,%rax
jbe    11f6 <scalecomplex+0xf6>
lea    (%rdi,%rdx,4),%rax
cmp    %rsi,%rax
jbe    11f6 <scalecomplex+0xf6>
xor    %r9d,%r9d
mov    %r9,%rax
not    %rax
test   $0x1,%dl
je     1171 <scalecomplex+0x71>
movss  (%rdi),%xmm2
movss  (%rsi),%xmm3
movaps %xmm3,%xmm4
mulss  %xmm1,%xmm4
movaps %xmm2,%xmm5
mulss  %xmm0,%xmm5
subss  %xmm4,%xmm5
mulss  %xmm0,%xmm3
mulss  %xmm1,%xmm2
addss  %xmm3,%xmm2
movss  %xmm5,(%rdi)
add    $0x4,%rdi
movss  %xmm2,(%rsi)
add    $0x4,%rsi
or     $0x1,%r9
add    %rdx,%rax
je     12c1 <scalecomplex+0x1c1>
sub    %r9,%rdx
xor    %eax,%eax
nop
movss  (%rdi,%rax,4),%xmm2
movss  (%rsi,%rax,4),%xmm3
movaps %xmm3,%xmm4
mulss  %xmm1,%xmm4
movaps %xmm2,%xmm5
mulss  %xmm0,%xmm5
subss  %xmm4,%xmm5
mulss  %xmm0,%xmm3
mulss  %xmm1,%xmm2
addss  %xmm3,%xmm2
movss  %xmm5,(%rdi,%rax,4)
movss  %xmm2,(%rsi,%rax,4)
movss  0x4(%rdi,%rax,4),%xmm2
movss  0x4(%rsi,%rax,4),%xmm3
movaps %xmm3,%xmm4
mulss  %xmm1,%xmm4
movaps %xmm2,%xmm5
mulss  %xmm0,%xmm5
subss  %xmm4,%xmm5
mulss  %xmm0,%xmm3
mulss  %xmm1,%xmm2
addss  %xmm3,%xmm2
movss  %xmm5,0x4(%rdi,%rax,4)
movss  %xmm2,0x4(%rsi,%rax,4)
add    $0x2,%rax
cmp    %rax,%rdx
jne    1180 <scalecomplex+0x80>
jmp    12c1 <scalecomplex+0x1c1>
mov    %rdx,%r9
and    $0xfffffffffffffffc,%r9
movaps %xmm1,%xmm2
shufps $0x0,%xmm1,%xmm2
movaps %xmm0,%xmm3
shufps $0x0,%xmm0,%xmm3
lea    -0x4(%r9),%rax
mov    %rax,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rax,%rax
je     12cf <scalecomplex+0x1cf>
mov    %r8,%rax
and    $0xfffffffffffffffe,%rax
xor    %ecx,%ecx
nopl   0x0(%rax)
movups (%rdi,%rcx,4),%xmm4
movups (%rsi,%rcx,4),%xmm5
movaps %xmm2,%xmm6
mulps  %xmm5,%xmm6
movaps %xmm4,%xmm7
mulps  %xmm3,%xmm7
subps  %xmm6,%xmm7
mulps  %xmm3,%xmm5
mulps  %xmm2,%xmm4
addps  %xmm5,%xmm4
movups %xmm7,(%rdi,%rcx,4)
movups %xmm4,(%rsi,%rcx,4)
movups 0x10(%rdi,%rcx,4),%xmm4
movups 0x10(%rsi,%rcx,4),%xmm5
movaps %xmm2,%xmm6
mulps  %xmm5,%xmm6
movaps %xmm4,%xmm7
mulps  %xmm3,%xmm7
subps  %xmm6,%xmm7
mulps  %xmm3,%xmm5
mulps  %xmm2,%xmm4
addps  %xmm5,%xmm4
movups %xmm7,0x10(%rdi,%rcx,4)
movups %xmm4,0x10(%rsi,%rcx,4)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rax
jne    1230 <scalecomplex+0x130>
test   $0x1,%r8b
je     12bc <scalecomplex+0x1bc>
movups (%rdi,%rcx,4),%xmm4
movups (%rsi,%rcx,4),%xmm5
movaps %xmm2,%xmm6
mulps  %xmm5,%xmm6
movaps %xmm4,%xmm7
mulps  %xmm3,%xmm7
subps  %xmm6,%xmm7
mulps  %xmm3,%xmm5
mulps  %xmm2,%xmm4
addps  %xmm5,%xmm4
movups %xmm7,(%rdi,%rcx,4)
movups %xmm4,(%rsi,%rcx,4)
cmp    %rdx,%r9
jne    12c2 <scalecomplex+0x1c2>
ret
lea    (%rdi,%r9,4),%rdi
lea    (%rsi,%r9,4),%rsi
jmp    112c <scalecomplex+0x2c>
xor    %ecx,%ecx
test   $0x1,%r8b
jne    1294 <scalecomplex+0x194>
jmp    12bc <scalecomplex+0x1bc>

LLVM IR (-O0):
define void @scalecomplex(float* %0, float* %1, float %2, float %3, i64 %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca float
%9 = alloca float
%10 = alloca i64
%11 = alloca i64
%12 = alloca float
%13 = alloca float
store float* %0, float** %6
store float* %1, float** %7
store float %2, float* %8
store float %3, float* %9
store i64 %4, i64* %10
store i64 0, i64* %11
br label %14
14: ; preds = %42, %5
%15 = load i64, i64* %11
%16 = load i64, i64* %10
%17 = icmp slt i64 %15, %16
br i1 %17, label %18, label %45
18: ; preds = %14
%19 = load float*, float** %6
%20 = load float, float* %19
%21 = load float, float* %8
%22 = load float*, float** %7
%23 = load float, float* %22
%24 = load float, float* %9
%25 = fmul float %23, %24
%26 = fneg float %25
%27 = call float @llvm.fmuladd.f32(float %20, float %21, float %26)
store float %27, float* %12
%28 = load float*, float** %6
%29 = load float, float* %28
%30 = load float, float* %9
%31 = load float*, float** %7
%32 = load float, float* %31
%33 = load float, float* %8
%34 = fmul float %32, %33
%35 = call float @llvm.fmuladd.f32(float %29, float %30, float %34)
store float %35, float* %13
%36 = load float, float* %12
%37 = load float*, float** %6
%38 = getelementptr inbounds float, float* %37, i32 1
store float* %38, float** %6
store float %36, float* %37
%39 = load float, float* %13
%40 = load float*, float** %7
%41 = getelementptr inbounds float, float* %40, i32 1
store float* %41, float** %7
store float %39, float* %40
br label %42
42: ; preds = %18
%43 = load i64, i64* %11
%44 = add nsw i64 %43, 1
store i64 %44, i64* %11
br label %14
45: ; preds = %14
ret void
}

LLVM IR (-O2):
define void @scalecomplex(float* nocapture %0, float* nocapture %1, float %2, float %3, i64 %4) {
%6 = icmp sgt i64 %4, 0
br i1 %6, label %7, label %130
7: ; preds = %5
%8 = icmp ult i64 %4, 4
br i1 %8, label %82, label %9
9: ; preds = %7
%10 = getelementptr float, float* %0, i64 %4
%11 = getelementptr float, float* %1, i64 %4
%12 = icmp ugt float* %11, %0
%13 = icmp ugt float* %10, %1
%14 = and i1 %12, %13
br i1 %14, label %82, label %15
15: ; preds = %9
%16 = and i64 %4, -4
%17 = getelementptr float, float* %0, i64 %16
%18 = getelementptr float, float* %1, i64 %16
%19 = insertelement <4 x float> poison, float %3, i64 0
%20 = shufflevector <4 x float> %19, <4 x float> poison, <4 x i32> zeroinitializer
%21 = insertelement <4 x float> poison, float %2, i64 0
%22 = shufflevector <4 x float> %21, <4 x float> poison, <4 x i32> zeroinitializer
%23 = add i64 %16, -4
%24 = lshr exact i64 %23, 2
%25 = add nuw nsw i64 %24, 1
%26 = and i64 %25, 1
%27 = icmp eq i64 %23, 0
br i1 %27, label %63, label %28
28: ; preds = %15
%29 = and i64 %25, 9223372036854775806
br label %30
30: ; preds = %30, %28
%31 = phi i64 [ 0, %28 ], [ %60, %30 ]
%32 = phi i64 [ 0, %28 ], [ %61, %30 ]
%33 = getelementptr float, float* %0, i64 %31
%34 = getelementptr float, float* %1, i64 %31
%35 = bitcast float* %33 to <4 x float>*
%36 = load <4 x float>, <4 x float>* %35
%37 = bitcast float* %34 to <4 x float>*
%38 = load <4 x float>, <4 x float>* %37
%39 = fneg <4 x float> %38
%40 = fmul <4 x float> %20, %39
%41 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %36, <4 x float> %22, <4 x float> %40)
%42 = fmul <4 x float> %38, %22
%43 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %36, <4 x float> %20, <4 x float> %42)
%44 = bitcast float* %33 to <4 x float>*
store <4 x float> %41, <4 x float>* %44
%45 = bitcast float* %34 to <4 x float>*
store <4 x float> %43, <4 x float>* %45
%46 = or i64 %31, 4
%47 = getelementptr float, float* %0, i64 %46
%48 = getelementptr float, float* %1, i64 %46
%49 = bitcast float* %47 to <4 x float>*
%50 = load <4 x float>, <4 x float>* %49
%51 = bitcast float* %48 to <4 x float>*
%52 = load <4 x float>, <4 x float>* %51
%53 = fneg <4 x float> %52
%54 = fmul <4 x float> %20, %53
%55 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %50, <4 x float> %22, <4 x float> %54)
%56 = fmul <4 x float> %52, %22
%57 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %50, <4 x float> %20, <4 x float> %56)
%58 = bitcast float* %47 to <4 x float>*
store <4 x float> %55, <4 x float>* %58
%59 = bitcast float* %48 to <4 x float>*
store <4 x float> %57, <4 x float>* %59
%60 = add nuw i64 %31, 8
%61 = add i64 %32, 2
%62 = icmp eq i64 %61, %29
br i1 %62, label %63, label %30
63: ; preds = %30, %15
%64 = phi i64 [ 0, %15 ], [ %60, %30 ]
%65 = icmp eq i64 %26, 0
br i1 %65, label %80, label %66
66: ; preds = %63
%67 = getelementptr float, float* %0, i64 %64
%68 = getelementptr float, float* %1, i64 %64
%69 = bitcast float* %67 to <4 x float>*
%70 = load <4 x float>, <4 x float>* %69
%71 = bitcast float* %68 to <4 x float>*
%72 = load <4 x float>, <4 x float>* %71
%73 = fneg <4 x float> %72
%74 = fmul <4 x float> %20, %73
%75 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %70, <4 x float> %22, <4 x float> %74)
%76 = fmul <4 x float> %72, %22
%77 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %70, <4 x float> %20, <4 x float> %76)
%78 = bitcast float* %67 to <4 x float>*
store <4 x float> %75, <4 x float>* %78
%79 = bitcast float* %68 to <4 x float>*
store <4 x float> %77, <4 x float>* %79
br label %80
80: ; preds = %63, %66
%81 = icmp eq i64 %16, %4
br i1 %81, label %130, label %82
82: ; preds = %9, %7, %80
%83 = phi float* [ %0, %9 ], [ %0, %7 ], [ %17, %80 ]
%84 = phi float* [ %1, %9 ], [ %1, %7 ], [ %18, %80 ]
%85 = phi i64 [ 0, %9 ], [ 0, %7 ], [ %16, %80 ]
%86 = xor i64 %85, -1
%87 = and i64 %4, 1
%88 = icmp eq i64 %87, 0
br i1 %88, label %100, label %89
89: ; preds = %82
%90 = load float, float* %83
%91 = load float, float* %84
%92 = fneg float %91
%93 = fmul float %92, %3
%94 = tail call float @llvm.fmuladd.f32(float %90, float %2, float %93)
%95 = fmul float %91, %2
%96 = tail call float @llvm.fmuladd.f32(float %90, float %3, float %95)
%97 = getelementptr inbounds float, float* %83, i64 1
store float %94, float* %83
%98 = getelementptr inbounds float, float* %84, i64 1
store float %96, float* %84
%99 = or i64 %85, 1
br label %100
100: ; preds = %89, %82
%101 = phi float* [ %83, %82 ], [ %97, %89 ]
%102 = phi float* [ %84, %82 ], [ %98, %89 ]
%103 = phi i64 [ %85, %82 ], [ %99, %89 ]
%104 = sub i64 0, %4
%105 = icmp eq i64 %86, %104
br i1 %105, label %130, label %106
106: ; preds = %100, %106
%107 = phi float* [ %126, %106 ], [ %101, %100 ]
%108 = phi float* [ %127, %106 ], [ %102, %100 ]
%109 = phi i64 [ %128, %106 ], [ %103, %100 ]
%110 = load float, float* %107
%111 = load float, float* %108
%112 = fneg float %111
%113 = fmul float %112, %3
%114 = tail call float @llvm.fmuladd.f32(float %110, float %2, float %113)
%115 = fmul float %111, %2
%116 = tail call float @llvm.fmuladd.f32(float %110, float %3, float %115)
%117 = getelementptr inbounds float, float* %107, i64 1
store float %114, float* %107
%118 = getelementptr inbounds float, float* %108, i64 1
store float %116, float* %108
%119 = load float, float* %117
%120 = load float, float* %118
%121 = fneg float %120
%122 = fmul float %121, %3
%123 = tail call float @llvm.fmuladd.f32(float %119, float %2, float %122)
%124 = fmul float %120, %2
%125 = tail call float @llvm.fmuladd.f32(float %119, float %3, float %124)
%126 = getelementptr inbounds float, float* %107, i64 2
store float %123, float* %117
%127 = getelementptr inbounds float, float* %108, i64 2
store float %125, float* %118
%128 = add nuw nsw i64 %109, 2
%129 = icmp eq i64 %128, %4
br i1 %129, label %130, label %106
130: ; preds = %100, %106, %80, %5
ret void
}