float detect_orbital_period(
    float timestamp[], float series[], int series_length, float min_period_days,
    float max_period_days, float increment_days, float min_dipped_density,
    float max_dipped_percent, float min_intermediate_percent,
    float max_intermediate_percent, float expected_dip_radius_percent,
    float peak_threshold, float max_vacancy_density, float dip_threshold) {
  float period_days = 0;
  float max_response = 0;
  int expected_width = (int)(256 * expected_dip_radius_percent / 100.0f);
  int max_dipped = (int)(256 * max_dipped_percent / 100.0f);
  int max_intermediates = (int)(256 * max_intermediate_percent / 100.0f);
  int min_intermediates = (int)(256 * min_intermediate_percent / 100.0f);
  int step = 0;
  int steps = (int)((max_period_days - min_period_days) / increment_days);
  float response[1000000];
  if (steps > 1000000) {
    printf("Maximum number of time steps exceeded\n");
    return 0;
  }
#pragma omp parallel for
  for (step = 0; step < steps; step++) {
    float curve[256];
    float density[256];
    float orbital_period_days = min_period_days + (step * increment_days);
    if (light_curve(timestamp, series, series_length, orbital_period_days,
                    curve, density, 256) != 0)
      continue;
    float variance = light_curve_variance(orbital_period_days, timestamp,
                                          series, series_length, curve, 256);
    float mean = 0;
    int hits = 0;
    for (int j = 0; j < 256; j++) {
      if (curve[j] > 0) {
        mean += curve[j];
        hits++;
      }
    }
    if (hits < 256) {
      response[step] = 0;
      continue;
    }
    mean /= hits;
    float mean_density = 0;
    hits = 0;
    for (int j = 0; j < 256; j++) {
      if (density[j] > 0) {
        mean_density += density[j];
        hits++;
      }
    }
    mean_density /= hits;
    float density_variance = 0;
    hits = 0;
    for (int j = 0; j < 256; j++) {
      if (density[j] > 0) {
        density_variance +=
            (density[j] - mean_density) * (density[j] - mean_density);
        hits++;
      }
    }
    density_variance = (float)(density_variance / hits);
    float minimum = 0;
    for (int j = 0; j < 256; j++) {
      float v = 0;
      hits = 0;
      for (int k = j - expected_width; k <= j + expected_width; k++) {
        int l = k;
        if (l < 0)
          l += 256;
        if (l >= 256)
          l -= 256;
        if (curve[l] > 0) {
          v += curve[l];
          hits++;
          if (k == j) {
            v += curve[l];
            hits++;
          }
        }
      }
      if (hits > 0) {
        v /= hits;
        if ((v < minimum) || (minimum == 0)) {
          minimum = v;
        }
      }
    }
    int start_index = -1;
    int end_index = -1;
    int dipped = 0;
    float dipped_density = 0;
    float threshold_dipped = minimum + ((mean - minimum) * dip_threshold);
    for (int j = 0; j < 256; j++) {
      if (curve[j] < threshold_dipped) {
        if (start_index == -1)
          start_index = j;
        end_index = j;
        dipped++;
        if (dipped > max_dipped) {
          break;
        }
        dipped_density += density[j];
      }
    }
    if ((start_index == -1) && (end_index == -1)) {
      response[step] = 0;
      continue;
    }
    if (end_index - start_index > 256 * 10 / 100) {
      response[step] = 0;
      continue;
    }
    if (dipped > max_dipped) {
      dipped = 0;
    }
    if (dipped == 0) {
      response[step] = 0;
      continue;
    }
    dipped_density /= dipped;
    if (dipped_density < min_dipped_density) {
      response[step] = 0;
      continue;
    }
    int peaked = 0;
    float threshold_peaked = mean + ((mean - minimum) * peak_threshold);
    for (int j = 0; j < 256; j++) {
      if (curve[j] > threshold_peaked) {
        peaked++;
        break;
      }
    }
    if (peaked > 0) {
      response[step] = 0;
      continue;
    }
    int nondipped = 0;
    float threshold_upper = mean - ((mean - minimum) * 0.2);
    for (int j = 0; j < 256; j++) {
      if ((curve[j] < threshold_upper) && (curve[j] > threshold_dipped)) {
        nondipped++;
        if (nondipped > max_intermediates) {
          break;
        }
      }
    }
    if ((nondipped < min_intermediates) || (nondipped > max_intermediates)) {
      response[step] = 0;
      continue;
    }
    variance = 0;
    hits = 0;
    for (int j = 0; j < 256; j++) {
      if (curve[j] > 0) {
        variance += (curve[j] - mean) * (curve[j] - mean);
        hits++;
      }
    }
    response[step] = 0;
    if (hits > 0) {
      response[step] = (mean - minimum) / (float)sqrt(variance / hits);
      response[step] =
          (mean - minimum) * dipped * 100 / (mean * (1 + nondipped));
      response[step] /= (density_variance * variance);
      float vacancy_density =
          dip_vacancy(start_index, end_index, timestamp, series, series_length,
                      orbital_period_days, curve, 256);
      if (vacancy_density > max_vacancy_density) {
        response[step] = 0;
      } else {
        response[step] /= (1.0f + (vacancy_density * 10));
      }
    }
  }
  for (int i = 0; i < steps; i++) {
    if (response[i] > max_response) {
      max_response = response[i];
      period_days = min_period_days + (i * increment_days);
    }
  }
  return period_days;
}


