#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ free(double *);
double *initAminoAcidMass();
void _computeFragmentIons(int *n_, double *W_, double *b_, double *y_) {
  int i;
  double b;
  double y;
  double C_term;
  double N_term;
  double Electron;
  double Hydrogen;
  double *M;
  M = initAminoAcidMass();
  C_term = 17.002740;
  N_term = 1.007825;
  Electron = 0.000549;
  Hydrogen = 1.007825;
  b = N_term - Electron;
  y = C_term + N_term + Hydrogen - Electron;
  for (i = 0; i < *n_; i++) {
    y += W_[i];
  }
  for (i = 0; i < *n_; i++) {
    b += W_[i];
    b_[i] = b;
    y_[*n_ - i - 1] = y;
    y -= W_[i];
  }
  free(M);
}


x86 assembly (clang, -O0):
<_computeFragmentIons>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    $0x0,%al
call   1040 <initAminoAcidMass@plt>
mov    %rax,-0x60(%rbp)
movsd  0xec5(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xea8(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xea3(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xe8e(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x8(%rbp),%rcx
cmp    (%rcx),%eax
jge    11da <_computeFragmentIons+0xba>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a6 <_computeFragmentIons+0x86>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x8(%rbp),%rcx
cmp    (%rcx),%eax
jge    125f <_computeFragmentIons+0x13f>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    (%rcx),%ecx
sub    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x38(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11e1 <_computeFragmentIons+0xc1>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<_computeFragmentIons>:
push   %r15
push   %r14
push   %r12
push   %rbx
push   %rax
mov    %rcx,%r14
mov    %rdx,%r15
mov    %rsi,%rbx
mov    %rdi,%r12
xor    %eax,%eax
call   1040 <initAminoAcidMass@plt>
mov    (%r12),%r8d
test   %r8d,%r8d
jle    1273 <_computeFragmentIons+0x153>
lea    -0x1(%r8),%rdx
mov    %r8d,%esi
and    $0x7,%esi
cmp    $0x7,%rdx
jae    1164 <_computeFragmentIons+0x44>
movsd  0xea0(%rip),%xmm0
xor    %edi,%edi
jmp    11b8 <_computeFragmentIons+0x98>
mov    %r8d,%ecx
and    $0xfffffff8,%ecx
movsd  0xe8e(%rip),%xmm0
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addsd  (%rbx,%rdi,8),%xmm0
addsd  0x8(%rbx,%rdi,8),%xmm0
addsd  0x10(%rbx,%rdi,8),%xmm0
addsd  0x18(%rbx,%rdi,8),%xmm0
addsd  0x20(%rbx,%rdi,8),%xmm0
addsd  0x28(%rbx,%rdi,8),%xmm0
addsd  0x30(%rbx,%rdi,8),%xmm0
addsd  0x38(%rbx,%rdi,8),%xmm0
add    $0x8,%rdi
cmp    %rdi,%rcx
jne    1180 <_computeFragmentIons+0x60>
test   %rsi,%rsi
je     11de <_computeFragmentIons+0xbe>
lea    (%rbx,%rdi,8),%rcx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rcx,%rdi,8),%xmm0
add    $0x1,%rdi
cmp    %rdi,%rsi
jne    11d0 <_computeFragmentIons+0xb0>
test   %r8d,%r8d
jle    1273 <_computeFragmentIons+0x153>
test   %rdx,%rdx
je     1286 <_computeFragmentIons+0x166>
mov    %r8d,%esi
and    $0xfffffffe,%esi
lea    -0x2(%r8),%edi
movsd  0xe06(%rip),%xmm1
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addsd  (%rbx,%rdx,8),%xmm1
movsd  %xmm1,(%r15,%rdx,8)
lea    0x1(%rdi),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%r14,%rcx,8)
subsd  (%rbx,%rdx,8),%xmm0
addsd  0x8(%rbx,%rdx,8),%xmm1
movsd  %xmm1,0x8(%r15,%rdx,8)
movslq %edi,%rdi
movsd  %xmm0,(%r14,%rdi,8)
subsd  0x8(%rbx,%rdx,8),%xmm0
add    $0x2,%rdx
add    $0xfffffffe,%edi
cmp    %rdx,%rsi
jne    1210 <_computeFragmentIons+0xf0>
test   $0x1,%r8b
je     1273 <_computeFragmentIons+0x153>
addsd  (%rbx,%rdx,8),%xmm1
movsd  %xmm1,(%r15,%rdx,8)
not    %edx
add    %r8d,%edx
movslq %edx,%rcx
movsd  %xmm0,(%r14,%rcx,8)
mov    %rax,%rdi
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
jmp    1030 <free@plt>
movsd  0xd7a(%rip),%xmm1
xor    %edx,%edx
test   $0x1,%r8b
jne    125a <_computeFragmentIons+0x13a>
jmp    1273 <_computeFragmentIons+0x153>

LLVM IR (-O0):
define void @_computeFragmentIons(i32* %0, double* %1, double* %2, double* %3) {
%5 = alloca i32*
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double*
store i32* %0, i32** %5
store double* %1, double** %6
store double* %2, double** %7
store double* %3, double** %8
%17 = call double* (...) @initAminoAcidMass()
store double* %17, double** %16
store double 0x403100B39192641B, double* %12
store double 0x3FF0200D1B71758E, double* %13
store double 5.490000e-04, double* %14
store double 0x3FF0200D1B71758E, double* %15
%18 = load double, double* %13
%19 = load double, double* %14
%20 = fsub double %18, %19
store double %20, double* %10
%21 = load double, double* %12
%22 = load double, double* %13
%23 = fadd double %21, %22
%24 = load double, double* %15
%25 = fadd double %23, %24
%26 = load double, double* %14
%27 = fsub double %25, %26
store double %27, double* %11
store i32 0, i32* %9
br label %28
28: ; preds = %41, %4
%29 = load i32, i32* %9
%30 = load i32*, i32** %5
%31 = load i32, i32* %30
%32 = icmp slt i32 %29, %31
br i1 %32, label %33, label %44
33: ; preds = %28
%34 = load double*, double** %6
%35 = load i32, i32* %9
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
%39 = load double, double* %11
%40 = fadd double %39, %38
store double %40, double* %11
br label %41
41: ; preds = %33
%42 = load i32, i32* %9
%43 = add nsw i32 %42, 1
store i32 %43, i32* %9
br label %28
44: ; preds = %28
store i32 0, i32* %9
br label %45
45: ; preds = %79, %44
%46 = load i32, i32* %9
%47 = load i32*, i32** %5
%48 = load i32, i32* %47
%49 = icmp slt i32 %46, %48
br i1 %49, label %50, label %82
50: ; preds = %45
%51 = load double*, double** %6
%52 = load i32, i32* %9
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %51, i64 %53
%55 = load double, double* %54
%56 = load double, double* %10
%57 = fadd double %56, %55
store double %57, double* %10
%58 = load double, double* %10
%59 = load double*, double** %7
%60 = load i32, i32* %9
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %59, i64 %61
store double %58, double* %62
%63 = load double, double* %11
%64 = load double*, double** %8
%65 = load i32*, i32** %5
%66 = load i32, i32* %65
%67 = load i32, i32* %9
%68 = sub nsw i32 %66, %67
%69 = sub nsw i32 %68, 1
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %64, i64 %70
store double %63, double* %71
%72 = load double*, double** %6
%73 = load i32, i32* %9
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
%76 = load double, double* %75
%77 = load double, double* %11
%78 = fsub double %77, %76
store double %78, double* %11
br label %79
79: ; preds = %50
%80 = load i32, i32* %9
%81 = add nsw i32 %80, 1
store i32 %81, i32* %9
br label %45
82: ; preds = %45
%83 = load double*, double** %16
%84 = call i32 @free(double* %83)
ret void
}

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