void normalise_fluxes(double *s, int Nreact, double **locked, int n_locked,
                      double *lock_value) {
  double *dummy, **dummy_l, Z = 0., Z_l = 0.;
  for (dummy = s; dummy < s + Nreact; dummy++)
    Z += *dummy;
  for (dummy_l = locked; dummy_l < locked + n_locked; dummy_l++) {
    Z -= **dummy_l;
    Z_l += **dummy_l;
  }
  for (dummy = s; dummy < s + Nreact; dummy++) {
    *dummy *= ((double)Nreact - Z_l) / Z;
  }
  dummy = lock_value;
  for (dummy_l = locked; dummy_l < locked + n_locked; dummy_l++) {
    **dummy_l = *dummy;
    dummy++;
  }
}


x86 assembly (clang, -O0):
<normalise_fluxes>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0xc(%rbp),%rdx
shl    $0x3,%rdx
add    %rdx,%rcx
cmp    %rcx,%rax
jae    116d <normalise_fluxes+0x6d>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x30(%rbp)
jmp    112e <normalise_fluxes+0x2e>
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
shl    $0x3,%rdx
add    %rdx,%rcx
cmp    %rcx,%rax
jae    11d0 <normalise_fluxes+0xd0>
mov    -0x38(%rbp),%rax
mov    (%rax),%rax
movsd  (%rax),%xmm1
movsd  -0x40(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%rax
mov    (%rax),%rax
movsd  (%rax),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x38(%rbp)
jmp    1175 <normalise_fluxes+0x75>
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0xc(%rbp),%rdx
shl    $0x3,%rdx
add    %rdx,%rcx
cmp    %rcx,%rax
jae    1220 <normalise_fluxes+0x120>
cvtsi2sdl -0xc(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x30(%rbp)
jmp    11d8 <normalise_fluxes+0xd8>
mov    -0x28(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
shl    $0x3,%rdx
add    %rdx,%rcx
cmp    %rcx,%rax
jae    127c <normalise_fluxes+0x17c>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    (%rax),%rax
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x30(%rbp)
mov    -0x38(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x38(%rbp)
jmp    1230 <normalise_fluxes+0x130>
pop    %rbp
ret

x86 assembly (clang, -O2):
<normalise_fluxes>:
push   %r15
push   %r14
push   %rbx
movslq %esi,%rax
lea    (%rdi,%rax,8),%r9
xorpd  %xmm2,%xmm2
xorpd  %xmm1,%xmm1
test   %eax,%eax
jle    112d <normalise_fluxes+0x2d>
xorpd  %xmm1,%xmm1
mov    %rdi,%rax
nop
addsd  (%rax),%xmm1
add    $0x8,%rax
cmp    %r9,%rax
jb     1120 <normalise_fluxes+0x20>
movslq %ecx,%rbx
lea    (%rdx,%rbx,8),%r15
test   %ebx,%ebx
jle    1158 <normalise_fluxes+0x58>
xorpd  %xmm2,%xmm2
mov    %rdx,%rbx
nop
mov    (%rbx),%rax
movsd  (%rax),%xmm0
subsd  %xmm0,%xmm1
addsd  %xmm0,%xmm2
add    $0x8,%rbx
cmp    %r15,%rbx
jb     1140 <normalise_fluxes+0x40>
test   %esi,%esi
jle    1265 <normalise_fluxes+0x165>
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
subsd  %xmm2,%xmm0
lea    0x8(%rdi),%rsi
cmp    %rsi,%r9
cmova  %r9,%rsi
divsd  %xmm1,%xmm0
mov    %rdi,%r11
not    %r11
add    %rsi,%r11
cmp    $0x18,%r11
jb     1250 <normalise_fluxes+0x150>
shr    $0x3,%r11
add    $0x1,%r11
mov    %r11,%r10
and    $0xfffffffffffffffc,%r10
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r10),%rsi
mov    %rsi,%r14
shr    $0x2,%r14
add    $0x1,%r14
test   %rsi,%rsi
je     128f <normalise_fluxes+0x18f>
mov    %r14,%rbx
and    $0xfffffffffffffffe,%rbx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movupd (%rdi,%rsi,8),%xmm2
movupd 0x10(%rdi,%rsi,8),%xmm3
movupd 0x20(%rdi,%rsi,8),%xmm4
movupd 0x30(%rdi,%rsi,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rdi,%rsi,8)
movupd %xmm3,0x10(%rdi,%rsi,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,0x20(%rdi,%rsi,8)
movupd %xmm5,0x30(%rdi,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rbx
jne    11d0 <normalise_fluxes+0xd0>
test   $0x1,%r14b
je     123c <normalise_fluxes+0x13c>
movupd (%rdi,%rsi,8),%xmm2
movupd 0x10(%rdi,%rsi,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rdi,%rsi,8)
movupd %xmm3,0x10(%rdi,%rsi,8)
cmp    %r10,%r11
je     1265 <normalise_fluxes+0x165>
lea    (%rdi,%r10,8),%rdi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdi),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi)
add    $0x8,%rdi
cmp    %r9,%rdi
jb     1250 <normalise_fluxes+0x150>
test   %ecx,%ecx
jle    1289 <normalise_fluxes+0x189>
nopl   0x0(%rax)
movsd  (%r8),%xmm0
mov    (%rdx),%rax
movsd  %xmm0,(%rax)
add    $0x8,%r8
add    $0x8,%rdx
cmp    %r15,%rdx
jb     1270 <normalise_fluxes+0x170>
pop    %rbx
pop    %r14
pop    %r15
ret
xor    %esi,%esi
test   $0x1,%r14b
jne    121e <normalise_fluxes+0x11e>
jmp    123c <normalise_fluxes+0x13c>

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

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