void updateSolutionEconomic_FP64(float *mat, double *x, double *r, int n,
                                 long long *Idxs) {
  int i, j, offset1, offset2;
  double t = r[n - 1];
  for (i = 0; i < n - 1; ++i) {
    t -= mat[i] * r[i];
  }
  t *= mat[n - 1];
  offset1 = n;
  offset2 = 2 * n - 1;
  x[Idxs[n - 1] - 1] += t;
  for (i = 0; i < n - 1; ++i) {
    x[Idxs[i] - 1] += mat[i + offset1] * (r[i] - mat[i + offset2] * t);
  }
}


x86 assembly (clang, -O0):
<updateSolutionEconomic_FP64>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    119b <updateSolutionEconomic_FP64+0x9b>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1135 <updateSolutionEconomic_FP64+0x35>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    -0x1c(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
mov    (%rcx,%rdx,8),%rcx
sub    $0x1,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    129e <updateSolutionEconomic_FP64+0x19e>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
movsd  -0x40(%rbp),%xmm3
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x2c(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11fc <updateSolutionEconomic_FP64+0xfc>
pop    %rbp
ret

x86 assembly (clang, -O2):
<updateSolutionEconomic_FP64>:
movslq %ecx,%rax
lea    -0x1(%rax),%r10
movsd  -0x8(%rdx,%rax,8),%xmm0
cmp    $0x2,%eax
jl     1177 <updateSolutionEconomic_FP64+0x77>
mov    %r10d,%r9d
cmp    $0x1,%r10d
jne    111f <updateSolutionEconomic_FP64+0x1f>
xor    %eax,%eax
jmp    115f <updateSolutionEconomic_FP64+0x5f>
mov    %r9d,%r11d
and    $0xfffffffe,%r11d
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movss  (%rdi,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  (%rdx,%rax,8),%xmm1
movss  0x4(%rdi,%rax,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  0x8(%rdx,%rax,8),%xmm2
subsd  %xmm1,%xmm0
subsd  %xmm2,%xmm0
add    $0x2,%rax
cmp    %rax,%r11
jne    1130 <updateSolutionEconomic_FP64+0x30>
test   $0x1,%r9b
je     1177 <updateSolutionEconomic_FP64+0x77>
movss  (%rdi,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  (%rdx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movss  (%rdi,%r10,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm0,%xmm1
mov    (%r8,%r10,8),%rax
movsd  -0x8(%rsi,%rax,8),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rsi,%rax,8)
cmp    $0x2,%ecx
jl     11fd <updateSolutionEconomic_FP64+0xfd>
mov    %ecx,%eax
add    %ecx,%ecx
movslq %ecx,%rcx
mov    %r10d,%r9d
lea    (%rdi,%rcx,4),%r10
add    $0xfffffffffffffffc,%r10
lea    (%rdi,%rax,4),%rdi
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
movss  (%rdi,%rax,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  (%rdx,%rax,8),%xmm2
movss  (%r10,%rax,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm1,%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm0,%xmm2
mov    (%r8,%rax,8),%rcx
addsd  -0x8(%rsi,%rcx,8),%xmm2
movsd  %xmm2,-0x8(%rsi,%rcx,8)
add    $0x1,%rax
cmp    %rax,%r9
jne    11c0 <updateSolutionEconomic_FP64+0xc0>
ret

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

LLVM IR (-O2):
define void @updateSolutionEconomic_FP64(float* nocapture %0, double* nocapture %1, double* nocapture %2, i32 %3, i64* nocapture %4) {
%6 = add nsw i32 %3, -1
%7 = sext i32 %6 to i64
%8 = getelementptr inbounds double, double* %2, i64 %7
%9 = load double, double* %8
%10 = icmp sgt i32 %3, 1
br i1 %10, label %11, label %52
11: ; preds = %5
%12 = zext i32 %6 to i64
%13 = and i64 %12, 1
%14 = icmp eq i32 %6, 1
br i1 %14, label %39, label %15
15: ; preds = %11
%16 = and i64 %12, 4294967294
br label %17
17: ; preds = %17, %15
%18 = phi i64 [ 0, %15 ], [ %36, %17 ]
%19 = phi double [ %9, %15 ], [ %35, %17 ]
%20 = phi i64 [ 0, %15 ], [ %37, %17 ]
%21 = getelementptr inbounds float, float* %0, i64 %18
%22 = load float, float* %21
%23 = fpext float %22 to double
%24 = getelementptr inbounds double, double* %2, i64 %18
%25 = load double, double* %24
%26 = fneg double %23
%27 = tail call double @llvm.fmuladd.f64(double %26, double %25, double %19)
%28 = or i64 %18, 1
%29 = getelementptr inbounds float, float* %0, i64 %28
%30 = load float, float* %29
%31 = fpext float %30 to double
%32 = getelementptr inbounds double, double* %2, i64 %28
%33 = load double, double* %32
%34 = fneg double %31
%35 = tail call double @llvm.fmuladd.f64(double %34, double %33, double %27)
%36 = add nuw nsw i64 %18, 2
%37 = add i64 %20, 2
%38 = icmp eq i64 %37, %16
br i1 %38, label %39, label %17
39: ; preds = %17, %11
%40 = phi double [ undef, %11 ], [ %35, %17 ]
%41 = phi i64 [ 0, %11 ], [ %36, %17 ]
%42 = phi double [ %9, %11 ], [ %35, %17 ]
%43 = icmp eq i64 %13, 0
br i1 %43, label %52, label %44
44: ; preds = %39
%45 = getelementptr inbounds float, float* %0, i64 %41
%46 = load float, float* %45
%47 = fpext float %46 to double
%48 = getelementptr inbounds double, double* %2, i64 %41
%49 = load double, double* %48
%50 = fneg double %47
%51 = tail call double @llvm.fmuladd.f64(double %50, double %49, double %42)
br label %52
52: ; preds = %44, %39, %5
%53 = phi double [ %9, %5 ], [ %40, %39 ], [ %51, %44 ]
%54 = getelementptr inbounds float, float* %0, i64 %7
%55 = load float, float* %54
%56 = fpext float %55 to double
%57 = fmul double %53, %56
%58 = getelementptr inbounds i64, i64* %4, i64 %7
%59 = load i64, i64* %58
%60 = add nsw i64 %59, -1
%61 = getelementptr inbounds double, double* %1, i64 %60
%62 = load double, double* %61
%63 = fadd double %62, %57
store double %63, double* %61
br i1 %10, label %64, label %92
64: ; preds = %52
%65 = shl nuw nsw i32 %3, 1
%66 = add nsw i32 %65, -1
%67 = zext i32 %3 to i64
%68 = sext i32 %66 to i64
%69 = zext i32 %6 to i64
br label %70
70: ; preds = %64, %70
%71 = phi i64 [ 0, %64 ], [ %90, %70 ]
%72 = add nuw nsw i64 %71, %67
%73 = getelementptr inbounds float, float* %0, i64 %72
%74 = load float, float* %73
%75 = fpext float %74 to double
%76 = getelementptr inbounds double, double* %2, i64 %71
%77 = load double, double* %76
%78 = add nsw i64 %71, %68
%79 = getelementptr inbounds float, float* %0, i64 %78
%80 = load float, float* %79
%81 = fpext float %80 to double
%82 = fneg double %81
%83 = tail call double @llvm.fmuladd.f64(double %82, double %57, double %77)
%84 = getelementptr inbounds i64, i64* %4, i64 %71
%85 = load i64, i64* %84
%86 = add nsw i64 %85, -1
%87 = getelementptr inbounds double, double* %1, i64 %86
%88 = load double, double* %87
%89 = tail call double @llvm.fmuladd.f64(double %75, double %83, double %88)
store double %89, double* %87
%90 = add nuw nsw i64 %71, 1
%91 = icmp eq i64 %90, %69
br i1 %91, label %92, label %70
92: ; preds = %70, %52
ret void
}