int LSF_check(float *lsf, int dim, int NoAn) {
  int k;
  int n;
  int m;
  int Nit = 2;
  int change = 0;
  int pos;
  float tmp;
  static const float eps = 0.039f;
  static const float eps2 = 0.0195f;
  static const float maxlsf = 3.14f;
  static const float minlsf = 0.01f;
  for (n = 0; n < Nit; n++) {
    for (m = 0; m < NoAn; m++) {
      for (k = 0; k < (dim - 1); k++) {
        pos = m * dim + k;
        if ((lsf[pos + 1] - lsf[pos]) < eps) {
          if (lsf[pos + 1] < lsf[pos]) {
            tmp = lsf[pos + 1];
            lsf[pos + 1] = lsf[pos] + eps2;
            lsf[pos] = lsf[pos + 1] - eps2;
          } else {
            lsf[pos] -= eps2;
            lsf[pos + 1] += eps2;
          }
          change = 1;
        }
        if (lsf[pos] < minlsf) {
          lsf[pos] = minlsf;
          change = 1;
        }
        if (lsf[pos] > maxlsf) {
          lsf[pos] = maxlsf;
          change = 1;
        }
      }
    }
  }
  return change;
}


x86 assembly (clang, -O0):
<LSF_check>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x2,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    131a <LSF_check+0x21a>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1307 <LSF_check+0x207>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12f4 <LSF_check+0x1f4>
mov    -0x1c(%rbp),%eax
imul   -0xc(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
movss  0xe72(%rip),%xmm0
ucomiss %xmm1,%xmm0
jbe    1271 <LSF_check+0x171>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss %xmm1,%xmm0
jbe    122d <LSF_check+0x12d>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  0xe1e(%rip),%xmm0
addss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  0xded(%rip),%xmm1
subss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    126a <LSF_check+0x16a>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  0xdc2(%rip),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  0xda4(%rip),%xmm0
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  0xd87(%rip),%xmm0
ucomiss (%rax,%rcx,4),%xmm0
jbe    12a7 <LSF_check+0x1a7>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  0xd6d(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  0xd50(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    12e1 <LSF_check+0x1e1>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  0xd37(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x24(%rbp)
jmp    12e6 <LSF_check+0x1e6>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1149 <LSF_check+0x49>
jmp    12f9 <LSF_check+0x1f9>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1136 <LSF_check+0x36>
jmp    130c <LSF_check+0x20c>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1123 <LSF_check+0x23>
mov    -0x24(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<LSF_check>:
push   %rbx
test   %edx,%edx
jle    12ac <LSF_check+0x1ac>
lea    -0x1(%rsi),%eax
movslq %esi,%r9
mov    %edx,%r8d
mov    %eax,%r10d
lea    0x4(%rdi),%rcx
shl    $0x2,%r9
xor    %r11d,%r11d
movss  0xed8(%rip),%xmm0
movss  0xed8(%rip),%xmm1
movss  0xecc(%rip),%xmm2
movss  0xecc(%rip),%xmm3
movss  0xec8(%rip),%xmm4
xor    %eax,%eax
jmp    1160 <LSF_check+0x60>
nopl   0x0(%rax)
add    $0x1,%r11
add    %r9,%rcx
cmp    %r8,%r11
je     11ea <LSF_check+0xea>
cmp    $0x2,%esi
jl     1150 <LSF_check+0x50>
xor    %ebx,%ebx
jmp    1184 <LSF_check+0x84>
nopl   0x0(%rax)
movss  %xmm6,-0x4(%rcx,%rbx,4)
mov    $0x1,%eax
add    $0x1,%rbx
cmp    %rbx,%r10
je     1150 <LSF_check+0x50>
movss  -0x4(%rcx,%rbx,4),%xmm5
movss  (%rcx,%rbx,4),%xmm6
movaps %xmm6,%xmm7
subss  %xmm5,%xmm7
ucomiss %xmm7,%xmm0
jbe    11d8 <LSF_check+0xd8>
ucomiss %xmm6,%xmm5
jbe    11c0 <LSF_check+0xc0>
addss  %xmm1,%xmm5
movss  %xmm5,(%rcx,%rbx,4)
addss  %xmm2,%xmm5
movss  %xmm5,-0x4(%rcx,%rbx,4)
jmp    11d3 <LSF_check+0xd3>
cs nopw 0x0(%rax,%rax,1)
nop
addss  %xmm2,%xmm5
movss  %xmm5,-0x4(%rcx,%rbx,4)
addss  %xmm1,%xmm6
movss  %xmm6,(%rcx,%rbx,4)
mov    $0x1,%eax
ucomiss %xmm5,%xmm3
movaps %xmm3,%xmm6
ja     1170 <LSF_check+0x70>
ucomiss %xmm4,%xmm5
movaps %xmm4,%xmm6
ja     1170 <LSF_check+0x70>
jmp    117b <LSF_check+0x7b>
test   %edx,%edx
jle    12aa <LSF_check+0x1aa>
add    $0x4,%rdi
xor    %ecx,%ecx
movss  0xe0c(%rip),%xmm3
movss  0xe08(%rip),%xmm4
jmp    1220 <LSF_check+0x120>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rcx
add    %r9,%rdi
cmp    %r8,%rcx
je     12aa <LSF_check+0x1aa>
cmp    $0x2,%esi
jl     1210 <LSF_check+0x110>
xor    %edx,%edx
jmp    1244 <LSF_check+0x144>
nopl   0x0(%rax)
movss  %xmm6,-0x4(%rdi,%rdx,4)
mov    $0x1,%eax
add    $0x1,%rdx
cmp    %rdx,%r10
je     1210 <LSF_check+0x110>
movss  -0x4(%rdi,%rdx,4),%xmm5
movss  (%rdi,%rdx,4),%xmm6
movaps %xmm6,%xmm7
subss  %xmm5,%xmm7
ucomiss %xmm7,%xmm0
jbe    1298 <LSF_check+0x198>
ucomiss %xmm6,%xmm5
jbe    1280 <LSF_check+0x180>
addss  %xmm1,%xmm5
movss  %xmm5,(%rdi,%rdx,4)
addss  %xmm2,%xmm5
movss  %xmm5,-0x4(%rdi,%rdx,4)
jmp    1293 <LSF_check+0x193>
cs nopw 0x0(%rax,%rax,1)
nop
addss  %xmm2,%xmm5
movss  %xmm5,-0x4(%rdi,%rdx,4)
addss  %xmm1,%xmm6
movss  %xmm6,(%rdi,%rdx,4)
mov    $0x1,%eax
ucomiss %xmm5,%xmm3
movaps %xmm3,%xmm6
ja     1230 <LSF_check+0x130>
ucomiss %xmm4,%xmm5
movaps %xmm4,%xmm6
ja     1230 <LSF_check+0x130>
jmp    123b <LSF_check+0x13b>
pop    %rbx
ret
xor    %eax,%eax
pop    %rbx
ret

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

LLVM IR (-O2):
define i32 @LSF_check(float* nocapture %0, i32 %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
%5 = add i32 %1, -1
%6 = icmp sgt i32 %1, 1
%7 = sext i32 %1 to i64
br i1 %4, label %8, label %95
8: ; preds = %3
%9 = zext i32 %2 to i64
%10 = zext i32 %5 to i64
br label %11
11: ; preds = %8, %47
%12 = phi i64 [ 0, %8 ], [ %49, %47 ]
%13 = phi i32 [ 0, %8 ], [ %48, %47 ]
br i1 %6, label %14, label %47
14: ; preds = %11
%15 = mul nsw i64 %12, %7
br label %16
16: ; preds = %14, %43
%17 = phi i64 [ 0, %14 ], [ %45, %43 ]
%18 = phi i32 [ %13, %14 ], [ %44, %43 ]
%19 = add nsw i64 %17, %15
%20 = add nsw i64 %19, 1
%21 = getelementptr inbounds float, float* %0, i64 %20
%22 = load float, float* %21
%23 = getelementptr inbounds float, float* %0, i64 %19
%24 = load float, float* %23
%25 = fsub float %22, %24
%26 = fcmp olt float %25, 0x3FA3F7CEE0000000
br i1 %26, label %27, label %35
27: ; preds = %16
%28 = fcmp olt float %22, %24
br i1 %28, label %29, label %32
29: ; preds = %27
%30 = fadd float %24, 0x3F93F7CEE0000000
store float %30, float* %21
%31 = fadd float %30, 0xBF93F7CEE0000000
store float %31, float* %23
br label %35
32: ; preds = %27
%33 = fadd float %24, 0xBF93F7CEE0000000
store float %33, float* %23
%34 = fadd float %22, 0x3F93F7CEE0000000
store float %34, float* %21
br label %35
35: ; preds = %29, %32, %16
%36 = phi float [ %31, %29 ], [ %33, %32 ], [ %24, %16 ]
%37 = phi i32 [ 1, %29 ], [ 1, %32 ], [ %18, %16 ]
%38 = fcmp olt float %36, 0x3F847AE140000000
br i1 %38, label %41, label %39
39: ; preds = %35
%40 = fcmp ogt float %36, 0x40091EB860000000
br i1 %40, label %41, label %43
41: ; preds = %39, %35
%42 = phi float [ 0x3F847AE140000000, %35 ], [ 0x40091EB860000000, %39 ]
store float %42, float* %23
br label %43
43: ; preds = %41, %39
%44 = phi i32 [ %37, %39 ], [ 1, %41 ]
%45 = add nuw nsw i64 %17, 1
%46 = icmp eq i64 %45, %10
br i1 %46, label %47, label %16
47: ; preds = %43, %11
%48 = phi i32 [ %13, %11 ], [ %44, %43 ]
%49 = add nuw nsw i64 %12, 1
%50 = icmp eq i64 %49, %9
br i1 %50, label %51, label %11
51: ; preds = %47
br i1 %4, label %52, label %95
52: ; preds = %51
%53 = zext i32 %2 to i64
%54 = zext i32 %5 to i64
br label %55
55: ; preds = %91, %52
%56 = phi i64 [ 0, %52 ], [ %93, %91 ]
%57 = phi i32 [ %48, %52 ], [ %92, %91 ]
br i1 %6, label %58, label %91
58: ; preds = %55
%59 = mul nsw i64 %56, %7
br label %60
60: ; preds = %87, %58
%61 = phi i64 [ 0, %58 ], [ %89, %87 ]
%62 = phi i32 [ %57, %58 ], [ %88, %87 ]
%63 = add nsw i64 %61, %59
%64 = add nsw i64 %63, 1
%65 = getelementptr inbounds float, float* %0, i64 %64
%66 = load float, float* %65
%67 = getelementptr inbounds float, float* %0, i64 %63
%68 = load float, float* %67
%69 = fsub float %66, %68
%70 = fcmp olt float %69, 0x3FA3F7CEE0000000
br i1 %70, label %71, label %79
71: ; preds = %60
%72 = fcmp olt float %66, %68
br i1 %72, label %76, label %73
73: ; preds = %71
%74 = fadd float %68, 0xBF93F7CEE0000000
store float %74, float* %67
%75 = fadd float %66, 0x3F93F7CEE0000000
store float %75, float* %65
br label %79
76: ; preds = %71
%77 = fadd float %68, 0x3F93F7CEE0000000
store float %77, float* %65
%78 = fadd float %77, 0xBF93F7CEE0000000
store float %78, float* %67
br label %79
79: ; preds = %76, %73, %60
%80 = phi float [ %78, %76 ], [ %74, %73 ], [ %68, %60 ]
%81 = phi i32 [ 1, %76 ], [ 1, %73 ], [ %62, %60 ]
%82 = fcmp olt float %80, 0x3F847AE140000000
br i1 %82, label %85, label %83
83: ; preds = %79
%84 = fcmp ogt float %80, 0x40091EB860000000
br i1 %84, label %85, label %87
85: ; preds = %79, %83
%86 = phi float [ 0x40091EB860000000, %83 ], [ 0x3F847AE140000000, %79 ]
store float %86, float* %67
br label %87
87: ; preds = %85, %83
%88 = phi i32 [ %81, %83 ], [ 1, %85 ]
%89 = add nuw nsw i64 %61, 1
%90 = icmp eq i64 %89, %54
br i1 %90, label %91, label %60
91: ; preds = %87, %55
%92 = phi i32 [ %57, %55 ], [ %88, %87 ]
%93 = add nuw nsw i64 %56, 1
%94 = icmp eq i64 %93, %53
br i1 %94, label %95, label %55
95: ; preds = %91, %3, %51
%96 = phi i32 [ %48, %51 ], [ 0, %3 ], [ %92, %91 ]
ret i32 %96
}