#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 */
double MIN(double, double);
void scaleGNdirCare(int ncp, int kcp, int *matInds, double *cp, double *gnDir) {
  int ii, numCP;
  double alpha = 1.0;
  numCP = ncp * kcp;
  for (ii = 1; ii <= (numCP - 2); ii++) {
    if (matInds[ii] < numCP) {
      if (gnDir[2 * matInds[ii] + 1] < -0.5 * cp[4 * ii + 3]) {
        alpha = MIN(alpha, -0.5 * cp[4 * ii + 3] / gnDir[2 * matInds[ii] + 1]);
      } else if (gnDir[2 * matInds[ii] + 1] > cp[4 * ii + 3]) {
        alpha = MIN(alpha, cp[4 * ii + 3] / gnDir[2 * matInds[ii] + 1]);
      }
    }
  }
  alpha = -alpha;
  for (ii = 0; ii < 2 * (numCP - 4); ii++) {
    gnDir[ii] = alpha * gnDir[ii];
  }
  return;
}


x86 assembly (clang, -O0):
<scaleGNdirCare>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movsd  0xece(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x4(%rbp),%eax
imul   -0x8(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jg     1295 <scaleGNdirCare+0x185>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x28(%rbp),%eax
jge    1282 <scaleGNdirCare+0x172>
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  0xe66(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    11ff <scaleGNdirCare+0xef>
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  0xe3a(%rip),%xmm1
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm1
call   1030 <MIN@plt>
movsd  %xmm0,-0x30(%rbp)
jmp    127d <scaleGNdirCare+0x16d>
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1278 <scaleGNdirCare+0x168>
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm1
call   1030 <MIN@plt>
movsd  %xmm0,-0x30(%rbp)
jmp    127d <scaleGNdirCare+0x16d>
jmp    1282 <scaleGNdirCare+0x172>
jmp    1287 <scaleGNdirCare+0x177>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1148 <scaleGNdirCare+0x38>
movsd  -0x30(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    $0x4,%ecx
shl    $0x1,%ecx
cmp    %ecx,%eax
jge    12fe <scaleGNdirCare+0x1ee>
movsd  -0x30(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12bd <scaleGNdirCare+0x1ad>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<scaleGNdirCare>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r8,%r13
mov    %edi,%r14d
imul   %esi,%r14d
cmp    $0x3,%r14d
jge    115c <scaleGNdirCare+0x4c>
movsd  0xecd(%rip),%xmm0
lea    (%r14,%r14,1),%ecx
add    $0xfffffff8,%ecx
test   %ecx,%ecx
jle    12bb <scaleGNdirCare+0x1ab>
xorpd  0xec6(%rip),%xmm0
mov    %ecx,%eax
cmp    $0x4,%ecx
jae    11d5 <scaleGNdirCare+0xc5>
xor    %ecx,%ecx
jmp    12a0 <scaleGNdirCare+0x190>
mov    %rcx,%r15
mov    %rdx,%r12
lea    -0x1(%r14),%ebx
add    $0xffffffffffffffff,%rbx
movsd  0xe8e(%rip),%xmm0
mov    $0x4,%ebp
jmp    1193 <scaleGNdirCare+0x83>
nopl   0x0(%rax)
divsd  %xmm2,%xmm1
call   1030 <MIN@plt>
add    $0x4,%rbp
add    $0xffffffffffffffff,%rbx
je     1133 <scaleGNdirCare+0x23>
movslq (%r12,%rbp,1),%rax
cmp    %r14d,%eax
jge    1189 <scaleGNdirCare+0x79>
shl    $0x4,%rax
movsd  0x8(%r13,%rax,1),%xmm2
mov    %ebp,%eax
shl    $0x3,%rax
or     $0x18,%rax
movsd  (%r15,%rax,1),%xmm3
movapd %xmm3,%xmm1
mulsd  0xe45(%rip),%xmm1
ucomisd %xmm2,%xmm1
ja     1180 <scaleGNdirCare+0x70>
ucomisd %xmm3,%xmm2
movapd %xmm3,%xmm1
ja     1180 <scaleGNdirCare+0x70>
jmp    1189 <scaleGNdirCare+0x79>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rcx),%rdx
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
test   %rdx,%rdx
je     12ca <scaleGNdirCare+0x1ba>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd 0x0(%r13,%rdx,8),%xmm2
movupd 0x10(%r13,%rdx,8),%xmm3
movupd 0x20(%r13,%rdx,8),%xmm4
movupd 0x30(%r13,%rdx,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x0(%r13,%rdx,8)
movupd %xmm3,0x10(%r13,%rdx,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,0x20(%r13,%rdx,8)
movupd %xmm5,0x30(%r13,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1210 <scaleGNdirCare+0x100>
test   $0x1,%sil
je     128c <scaleGNdirCare+0x17c>
movupd 0x0(%r13,%rdx,8),%xmm2
movupd 0x10(%r13,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,0x0(%r13,%rdx,8)
movupd %xmm3,0x10(%r13,%rdx,8)
cmp    %rax,%rcx
je     12bb <scaleGNdirCare+0x1ab>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    12a0 <scaleGNdirCare+0x190>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%sil
jne    1268 <scaleGNdirCare+0x158>
jmp    128c <scaleGNdirCare+0x17c>

LLVM IR (-O0):
define void @scaleGNdirCare(i32 %0, i32 %1, i32* %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32*
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32* %2, i32** %8
store double* %3, double** %9
store double* %4, double** %10
store double 1.000000e+00, double* %13
%14 = load i32, i32* %6
%15 = load i32, i32* %7
%16 = mul nsw i32 %14, %15
store i32 %16, i32* %12
store i32 1, i32* %11
br label %17
17: ; preds = %119, %5
%18 = load i32, i32* %11
%19 = load i32, i32* %12
%20 = sub nsw i32 %19, 2
%21 = icmp sle i32 %18, %20
br i1 %21, label %22, label %122
22: ; preds = %17
%23 = load i32*, i32** %8
%24 = load i32, i32* %11
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds i32, i32* %23, i64 %25
%27 = load i32, i32* %26
%28 = load i32, i32* %12
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %118
30: ; preds = %22
%31 = load double*, double** %10
%32 = load i32*, i32** %8
%33 = load i32, i32* %11
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds i32, i32* %32, i64 %34
%36 = load i32, i32* %35
%37 = mul nsw i32 2, %36
%38 = add nsw i32 %37, 1
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %31, i64 %39
%41 = load double, double* %40
%42 = load double*, double** %9
%43 = load i32, i32* %11
%44 = mul nsw i32 4, %43
%45 = add nsw i32 %44, 3
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %42, i64 %46
%48 = load double, double* %47
%49 = fmul double -5.000000e-01, %48
%50 = fcmp olt double %41, %49
br i1 %50, label %51, label %74
51: ; preds = %30
%52 = load double, double* %13
%53 = load double*, double** %9
%54 = load i32, i32* %11
%55 = mul nsw i32 4, %54
%56 = add nsw i32 %55, 3
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %53, i64 %57
%59 = load double, double* %58
%60 = fmul double -5.000000e-01, %59
%61 = load double*, double** %10
%62 = load i32*, i32** %8
%63 = load i32, i32* %11
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds i32, i32* %62, i64 %64
%66 = load i32, i32* %65
%67 = mul nsw i32 2, %66
%68 = add nsw i32 %67, 1
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %61, i64 %69
%71 = load double, double* %70
%72 = fdiv double %60, %71
%73 = call double @MIN(double %52, double %72)
store double %73, double* %13
br label %117
74: ; preds = %30
%75 = load double*, double** %10
%76 = load i32*, i32** %8
%77 = load i32, i32* %11
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds i32, i32* %76, i64 %78
%80 = load i32, i32* %79
%81 = mul nsw i32 2, %80
%82 = add nsw i32 %81, 1
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds double, double* %75, i64 %83
%85 = load double, double* %84
%86 = load double*, double** %9
%87 = load i32, i32* %11
%88 = mul nsw i32 4, %87
%89 = add nsw i32 %88, 3
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %86, i64 %90
%92 = load double, double* %91
%93 = fcmp ogt double %85, %92
br i1 %93, label %94, label %116
94: ; preds = %74
%95 = load double, double* %13
%96 = load double*, double** %9
%97 = load i32, i32* %11
%98 = mul nsw i32 4, %97
%99 = add nsw i32 %98, 3
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %96, i64 %100
%102 = load double, double* %101
%103 = load double*, double** %10
%104 = load i32*, i32** %8
%105 = load i32, i32* %11
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds i32, i32* %104, i64 %106
%108 = load i32, i32* %107
%109 = mul nsw i32 2, %108
%110 = add nsw i32 %109, 1
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds double, double* %103, i64 %111
%113 = load double, double* %112
%114 = fdiv double %102, %113
%115 = call double @MIN(double %95, double %114)
store double %115, double* %13
br label %116
116: ; preds = %94, %74
br label %117
117: ; preds = %116, %51
br label %118
118: ; preds = %117, %22
br label %119
119: ; preds = %118
%120 = load i32, i32* %11
%121 = add nsw i32 %120, 1
store i32 %121, i32* %11
br label %17
122: ; preds = %17
%123 = load double, double* %13
%124 = fneg double %123
store double %124, double* %13
store i32 0, i32* %11
br label %125
125: ; preds = %143, %122
%126 = load i32, i32* %11
%127 = load i32, i32* %12
%128 = sub nsw i32 %127, 4
%129 = mul nsw i32 2, %128
%130 = icmp slt i32 %126, %129
br i1 %130, label %131, label %146
131: ; preds = %125
%132 = load double, double* %13
%133 = load double*, double** %10
%134 = load i32, i32* %11
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds double, double* %133, i64 %135
%137 = load double, double* %136
%138 = fmul double %132, %137
%139 = load double*, double** %10
%140 = load i32, i32* %11
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %139, i64 %141
store double %138, double* %142
br label %143
143: ; preds = %131
%144 = load i32, i32* %11
%145 = add nsw i32 %144, 1
store i32 %145, i32* %11
br label %125
146: ; preds = %125
ret void
}

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