x86 assembly (clang, -O0):
<detect_orbital_period>:
push   %rbp
mov    %rsp,%rbp
sub    $0x3d11e0,%rsp
movss  0x20(%rbp),%xmm8
movss  0x18(%rbp),%xmm8
movss  0x10(%rbp),%xmm8
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movss  %xmm0,-0x20(%rbp)
movss  %xmm1,-0x24(%rbp)
movss  %xmm2,-0x28(%rbp)
movss  %xmm3,-0x2c(%rbp)
movss  %xmm4,-0x30(%rbp)
movss  %xmm5,-0x34(%rbp)
movss  %xmm6,-0x38(%rbp)
movss  %xmm7,-0x3c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  0xe4c(%rip),%xmm0
mulss  -0x3c(%rbp),%xmm0
movss  0xe3b(%rip),%xmm1
divss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x48(%rbp)
movss  0xe2c(%rip),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  0xe1b(%rip),%xmm1
divss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x4c(%rbp)
movss  0xe0c(%rip),%xmm0
mulss  -0x38(%rbp),%xmm0
movss  0xdfb(%rip),%xmm1
divss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x50(%rbp)
movss  0xdec(%rip),%xmm0
mulss  -0x34(%rbp),%xmm0
movss  0xddb(%rip),%xmm1
divss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x54(%rbp)
movl   $0x0,-0x58(%rbp)
movss  -0x24(%rbp),%xmm0
subss  -0x20(%rbp),%xmm0
divss  -0x28(%rbp),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x5c(%rbp)
cmpl   $0xf4240,-0x5c(%rbp)
jle    1275 <detect_orbital_period+0x125>
lea    0xdb7(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1e1b <detect_orbital_period+0xccb>
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jge    1d91 <detect_orbital_period+0xc41>
movss  -0x20(%rbp),%xmm1
cvtsi2ssl -0x58(%rbp),%xmm0
movss  -0x28(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3d1164(%rbp)
ff
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x1c(%rbp),%edx
movss  -0x3d1164(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
lea    -0x3d0d60(%rbp),%rcx
lea    -0x3d1160(%rbp),%r8
mov    $0x100,%r9d
mov    $0x1,%al
call   1060 <light_curve@plt>
cmp    $0x0,%eax
je     12e7 <detect_orbital_period+0x197>
jmp    1d83 <detect_orbital_period+0xc33>
movss  -0x3d1164(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x1c(%rbp),%edx
lea    -0x3d0d60(%rbp),%rcx
mov    $0x100,%r8d
mov    $0x1,%al
call   1070 <light_curve_variance@plt>
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x3d1168(%rbp)
ff
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d116c(%rbp)
ff
movl   $0x0,-0x3d1170(%rbp)
movl   $0x0,-0x3d1174(%rbp)
cmpl   $0x100,-0x3d1174(%rbp)
jge    13b1 <detect_orbital_period+0x261>
movslq -0x3d1174(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    1398 <detect_orbital_period+0x248>
movslq -0x3d1174(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
addss  -0x3d116c(%rbp),%xmm0
ff
movss  %xmm0,-0x3d116c(%rbp)
ff
mov    -0x3d1170(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1170(%rbp)
jmp    139d <detect_orbital_period+0x24d>
mov    -0x3d1174(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1174(%rbp)
jmp    133d <detect_orbital_period+0x1ed>
cmpl   $0x100,-0x3d1170(%rbp)
jge    13d6 <detect_orbital_period+0x286>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d83 <detect_orbital_period+0xc33>
cvtsi2ssl -0x3d1170(%rbp),%xmm1
ff
movss  -0x3d116c(%rbp),%xmm0
ff
divss  %xmm1,%xmm0
movss  %xmm0,-0x3d116c(%rbp)
ff
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d1178(%rbp)
ff
movl   $0x0,-0x3d1170(%rbp)
movl   $0x0,-0x3d117c(%rbp)
cmpl   $0x100,-0x3d117c(%rbp)
jge    1485 <detect_orbital_period+0x335>
movslq -0x3d117c(%rbp),%rax
movss  -0x3d1160(%rbp,%rax,4),%xmm0
c2 ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    146c <detect_orbital_period+0x31c>
movslq -0x3d117c(%rbp),%rax
movss  -0x3d1160(%rbp,%rax,4),%xmm0
c2 ff
addss  -0x3d1178(%rbp),%xmm0
ff
movss  %xmm0,-0x3d1178(%rbp)
ff
mov    -0x3d1170(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1170(%rbp)
jmp    1471 <detect_orbital_period+0x321>
mov    -0x3d117c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d117c(%rbp)
jmp    1411 <detect_orbital_period+0x2c1>
cvtsi2ssl -0x3d1170(%rbp),%xmm1
ff
movss  -0x3d1178(%rbp),%xmm0
ff
divss  %xmm1,%xmm0
movss  %xmm0,-0x3d1178(%rbp)
ff
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d1180(%rbp)
ff
movl   $0x0,-0x3d1170(%rbp)
movl   $0x0,-0x3d1184(%rbp)
cmpl   $0x100,-0x3d1184(%rbp)
jge    155c <detect_orbital_period+0x40c>
movslq -0x3d1184(%rbp),%rax
movss  -0x3d1160(%rbp,%rax,4),%xmm0
c2 ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    1543 <detect_orbital_period+0x3f3>
movslq -0x3d1184(%rbp),%rax
movss  -0x3d1160(%rbp,%rax,4),%xmm0
c2 ff
subss  -0x3d1178(%rbp),%xmm0
ff
movslq -0x3d1184(%rbp),%rax
movss  -0x3d1160(%rbp,%rax,4),%xmm2
c2 ff
subss  -0x3d1178(%rbp),%xmm2
ff
movss  -0x3d1180(%rbp),%xmm1
ff
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3d1180(%rbp)
ff
mov    -0x3d1170(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1170(%rbp)
jmp    1548 <detect_orbital_period+0x3f8>
mov    -0x3d1184(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1184(%rbp)
jmp    14c0 <detect_orbital_period+0x370>
movss  -0x3d1180(%rbp),%xmm0
ff
cvtsi2ssl -0x3d1170(%rbp),%xmm1
ff
divss  %xmm1,%xmm0
movss  %xmm0,-0x3d1180(%rbp)
ff
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d1188(%rbp)
ff
movl   $0x0,-0x3d118c(%rbp)
cmpl   $0x100,-0x3d118c(%rbp)
jge    1757 <detect_orbital_period+0x607>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d1190(%rbp)
ff
movl   $0x0,-0x3d1170(%rbp)
mov    -0x3d118c(%rbp),%eax
sub    -0x48(%rbp),%eax
mov    %eax,-0x3d1194(%rbp)
mov    -0x3d1194(%rbp),%eax
mov    -0x3d118c(%rbp),%ecx
add    -0x48(%rbp),%ecx
cmp    %ecx,%eax
jg     16cd <detect_orbital_period+0x57d>
mov    -0x3d1194(%rbp),%eax
mov    %eax,-0x3d1198(%rbp)
cmpl   $0x0,-0x3d1198(%rbp)
jge    1602 <detect_orbital_period+0x4b2>
mov    -0x3d1198(%rbp),%eax
add    $0x100,%eax
mov    %eax,-0x3d1198(%rbp)
cmpl   $0x100,-0x3d1198(%rbp)
jl     1623 <detect_orbital_period+0x4d3>
mov    -0x3d1198(%rbp),%eax
sub    $0x100,%eax
mov    %eax,-0x3d1198(%rbp)
movslq -0x3d1198(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    16b4 <detect_orbital_period+0x564>
movslq -0x3d1198(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
addss  -0x3d1190(%rbp),%xmm0
ff
movss  %xmm0,-0x3d1190(%rbp)
ff
mov    -0x3d1170(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1170(%rbp)
mov    -0x3d1194(%rbp),%eax
cmp    -0x3d118c(%rbp),%eax
jne    16af <detect_orbital_period+0x55f>
movslq -0x3d1198(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
addss  -0x3d1190(%rbp),%xmm0
ff
movss  %xmm0,-0x3d1190(%rbp)
ff
mov    -0x3d1170(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1170(%rbp)
jmp    16b4 <detect_orbital_period+0x564>
jmp    16b9 <detect_orbital_period+0x569>
mov    -0x3d1194(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1194(%rbp)
jmp    15c1 <detect_orbital_period+0x471>
cmpl   $0x0,-0x3d1170(%rbp)
jle    173e <detect_orbital_period+0x5ee>
cvtsi2ssl -0x3d1170(%rbp),%xmm1
ff
movss  -0x3d1190(%rbp),%xmm0
ff
divss  %xmm1,%xmm0
movss  %xmm0,-0x3d1190(%rbp)
ff
movss  -0x3d1190(%rbp),%xmm1
ff
movss  -0x3d1188(%rbp),%xmm0
ff
ucomiss %xmm1,%xmm0
ja     1729 <detect_orbital_period+0x5d9>
movss  -0x3d1188(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1739 <detect_orbital_period+0x5e9>
jp     1739 <detect_orbital_period+0x5e9>
movss  -0x3d1190(%rbp),%xmm0
ff
movss  %xmm0,-0x3d1188(%rbp)
ff
jmp    173e <detect_orbital_period+0x5ee>
jmp    1743 <detect_orbital_period+0x5f3>
mov    -0x3d118c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d118c(%rbp)
jmp    158d <detect_orbital_period+0x43d>
movl   $0xffffffff,-0x3d119c(%rbp)
ff ff ff
movl   $0xffffffff,-0x3d11a0(%rbp)
ff ff ff
movl   $0x0,-0x3d11a4(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d11a8(%rbp)
ff
movss  -0x3d1188(%rbp),%xmm1
ff
movss  -0x3d116c(%rbp),%xmm0
ff
subss  -0x3d1188(%rbp),%xmm0
ff
movss  0x20(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3d11ac(%rbp)
ff
movl   $0x0,-0x3d11b0(%rbp)
cmpl   $0x100,-0x3d11b0(%rbp)
jge    1869 <detect_orbital_period+0x719>
movslq -0x3d11b0(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm1
c2 ff
movss  -0x3d11ac(%rbp),%xmm0
ff
ucomiss %xmm1,%xmm0
jbe    1850 <detect_orbital_period+0x700>
cmpl   $0xffffffff,-0x3d119c(%rbp)
jne    1801 <detect_orbital_period+0x6b1>
mov    -0x3d11b0(%rbp),%eax
mov    %eax,-0x3d119c(%rbp)
mov    -0x3d11b0(%rbp),%eax
mov    %eax,-0x3d11a0(%rbp)
mov    -0x3d11a4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11a4(%rbp)
mov    -0x3d11a4(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jle    1830 <detect_orbital_period+0x6e0>
jmp    1869 <detect_orbital_period+0x719>
movslq -0x3d11b0(%rbp),%rax
movss  -0x3d1160(%rbp,%rax,4),%xmm0
c2 ff
addss  -0x3d11a8(%rbp),%xmm0
ff
movss  %xmm0,-0x3d11a8(%rbp)
ff
jmp    1855 <detect_orbital_period+0x705>
mov    -0x3d11b0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11b0(%rbp)
jmp    17b7 <detect_orbital_period+0x667>
cmpl   $0xffffffff,-0x3d119c(%rbp)
jne    1898 <detect_orbital_period+0x748>
cmpl   $0xffffffff,-0x3d11a0(%rbp)
jne    1898 <detect_orbital_period+0x748>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d83 <detect_orbital_period+0xc33>
mov    -0x3d11a0(%rbp),%eax
sub    -0x3d119c(%rbp),%eax
cmp    $0x19,%eax
jle    18c2 <detect_orbital_period+0x772>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d83 <detect_orbital_period+0xc33>
mov    -0x3d11a4(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jle    18db <detect_orbital_period+0x78b>
movl   $0x0,-0x3d11a4(%rbp)
cmpl   $0x0,-0x3d11a4(%rbp)
jne    18fd <detect_orbital_period+0x7ad>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d83 <detect_orbital_period+0xc33>
cvtsi2ssl -0x3d11a4(%rbp),%xmm1
ff
movss  -0x3d11a8(%rbp),%xmm0
ff
divss  %xmm1,%xmm0
movss  %xmm0,-0x3d11a8(%rbp)
ff
movss  -0x3d11a8(%rbp),%xmm1
ff
movss  -0x2c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1944 <detect_orbital_period+0x7f4>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d83 <detect_orbital_period+0xc33>
movl   $0x0,-0x3d11b4(%rbp)
movss  -0x3d116c(%rbp),%xmm1
ff
movss  -0x3d116c(%rbp),%xmm0
ff
subss  -0x3d1188(%rbp),%xmm0
ff
movss  0x10(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3d11b8(%rbp)
ff
movl   $0x0,-0x3d11bc(%rbp)
cmpl   $0x100,-0x3d11bc(%rbp)
jge    19df <detect_orbital_period+0x88f>
movslq -0x3d11bc(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
ucomiss -0x3d11b8(%rbp),%xmm0
jbe    19c6 <detect_orbital_period+0x876>
mov    -0x3d11b4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11b4(%rbp)
jmp    19df <detect_orbital_period+0x88f>
jmp    19cb <detect_orbital_period+0x87b>
mov    -0x3d11bc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11bc(%rbp)
jmp    1985 <detect_orbital_period+0x835>
cmpl   $0x0,-0x3d11b4(%rbp)
jle    1a01 <detect_orbital_period+0x8b1>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d83 <detect_orbital_period+0xc33>
movl   $0x0,-0x3d11c0(%rbp)
movss  -0x3d116c(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movss  -0x3d116c(%rbp),%xmm0
ff
subss  -0x3d1188(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0x5c6(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3d11c4(%rbp)
ff
movl   $0x0,-0x3d11c8(%rbp)
cmpl   $0x100,-0x3d11c8(%rbp)
jge    1af7 <detect_orbital_period+0x9a7>
movslq -0x3d11c8(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm1
c2 ff
movss  -0x3d11c4(%rbp),%xmm0
ff
ucomiss %xmm1,%xmm0
jbe    1ade <detect_orbital_period+0x98e>
movslq -0x3d11c8(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
ucomiss -0x3d11ac(%rbp),%xmm0
jbe    1ade <detect_orbital_period+0x98e>
mov    -0x3d11c0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11c0(%rbp)
mov    -0x3d11c0(%rbp),%eax
cmp    -0x50(%rbp),%eax
jle    1ad9 <detect_orbital_period+0x989>
jmp    1af7 <detect_orbital_period+0x9a7>
jmp    1ade <detect_orbital_period+0x98e>
jmp    1ae3 <detect_orbital_period+0x993>
mov    -0x3d11c8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11c8(%rbp)
jmp    1a68 <detect_orbital_period+0x918>
mov    -0x3d11c0(%rbp),%eax
cmp    -0x54(%rbp),%eax
jl     1b15 <detect_orbital_period+0x9c5>
mov    -0x3d11c0(%rbp),%eax
cmp    -0x50(%rbp),%eax
jle    1b2a <detect_orbital_period+0x9da>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d83 <detect_orbital_period+0xc33>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d1168(%rbp)
ff
movl   $0x0,-0x3d1170(%rbp)
movl   $0x0,-0x3d11cc(%rbp)
cmpl   $0x100,-0x3d11cc(%rbp)
jge    1be5 <detect_orbital_period+0xa95>
movslq -0x3d11cc(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    1bcc <detect_orbital_period+0xa7c>
movslq -0x3d11cc(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm0
c2 ff
subss  -0x3d116c(%rbp),%xmm0
ff
movslq -0x3d11cc(%rbp),%rax
movss  -0x3d0d60(%rbp,%rax,4),%xmm2
c2 ff
subss  -0x3d116c(%rbp),%xmm2
ff
movss  -0x3d1168(%rbp),%xmm1
ff
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3d1168(%rbp)
ff
mov    -0x3d1170(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d1170(%rbp)
jmp    1bd1 <detect_orbital_period+0xa81>
mov    -0x3d11cc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11cc(%rbp)
jmp    1b49 <detect_orbital_period+0x9f9>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
cmpl   $0x0,-0x3d1170(%rbp)
jle    1d7e <detect_orbital_period+0xc2e>
movss  -0x3d116c(%rbp),%xmm0
ff
subss  -0x3d1188(%rbp),%xmm0
ff
movss  %xmm0,-0x3d11d8(%rbp)
ff
movss  -0x3d1168(%rbp),%xmm0
ff
cvtsi2ssl -0x3d1170(%rbp),%xmm1
ff
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x3d11d8(%rbp),%xmm0
ff
cvtsd2ss %xmm1,%xmm1
divss  %xmm1,%xmm0
movslq -0x58(%rbp),%rax
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
movss  -0x3d116c(%rbp),%xmm0
ff
subss  -0x3d1188(%rbp),%xmm0
ff
cvtsi2ssl -0x3d11a4(%rbp),%xmm1
ff
mulss  %xmm1,%xmm0
movss  0x385(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  -0x3d116c(%rbp),%xmm1
ff
mov    -0x3d11c0(%rbp),%eax
add    $0x1,%eax
cvtsi2ss %eax,%xmm2
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
movslq -0x58(%rbp),%rax
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
movss  -0x3d1180(%rbp),%xmm1
ff
mulss  -0x3d1168(%rbp),%xmm1
ff
movslq -0x58(%rbp),%rax
movss  -0x3d0960(%rbp,%rax,4),%xmm0
c2 ff
divss  %xmm1,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
mov    -0x3d119c(%rbp),%edi
mov    -0x3d11a0(%rbp),%esi
mov    -0x10(%rbp),%rdx
mov    -0x18(%rbp),%rcx
mov    -0x1c(%rbp),%r8d
movss  -0x3d1164(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
lea    -0x3d0d60(%rbp),%r9
movl   $0x100,(%rsp)
mov    $0x1,%al
call   1040 <dip_vacancy@plt>
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x3d11d0(%rbp)
ff
movss  -0x3d11d0(%rbp),%xmm0
ff
ucomiss 0x18(%rbp),%xmm0
jbe    1d3f <detect_orbital_period+0xbef>
movslq -0x58(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d79 <detect_orbital_period+0xc29>
movss  -0x3d11d0(%rbp),%xmm1
ff
movss  0x2b9(%rip),%xmm0
mulss  %xmm0,%xmm1
movss  0x2b1(%rip),%xmm0
addss  %xmm0,%xmm1
movslq -0x58(%rbp),%rax
movss  -0x3d0960(%rbp,%rax,4),%xmm0
c2 ff
divss  %xmm1,%xmm0
movss  %xmm0,-0x3d0960(%rbp,%rax,4)
c2 ff
jmp    1d7e <detect_orbital_period+0xc2e>
jmp    1d83 <detect_orbital_period+0xc33>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    127c <detect_orbital_period+0x12c>
movl   $0x0,-0x3d11d4(%rbp)
mov    -0x3d11d4(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jge    1e11 <detect_orbital_period+0xcc1>
movslq -0x3d11d4(%rbp),%rax
movss  -0x3d0960(%rbp,%rax,4),%xmm0
c2 ff
ucomiss -0x44(%rbp),%xmm0
jbe    1df8 <detect_orbital_period+0xca8>
movslq -0x3d11d4(%rbp),%rax
movss  -0x3d0960(%rbp,%rax,4),%xmm0
c2 ff
movss  %xmm0,-0x44(%rbp)
movss  -0x20(%rbp),%xmm1
cvtsi2ssl -0x3d11d4(%rbp),%xmm0
ff
movss  -0x28(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
jmp    1dfd <detect_orbital_period+0xcad>
mov    -0x3d11d4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3d11d4(%rbp)
jmp    1d9b <detect_orbital_period+0xc4b>
movss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x3d11e0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<detect_orbital_period>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x3d1198,%rsp
movss  %xmm3,0x44(%rsp)
mov    %edx,0x14(%rsp)
mov    %rsi,0x28(%rsp)
mov    %rdi,0x20(%rsp)
subss  %xmm0,%xmm1
divss  %xmm2,%xmm1
cvttss2si %xmm1,%eax
cmp    $0xf4240,%eax
jle    119c <detect_orbital_period+0x4c>
lea    0xe91(%rip),%rdi
call   1030 <puts@plt>
xorps  %xmm0,%xmm0
jmp    191b <detect_orbital_period+0x7cb>
test   %eax,%eax
jle    185a <detect_orbital_period+0x70a>
movaps %xmm0,%xmm8
movss  0xe50(%rip),%xmm0
mulss  %xmm0,%xmm7
movss  0xe48(%rip),%xmm1
divss  %xmm1,%xmm7
cvttss2si %xmm7,%r15d
mulss  %xmm0,%xmm4
divss  %xmm1,%xmm4
cvttss2si %xmm4,%r14d
mulss  %xmm0,%xmm6
divss  %xmm1,%xmm6
cvttss2si %xmm6,%ecx
mov    %ecx,0x3c(%rsp)
mulss  %xmm0,%xmm5
divss  %xmm1,%xmm5
cvttss2si %xmm5,%ecx
mov    %ecx,0x38(%rsp)
mov    %eax,0x40(%rsp)
mov    %eax,%eax
mov    %rax,0x18(%rsp)
lea    (%r15,%r15,1),%r13d
not    %r13d
xor    %ebp,%ebp
movss  %xmm8,0x4c(%rsp)
movss  %xmm2,0x48(%rsp)
jmp    1247 <detect_orbital_period+0xf7>
cs nopw 0x0(%rax,%rax,1)
nop
movl   $0x0,0x890(%rsp,%rbp,4)
add    $0x1,%rbp
cmp    0x18(%rsp),%rbp
movss  0x4c(%rsp),%xmm8
movss  0x48(%rsp),%xmm2
je     1840 <detect_orbital_period+0x6f0>
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm8,%xmm0
cvtss2sd %xmm0,%xmm0
mov    0x20(%rsp),%rdi
mov    0x28(%rsp),%rsi
mov    0x14(%rsp),%edx
movsd  %xmm0,0x50(%rsp)
lea    0x90(%rsp),%rcx
lea    0x490(%rsp),%r8
mov    $0x100,%r9d
mov    $0x1,%al
call   1060 <light_curve@plt>
test   %eax,%eax
jne    122b <detect_orbital_period+0xdb>
movsd  0x50(%rsp),%xmm0
mov    0x20(%rsp),%rdi
mov    0x28(%rsp),%rsi
mov    0x14(%rsp),%edx
lea    0x90(%rsp),%rcx
mov    $0x100,%r8d
mov    $0x1,%al
call   1070 <light_curve_variance@plt>
xorps  %xmm8,%xmm8
xor    %ecx,%ecx
xor    %eax,%eax
xorps  %xmm7,%xmm7
movss  0xd3c(%rip),%xmm4
nopl   0x0(%rax)
movss  0x90(%rsp,%rcx,4),%xmm0
xor    %edx,%edx
ucomiss %xmm7,%xmm0
movss  0x94(%rsp,%rcx,4),%xmm1
seta   %dl
maxss  %xmm4,%xmm0
add    %eax,%edx
xor    %eax,%eax
ucomiss %xmm7,%xmm1
addss  %xmm8,%xmm0
seta   %al
maxss  %xmm4,%xmm1
movaps %xmm1,%xmm8
addss  %xmm0,%xmm8
add    %edx,%eax
add    $0x2,%rcx
cmp    $0x100,%rcx
jne    12d0 <detect_orbital_period+0x180>
cmp    $0x100,%eax
jb     1220 <detect_orbital_period+0xd0>
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
xor    %eax,%eax
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  0x490(%rsp,%rax,4),%xmm1
xor    %edx,%edx
ucomiss %xmm7,%xmm1
movss  0x494(%rsp,%rax,4),%xmm3
seta   %dl
add    %ecx,%edx
maxss  %xmm4,%xmm1
xor    %ecx,%ecx
ucomiss %xmm7,%xmm3
addss  %xmm2,%xmm1
seta   %cl
add    %edx,%ecx
maxss  %xmm4,%xmm3
movaps %xmm3,%xmm2
addss  %xmm1,%xmm2
add    $0x2,%rax
cmp    $0x100,%rax
jne    1340 <detect_orbital_period+0x1f0>
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
divss  %xmm1,%xmm2
xor    %eax,%eax
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
jmp    13b0 <detect_orbital_period+0x260>
nopl   0x0(%rax)
mov    %bl,%dl
add    %edx,%ecx
add    $0x2,%rax
cmp    $0x100,%rax
je     13f4 <detect_orbital_period+0x2a4>
movss  0x490(%rsp,%rax,4),%xmm3
xor    %edx,%edx
ucomiss %xmm7,%xmm3
seta   %bl
jbe    13cf <detect_orbital_period+0x27f>
subss  %xmm2,%xmm3
mulss  %xmm3,%xmm3
addss  %xmm3,%xmm1
movss  0x494(%rsp,%rax,4),%xmm3
mov    %bl,%dl
add    %edx,%ecx
xor    %edx,%edx
ucomiss %xmm7,%xmm3
seta   %bl
jbe    13a0 <detect_orbital_period+0x250>
subss  %xmm2,%xmm3
mulss  %xmm3,%xmm3
addss  %xmm3,%xmm1
jmp    13a0 <detect_orbital_period+0x250>
xorps  %xmm3,%xmm3
cvtsi2ss %ecx,%xmm3
xorps  %xmm2,%xmm2
xor    %eax,%eax
jmp    141e <detect_orbital_period+0x2ce>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%eax
cmp    $0x100,%eax
je     14c3 <detect_orbital_period+0x373>
mov    %eax,%ecx
sub    %r15d,%ecx
lea    (%rax,%r15,1),%edx
cmp    %edx,%ecx
jg     1410 <detect_orbital_period+0x2c0>
xorps  %xmm4,%xmm4
xor    %esi,%esi
xor    %edx,%edx
jmp    144e <detect_orbital_period+0x2fe>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%edx
add    $0x1,%ecx
add    $0xffffffff,%esi
cmp    %esi,%r13d
je     1490 <detect_orbital_period+0x340>
lea    0x100(%rcx),%edi
test   %ecx,%ecx
cmovns %ecx,%edi
lea    -0x100(%rdi),%ebx
cmp    $0x100,%edi
cmovl  %edi,%ebx
movslq %ebx,%rdi
movss  0x90(%rsp,%rdi,4),%xmm5
ucomiss %xmm7,%xmm5
jbe    1443 <detect_orbital_period+0x2f3>
mov    %r15d,%edi
addss  %xmm5,%xmm4
add    %esi,%edi
jne    1440 <detect_orbital_period+0x2f0>
addss  %xmm5,%xmm4
add    $0x2,%edx
jmp    1443 <detect_orbital_period+0x2f3>
nopl   (%rax)
test   %edx,%edx
jle    1410 <detect_orbital_period+0x2c0>
xorps  %xmm5,%xmm5
cvtsi2ss %edx,%xmm5
divss  %xmm5,%xmm4
movaps %xmm2,%xmm5
cmpeqss %xmm7,%xmm5
movaps %xmm5,%xmm6
andps  %xmm4,%xmm6
minss  %xmm2,%xmm4
andnps %xmm4,%xmm5
orps   %xmm6,%xmm5
movaps %xmm5,%xmm2
jmp    1410 <detect_orbital_period+0x2c0>
unpcklps %xmm1,%xmm8
unpcklps %xmm3,%xmm0
divps  %xmm0,%xmm8
movaps %xmm8,%xmm3
subss  %xmm2,%xmm3
movaps %xmm3,%xmm0
mulss  0x3d11e0(%rsp),%xmm0
addss  %xmm2,%xmm0
mov    $0xffffffff,%ebx
xor    %ecx,%ecx
xorps  %xmm1,%xmm1
xor    %eax,%eax
mov    $0xffffffff,%r12d
jmp    1512 <detect_orbital_period+0x3c2>
nopw   0x0(%rax,%rax,1)
mov    %esi,%r12d
mov    %edx,%eax
add    $0x2,%rcx
cmp    $0x100,%rcx
je     1580 <detect_orbital_period+0x430>
ucomiss 0x90(%rsp,%rcx,4),%xmm0
jbe    1550 <detect_orbital_period+0x400>
cmp    $0xffffffff,%r12d
mov    %ecx,%esi
cmovne %r12d,%esi
lea    0x1(%rax),%edx
cmp    %r14d,%eax
jge    1591 <detect_orbital_period+0x441>
addss  0x490(%rsp,%rcx,4),%xmm1
mov    %ecx,%ebx
ucomiss 0x94(%rsp,%rcx,4),%xmm0
jbe    1500 <detect_orbital_period+0x3b0>
jmp    155f <detect_orbital_period+0x40f>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r12d,%esi
mov    %eax,%edx
ucomiss 0x94(%rsp,%rcx,4),%xmm0
jbe    1500 <detect_orbital_period+0x3b0>
lea    0x1(%rcx),%rbx
cmp    $0xffffffff,%esi
mov    %ebx,%r12d
cmovne %esi,%r12d
lea    0x1(%rdx),%eax
cmp    %r14d,%edx
jge    1599 <detect_orbital_period+0x449>
addss  0x494(%rsp,%rcx,4),%xmm1
jmp    1505 <detect_orbital_period+0x3b5>
cmp    $0xffffffff,%r12d
jne    1599 <detect_orbital_period+0x449>
cmp    $0xffffffff,%ebx
je     1220 <detect_orbital_period+0xd0>
jmp    1599 <detect_orbital_period+0x449>
mov    %rcx,%rbx
mov    %esi,%r12d
mov    %edx,%eax
mov    %ebx,%ecx
sub    %r12d,%ecx
cmp    $0x1a,%ecx
jge    1220 <detect_orbital_period+0xd0>
cmp    %r14d,%eax
jg     1220 <detect_orbital_period+0xd0>
test   %eax,%eax
je     1220 <detect_orbital_period+0xd0>
xorps  %xmm4,%xmm4
cvtsi2ss %eax,%xmm4
divss  %xmm4,%xmm1
movss  0x44(%rsp),%xmm2
ucomiss %xmm1,%xmm2
ja     1220 <detect_orbital_period+0xd0>
movaps %xmm3,%xmm1
mulss  0x3d11d0(%rsp),%xmm1
addss  %xmm8,%xmm1
xor    %eax,%eax
mov    0x3c(%rsp),%edx
nopl   0x0(%rax)
movss  0x90(%rsp,%rax,4),%xmm2
ucomiss %xmm1,%xmm2
ja     1220 <detect_orbital_period+0xd0>
movss  0x94(%rsp,%rax,4),%xmm2
ucomiss %xmm1,%xmm2
ja     1220 <detect_orbital_period+0xd0>
movss  0x98(%rsp,%rax,4),%xmm2
ucomiss %xmm1,%xmm2
ja     1220 <detect_orbital_period+0xd0>
movss  0x9c(%rsp,%rax,4),%xmm2
ucomiss %xmm1,%xmm2
ja     1220 <detect_orbital_period+0xd0>
add    $0x4,%rax
cmp    $0x100,%rax
jne    15f0 <detect_orbital_period+0x4a0>
xorps  %xmm1,%xmm1
cvtss2sd %xmm8,%xmm1
xorps  %xmm2,%xmm2
cvtss2sd %xmm3,%xmm2
mulsd  0x9bd(%rip),%xmm2
addsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
xor    %eax,%eax
xor    %edi,%edi
jmp    167c <detect_orbital_period+0x52c>
nopl   0x0(%rax)
add    $0x2,%rax
cmp    $0x100,%rax
je     16b4 <detect_orbital_period+0x564>
movss  0x90(%rsp,%rax,4),%xmm2
ucomiss %xmm2,%xmm1
jbe    1698 <detect_orbital_period+0x548>
ucomiss %xmm0,%xmm2
jbe    1698 <detect_orbital_period+0x548>
lea    0x1(%rdi),%ecx
cmp    %edx,%edi
mov    %ecx,%edi
jge    16b4 <detect_orbital_period+0x564>
movss  0x94(%rsp,%rax,4),%xmm2
ucomiss %xmm2,%xmm1
jbe    1670 <detect_orbital_period+0x520>
ucomiss %xmm0,%xmm2
jbe    1670 <detect_orbital_period+0x520>
lea    0x1(%rdi),%ecx
cmp    %edx,%edi
mov    %ecx,%edi
jl     1670 <detect_orbital_period+0x520>
cmp    0x38(%rsp),%edi
jl     1220 <detect_orbital_period+0xd0>
cmp    %edx,%edi
jg     1220 <detect_orbital_period+0xd0>
xorps  %xmm5,%xmm5
xor    %ecx,%ecx
xor    %eax,%eax
xorps  %xmm1,%xmm1
jmp    16f2 <detect_orbital_period+0x5a2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %sil,%dl
add    %edx,%eax
add    $0x2,%rcx
cmp    $0x100,%rcx
je     173b <detect_orbital_period+0x5eb>
movss  0x90(%rsp,%rcx,4),%xmm0
xor    %edx,%edx
ucomiss %xmm1,%xmm0
seta   %sil
jbe    1713 <detect_orbital_period+0x5c3>
subss  %xmm8,%xmm0
mulss  %xmm0,%xmm0
addss  %xmm0,%xmm5
movss  0x94(%rsp,%rcx,4),%xmm0
mov    %sil,%dl
add    %edx,%eax
xor    %edx,%edx
ucomiss %xmm1,%xmm0
seta   %sil
jbe    16e0 <detect_orbital_period+0x590>
subss  %xmm8,%xmm0
mulss  %xmm0,%xmm0
addss  %xmm0,%xmm5
jmp    16e0 <detect_orbital_period+0x590>
xorps  %xmm2,%xmm2
movl   $0x0,0x890(%rsp,%rbp,4)
test   %eax,%eax
je     122b <detect_orbital_period+0xdb>
mov    %rdi,0x58(%rsp)
movss  %xmm4,0x34(%rsp)
movaps %xmm3,0x70(%rsp)
movaps %xmm8,0x80(%rsp)
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
movaps %xmm5,0x60(%rsp)
movaps %xmm5,%xmm0
divss  %xmm1,%xmm0
ucomiss %xmm0,%xmm2
ja     182d <detect_orbital_period+0x6dd>
movl   $0x100,(%rsp)
mov    %r12d,%edi
mov    %ebx,%esi
mov    0x20(%rsp),%rdx
mov    0x28(%rsp),%rcx
mov    0x14(%rsp),%r8d
movsd  0x50(%rsp),%xmm0
lea    0x90(%rsp),%r9
mov    $0x1,%al
call   1050 <dip_vacancy@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
xorps  %xmm1,%xmm1
ucomiss 0x3d11d8(%rsp),%xmm0
ja     181f <detect_orbital_period+0x6cf>
movaps 0x70(%rsp),%xmm3
mulss  0x34(%rsp),%xmm3
mulss  0x827(%rip),%xmm3
mov    0x58(%rsp),%rax
add    $0x1,%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
unpcklps 0x60(%rsp),%xmm1
movaps 0x80(%rsp),%xmm2
mulps  %xmm1,%xmm2
divss  %xmm2,%xmm3
shufps $0x55,%xmm2,%xmm2
divss  %xmm2,%xmm3
mulss  0x7fc(%rip),%xmm0
addss  0x7f8(%rip),%xmm0
divss  %xmm0,%xmm3
movaps %xmm3,%xmm1
movss  %xmm1,0x890(%rsp,%rbp,4)
jmp    122b <detect_orbital_period+0xdb>
ucomiss %xmm2,%xmm0
jae    1786 <detect_orbital_period+0x636>
call   1040 <sqrtf@plt>
jmp    1786 <detect_orbital_period+0x636>
mov    0x40(%rsp),%eax
test   %eax,%eax
jle    185a <detect_orbital_period+0x70a>
cmp    $0x1,%eax
jne    1862 <detect_orbital_period+0x712>
xorps  %xmm0,%xmm0
xor    %eax,%eax
xorps  %xmm1,%xmm1
jmp    18ee <detect_orbital_period+0x79e>
xorps  %xmm0,%xmm0
jmp    191b <detect_orbital_period+0x7cb>
mov    0x18(%rsp),%rcx
and    $0xfffffffe,%ecx
xorps  %xmm0,%xmm0
xor    %eax,%eax
xorps  %xmm1,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorps  %xmm6,%xmm6
cvtsi2ss %eax,%xmm6
mulss  %xmm2,%xmm6
addss  %xmm8,%xmm6
movss  0x890(%rsp,%rax,4),%xmm3
movss  0x894(%rsp,%rax,4),%xmm4
movaps %xmm3,%xmm5
maxss  %xmm1,%xmm5
cmpltss %xmm3,%xmm1
andps  %xmm1,%xmm6
andnps %xmm0,%xmm1
movaps %xmm1,%xmm0
orps   %xmm6,%xmm0
lea    0x1(%rax),%edx
xorps  %xmm3,%xmm3
cvtsi2ss %edx,%xmm3
mulss  %xmm2,%xmm3
addss  %xmm8,%xmm3
movaps %xmm4,%xmm1
maxss  %xmm5,%xmm1
cmpltss %xmm4,%xmm5
andps  %xmm5,%xmm3
andnps %xmm0,%xmm5
orps   %xmm3,%xmm5
movaps %xmm5,%xmm0
add    $0x2,%rax
cmp    %rcx,%rax
jne    1880 <detect_orbital_period+0x730>
testb  $0x1,0x18(%rsp)
je     191b <detect_orbital_period+0x7cb>
xorps  %xmm3,%xmm3
cvtsi2ss %eax,%xmm3
mulss  %xmm2,%xmm3
cmpltss 0x890(%rsp,%rax,4),%xmm1
addss  %xmm8,%xmm3
andps  %xmm1,%xmm3
andnps %xmm0,%xmm1
orps   %xmm3,%xmm1
movaps %xmm1,%xmm0
add    $0x3d1198,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [39 x i8] c"Maximum number of time steps exceeded\0A\00"
define float @detect_orbital_period(float* %0, float* %1, i32 %2, float %3, float %4, float %5, float %6, float %7, float %8, float %9, float %10, float %11, float %12, float %13) {
%15 = alloca float
%16 = alloca float*
%17 = alloca float*
%18 = alloca i32
%19 = alloca float
%20 = alloca float
%21 = alloca float
%22 = alloca float
%23 = alloca float
%24 = alloca float
%25 = alloca float
%26 = alloca float
%27 = alloca float
%28 = alloca float
%29 = alloca float
%30 = alloca float
%31 = alloca float
%32 = alloca i32
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
%36 = alloca i32
%37 = alloca i32
%38 = alloca [1000000 x float]
%39 = alloca [256 x float]
%40 = alloca [256 x float]
%41 = alloca float
%42 = alloca float
%43 = alloca float
%44 = alloca i32
%45 = alloca i32
%46 = alloca float
%47 = alloca i32
%48 = alloca float
%49 = alloca i32
%50 = alloca float
%51 = alloca i32
%52 = alloca float
%53 = alloca i32
%54 = alloca i32
%55 = alloca i32
%56 = alloca i32
%57 = alloca i32
%58 = alloca float
%59 = alloca float
%60 = alloca i32
%61 = alloca i32
%62 = alloca float
%63 = alloca i32
%64 = alloca i32
%65 = alloca float
%66 = alloca i32
%67 = alloca i32
%68 = alloca float
%69 = alloca i32
store float* %0, float** %16
store float* %1, float** %17
store i32 %2, i32* %18
store float %3, float* %19
store float %4, float* %20
store float %5, float* %21
store float %6, float* %22
store float %7, float* %23
store float %8, float* %24
store float %9, float* %25
store float %10, float* %26
store float %11, float* %27
store float %12, float* %28
store float %13, float* %29
store float 0.000000e+00, float* %30
store float 0.000000e+00, float* %31
%70 = load float, float* %26
%71 = fmul float 2.560000e+02, %70
%72 = fdiv float %71, 1.000000e+02
%73 = fptosi float %72 to i32
store i32 %73, i32* %32
%74 = load float, float* %23
%75 = fmul float 2.560000e+02, %74
%76 = fdiv float %75, 1.000000e+02
%77 = fptosi float %76 to i32
store i32 %77, i32* %33
%78 = load float, float* %25
%79 = fmul float 2.560000e+02, %78
%80 = fdiv float %79, 1.000000e+02
%81 = fptosi float %80 to i32
store i32 %81, i32* %34
%82 = load float, float* %24
%83 = fmul float 2.560000e+02, %82
%84 = fdiv float %83, 1.000000e+02
%85 = fptosi float %84 to i32
store i32 %85, i32* %35
store i32 0, i32* %36
%86 = load float, float* %20
%87 = load float, float* %19
%88 = fsub float %86, %87
%89 = load float, float* %21
%90 = fdiv float %88, %89
%91 = fptosi float %90 to i32
store i32 %91, i32* %37
%92 = load i32, i32* %37
%93 = icmp sgt i32 %92, 1000000
br i1 %93, label %94, label %96
94: ; preds = %14
%95 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str, i64 0, i64 0))
store float 0.000000e+00, float* %15
br label %601
96: ; preds = %14
store i32 0, i32* %36
br label %97
97: ; preds = %570, %96
%98 = load i32, i32* %36
%99 = load i32, i32* %37
%100 = icmp slt i32 %98, %99
br i1 %100, label %101, label %573
101: ; preds = %97
%102 = load float, float* %19
%103 = load i32, i32* %36
%104 = sitofp i32 %103 to float
%105 = load float, float* %21
%106 = call float @llvm.fmuladd.f32(float %104, float %105, float %102)
store float %106, float* %41
%107 = load float*, float** %16
%108 = load float*, float** %17
%109 = load i32, i32* %18
%110 = load float, float* %41
%111 = fpext float %110 to double
%112 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 0
%113 = getelementptr inbounds [256 x float], [256 x float]* %40, i64 0, i64 0
%114 = call i32 (float*, float*, i32, double, float*, float*, i32, ...) bitcast (i32 (...)* @light_curve to i32 (float*, float*, i32, double, float*, float*, i32, ...)*)(float* %107, float* %108, i32 %109, double %111, float* %112, float* %113, i32 256)
%115 = icmp ne i32 %114, 0
br i1 %115, label %116, label %117
116: ; preds = %101
br label %570
117: ; preds = %101
%118 = load float, float* %41
%119 = fpext float %118 to double
%120 = load float*, float** %16
%121 = load float*, float** %17
%122 = load i32, i32* %18
%123 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 0
%124 = call i32 (double, float*, float*, i32, float*, i32, ...) bitcast (i32 (...)* @light_curve_variance to i32 (double, float*, float*, i32, float*, i32, ...)*)(double %119, float* %120, float* %121, i32 %122, float* %123, i32 256)
%125 = sitofp i32 %124 to float
store float %125, float* %42
store float 0.000000e+00, float* %43
store i32 0, i32* %44
store i32 0, i32* %45
br label %126
126: ; preds = %145, %117
%127 = load i32, i32* %45
%128 = icmp slt i32 %127, 256
br i1 %128, label %129, label %148
129: ; preds = %126
%130 = load i32, i32* %45
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %131
%133 = load float, float* %132
%134 = fcmp ogt float %133, 0.000000e+00
br i1 %134, label %135, label %144
135: ; preds = %129
%136 = load i32, i32* %45
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %137
%139 = load float, float* %138
%140 = load float, float* %43
%141 = fadd float %140, %139
store float %141, float* %43
%142 = load i32, i32* %44
%143 = add nsw i32 %142, 1
store i32 %143, i32* %44
br label %144
144: ; preds = %135, %129
br label %145
145: ; preds = %144
%146 = load i32, i32* %45
%147 = add nsw i32 %146, 1
store i32 %147, i32* %45
br label %126
148: ; preds = %126
%149 = load i32, i32* %44
%150 = icmp slt i32 %149, 256
br i1 %150, label %151, label %155
151: ; preds = %148
%152 = load i32, i32* %36
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %153
store float 0.000000e+00, float* %154
br label %570
155: ; preds = %148
%156 = load i32, i32* %44
%157 = sitofp i32 %156 to float
%158 = load float, float* %43
%159 = fdiv float %158, %157
store float %159, float* %43
store float 0.000000e+00, float* %46
store i32 0, i32* %44
store i32 0, i32* %47
br label %160
160: ; preds = %179, %155
%161 = load i32, i32* %47
%162 = icmp slt i32 %161, 256
br i1 %162, label %163, label %182
163: ; preds = %160
%164 = load i32, i32* %47
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds [256 x float], [256 x float]* %40, i64 0, i64 %165
%167 = load float, float* %166
%168 = fcmp ogt float %167, 0.000000e+00
br i1 %168, label %169, label %178
169: ; preds = %163
%170 = load i32, i32* %47
%171 = sext i32 %170 to i64
%172 = getelementptr inbounds [256 x float], [256 x float]* %40, i64 0, i64 %171
%173 = load float, float* %172
%174 = load float, float* %46
%175 = fadd float %174, %173
store float %175, float* %46
%176 = load i32, i32* %44
%177 = add nsw i32 %176, 1
store i32 %177, i32* %44
br label %178
178: ; preds = %169, %163
br label %179
179: ; preds = %178
%180 = load i32, i32* %47
%181 = add nsw i32 %180, 1
store i32 %181, i32* %47
br label %160
182: ; preds = %160
%183 = load i32, i32* %44
%184 = sitofp i32 %183 to float
%185 = load float, float* %46
%186 = fdiv float %185, %184
store float %186, float* %46
store float 0.000000e+00, float* %48
store i32 0, i32* %44
store i32 0, i32* %49
br label %187
187: ; preds = %214, %182
%188 = load i32, i32* %49
%189 = icmp slt i32 %188, 256
br i1 %189, label %190, label %217
190: ; preds = %187
%191 = load i32, i32* %49
%192 = sext i32 %191 to i64
%193 = getelementptr inbounds [256 x float], [256 x float]* %40, i64 0, i64 %192
%194 = load float, float* %193
%195 = fcmp ogt float %194, 0.000000e+00
br i1 %195, label %196, label %213
196: ; preds = %190
%197 = load i32, i32* %49
%198 = sext i32 %197 to i64
%199 = getelementptr inbounds [256 x float], [256 x float]* %40, i64 0, i64 %198
%200 = load float, float* %199
%201 = load float, float* %46
%202 = fsub float %200, %201
%203 = load i32, i32* %49
%204 = sext i32 %203 to i64
%205 = getelementptr inbounds [256 x float], [256 x float]* %40, i64 0, i64 %204
%206 = load float, float* %205
%207 = load float, float* %46
%208 = fsub float %206, %207
%209 = load float, float* %48
%210 = call float @llvm.fmuladd.f32(float %202, float %208, float %209)
store float %210, float* %48
%211 = load i32, i32* %44
%212 = add nsw i32 %211, 1
store i32 %212, i32* %44
br label %213
213: ; preds = %196, %190
br label %214
214: ; preds = %213
%215 = load i32, i32* %49
%216 = add nsw i32 %215, 1
store i32 %216, i32* %49
br label %187
217: ; preds = %187
%218 = load float, float* %48
%219 = load i32, i32* %44
%220 = sitofp i32 %219 to float
%221 = fdiv float %218, %220
store float %221, float* %48
store float 0.000000e+00, float* %50
store i32 0, i32* %51
br label %222
222: ; preds = %298, %217
%223 = load i32, i32* %51
%224 = icmp slt i32 %223, 256
br i1 %224, label %225, label %301
225: ; preds = %222
store float 0.000000e+00, float* %52
store i32 0, i32* %44
%226 = load i32, i32* %51
%227 = load i32, i32* %32
%228 = sub nsw i32 %226, %227
store i32 %228, i32* %53
br label %229
229: ; preds = %277, %225
%230 = load i32, i32* %53
%231 = load i32, i32* %51
%232 = load i32, i32* %32
%233 = add nsw i32 %231, %232
%234 = icmp sle i32 %230, %233
br i1 %234, label %235, label %280
235: ; preds = %229
%236 = load i32, i32* %53
store i32 %236, i32* %54
%237 = load i32, i32* %54
%238 = icmp slt i32 %237, 0
br i1 %238, label %239, label %242
239: ; preds = %235
%240 = load i32, i32* %54
%241 = add nsw i32 %240, 256
store i32 %241, i32* %54
br label %242
242: ; preds = %239, %235
%243 = load i32, i32* %54
%244 = icmp sge i32 %243, 256
br i1 %244, label %245, label %248
245: ; preds = %242
%246 = load i32, i32* %54
%247 = sub nsw i32 %246, 256
store i32 %247, i32* %54
br label %248
248: ; preds = %245, %242
%249 = load i32, i32* %54
%250 = sext i32 %249 to i64
%251 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %250
%252 = load float, float* %251
%253 = fcmp ogt float %252, 0.000000e+00
br i1 %253, label %254, label %276
254: ; preds = %248
%255 = load i32, i32* %54
%256 = sext i32 %255 to i64
%257 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %256
%258 = load float, float* %257
%259 = load float, float* %52
%260 = fadd float %259, %258
store float %260, float* %52
%261 = load i32, i32* %44
%262 = add nsw i32 %261, 1
store i32 %262, i32* %44
%263 = load i32, i32* %53
%264 = load i32, i32* %51
%265 = icmp eq i32 %263, %264
br i1 %265, label %266, label %275
266: ; preds = %254
%267 = load i32, i32* %54
%268 = sext i32 %267 to i64
%269 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %268
%270 = load float, float* %269
%271 = load float, float* %52
%272 = fadd float %271, %270
store float %272, float* %52
%273 = load i32, i32* %44
%274 = add nsw i32 %273, 1
store i32 %274, i32* %44
br label %275
275: ; preds = %266, %254
br label %276
276: ; preds = %275, %248
br label %277
277: ; preds = %276
%278 = load i32, i32* %53
%279 = add nsw i32 %278, 1
store i32 %279, i32* %53
br label %229
280: ; preds = %229
%281 = load i32, i32* %44
%282 = icmp sgt i32 %281, 0
br i1 %282, label %283, label %297
283: ; preds = %280
%284 = load i32, i32* %44
%285 = sitofp i32 %284 to float
%286 = load float, float* %52
%287 = fdiv float %286, %285
store float %287, float* %52
%288 = load float, float* %52
%289 = load float, float* %50
%290 = fcmp olt float %288, %289
br i1 %290, label %294, label %291
291: ; preds = %283
%292 = load float, float* %50
%293 = fcmp oeq float %292, 0.000000e+00
br i1 %293, label %294, label %296
294: ; preds = %291, %283
%295 = load float, float* %52
store float %295, float* %50
br label %296
296: ; preds = %294, %291
br label %297
297: ; preds = %296, %280
br label %298
298: ; preds = %297
%299 = load i32, i32* %51
%300 = add nsw i32 %299, 1
store i32 %300, i32* %51
br label %222
301: ; preds = %222
store i32 -1, i32* %55
store i32 -1, i32* %56
store i32 0, i32* %57
store float 0.000000e+00, float* %58
%302 = load float, float* %50
%303 = load float, float* %43
%304 = load float, float* %50
%305 = fsub float %303, %304
%306 = load float, float* %29
%307 = call float @llvm.fmuladd.f32(float %305, float %306, float %302)
store float %307, float* %59
store i32 0, i32* %60
br label %308
308: ; preds = %339, %301
%309 = load i32, i32* %60
%310 = icmp slt i32 %309, 256
br i1 %310, label %311, label %342
311: ; preds = %308
%312 = load i32, i32* %60
%313 = sext i32 %312 to i64
%314 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %313
%315 = load float, float* %314
%316 = load float, float* %59
%317 = fcmp olt float %315, %316
br i1 %317, label %318, label %338
318: ; preds = %311
%319 = load i32, i32* %55
%320 = icmp eq i32 %319, -1
br i1 %320, label %321, label %323
321: ; preds = %318
%322 = load i32, i32* %60
store i32 %322, i32* %55
br label %323
323: ; preds = %321, %318
%324 = load i32, i32* %60
store i32 %324, i32* %56
%325 = load i32, i32* %57
%326 = add nsw i32 %325, 1
store i32 %326, i32* %57
%327 = load i32, i32* %57
%328 = load i32, i32* %33
%329 = icmp sgt i32 %327, %328
br i1 %329, label %330, label %331
330: ; preds = %323
br label %342
331: ; preds = %323
%332 = load i32, i32* %60
%333 = sext i32 %332 to i64
%334 = getelementptr inbounds [256 x float], [256 x float]* %40, i64 0, i64 %333
%335 = load float, float* %334
%336 = load float, float* %58
%337 = fadd float %336, %335
store float %337, float* %58
br label %338
338: ; preds = %331, %311
br label %339
339: ; preds = %338
%340 = load i32, i32* %60
%341 = add nsw i32 %340, 1
store i32 %341, i32* %60
br label %308
342: ; preds = %330, %308
%343 = load i32, i32* %55
%344 = icmp eq i32 %343, -1
br i1 %344, label %345, label %352
345: ; preds = %342
%346 = load i32, i32* %56
%347 = icmp eq i32 %346, -1
br i1 %347, label %348, label %352
348: ; preds = %345
%349 = load i32, i32* %36
%350 = sext i32 %349 to i64
%351 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %350
store float 0.000000e+00, float* %351
br label %570
352: ; preds = %345, %342
%353 = load i32, i32* %56
%354 = load i32, i32* %55
%355 = sub nsw i32 %353, %354
%356 = icmp sgt i32 %355, 25
br i1 %356, label %357, label %361
357: ; preds = %352
%358 = load i32, i32* %36
%359 = sext i32 %358 to i64
%360 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %359
store float 0.000000e+00, float* %360
br label %570
361: ; preds = %352
%362 = load i32, i32* %57
%363 = load i32, i32* %33
%364 = icmp sgt i32 %362, %363
br i1 %364, label %365, label %366
365: ; preds = %361
store i32 0, i32* %57
br label %366
366: ; preds = %365, %361
%367 = load i32, i32* %57
%368 = icmp eq i32 %367, 0
br i1 %368, label %369, label %373
369: ; preds = %366
%370 = load i32, i32* %36
%371 = sext i32 %370 to i64
%372 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %371
store float 0.000000e+00, float* %372
br label %570
373: ; preds = %366
%374 = load i32, i32* %57
%375 = sitofp i32 %374 to float
%376 = load float, float* %58
%377 = fdiv float %376, %375
store float %377, float* %58
%378 = load float, float* %58
%379 = load float, float* %22
%380 = fcmp olt float %378, %379
br i1 %380, label %381, label %385
381: ; preds = %373
%382 = load i32, i32* %36
%383 = sext i32 %382 to i64
%384 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %383
store float 0.000000e+00, float* %384
br label %570
385: ; preds = %373
store i32 0, i32* %61
%386 = load float, float* %43
%387 = load float, float* %43
%388 = load float, float* %50
%389 = fsub float %387, %388
%390 = load float, float* %27
%391 = call float @llvm.fmuladd.f32(float %389, float %390, float %386)
store float %391, float* %62
store i32 0, i32* %63
br label %392
392: ; preds = %406, %385
%393 = load i32, i32* %63
%394 = icmp slt i32 %393, 256
br i1 %394, label %395, label %409
395: ; preds = %392
%396 = load i32, i32* %63
%397 = sext i32 %396 to i64
%398 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %397
%399 = load float, float* %398
%400 = load float, float* %62
%401 = fcmp ogt float %399, %400
br i1 %401, label %402, label %405
402: ; preds = %395
%403 = load i32, i32* %61
%404 = add nsw i32 %403, 1
store i32 %404, i32* %61
br label %409
405: ; preds = %395
br label %406
406: ; preds = %405
%407 = load i32, i32* %63
%408 = add nsw i32 %407, 1
store i32 %408, i32* %63
br label %392
409: ; preds = %402, %392
%410 = load i32, i32* %61
%411 = icmp sgt i32 %410, 0
br i1 %411, label %412, label %416
412: ; preds = %409
%413 = load i32, i32* %36
%414 = sext i32 %413 to i64
%415 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %414
store float 0.000000e+00, float* %415
br label %570
416: ; preds = %409
store i32 0, i32* %64
%417 = load float, float* %43
%418 = fpext float %417 to double
%419 = load float, float* %43
%420 = load float, float* %50
%421 = fsub float %419, %420
%422 = fpext float %421 to double
%423 = fneg double %422
%424 = call double @llvm.fmuladd.f64(double %423, double 2.000000e-01, double %418)
%425 = fptrunc double %424 to float
store float %425, float* %65
store i32 0, i32* %66
br label %426
426: ; preds = %452, %416
%427 = load i32, i32* %66
%428 = icmp slt i32 %427, 256
br i1 %428, label %429, label %455
429: ; preds = %426
%430 = load i32, i32* %66
%431 = sext i32 %430 to i64
%432 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %431
%433 = load float, float* %432
%434 = load float, float* %65
%435 = fcmp olt float %433, %434
br i1 %435, label %436, label %451
436: ; preds = %429
%437 = load i32, i32* %66
%438 = sext i32 %437 to i64
%439 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %438
%440 = load float, float* %439
%441 = load float, float* %59
%442 = fcmp ogt float %440, %441
br i1 %442, label %443, label %451
443: ; preds = %436
%444 = load i32, i32* %64
%445 = add nsw i32 %444, 1
store i32 %445, i32* %64
%446 = load i32, i32* %64
%447 = load i32, i32* %34
%448 = icmp sgt i32 %446, %447
br i1 %448, label %449, label %450
449: ; preds = %443
br label %455
450: ; preds = %443
br label %451
451: ; preds = %450, %436, %429
br label %452
452: ; preds = %451
%453 = load i32, i32* %66
%454 = add nsw i32 %453, 1
store i32 %454, i32* %66
br label %426
455: ; preds = %449, %426
%456 = load i32, i32* %64
%457 = load i32, i32* %35
%458 = icmp slt i32 %456, %457
br i1 %458, label %463, label %459
459: ; preds = %455
%460 = load i32, i32* %64
%461 = load i32, i32* %34
%462 = icmp sgt i32 %460, %461
br i1 %462, label %463, label %467
463: ; preds = %459, %455
%464 = load i32, i32* %36
%465 = sext i32 %464 to i64
%466 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %465
store float 0.000000e+00, float* %466
br label %570
467: ; preds = %459
store float 0.000000e+00, float* %42
store i32 0, i32* %44
store i32 0, i32* %67
br label %468
468: ; preds = %495, %467
%469 = load i32, i32* %67
%470 = icmp slt i32 %469, 256
br i1 %470, label %471, label %498
471: ; preds = %468
%472 = load i32, i32* %67
%473 = sext i32 %472 to i64
%474 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %473
%475 = load float, float* %474
%476 = fcmp ogt float %475, 0.000000e+00
br i1 %476, label %477, label %494
477: ; preds = %471
%478 = load i32, i32* %67
%479 = sext i32 %478 to i64
%480 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %479
%481 = load float, float* %480
%482 = load float, float* %43
%483 = fsub float %481, %482
%484 = load i32, i32* %67
%485 = sext i32 %484 to i64
%486 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 %485
%487 = load float, float* %486
%488 = load float, float* %43
%489 = fsub float %487, %488
%490 = load float, float* %42
%491 = call float @llvm.fmuladd.f32(float %483, float %489, float %490)
store float %491, float* %42
%492 = load i32, i32* %44
%493 = add nsw i32 %492, 1
store i32 %493, i32* %44
br label %494
494: ; preds = %477, %471
br label %495
495: ; preds = %494
%496 = load i32, i32* %67
%497 = add nsw i32 %496, 1
store i32 %497, i32* %67
br label %468
498: ; preds = %468
%499 = load i32, i32* %36
%500 = sext i32 %499 to i64
%501 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %500
store float 0.000000e+00, float* %501
%502 = load i32, i32* %44
%503 = icmp sgt i32 %502, 0
br i1 %503, label %504, label %569
504: ; preds = %498
%505 = load float, float* %43
%506 = load float, float* %50
%507 = fsub float %505, %506
%508 = load float, float* %42
%509 = load i32, i32* %44
%510 = sitofp i32 %509 to float
%511 = fdiv float %508, %510
%512 = fpext float %511 to double
%513 = call double @sqrt(double %512)
%514 = fptrunc double %513 to float
%515 = fdiv float %507, %514
%516 = load i32, i32* %36
%517 = sext i32 %516 to i64
%518 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %517
store float %515, float* %518
%519 = load float, float* %43
%520 = load float, float* %50
%521 = fsub float %519, %520
%522 = load i32, i32* %57
%523 = sitofp i32 %522 to float
%524 = fmul float %521, %523
%525 = fmul float %524, 1.000000e+02
%526 = load float, float* %43
%527 = load i32, i32* %64
%528 = add nsw i32 1, %527
%529 = sitofp i32 %528 to float
%530 = fmul float %526, %529
%531 = fdiv float %525, %530
%532 = load i32, i32* %36
%533 = sext i32 %532 to i64
%534 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %533
store float %531, float* %534
%535 = load float, float* %48
%536 = load float, float* %42
%537 = fmul float %535, %536
%538 = load i32, i32* %36
%539 = sext i32 %538 to i64
%540 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %539
%541 = load float, float* %540
%542 = fdiv float %541, %537
store float %542, float* %540
%543 = load i32, i32* %55
%544 = load i32, i32* %56
%545 = load float*, float** %16
%546 = load float*, float** %17
%547 = load i32, i32* %18
%548 = load float, float* %41
%549 = fpext float %548 to double
%550 = getelementptr inbounds [256 x float], [256 x float]* %39, i64 0, i64 0
%551 = call i32 (i32, i32, float*, float*, i32, double, float*, i32, ...) bitcast (i32 (...)* @dip_vacancy to i32 (i32, i32, float*, float*, i32, double, float*, i32, ...)*)(i32 %543, i32 %544, float* %545, float* %546, i32 %547, double %549, float* %550, i32 256)
%552 = sitofp i32 %551 to float
store float %552, float* %68
%553 = load float, float* %68
%554 = load float, float* %28
%555 = fcmp ogt float %553, %554
br i1 %555, label %556, label %560
556: ; preds = %504
%557 = load i32, i32* %36
%558 = sext i32 %557 to i64
%559 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %558
store float 0.000000e+00, float* %559
br label %568
560: ; preds = %504
%561 = load float, float* %68
%562 = call float @llvm.fmuladd.f32(float %561, float 1.000000e+01, float 1.000000e+00)
%563 = load i32, i32* %36
%564 = sext i32 %563 to i64
%565 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %564
%566 = load float, float* %565
%567 = fdiv float %566, %562
store float %567, float* %565
br label %568
568: ; preds = %560, %556
br label %569
569: ; preds = %568, %498
br label %570
570: ; preds = %569, %463, %412, %381, %369, %357, %348, %151, %116
%571 = load i32, i32* %36
%572 = add nsw i32 %571, 1
store i32 %572, i32* %36
br label %97
573: ; preds = %97
store i32 0, i32* %69
br label %574
574: ; preds = %596, %573
%575 = load i32, i32* %69
%576 = load i32, i32* %37
%577 = icmp slt i32 %575, %576
br i1 %577, label %578, label %599
578: ; preds = %574
%579 = load i32, i32* %69
%580 = sext i32 %579 to i64
%581 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %580
%582 = load float, float* %581
%583 = load float, float* %31
%584 = fcmp ogt float %582, %583
br i1 %584, label %585, label %595
585: ; preds = %578
%586 = load i32, i32* %69
%587 = sext i32 %586 to i64
%588 = getelementptr inbounds [1000000 x float], [1000000 x float]* %38, i64 0, i64 %587
%589 = load float, float* %588
store float %589, float* %31
%590 = load float, float* %19
%591 = load i32, i32* %69
%592 = sitofp i32 %591 to float
%593 = load float, float* %21
%594 = call float @llvm.fmuladd.f32(float %592, float %593, float %590)
store float %594, float* %30
br label %595
595: ; preds = %585, %578
br label %596
596: ; preds = %595
%597 = load i32, i32* %69
%598 = add nsw i32 %597, 1
store i32 %598, i32* %69
br label %574
599: ; preds = %574
%600 = load float, float* %30
store float %600, float* %15
br label %601
601: ; preds = %599, %94
%602 = load float, float* %15
ret float %602
}

LLVM IR (-O2):
@str = [38 x i8] c"Maximum number of time steps exceeded\00"
define float @detect_orbital_period(float* %0, float* %1, i32 %2, float %3, float %4, float %5, float %6, float %7, float %8, float %9, float %10, float %11, float %12, float %13) {
%15 = alloca [1000000 x float]
%16 = alloca [256 x float]
%17 = alloca [256 x float]
%18 = fmul float %10, 2.560000e+02
%19 = fdiv float %18, 1.000000e+02
%20 = fptosi float %19 to i32
%21 = fmul float %7, 2.560000e+02
%22 = fdiv float %21, 1.000000e+02
%23 = fptosi float %22 to i32
%24 = fmul float %9, 2.560000e+02
%25 = fdiv float %24, 1.000000e+02
%26 = fptosi float %25 to i32
%27 = fmul float %8, 2.560000e+02
%28 = fdiv float %27, 1.000000e+02
%29 = fptosi float %28 to i32
%30 = fsub float %4, %3
%31 = fdiv float %30, %5
%32 = fptosi float %31 to i32
%33 = bitcast [1000000 x float]* %15 to i8*
%34 = icmp sgt i32 %32, 1000000
br i1 %34, label %43, label %35
35: ; preds = %14
%36 = icmp sgt i32 %32, 0
br i1 %36, label %37, label %435
37: ; preds = %35
%38 = bitcast [256 x float]* %16 to i8*
%39 = bitcast [256 x float]* %17 to i8*
%40 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 0
%41 = getelementptr inbounds [256 x float], [256 x float]* %17, i64 0, i64 0
%42 = zext i32 %32 to i64
br label %51
43: ; preds = %14
%44 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([38 x i8], [38 x i8]* @str, i64 0, i64 0))
br label %435
45: ; preds = %393
br i1 %36, label %46, label %435
46: ; preds = %45
%47 = and i64 %42, 1
%48 = icmp eq i32 %32, 1
br i1 %48, label %421, label %49
49: ; preds = %46
%50 = and i64 %42, 4294967294
br label %396
51: ; preds = %37, %393
%52 = phi i64 [ 0, %37 ], [ %394, %393 ]
%53 = trunc i64 %52 to i32
%54 = sitofp i32 %53 to float
%55 = call float @llvm.fmuladd.f32(float %54, float %5, float %3)
%56 = fpext float %55 to double
%57 = call i32 (float*, float*, i32, double, float*, float*, i32, ...) bitcast (i32 (...)* @light_curve to i32 (float*, float*, i32, double, float*, float*, i32, ...)*)(float* %0, float* %1, i32 %2, double %56, float* %40, float* %41, i32 256)
%58 = icmp eq i32 %57, 0
br i1 %58, label %59, label %393
59: ; preds = %51
%60 = call i32 (double, float*, float*, i32, float*, i32, ...) bitcast (i32 (...)* @light_curve_variance to i32 (double, float*, float*, i32, float*, i32, ...)*)(double %56, float* %0, float* %1, i32 %2, float* %40, i32 256)
br label %63
61: ; preds = %63
%62 = icmp ult i32 %83, 256
br i1 %62, label %86, label %88
63: ; preds = %63, %59
%64 = phi i64 [ 0, %59 ], [ %84, %63 ]
%65 = phi i32 [ 0, %59 ], [ %83, %63 ]
%66 = phi float [ 0.000000e+00, %59 ], [ %81, %63 ]
%67 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %64
%68 = load float, float* %67
%69 = fcmp ogt float %68, 0.000000e+00
%70 = fcmp ogt float %68, -0.000000e+00
%71 = select i1 %70, float %68, float -0.000000e+00
%72 = fadd float %66, %71
%73 = zext i1 %69 to i32
%74 = add nuw nsw i32 %65, %73
%75 = or i64 %64, 1
%76 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %75
%77 = load float, float* %76
%78 = fcmp ogt float %77, 0.000000e+00
%79 = fcmp ogt float %77, -0.000000e+00
%80 = select i1 %79, float %77, float -0.000000e+00
%81 = fadd float %72, %80
%82 = zext i1 %78 to i32
%83 = add nuw nsw i32 %74, %82
%84 = add nuw nsw i64 %64, 2
%85 = icmp eq i64 %84, 256
br i1 %85, label %61, label %63
86: ; preds = %61
%87 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %87
br label %393
88: ; preds = %61
%89 = sitofp i32 %83 to float
br label %93
90: ; preds = %93
%91 = sitofp i32 %110 to float
%92 = fdiv float %113, %91
br label %118
93: ; preds = %93, %88
%94 = phi i64 [ 0, %88 ], [ %114, %93 ]
%95 = phi float [ 0.000000e+00, %88 ], [ %113, %93 ]
%96 = phi i32 [ 0, %88 ], [ %110, %93 ]
%97 = getelementptr inbounds [256 x float], [256 x float]* %17, i64 0, i64 %94
%98 = load float, float* %97
%99 = fcmp ogt float %98, 0.000000e+00
%100 = zext i1 %99 to i32
%101 = add nuw nsw i32 %96, %100
%102 = fcmp ogt float %98, -0.000000e+00
%103 = select i1 %102, float %98, float -0.000000e+00
%104 = fadd float %95, %103
%105 = or i64 %94, 1
%106 = getelementptr inbounds [256 x float], [256 x float]* %17, i64 0, i64 %105
%107 = load float, float* %106
%108 = fcmp ogt float %107, 0.000000e+00
%109 = zext i1 %108 to i32
%110 = add nuw nsw i32 %101, %109
%111 = fcmp ogt float %107, -0.000000e+00
%112 = select i1 %111, float %107, float -0.000000e+00
%113 = fadd float %104, %112
%114 = add nuw nsw i64 %94, 2
%115 = icmp eq i64 %114, 256
br i1 %115, label %90, label %93
116: ; preds = %118
%117 = sitofp i32 %137 to float
br label %150
118: ; preds = %118, %90
%119 = phi i64 [ 0, %90 ], [ %139, %118 ]
%120 = phi float [ 0.000000e+00, %90 ], [ %138, %118 ]
%121 = phi i32 [ 0, %90 ], [ %137, %118 ]
%122 = getelementptr inbounds [256 x float], [256 x float]* %17, i64 0, i64 %119
%123 = load float, float* %122
%124 = fcmp ogt float %123, 0.000000e+00
%125 = fsub float %123, %92
%126 = call float @llvm.fmuladd.f32(float %125, float %125, float %120)
%127 = zext i1 %124 to i32
%128 = add nuw nsw i32 %121, %127
%129 = select i1 %124, float %126, float %120
%130 = or i64 %119, 1
%131 = getelementptr inbounds [256 x float], [256 x float]* %17, i64 0, i64 %130
%132 = load float, float* %131
%133 = fcmp ogt float %132, 0.000000e+00
%134 = fsub float %132, %92
%135 = call float @llvm.fmuladd.f32(float %134, float %134, float %129)
%136 = zext i1 %133 to i32
%137 = add nuw nsw i32 %128, %136
%138 = select i1 %133, float %135, float %129
%139 = add nuw nsw i64 %119, 2
%140 = icmp eq i64 %139, 256
br i1 %140, label %116, label %118
141: ; preds = %193
%142 = insertelement <2 x float> poison, float %81, i64 0
%143 = insertelement <2 x float> %142, float %138, i64 1
%144 = insertelement <2 x float> poison, float %89, i64 0
%145 = insertelement <2 x float> %144, float %117, i64 1
%146 = fdiv <2 x float> %143, %145
%147 = extractelement <2 x float> %146, i64 0
%148 = fsub float %147, %194
%149 = call float @llvm.fmuladd.f32(float %148, float %13, float %194)
br label %197
150: ; preds = %116, %193
%151 = phi i32 [ %20, %116 ], [ %154, %193 ]
%152 = phi i32 [ 0, %116 ], [ %195, %193 ]
%153 = phi float [ 0.000000e+00, %116 ], [ %194, %193 ]
%154 = add i32 %151, 1
%155 = sub nsw i32 %152, %20
%156 = add nsw i32 %152, %20
%157 = icmp sgt i32 %155, %156
br i1 %157, label %193, label %160
158: ; preds = %181
%159 = icmp sgt i32 %182, 0
br i1 %159, label %186, label %193
160: ; preds = %150, %181
%161 = phi i32 [ %184, %181 ], [ %155, %150 ]
%162 = phi float [ %183, %181 ], [ 0.000000e+00, %150 ]
%163 = phi i32 [ %182, %181 ], [ 0, %150 ]
%164 = icmp slt i32 %161, 0
%165 = add nsw i32 %161, 256
%166 = select i1 %164, i32 %165, i32 %161
%167 = icmp sgt i32 %166, 255
%168 = add nsw i32 %166, -256
%169 = select i1 %167, i32 %168, i32 %166
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %170
%172 = load float, float* %171
%173 = fcmp ogt float %172, 0.000000e+00
br i1 %173, label %174, label %181
174: ; preds = %160
%175 = fadd float %162, %172
%176 = add nsw i32 %163, 1
%177 = icmp eq i32 %161, %152
br i1 %177, label %178, label %181
178: ; preds = %174
%179 = fadd float %172, %175
%180 = add nsw i32 %163, 2
br label %181
181: ; preds = %174, %178, %160
%182 = phi i32 [ %180, %178 ], [ %176, %174 ], [ %163, %160 ]
%183 = phi float [ %179, %178 ], [ %175, %174 ], [ %162, %160 ]
%184 = add i32 %161, 1
%185 = icmp eq i32 %161, %151
br i1 %185, label %158, label %160
186: ; preds = %158
%187 = sitofp i32 %182 to float
%188 = fdiv float %183, %187
%189 = fcmp olt float %188, %153
%190 = fcmp oeq float %153, 0.000000e+00
%191 = or i1 %190, %189
%192 = select i1 %191, float %188, float %153
br label %193
193: ; preds = %150, %186, %158
%194 = phi float [ %153, %158 ], [ %192, %186 ], [ %153, %150 ]
%195 = add nuw nsw i32 %152, 1
%196 = icmp eq i32 %195, 256
br i1 %196, label %141, label %150
197: ; preds = %235, %141
%198 = phi i64 [ 0, %141 ], [ %240, %235 ]
%199 = phi float [ 0.000000e+00, %141 ], [ %239, %235 ]
%200 = phi i32 [ 0, %141 ], [ %238, %235 ]
%201 = phi i32 [ -1, %141 ], [ %237, %235 ]
%202 = phi i32 [ -1, %141 ], [ %236, %235 ]
%203 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %198
%204 = load float, float* %203
%205 = fcmp olt float %204, %149
br i1 %205, label %206, label %216
206: ; preds = %197
%207 = icmp eq i32 %202, -1
%208 = trunc i64 %198 to i32
%209 = select i1 %207, i32 %208, i32 %202
%210 = add nsw i32 %200, 1
%211 = icmp slt i32 %200, %23
br i1 %211, label %212, label %248
212: ; preds = %206
%213 = getelementptr inbounds [256 x float], [256 x float]* %17, i64 0, i64 %198
%214 = load float, float* %213
%215 = fadd float %199, %214
br label %216
216: ; preds = %197, %212
%217 = phi i32 [ %209, %212 ], [ %202, %197 ]
%218 = phi i32 [ %208, %212 ], [ %201, %197 ]
%219 = phi i32 [ %210, %212 ], [ %200, %197 ]
%220 = phi float [ %215, %212 ], [ %199, %197 ]
%221 = or i64 %198, 1
%222 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %221
%223 = load float, float* %222
%224 = fcmp olt float %223, %149
br i1 %224, label %225, label %235
225: ; preds = %216
%226 = icmp eq i32 %217, -1
%227 = trunc i64 %221 to i32
%228 = select i1 %226, i32 %227, i32 %217
%229 = add nsw i32 %219, 1
%230 = icmp slt i32 %219, %23
br i1 %230, label %231, label %248
231: ; preds = %225
%232 = getelementptr inbounds [256 x float], [256 x float]* %17, i64 0, i64 %221
%233 = load float, float* %232
%234 = fadd float %220, %233
br label %235
235: ; preds = %231, %216
%236 = phi i32 [ %228, %231 ], [ %217, %216 ]
%237 = phi i32 [ %227, %231 ], [ %218, %216 ]
%238 = phi i32 [ %229, %231 ], [ %219, %216 ]
%239 = phi float [ %234, %231 ], [ %220, %216 ]
%240 = add nuw nsw i64 %198, 2
%241 = icmp eq i64 %240, 256
br i1 %241, label %242, label %197
242: ; preds = %235
%243 = icmp eq i32 %236, -1
%244 = icmp eq i32 %237, -1
%245 = select i1 %243, i1 %244, i1 false
br i1 %245, label %246, label %254
246: ; preds = %242
%247 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %247
br label %393
248: ; preds = %225, %206
%249 = phi i64 [ %198, %206 ], [ %221, %225 ]
%250 = phi float [ %199, %206 ], [ %220, %225 ]
%251 = phi i32 [ %209, %206 ], [ %228, %225 ]
%252 = phi i32 [ %210, %206 ], [ %229, %225 ]
%253 = trunc i64 %249 to i32
br label %254
254: ; preds = %248, %242
%255 = phi float [ %239, %242 ], [ %250, %248 ]
%256 = phi i32 [ %238, %242 ], [ %252, %248 ]
%257 = phi i32 [ %237, %242 ], [ %253, %248 ]
%258 = phi i32 [ %236, %242 ], [ %251, %248 ]
%259 = sub nsw i32 %257, %258
%260 = icmp sgt i32 %259, 25
br i1 %260, label %261, label %263
261: ; preds = %254
%262 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %262
br label %393
263: ; preds = %254
%264 = icmp sgt i32 %256, %23
%265 = icmp eq i32 %256, 0
%266 = or i1 %264, %265
br i1 %266, label %267, label %269
267: ; preds = %263
%268 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %268
br label %393
269: ; preds = %263
%270 = sitofp i32 %256 to float
%271 = fdiv float %255, %270
%272 = fcmp olt float %271, %6
br i1 %272, label %273, label %275
273: ; preds = %269
%274 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %274
br label %393
275: ; preds = %269
%276 = call float @llvm.fmuladd.f32(float %148, float %11, float %147)
br label %295
277: ; preds = %295
%278 = or i64 %296, 1
%279 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %278
%280 = load float, float* %279
%281 = fcmp ogt float %280, %276
br i1 %281, label %300, label %282
282: ; preds = %277
%283 = or i64 %296, 2
%284 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %283
%285 = load float, float* %284
%286 = fcmp ogt float %285, %276
br i1 %286, label %300, label %287
287: ; preds = %282
%288 = or i64 %296, 3
%289 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %288
%290 = load float, float* %289
%291 = fcmp ogt float %290, %276
br i1 %291, label %300, label %292
292: ; preds = %287
%293 = add nuw nsw i64 %296, 4
%294 = icmp eq i64 %293, 256
br i1 %294, label %302, label %295
295: ; preds = %292, %275
%296 = phi i64 [ 0, %275 ], [ %293, %292 ]
%297 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %296
%298 = load float, float* %297
%299 = fcmp ogt float %298, %276
br i1 %299, label %300, label %277
300: ; preds = %287, %282, %277, %295
%301 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %301
br label %393
302: ; preds = %292
%303 = fpext float %147 to double
%304 = fpext float %148 to double
%305 = fneg double %304
%306 = call double @llvm.fmuladd.f64(double %305, double 2.000000e-01, double %303)
%307 = fptrunc double %306 to float
br label %308
308: ; preds = %330, %302
%309 = phi i64 [ 0, %302 ], [ %332, %330 ]
%310 = phi i32 [ 0, %302 ], [ %331, %330 ]
%311 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %309
%312 = load float, float* %311
%313 = fcmp olt float %312, %307
%314 = fcmp ogt float %312, %149
%315 = select i1 %313, i1 %314, i1 false
br i1 %315, label %316, label %319
316: ; preds = %308
%317 = add nsw i32 %310, 1
%318 = icmp slt i32 %310, %26
br i1 %318, label %319, label %334
319: ; preds = %308, %316
%320 = phi i32 [ %317, %316 ], [ %310, %308 ]
%321 = or i64 %309, 1
%322 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %321
%323 = load float, float* %322
%324 = fcmp olt float %323, %307
%325 = fcmp ogt float %323, %149
%326 = select i1 %324, i1 %325, i1 false
br i1 %326, label %327, label %330
327: ; preds = %319
%328 = add nsw i32 %320, 1
%329 = icmp slt i32 %320, %26
br i1 %329, label %330, label %334
330: ; preds = %327, %319
%331 = phi i32 [ %328, %327 ], [ %320, %319 ]
%332 = add nuw nsw i64 %309, 2
%333 = icmp eq i64 %332, 256
br i1 %333, label %334, label %308
334: ; preds = %330, %327, %316
%335 = phi i32 [ %317, %316 ], [ %328, %327 ], [ %331, %330 ]
%336 = icmp slt i32 %335, %29
%337 = icmp sgt i32 %335, %26
%338 = select i1 %336, i1 true, i1 %337
br i1 %338, label %339, label %344
339: ; preds = %334
%340 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %340
br label %393
341: ; preds = %344
%342 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %52
store float 0.000000e+00, float* %342
%343 = icmp eq i32 %364, 0
br i1 %343, label %393, label %367
344: ; preds = %334, %344
%345 = phi i64 [ %365, %344 ], [ 0, %334 ]
%346 = phi i32 [ %364, %344 ], [ 0, %334 ]
%347 = phi float [ %362, %344 ], [ 0.000000e+00, %334 ]
%348 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %345
%349 = load float, float* %348
%350 = fcmp ogt float %349, 0.000000e+00
%351 = fsub float %349, %147
%352 = call float @llvm.fmuladd.f32(float %351, float %351, float %347)
%353 = select i1 %350, float %352, float %347
%354 = zext i1 %350 to i32
%355 = add nuw nsw i32 %346, %354
%356 = or i64 %345, 1
%357 = getelementptr inbounds [256 x float], [256 x float]* %16, i64 0, i64 %356
%358 = load float, float* %357
%359 = fcmp ogt float %358, 0.000000e+00
%360 = fsub float %358, %147
%361 = call float @llvm.fmuladd.f32(float %360, float %360, float %353)
%362 = select i1 %359, float %361, float %353
%363 = zext i1 %359 to i32
%364 = add nuw nsw i32 %355, %363
%365 = add nuw nsw i64 %345, 2
%366 = icmp eq i64 %365, 256
br i1 %366, label %341, label %344
367: ; preds = %341
%368 = sitofp i32 %364 to float
%369 = fdiv float %362, %368
%370 = fcmp olt float %369, 0.000000e+00
br i1 %370, label %371, label %373
371: ; preds = %367
%372 = call float @sqrtf(float %369)
br label %373
373: ; preds = %367, %371
%374 = call i32 (i32, i32, float*, float*, i32, double, float*, i32, ...) bitcast (i32 (...)* @dip_vacancy to i32 (i32, i32, float*, float*, i32, double, float*, i32, ...)*)(i32 %258, i32 %257, float* %0, float* %1, i32 %2, double %56, float* %40, i32 256)
%375 = sitofp i32 %374 to float
%376 = fcmp ogt float %375, %12
br i1 %376, label %391, label %377
377: ; preds = %373
%378 = fmul float %148, %270
%379 = fmul float %378, 1.000000e+02
%380 = add nsw i32 %335, 1
%381 = sitofp i32 %380 to float
%382 = insertelement <2 x float> poison, float %381, i64 0
%383 = insertelement <2 x float> %382, float %362, i64 1
%384 = fmul <2 x float> %146, %383
%385 = extractelement <2 x float> %384, i64 0
%386 = fdiv float %379, %385
%387 = extractelement <2 x float> %384, i64 1
%388 = fdiv float %386, %387
%389 = call float @llvm.fmuladd.f32(float %375, float 1.000000e+01, float 1.000000e+00)
%390 = fdiv float %388, %389
br label %391
391: ; preds = %373, %377
%392 = phi float [ %390, %377 ], [ 0.000000e+00, %373 ]
store float %392, float* %342
br label %393
393: ; preds = %86, %300, %341, %391, %339, %273, %267, %261, %246, %51
%394 = add nuw nsw i64 %52, 1
%395 = icmp eq i64 %394, %42
br i1 %395, label %45, label %51
396: ; preds = %396, %49
%397 = phi i64 [ 0, %49 ], [ %418, %396 ]
%398 = phi float [ 0.000000e+00, %49 ], [ %417, %396 ]
%399 = phi float [ 0.000000e+00, %49 ], [ %416, %396 ]
%400 = phi i64 [ 0, %49 ], [ %419, %396 ]
%401 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %397
%402 = load float, float* %401
%403 = fcmp ogt float %402, %399
%404 = trunc i64 %397 to i32
%405 = sitofp i32 %404 to float
%406 = call float @llvm.fmuladd.f32(float %405, float %5, float %3)
%407 = select i1 %403, float %402, float %399
%408 = select i1 %403, float %406, float %398
%409 = or i64 %397, 1
%410 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %409
%411 = load float, float* %410
%412 = fcmp ogt float %411, %407
%413 = trunc i64 %409 to i32
%414 = sitofp i32 %413 to float
%415 = call float @llvm.fmuladd.f32(float %414, float %5, float %3)
%416 = select i1 %412, float %411, float %407
%417 = select i1 %412, float %415, float %408
%418 = add nuw nsw i64 %397, 2
%419 = add i64 %400, 2
%420 = icmp eq i64 %419, %50
br i1 %420, label %421, label %396
421: ; preds = %396, %46
%422 = phi float [ undef, %46 ], [ %417, %396 ]
%423 = phi i64 [ 0, %46 ], [ %418, %396 ]
%424 = phi float [ 0.000000e+00, %46 ], [ %417, %396 ]
%425 = phi float [ 0.000000e+00, %46 ], [ %416, %396 ]
%426 = icmp eq i64 %47, 0
br i1 %426, label %435, label %427
427: ; preds = %421
%428 = getelementptr inbounds [1000000 x float], [1000000 x float]* %15, i64 0, i64 %423
%429 = load float, float* %428
%430 = fcmp ogt float %429, %425
%431 = trunc i64 %423 to i32
%432 = sitofp i32 %431 to float
%433 = call float @llvm.fmuladd.f32(float %432, float %5, float %3)
%434 = select i1 %430, float %433, float %424
br label %435
435: ; preds = %427, %421, %35, %45, %43
%436 = phi float [ 0.000000e+00, %43 ], [ 0.000000e+00, %45 ], [ 0.000000e+00, %35 ], [ %422, %421 ], [ %434, %427 ]
ret float %436
}