#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 scaleGNdir(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.9 * cp[4 * ii + 3]) {
        alpha = MIN(alpha, -0.9 * cp[4 * ii + 3] / gnDir[2 * matInds[ii] + 1]);
      } else if (gnDir[2 * matInds[ii] + 1] > 9.0 * cp[4 * ii + 3]) {
        alpha = MIN(alpha, 9.0 * 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):
<scaleGNdir>:
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     12a9 <scaleGNdir+0x199>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x28(%rbp),%eax
jge    1296 <scaleGNdir+0x186>
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 <scaleGNdir+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    1291 <scaleGNdir+0x181>
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
movsd  0xddc(%rip),%xmm1
mulsd  (%rax,%rcx,8),%xmm1
ucomisd %xmm1,%xmm0
jbe    128c <scaleGNdir+0x17c>
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  0xdb0(%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    1291 <scaleGNdir+0x181>
jmp    1296 <scaleGNdir+0x186>
jmp    129b <scaleGNdir+0x18b>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1148 <scaleGNdir+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    1312 <scaleGNdir+0x202>
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    12d1 <scaleGNdir+0x1c1>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<scaleGNdir>:
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 <scaleGNdir+0x4c>
movsd  0xecd(%rip),%xmm0
lea    (%r14,%r14,1),%ecx
add    $0xfffffff8,%ecx
test   %ecx,%ecx
jle    12bb <scaleGNdir+0x1ab>
xorpd  0xed6(%rip),%xmm0
mov    %ecx,%eax
cmp    $0x4,%ecx
jae    11dd <scaleGNdir+0xcd>
xor    %ecx,%ecx
jmp    12a0 <scaleGNdir+0x190>
mov    %rcx,%r15
mov    %rdx,%r12
lea    -0x1(%r14),%ebx
add    $0xffffffffffffffff,%rbx
movsd  0xe8e(%rip),%xmm0
mov    $0x4,%ebp
jmp    1193 <scaleGNdir+0x83>
nopl   0x0(%rax)
divsd  %xmm2,%xmm1
call   1030 <MIN@plt>
add    $0x4,%rbp
add    $0xffffffffffffffff,%rbx
je     1133 <scaleGNdir+0x23>
movslq (%r12,%rbp,1),%rax
cmp    %r14d,%eax
jge    1189 <scaleGNdir+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 <scaleGNdir+0x70>
mulsd  0xe3f(%rip),%xmm3
ucomisd %xmm3,%xmm2
movapd %xmm3,%xmm1
ja     1180 <scaleGNdir+0x70>
jmp    1189 <scaleGNdir+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 <scaleGNdir+0x1ba>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
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 <scaleGNdir+0x100>
test   $0x1,%sil
je     128c <scaleGNdir+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 <scaleGNdir+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 <scaleGNdir+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 <scaleGNdir+0x158>
jmp    128c <scaleGNdir+0x17c>

LLVM IR (-O0):
define void @scaleGNdir(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 = %121, %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 %124
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 %120
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 -9.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 -9.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 %119
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 = fmul double 9.000000e+00, %92
%94 = fcmp ogt double %85, %93
br i1 %94, label %95, label %118
95: ; preds = %74
%96 = load double, double* %13
%97 = load double*, double** %9
%98 = load i32, i32* %11
%99 = mul nsw i32 4, %98
%100 = add nsw i32 %99, 3
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds double, double* %97, i64 %101
%103 = load double, double* %102
%104 = fmul double 9.000000e+00, %103
%105 = load double*, double** %10
%106 = load i32*, i32** %8
%107 = load i32, i32* %11
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds i32, i32* %106, i64 %108
%110 = load i32, i32* %109
%111 = mul nsw i32 2, %110
%112 = add nsw i32 %111, 1
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %105, i64 %113
%115 = load double, double* %114
%116 = fdiv double %104, %115
%117 = call double @MIN(double %96, double %116)
store double %117, double* %13
br label %118
118: ; preds = %95, %74
br label %119
119: ; preds = %118, %51
br label %120
120: ; preds = %119, %22
br label %121
121: ; preds = %120
%122 = load i32, i32* %11
%123 = add nsw i32 %122, 1
store i32 %123, i32* %11
br label %17
124: ; preds = %17
%125 = load double, double* %13
%126 = fneg double %125
store double %126, double* %13
store i32 0, i32* %11
br label %127
127: ; preds = %145, %124
%128 = load i32, i32* %11
%129 = load i32, i32* %12
%130 = sub nsw i32 %129, 4
%131 = mul nsw i32 2, %130
%132 = icmp slt i32 %128, %131
br i1 %132, label %133, label %148
133: ; preds = %127
%134 = load double, double* %13
%135 = load double*, double** %10
%136 = load i32, i32* %11
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %135, i64 %137
%139 = load double, double* %138
%140 = fmul double %134, %139
%141 = load double*, double** %10
%142 = load i32, i32* %11
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %141, i64 %143
store double %140, double* %144
br label %145
145: ; preds = %133
%146 = load i32, i32* %11
%147 = add nsw i32 %146, 1
store i32 %147, i32* %11
br label %127
148: ; preds = %127
ret void
}

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