#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 INFINITY;
void getD_2(double *KmatOut, double *KmatDiag, int *K, int *Nin, int *Nout,
            double *g, double *amount_points, double *connection_strength,
            double *d, double *err, double *summed_diameter) {
  int i, j, k;
  *err = 0;
  for (i = 0; i < *Nout; i++) {
    double minval = INFINITY;
    int kstar = 0;
    bool first = true;
    for (k = 0; k < *K; k++) {
      double summa = 0;
      for (j = 0; j < *Nin; j++)
        summa += connection_strength[k + *K * j] * KmatOut[i + *Nout * j];
      double h = -2.0 * summa / amount_points[k];
      d[k + *K * i] = KmatDiag[i] + h + g[k];
      double val = d[k + *K * i];
      if (val < minval || first) {
        first = false;
        kstar = k;
        minval = val;
      }
    }
    *err += minval;
    summed_diameter[kstar] += minval;
  }
}


x86 assembly (clang, -O0):
<getD_2>:
push   %rbp
mov    %rsp,%rbp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x28(%rbp),%rcx
cmp    (%rcx),%eax
jge    12f0 <getD_2+0x1f0>
mov    0x2e90(%rip),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x4c(%rbp)
movl   $0x1,-0x50(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x18(%rbp),%rcx
cmp    (%rcx),%eax
jge    12ba <getD_2+0x1ba>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
jge    11f9 <getD_2+0xf9>
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%rdx
mov    (%rdx),%edx
imul   -0x38(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x28(%rbp),%rdx
mov    (%rdx),%edx
imul   -0x38(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1194 <getD_2+0x94>
movsd  0xdff(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x60(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%rdx
mov    (%rdx),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%rdx
mov    (%rdx),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
ja     1290 <getD_2+0x190>
cmpl   $0x0,-0x50(%rbp)
je     12a7 <getD_2+0x1a7>
movl   $0x0,-0x50(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    12ac <getD_2+0x1ac>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1176 <getD_2+0x76>
movsd  -0x48(%rbp),%xmm0
mov    0x28(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
movsd  -0x48(%rbp),%xmm0
mov    0x30(%rbp),%rax
movslq -0x4c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1142 <getD_2+0x42>
pop    %rbp
ret

x86 assembly (clang, -O2):
<getD_2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9,-0x10(%rsp)
mov    %rcx,-0x28(%rsp)
mov    %rsi,-0x18(%rsp)
mov    %rdi,-0x40(%rsp)
mov    0x50(%rsp),%rax
movq   $0x0,(%rax)
mov    (%r8),%r13d
test   %r13d,%r13d
jle    12eb <getD_2+0x1eb>
mov    (%rdx),%eax
mov    %rax,-0x30(%rsp)
movslq %eax,%rbx
mov    %rbx,%r8
shl    $0x4,%r8
mov    %r13,%rax
shl    $0x4,%rax
xor    %r9d,%r9d
mov    -0x40(%rsp),%rcx
mov    %rcx,-0x38(%rsp)
jmp    119d <getD_2+0x9d>
nopl   (%rax)
movapd %xmm0,%xmm2
xor    %r10d,%r10d
mov    0x50(%rsp),%rcx
movsd  (%rcx),%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,(%rcx)
movslq %r10d,%rcx
mov    0x58(%rsp),%rdx
addsd  (%rdx,%rcx,8),%xmm2
movsd  %xmm2,(%rdx,%rcx,8)
add    $0x1,%r9
addq   $0x8,-0x38(%rsp)
cmp    %r13,%r9
je     12eb <getD_2+0x1eb>
mov    0x2e44(%rip),%rcx
movsd  (%rcx),%xmm0
cmpl   $0x0,-0x30(%rsp)
jle    1160 <getD_2+0x60>
mov    -0x28(%rsp),%rcx
mov    (%rcx),%r14d
mov    %r9,%rcx
imul   %rbx,%rcx
mov    %rcx,-0x8(%rsp)
mov    %r14d,%ecx
and    $0xfffffffe,%ecx
xor    %r11d,%r11d
mov    $0x1,%dil
mov    0x40(%rsp),%rsi
xor    %r10d,%r10d
mov    %r9,-0x20(%rsp)
jmp    11fd <getD_2+0xfd>
xchg   %ax,%ax
cmovne %r11d,%r10d
add    $0x1,%r11
add    $0x8,%rsi
xor    %edi,%edi
movapd %xmm2,%xmm0
cmp    -0x30(%rsp),%r11
je     1167 <getD_2+0x67>
xorpd  %xmm1,%xmm1
test   %r14d,%r14d
jle    1292 <getD_2+0x192>
cmp    $0x1,%r14d
jne    1220 <getD_2+0x120>
xor    %r12d,%r12d
jmp    125d <getD_2+0x15d>
cs nopw 0x0(%rax,%rax,1)
nop
mov    -0x38(%rsp),%r9
mov    %rsi,%r15
xor    %r12d,%r12d
nopl   0x0(%rax,%rax,1)
movsd  (%r15),%xmm2
mulsd  (%r9),%xmm2
addsd  %xmm1,%xmm2
movsd  (%r15,%rbx,8),%xmm1
mulsd  (%r9,%r13,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x2,%r12
add    %r8,%r15
add    %rax,%r9
cmp    %r12,%rcx
jne    1230 <getD_2+0x130>
test   $0x1,%r14b
mov    -0x20(%rsp),%r9
je     1292 <getD_2+0x192>
mov    %r12,%rbp
imul   %rbx,%rbp
add    %r11,%rbp
mov    0x40(%rsp),%rdx
movsd  (%rdx,%rbp,8),%xmm2
imul   %r13,%r12
add    %r9,%r12
mov    -0x40(%rsp),%rdx
mulsd  (%rdx,%r12,8),%xmm2
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm1
mov    0x38(%rsp),%rdx
divsd  (%rdx,%r11,8),%xmm1
mov    -0x18(%rsp),%rdx
movsd  (%rdx,%r9,8),%xmm2
subsd  %xmm1,%xmm2
mov    -0x10(%rsp),%rdx
addsd  (%rdx,%r11,8),%xmm2
mov    -0x8(%rsp),%rdx
lea    (%rdx,%r11,1),%rbp
mov    0x48(%rsp),%rdx
movsd  %xmm2,(%rdx,%rbp,8)
ucomisd %xmm2,%xmm0
seta   %dl
or     %dl,%dil
test   $0x1,%dil
jne    11e0 <getD_2+0xe0>
movapd %xmm0,%xmm2
jmp    11e0 <getD_2+0xe0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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