void updateSolution_CFP32(float _Complex *mat, float _Complex *x,
                          float _Complex *r, int n, long long *Idxs) {
  int i, j, offset;
  float _Complex t;
  for (i = 0; i < n; ++i) {
    t = 0.0;
    offset = i * n;
    for (j = 0; j < n; ++j) {
      t += conj(mat[offset + j]) * r[j];
    }
    x[Idxs[i] - 1] += t;
  }
}


x86 assembly (clang, -O0):
<updateSolution_CFP32>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12cf <updateSolution_CFP32+0x1cf>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x2c(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1285 <updateSolution_CFP32+0x185>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,8),%xmm1
movss  0x4(%rax,%rcx,8),%xmm0
cvtss2sd %xmm1,%xmm1
movsd  %xmm1,-0x80(%rbp)
cvtss2sd %xmm0,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
movsd  %xmm2,-0x78(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,8),%xmm3
movss  0x4(%rax,%rcx,8),%xmm0
cvtss2sd %xmm3,%xmm4
movsd  %xmm4,-0x70(%rbp)
cvtss2sd %xmm0,%xmm5
movsd  %xmm5,-0x68(%rbp)
movaps %xmm1,%xmm0
mulsd  %xmm4,%xmm0
movaps %xmm2,%xmm3
mulsd  %xmm5,%xmm3
mulsd  %xmm5,%xmm1
mulsd  %xmm4,%xmm2
subsd  %xmm3,%xmm0
movsd  %xmm0,-0x60(%rbp)
addsd  %xmm2,%xmm1
movsd  %xmm1,-0x58(%rbp)
ucomisd %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  %xmm1,-0x48(%rbp)
jnp    1241 <updateSolution_CFP32+0x141>
movsd  -0x58(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm0
ucomisd %xmm1,%xmm1
movsd  %xmm0,-0x50(%rbp)
movsd  %xmm1,-0x48(%rbp)
jnp    1241 <updateSolution_CFP32+0x141>
movsd  -0x68(%rbp),%xmm3
movsd  -0x70(%rbp),%xmm2
movsd  -0x78(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm0
call   12e0 <__muldc3>
movsd  %xmm0,-0x50(%rbp)
movsd  %xmm1,-0x48(%rbp)
movsd  -0x50(%rbp),%xmm3
movsd  -0x48(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
movss  -0x3c(%rbp),%xmm0
cvtss2sd %xmm1,%xmm1
cvtss2sd %xmm0,%xmm0
addsd  %xmm3,%xmm1
addsd  %xmm2,%xmm0
cvtsd2ss %xmm1,%xmm1
cvtsd2ss %xmm0,%xmm0
movss  %xmm1,-0x40(%rbp)
movss  %xmm0,-0x3c(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1152 <updateSolution_CFP32+0x52>
movss  -0x40(%rbp),%xmm3
movss  -0x3c(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x2c(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,8),%xmm1
movss  0x4(%rax,%rcx,8),%xmm0
addss  %xmm3,%xmm1
addss  %xmm2,%xmm0
movss  %xmm1,(%rax,%rcx,8)
movss  %xmm0,0x4(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1125 <updateSolution_CFP32+0x25>
add    $0x80,%rsp
pop    %rbp
ret
nopl   0x0(%rax,%rax,1)

x86 assembly (clang, -O2):
<updateSolution_CFP32>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
test   %ecx,%ecx
jle    1248 <updateSolution_CFP32+0x148>
mov    %rdx,%r12
mov    %rsi,%r15
mov    %rdi,%rbx
mov    %ecx,%r13d
lea    0x0(,%r13,8),%rcx
xor    %r14d,%r14d
movapd 0xecb(%rip),%xmm5
jmp    1170 <updateSolution_CFP32+0x70>
nopw   0x0(%rax,%rax,1)
mov    (%r8,%r14,8),%rax
addss  -0x8(%r15,%rax,8),%xmm6
addss  -0x4(%r15,%rax,8),%xmm7
movss  %xmm6,-0x8(%r15,%rax,8)
movss  %xmm7,-0x4(%r15,%rax,8)
add    $0x1,%r14
add    %rcx,%rbx
cmp    %r13,%r14
je     1248 <updateSolution_CFP32+0x148>
xorps  %xmm6,%xmm6
xor    %ebp,%ebp
xorps  %xmm7,%xmm7
nopl   0x0(%rax,%rax,1)
cvtps2pd (%rbx,%rbp,8),%xmm0
cvtps2pd (%r12,%rbp,8),%xmm2
movapd %xmm0,%xmm1
mulpd  %xmm2,%xmm1
movapd %xmm1,%xmm3
unpckhpd %xmm1,%xmm3
addsd  %xmm1,%xmm3
movapd %xmm2,%xmm1
shufpd $0x1,%xmm2,%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,%xmm4
unpckhpd %xmm1,%xmm4
subsd  %xmm4,%xmm1
ucomisd %xmm3,%xmm3
jp     11ee <updateSolution_CFP32+0xee>
xorps  %xmm0,%xmm0
cvtss2sd %xmm6,%xmm0
addsd  %xmm3,%xmm0
xorps  %xmm2,%xmm2
cvtss2sd %xmm7,%xmm2
addsd  %xmm1,%xmm2
xorps  %xmm6,%xmm6
cvtsd2ss %xmm0,%xmm6
xorps  %xmm7,%xmm7
cvtsd2ss %xmm2,%xmm7
add    $0x1,%rbp
cmp    %rbp,%r13
jne    1180 <updateSolution_CFP32+0x80>
jmp    1140 <updateSolution_CFP32+0x40>
ucomisd %xmm1,%xmm1
jnp    11bc <updateSolution_CFP32+0xbc>
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
xorpd  %xmm5,%xmm1
movapd %xmm2,%xmm3
unpckhpd %xmm2,%xmm3
mov    %r8,0x10(%rsp)
mov    %rcx,0x8(%rsp)
movss  %xmm6,0x4(%rsp)
movss  %xmm7,(%rsp)
call   1260 <__muldc3>
movss  (%rsp),%xmm7
movss  0x4(%rsp),%xmm6
movapd 0xdcb(%rip),%xmm5
mov    0x8(%rsp),%rcx
mov    0x10(%rsp),%r8
movapd %xmm0,%xmm3
jmp    11bc <updateSolution_CFP32+0xbc>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
nopw   0x0(%rax,%rax,1)

LLVM IR (-O0):
define void @updateSolution_CFP32({ float, float }* %0, { float, float }* %1, { float, float }* %2, i32 %3, i64* %4) {
%6 = alloca { float, float }*
%7 = alloca { float, float }*
%8 = alloca { float, float }*
%9 = alloca i32
%10 = alloca i64*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca { float, float }
store { float, float }* %0, { float, float }** %6
store { float, float }* %1, { float, float }** %7
store { float, float }* %2, { float, float }** %8
store i32 %3, i32* %9
store i64* %4, i64** %10
store i32 0, i32* %11
br label %15
15: ; preds = %105, %5
%16 = load i32, i32* %11
%17 = load i32, i32* %9
%18 = icmp slt i32 %16, %17
br i1 %18, label %19, label %108
19: ; preds = %15
%20 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 0
%21 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 1
store float 0.000000e+00, float* %20
store float 0.000000e+00, float* %21
%22 = load i32, i32* %11
%23 = load i32, i32* %9
%24 = mul nsw i32 %22, %23
store i32 %24, i32* %13
store i32 0, i32* %12
br label %25
25: ; preds = %81, %19
%26 = load i32, i32* %12
%27 = load i32, i32* %9
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %84
29: ; preds = %25
%30 = load { float, float }*, { float, float }** %6
%31 = load i32, i32* %13
%32 = load i32, i32* %12
%33 = add nsw i32 %31, %32
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds { float, float }, { float, float }* %30, i64 %34
%36 = getelementptr inbounds { float, float }, { float, float }* %35, i32 0, i32 0
%37 = load float, float* %36
%38 = getelementptr inbounds { float, float }, { float, float }* %35, i32 0, i32 1
%39 = load float, float* %38
%40 = fpext float %37 to double
%41 = fpext float %39 to double
%42 = fneg double %41
%43 = load { float, float }*, { float, float }** %8
%44 = load i32, i32* %12
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds { float, float }, { float, float }* %43, i64 %45
%47 = getelementptr inbounds { float, float }, { float, float }* %46, i32 0, i32 0
%48 = load float, float* %47
%49 = getelementptr inbounds { float, float }, { float, float }* %46, i32 0, i32 1
%50 = load float, float* %49
%51 = fpext float %48 to double
%52 = fpext float %50 to double
%53 = fmul double %40, %51
%54 = fmul double %42, %52
%55 = fmul double %40, %52
%56 = fmul double %42, %51
%57 = fsub double %53, %54
%58 = fadd double %55, %56
%59 = fcmp uno double %57, %57
br i1 %59, label %60, label %66
60: ; preds = %29
%61 = fcmp uno double %58, %58
br i1 %61, label %62, label %66
62: ; preds = %60
%63 = call { double, double } @__muldc3(double %40, double %42, double %51, double %52)
%64 = extractvalue { double, double } %63, 0
%65 = extractvalue { double, double } %63, 1
br label %66
66: ; preds = %62, %60, %29
%67 = phi double [ %57, %29 ], [ %57, %60 ], [ %64, %62 ]
%68 = phi double [ %58, %29 ], [ %58, %60 ], [ %65, %62 ]
%69 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 0
%70 = load float, float* %69
%71 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 1
%72 = load float, float* %71
%73 = fpext float %70 to double
%74 = fpext float %72 to double
%75 = fadd double %73, %67
%76 = fadd double %74, %68
%77 = fptrunc double %75 to float
%78 = fptrunc double %76 to float
%79 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 0
%80 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 1
store float %77, float* %79
store float %78, float* %80
br label %81
81: ; preds = %66
%82 = load i32, i32* %12
%83 = add nsw i32 %82, 1
store i32 %83, i32* %12
br label %25
84: ; preds = %25
%85 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 0
%86 = load float, float* %85
%87 = getelementptr inbounds { float, float }, { float, float }* %14, i32 0, i32 1
%88 = load float, float* %87
%89 = load { float, float }*, { float, float }** %7
%90 = load i64*, i64** %10
%91 = load i32, i32* %11
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds i64, i64* %90, i64 %92
%94 = load i64, i64* %93
%95 = sub nsw i64 %94, 1
%96 = getelementptr inbounds { float, float }, { float, float }* %89, i64 %95
%97 = getelementptr inbounds { float, float }, { float, float }* %96, i32 0, i32 0
%98 = load float, float* %97
%99 = getelementptr inbounds { float, float }, { float, float }* %96, i32 0, i32 1
%100 = load float, float* %99
%101 = fadd float %98, %86
%102 = fadd float %100, %88
%103 = getelementptr inbounds { float, float }, { float, float }* %96, i32 0, i32 0
%104 = getelementptr inbounds { float, float }, { float, float }* %96, i32 0, i32 1
store float %101, float* %103
store float %102, float* %104
br label %105
105: ; preds = %84
%106 = load i32, i32* %11
%107 = add nsw i32 %106, 1
store i32 %107, i32* %11
br label %15
108: ; preds = %15
ret void
}

LLVM IR (-O2):
define void @updateSolution_CFP32({ float, float }* nocapture %0, { float, float }* nocapture %1, { float, float }* nocapture %2, i32 %3, i64* nocapture %4) {
%6 = icmp sgt i32 %3, 0
br i1 %6, label %7, label %70
7: ; preds = %5
%8 = zext i32 %3 to i64
%9 = zext i32 %3 to i64
br label %10
10: ; preds = %58, %7
%11 = phi i64 [ 0, %7 ], [ %68, %58 ]
%12 = mul nsw i64 %11, %8
br label %13
13: ; preds = %10, %47
%14 = phi i64 [ 0, %10 ], [ %56, %47 ]
%15 = phi float [ 0.000000e+00, %10 ], [ %54, %47 ]
%16 = phi float [ 0.000000e+00, %10 ], [ %55, %47 ]
%17 = add nuw nsw i64 %14, %12
%18 = getelementptr inbounds { float, float }, { float, float }* %0, i64 %17, i32 0
%19 = bitcast float* %18 to <2 x float>*
%20 = load <2 x float>, <2 x float>* %19
%21 = fpext <2 x float> %20 to <2 x double>
%22 = extractelement <2 x double> %21, i64 1
%23 = fneg double %22
%24 = getelementptr inbounds { float, float }, { float, float }* %2, i64 %14, i32 0
%25 = bitcast float* %24 to <2 x float>*
%26 = load <2 x float>, <2 x float>* %25
%27 = fpext <2 x float> %26 to <2 x double>
%28 = fmul <2 x double> %21, %27
%29 = shufflevector <2 x double> %28, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%30 = fadd <2 x double> %28, %29
%31 = extractelement <2 x double> %30, i64 0
%32 = shufflevector <2 x double> %27, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%33 = fmul <2 x double> %32, %21
%34 = shufflevector <2 x double> %33, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%35 = fsub <2 x double> %33, %34
%36 = extractelement <2 x double> %35, i64 0
%37 = fcmp uno double %31, 0.000000e+00
br i1 %37, label %38, label %47
38: ; preds = %13
%39 = fcmp uno double %36, 0.000000e+00
br i1 %39, label %40, label %47
40: ; preds = %38
%41 = extractelement <2 x double> %21, i64 0
%42 = extractelement <2 x double> %27, i64 0
%43 = extractelement <2 x double> %27, i64 1
%44 = tail call { double, double } @__muldc3(double %41, double %23, double %42, double %43)
%45 = extractvalue { double, double } %44, 0
%46 = extractvalue { double, double } %44, 1
br label %47
47: ; preds = %40, %38, %13
%48 = phi double [ %31, %13 ], [ %31, %38 ], [ %45, %40 ]
%49 = phi double [ %36, %13 ], [ %36, %38 ], [ %46, %40 ]
%50 = fpext float %15 to double
%51 = fpext float %16 to double
%52 = fadd double %48, %50
%53 = fadd double %49, %51
%54 = fptrunc double %52 to float
%55 = fptrunc double %53 to float
%56 = add nuw nsw i64 %14, 1
%57 = icmp eq i64 %56, %9
br i1 %57, label %58, label %13
58: ; preds = %47
%59 = getelementptr inbounds i64, i64* %4, i64 %11
%60 = load i64, i64* %59
%61 = add nsw i64 %60, -1
%62 = getelementptr inbounds { float, float }, { float, float }* %1, i64 %61, i32 0
%63 = load float, float* %62
%64 = getelementptr inbounds { float, float }, { float, float }* %1, i64 %61, i32 1
%65 = load float, float* %64
%66 = fadd float %63, %54
%67 = fadd float %65, %55
store float %66, float* %62
store float %67, float* %64
%68 = add nuw nsw i64 %11, 1
%69 = icmp eq i64 %68, %9
br i1 %69, label %70, label %10
70: ; preds = %58, %5
ret void
}