int LSF_check(float *lsf, int dim, int NoAn) {
  int k, n, m, Nit = 2, change = 0, pos;
  float tmp;
  static float eps = (float)0.039;
  static float eps2 = (float)0.0195;
  static float maxlsf = (float)3.14;
  static float minlsf = (float)0.01;
  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    1316 <LSF_check+0x216>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1303 <LSF_check+0x203>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12f0 <LSF_check+0x1f0>
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  0x2e92(%rip),%xmm0
ucomiss %xmm1,%xmm0
jbe    126d <LSF_check+0x16d>
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    1229 <LSF_check+0x129>
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  (%rax,%rcx,4),%xmm0
addss  0x2e39(%rip),%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
subss  0x2e0d(%rip),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1266 <LSF_check+0x166>
movss  0x2df3(%rip),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  0x2dd5(%rip),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0x2daa(%rip),%xmm0
ucomiss %xmm1,%xmm0
jbe    12a7 <LSF_check+0x1a7>
movss  0x2d99(%rip),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss 0x2d6d(%rip),%xmm0
jbe    12dd <LSF_check+0x1dd>
movss  0x2d5f(%rip),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x24(%rbp)
jmp    12e2 <LSF_check+0x1e2>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1149 <LSF_check+0x49>
jmp    12f5 <LSF_check+0x1f5>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1136 <LSF_check+0x36>
jmp    1308 <LSF_check+0x208>
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 global float 0x3FA3F7CEE0000000
@LSF_check.eps2 = internal global float 0x3F93F7CEE0000000
@LSF_check.maxlsf = internal global float 0x40091EB860000000
@LSF_check.minlsf = internal global 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 = %147, %3
%15 = load i32, i32* %8
%16 = load i32, i32* %10
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %150
18: ; preds = %14
store i32 0, i32* %9
br label %19
19: ; preds = %143, %18
%20 = load i32, i32* %9
%21 = load i32, i32* %6
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %146
23: ; preds = %19
store i32 0, i32* %7
br label %24
24: ; preds = %139, %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 %142
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 = load float, float* @LSF_check.eps
%48 = fcmp olt float %46, %47
br i1 %48, label %49, label %110
49: ; preds = %29
%50 = load float*, float** %4
%51 = load i32, i32* %12
%52 = add nsw i32 %51, 1
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds float, float* %50, i64 %53
%55 = load float, float* %54
%56 = load float*, float** %4
%57 = load i32, i32* %12
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %56, i64 %58
%60 = load float, float* %59
%61 = fcmp olt float %55, %60
br i1 %61, label %62, label %93
62: ; preds = %49
%63 = load float*, float** %4
%64 = load i32, i32* %12
%65 = add nsw i32 %64, 1
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds float, float* %63, i64 %66
%68 = load float, float* %67
store float %68, float* %13
%69 = load float*, float** %4
%70 = load i32, i32* %12
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds float, float* %69, i64 %71
%73 = load float, float* %72
%74 = load float, float* @LSF_check.eps2
%75 = fadd float %73, %74
%76 = load float*, float** %4
%77 = load i32, i32* %12
%78 = add nsw i32 %77, 1
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %76, i64 %79
store float %75, float* %80
%81 = load float*, float** %4
%82 = load i32, i32* %12
%83 = add nsw i32 %82, 1
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %81, i64 %84
%86 = load float, float* %85
%87 = load float, float* @LSF_check.eps2
%88 = fsub float %86, %87
%89 = load float*, float** %4
%90 = load i32, i32* %12
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds float, float* %89, i64 %91
store float %88, float* %92
br label %109
93: ; preds = %49
%94 = load float, float* @LSF_check.eps2
%95 = load float*, float** %4
%96 = load i32, i32* %12
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds float, float* %95, i64 %97
%99 = load float, float* %98
%100 = fsub float %99, %94
store float %100, float* %98
%101 = load float, float* @LSF_check.eps2
%102 = load float*, float** %4
%103 = load i32, i32* %12
%104 = add nsw i32 %103, 1
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds float, float* %102, i64 %105
%107 = load float, float* %106
%108 = fadd float %107, %101
store float %108, float* %106
br label %109
109: ; preds = %93, %62
store i32 1, i32* %11
br label %110
110: ; preds = %109, %29
%111 = load float*, float** %4
%112 = load i32, i32* %12
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds float, float* %111, i64 %113
%115 = load float, float* %114
%116 = load float, float* @LSF_check.minlsf
%117 = fcmp olt float %115, %116
br i1 %117, label %118, label %124
118: ; preds = %110
%119 = load float, float* @LSF_check.minlsf
%120 = load float*, float** %4
%121 = load i32, i32* %12
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds float, float* %120, i64 %122
store float %119, float* %123
store i32 1, i32* %11
br label %124
124: ; preds = %118, %110
%125 = load float*, float** %4
%126 = load i32, i32* %12
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds float, float* %125, i64 %127
%129 = load float, float* %128
%130 = load float, float* @LSF_check.maxlsf
%131 = fcmp ogt float %129, %130
br i1 %131, label %132, label %138
132: ; preds = %124
%133 = load float, float* @LSF_check.maxlsf
%134 = load float*, float** %4
%135 = load i32, i32* %12
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds float, float* %134, i64 %136
store float %133, float* %137
store i32 1, i32* %11
br label %138
138: ; preds = %132, %124
br label %139
139: ; preds = %138
%140 = load i32, i32* %7
%141 = add nsw i32 %140, 1
store i32 %141, i32* %7
br label %24
142: ; preds = %24
br label %143
143: ; preds = %142
%144 = load i32, i32* %9
%145 = add nsw i32 %144, 1
store i32 %145, i32* %9
br label %19
146: ; preds = %19
br label %147
147: ; preds = %146
%148 = load i32, i32* %8
%149 = add nsw i32 %148, 1
store i32 %149, i32* %8
br label %14
150: ; preds = %14
%151 = load i32, i32* %11
ret i32 %151
}

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
}