float SymhExtension(const float *Src, int Stride, int N, int n) {
  while (1) {
    if (n < 0)
      n = -1 - n;
    else if (n >= N)
      n = 2 * N - 1 - n;
    else
      break;
  }
  return Src[Stride * n];
}


x86 assembly (clang, -O0):
<SymhExtension>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jge    112b <SymhExtension+0x2b>
mov    $0xffffffff,%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    1155 <SymhExtension+0x55>
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jl     114b <SymhExtension+0x4b>
mov    -0x10(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    1150 <SymhExtension+0x50>
jmp    115a <SymhExtension+0x5a>
jmp    1155 <SymhExtension+0x55>
jmp    1111 <SymhExtension+0x11>
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<SymhExtension>:
lea    (%rdx,%rdx,1),%eax
test   %ecx,%ecx
jns    1116 <SymhExtension+0x16>
nopw   0x0(%rax,%rax,1)
not    %ecx
test   %ecx,%ecx
js     1110 <SymhExtension+0x10>
cmp    %edx,%ecx
jl     1124 <SymhExtension+0x24>
not    %ecx
add    %eax,%ecx
test   %ecx,%ecx
jns    1116 <SymhExtension+0x16>
jmp    1110 <SymhExtension+0x10>
imul   %esi,%ecx
movslq %ecx,%rax
movss  (%rdi,%rax,4),%xmm0
ret

LLVM IR (-O0):
define float @SymhExtension(float* %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca float*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
store float* %0, float** %5
store i32 %1, i32* %6
store i32 %2, i32* %7
store i32 %3, i32* %8
br label %9
9: ; preds = %4, %27
%10 = load i32, i32* %8
%11 = icmp slt i32 %10, 0
br i1 %11, label %12, label %15
12: ; preds = %9
%13 = load i32, i32* %8
%14 = sub nsw i32 -1, %13
store i32 %14, i32* %8
br label %27
15: ; preds = %9
%16 = load i32, i32* %8
%17 = load i32, i32* %7
%18 = icmp sge i32 %16, %17
br i1 %18, label %19, label %25
19: ; preds = %15
%20 = load i32, i32* %7
%21 = mul nsw i32 2, %20
%22 = sub nsw i32 %21, 1
%23 = load i32, i32* %8
%24 = sub nsw i32 %22, %23
store i32 %24, i32* %8
br label %26
25: ; preds = %15
br label %28
26: ; preds = %19
br label %27
27: ; preds = %26, %12
br label %9
28: ; preds = %25
%29 = load float*, float** %5
%30 = load i32, i32* %6
%31 = load i32, i32* %8
%32 = mul nsw i32 %30, %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %29, i64 %33
%35 = load float, float* %34
ret float %35
}

LLVM IR (-O2):
define float @SymhExtension(float* nocapture %0, i32 %1, i32 %2, i32 %3) {
%5 = shl nsw i32 %2, 1
br label %6
6: ; preds = %16, %4
%7 = phi i32 [ %3, %4 ], [ %17, %16 ]
%8 = icmp slt i32 %7, 0
br i1 %8, label %9, label %11
9: ; preds = %6
%10 = xor i32 %7, -1
br label %16
11: ; preds = %6
%12 = icmp slt i32 %7, %2
br i1 %12, label %18, label %13
13: ; preds = %11
%14 = xor i32 %7, -1
%15 = add i32 %5, %14
br label %16
16: ; preds = %13, %9
%17 = phi i32 [ %10, %9 ], [ %15, %13 ]
br label %6
18: ; preds = %11
%19 = mul nsw i32 %7, %1
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds float, float* %0, i64 %20
%22 = load float, float* %21
ret float %22
}