void detect_lobe_spike_creation(int x, int y, int **stim_matrix,
                                float frequency, float phase, int t) {
  if (frequency == 0) {
    stim_matrix[x][y] = 0;
    return;
  }
  int T = (int)(20000.0 / frequency);
  int offset_steps = (int)((phase / 360.0) * T);
  if ((t >= offset_steps) && (((t - offset_steps) % T) == 0)) {
    stim_matrix[x][y] = 500;
  } else
    stim_matrix[x][y] = 0;
  return;
}


x86 assembly (clang, -O0):
<detect_lobe_spike_creation>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
movss  %xmm1,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movss  -0x14(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    114e <detect_lobe_spike_creation+0x4e>
jp     114e <detect_lobe_spike_creation+0x4e>
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
jmp    11e6 <detect_lobe_spike_creation+0xe6>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xea9(%rip),%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x20(%rbp)
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe85(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsi2sdl -0x20(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x24(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jl     11ca <detect_lobe_spike_creation+0xca>
mov    -0x1c(%rbp),%eax
sub    -0x24(%rbp),%eax
cltd
idivl  -0x20(%rbp)
cmp    $0x0,%edx
jne    11ca <detect_lobe_spike_creation+0xca>
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8(%rbp),%rcx
movl   $0x1f4,(%rax,%rcx,4)
jmp    11e1 <detect_lobe_spike_creation+0xe1>
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
jmp    11e6 <detect_lobe_spike_creation+0xe6>
pop    %rbp
ret

x86 assembly (clang, -O2):
<detect_lobe_spike_creation>:
mov    %ecx,%eax
mov    %rdx,%r8
xor    %ecx,%ecx
xorps  %xmm2,%xmm2
ucomiss %xmm2,%xmm0
jne    1111 <detect_lobe_spike_creation+0x11>
jnp    1158 <detect_lobe_spike_creation+0x58>
xorps  %xmm2,%xmm2
cvtss2sd %xmm0,%xmm2
movsd  0xee0(%rip),%xmm0
divsd  %xmm2,%xmm0
cvtss2sd %xmm1,%xmm1
divsd  0xed8(%rip),%xmm1
cvttpd2dq %xmm0,%xmm2
cvtdq2pd %xmm2,%xmm2
mulsd  %xmm1,%xmm2
cvttsd2si %xmm2,%edx
cmp    %eax,%edx
jg     1158 <detect_lobe_spike_creation+0x58>
cvttsd2si %xmm0,%r9d
sub    %edx,%eax
cltd
idiv   %r9d
test   %edx,%edx
jne    1158 <detect_lobe_spike_creation+0x58>
mov    $0x1f4,%ecx
movslq %edi,%rax
mov    (%r8,%rax,8),%rax
movslq %esi,%rdx
mov    %ecx,(%rax,%rdx,4)
ret

LLVM IR (-O0):
define void @detect_lobe_spike_creation(i32 %0, i32 %1, i32** %2, float %3, float %4, i32 %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32**
%10 = alloca float
%11 = alloca float
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
store i32 %0, i32* %7
store i32 %1, i32* %8
store i32** %2, i32*** %9
store float %3, float* %10
store float %4, float* %11
store i32 %5, i32* %12
%15 = load float, float* %10
%16 = fcmp oeq float %15, 0.000000e+00
br i1 %16, label %17, label %26
17: ; preds = %6
%18 = load i32**, i32*** %9
%19 = load i32, i32* %7
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds i32*, i32** %18, i64 %20
%22 = load i32*, i32** %21
%23 = load i32, i32* %8
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds i32, i32* %22, i64 %24
store i32 0, i32* %25
br label %67
26: ; preds = %6
%27 = load float, float* %10
%28 = fpext float %27 to double
%29 = fdiv double 2.000000e+04, %28
%30 = fptosi double %29 to i32
store i32 %30, i32* %13
%31 = load float, float* %11
%32 = fpext float %31 to double
%33 = fdiv double %32, 3.600000e+02
%34 = load i32, i32* %13
%35 = sitofp i32 %34 to double
%36 = fmul double %33, %35
%37 = fptosi double %36 to i32
store i32 %37, i32* %14
%38 = load i32, i32* %12
%39 = load i32, i32* %14
%40 = icmp sge i32 %38, %39
br i1 %40, label %41, label %57
41: ; preds = %26
%42 = load i32, i32* %12
%43 = load i32, i32* %14
%44 = sub nsw i32 %42, %43
%45 = load i32, i32* %13
%46 = srem i32 %44, %45
%47 = icmp eq i32 %46, 0
br i1 %47, label %48, label %57
48: ; preds = %41
%49 = load i32**, i32*** %9
%50 = load i32, i32* %7
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds i32*, i32** %49, i64 %51
%53 = load i32*, i32** %52
%54 = load i32, i32* %8
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds i32, i32* %53, i64 %55
store i32 500, i32* %56
br label %66
57: ; preds = %41, %26
%58 = load i32**, i32*** %9
%59 = load i32, i32* %7
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds i32*, i32** %58, i64 %60
%62 = load i32*, i32** %61
%63 = load i32, i32* %8
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds i32, i32* %62, i64 %64
store i32 0, i32* %65
br label %66
66: ; preds = %57, %48
br label %67
67: ; preds = %66, %17
ret void
}

LLVM IR (-O2):
define void @detect_lobe_spike_creation(i32 %0, i32 %1, i32** nocapture %2, float %3, float %4, i32 %5) {
%7 = fcmp oeq float %3, 0.000000e+00
br i1 %7, label %23, label %8
8: ; preds = %6
%9 = fpext float %3 to double
%10 = fdiv double 2.000000e+04, %9
%11 = fptosi double %10 to i32
%12 = fpext float %4 to double
%13 = fdiv double %12, 3.600000e+02
%14 = sitofp i32 %11 to double
%15 = fmul double %13, %14
%16 = fptosi double %15 to i32
%17 = icmp sgt i32 %16, %5
br i1 %17, label %22, label %18
18: ; preds = %8
%19 = sub nsw i32 %5, %16
%20 = srem i32 %19, %11
%21 = icmp eq i32 %20, 0
br i1 %21, label %23, label %22
22: ; preds = %18, %8
br label %23
23: ; preds = %18, %6, %22
%24 = phi i32 [ 0, %22 ], [ 0, %6 ], [ 500, %18 ]
%25 = sext i32 %0 to i64
%26 = getelementptr inbounds i32*, i32** %2, i64 %25
%27 = load i32*, i32** %26
%28 = sext i32 %1 to i64
%29 = getelementptr inbounds i32, i32* %27, i64 %28
store i32 %24, i32* %29
ret void
}