void demean(s, nt, dt, fdisp) float *s, *dt, *fdisp;
int nt;
{
  int it, j, test;
  float old, new, oabs, nabs, nsgn;
  float difp, tolp, tolm;
  float tol = 0.1;
  float sum = 0.0;
  tolp = 1.0 + tol;
  tolm = 1.0 - tol;
  j = 0;
  test = 1;
  while (test) {
    sum = 0.0;
    for (it = 0; it < j; it++)
      sum = sum + (nt - it) * s[it];
    for (it = j + 1; it < nt; it++)
      sum = sum + (nt - it) * s[it];
    new = *fdisp - sum / (nt - j);
    nabs = new;
    nsgn = 1.0;
    if (nabs < 0.0) {
      nsgn = -1.0;
      nabs = -new;
    }
    old = s[j];
    oabs = old;
    if (oabs < 0.0)
      oabs = -old;
    difp = nabs / oabs;
    if (difp > tolp)
      new = tolp *oabs *nsgn;
    else if (difp < tolm)
      new = tolm *oabs *nsgn;
    else
      test = 0;
    s[j] = new;
    j++;
  }
  printf("demean(): j=%d\n", j);
}


x86 assembly (clang, -O0):
<demean>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movss  0xed9(%rip),%xmm0
movss  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xeb3(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe95(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
je     137b <demean+0x26b>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    11e2 <demean+0xd2>
movss  -0x54(%rbp),%xmm1
mov    -0xc(%rbp),%eax
sub    -0x24(%rbp),%eax
cvtsi2ss %eax,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    119f <demean+0x8f>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    122e <demean+0x11e>
movss  -0x54(%rbp),%xmm1
mov    -0xc(%rbp),%eax
sub    -0x24(%rbp),%eax
cvtsi2ss %eax,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11eb <demean+0xdb>
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  -0x54(%rbp),%xmm1
mov    -0xc(%rbp),%eax
sub    -0x28(%rbp),%eax
cvtsi2ss %eax,%xmm2
divss  %xmm2,%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  0xda8(%rip),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x3c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    12a3 <demean+0x193>
movss  0xd89(%rip),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x34(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    12ec <demean+0x1dc>
movss  -0x30(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x3c(%rbp),%xmm0
divss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
ucomiss -0x48(%rbp),%xmm0
jbe    1323 <demean+0x213>
movss  -0x48(%rbp),%xmm0
mulss  -0x38(%rbp),%xmm0
mulss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
jmp    135b <demean+0x24b>
movss  -0x44(%rbp),%xmm1
movss  -0x4c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    134f <demean+0x23f>
movss  -0x4c(%rbp),%xmm0
mulss  -0x38(%rbp),%xmm0
mulss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
jmp    1356 <demean+0x246>
movl   $0x0,-0x2c(%rbp)
jmp    135b <demean+0x24b>
movss  -0x34(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1186 <demean+0x76>
mov    -0x28(%rbp),%esi
lea    0xc8f(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<demean>:
push   %r14
push   %rbx
mov    %esi,%r11d
lea    -0x1(%rsi),%r8d
xor    %r10d,%r10d
movaps 0xedc(%rip),%xmm0
xorps  %xmm9,%xmm9
movss  0xee0(%rip),%xmm2
movss  0xedb(%rip),%xmm8
mov    %rdi,%r14
jmp    1168 <demean+0x58>
xchg   %ax,%ax
mulss  %xmm3,%xmm6
movaps %xmm4,%xmm1
andnps %xmm6,%xmm1
xorps  %xmm0,%xmm6
andps  %xmm4,%xmm6
orps   %xmm1,%xmm6
movss  %xmm6,(%rdi,%r10,4)
add    $0xffffffff,%r8d
add    $0x4,%r14
add    $0xffffffffffffffff,%r11
mov    %r9,%r10
xorps  %xmm4,%xmm4
test   %r10,%r10
je     122d <demean+0x11d>
lea    -0x1(%r10),%rax
cmp    $0x3,%rax
jae    1190 <demean+0x80>
xorps  %xmm4,%xmm4
xor    %edx,%edx
jmp    11f8 <demean+0xe8>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r10,%r9
and    $0xfffffffffffffffc,%r9
xorps  %xmm4,%xmm4
mov    %esi,%ebx
xor    %edx,%edx
xchg   %ax,%ax
xorps  %xmm3,%xmm3
cvtsi2ss %ebx,%xmm3
mulss  (%rdi,%rdx,4),%xmm3
addss  %xmm4,%xmm3
lea    -0x1(%rbx),%eax
xorps  %xmm5,%xmm5
cvtsi2ss %eax,%xmm5
lea    -0x2(%rbx),%eax
xorps  %xmm6,%xmm6
cvtsi2ss %eax,%xmm6
mulss  0x4(%rdi,%rdx,4),%xmm5
mulss  0x8(%rdi,%rdx,4),%xmm6
lea    -0x3(%rbx),%eax
xorps  %xmm4,%xmm4
cvtsi2ss %eax,%xmm4
addss  %xmm3,%xmm5
addss  %xmm5,%xmm6
mulss  0xc(%rdi,%rdx,4),%xmm4
addss  %xmm6,%xmm4
add    $0x4,%rdx
add    $0xfffffffc,%ebx
cmp    %rdx,%r9
jne    11a0 <demean+0x90>
test   $0x3,%r10b
je     122d <demean+0x11d>
mov    %r10d,%eax
and    $0x3,%eax
lea    (%rdi,%rdx,4),%r9
mov    %esi,%ebx
sub    %edx,%ebx
xor    %edx,%edx
xchg   %ax,%ax
xorps  %xmm3,%xmm3
cvtsi2ss %ebx,%xmm3
mulss  (%r9,%rdx,4),%xmm3
addss  %xmm3,%xmm4
add    $0x1,%rdx
add    $0xffffffff,%ebx
cmp    %rdx,%rax
jne    1210 <demean+0x100>
lea    0x1(%r10),%r9
cmp    %esi,%r9d
jge    125d <demean+0x14d>
mov    $0x1,%eax
mov    %r8d,%edx
xchg   %ax,%ax
xorps  %xmm3,%xmm3
cvtsi2ss %edx,%xmm3
mulss  (%r14,%rax,4),%xmm3
addss  %xmm3,%xmm4
add    $0xffffffff,%edx
add    $0x1,%rax
cmp    %rax,%r11
jne    1240 <demean+0x130>
movss  (%rcx),%xmm5
mov    %esi,%eax
sub    %r10d,%eax
xorps  %xmm3,%xmm3
cvtsi2ss %eax,%xmm3
divss  %xmm3,%xmm4
subss  %xmm4,%xmm5
movaps %xmm5,%xmm3
xorps  %xmm0,%xmm3
movaps %xmm5,%xmm4
cmpltss %xmm9,%xmm4
andps  %xmm4,%xmm3
movaps %xmm4,%xmm7
andnps %xmm5,%xmm7
orps   %xmm3,%xmm7
movss  (%rdi,%r10,4),%xmm1
movaps %xmm1,%xmm3
xorps  %xmm0,%xmm3
movaps %xmm1,%xmm6
cmpltss %xmm9,%xmm6
andps  %xmm6,%xmm3
andnps %xmm1,%xmm6
orps   %xmm3,%xmm6
divss  %xmm6,%xmm7
ucomiss %xmm2,%xmm7
movaps %xmm2,%xmm3
ja     1140 <demean+0x30>
ucomiss %xmm7,%xmm8
movaps %xmm8,%xmm3
ja     1140 <demean+0x30>
movss  %xmm5,(%rdi,%r10,4)
lea    0xd3f(%rip),%rdi
mov    %r9d,%esi
xor    %eax,%eax
pop    %rbx
pop    %r14
jmp    1030 <printf@plt>

LLVM IR (-O0):
@.str = [16 x i8] c"demean(): j=%d\0A\00"
define void @demean(float* %0, i32 %1, float* %2, float* %3) {
%5 = alloca float*
%6 = alloca i32
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = alloca float
%20 = alloca float
%21 = alloca float
store float* %0, float** %5
store i32 %1, i32* %6
store float* %2, float** %7
store float* %3, float** %8
store float 0x3FB99999A0000000, float* %20
store float 0.000000e+00, float* %21
%22 = load float, float* %20
%23 = fpext float %22 to double
%24 = fadd double 1.000000e+00, %23
%25 = fptrunc double %24 to float
store float %25, float* %18
%26 = load float, float* %20
%27 = fpext float %26 to double
%28 = fsub double 1.000000e+00, %27
%29 = fptrunc double %28 to float
store float %29, float* %19
store i32 0, i32* %10
store i32 1, i32* %11
br label %30
30: ; preds = %130, %4
%31 = load i32, i32* %11
%32 = icmp ne i32 %31, 0
br i1 %32, label %33, label %138
33: ; preds = %30
store float 0.000000e+00, float* %21
store i32 0, i32* %9
br label %34
34: ; preds = %50, %33
%35 = load i32, i32* %9
%36 = load i32, i32* %10
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %53
38: ; preds = %34
%39 = load float, float* %21
%40 = load i32, i32* %6
%41 = load i32, i32* %9
%42 = sub nsw i32 %40, %41
%43 = sitofp i32 %42 to float
%44 = load float*, float** %5
%45 = load i32, i32* %9
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds float, float* %44, i64 %46
%48 = load float, float* %47
%49 = call float @llvm.fmuladd.f32(float %43, float %48, float %39)
store float %49, float* %21
br label %50
50: ; preds = %38
%51 = load i32, i32* %9
%52 = add nsw i32 %51, 1
store i32 %52, i32* %9
br label %34
53: ; preds = %34
%54 = load i32, i32* %10
%55 = add nsw i32 %54, 1
store i32 %55, i32* %9
br label %56
56: ; preds = %72, %53
%57 = load i32, i32* %9
%58 = load i32, i32* %6
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %75
60: ; preds = %56
%61 = load float, float* %21
%62 = load i32, i32* %6
%63 = load i32, i32* %9
%64 = sub nsw i32 %62, %63
%65 = sitofp i32 %64 to float
%66 = load float*, float** %5
%67 = load i32, i32* %9
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds float, float* %66, i64 %68
%70 = load float, float* %69
%71 = call float @llvm.fmuladd.f32(float %65, float %70, float %61)
store float %71, float* %21
br label %72
72: ; preds = %60
%73 = load i32, i32* %9
%74 = add nsw i32 %73, 1
store i32 %74, i32* %9
br label %56
75: ; preds = %56
%76 = load float*, float** %8
%77 = load float, float* %76
%78 = load float, float* %21
%79 = load i32, i32* %6
%80 = load i32, i32* %10
%81 = sub nsw i32 %79, %80
%82 = sitofp i32 %81 to float
%83 = fdiv float %78, %82
%84 = fsub float %77, %83
store float %84, float* %13
%85 = load float, float* %13
store float %85, float* %15
store float 1.000000e+00, float* %16
%86 = load float, float* %15
%87 = fpext float %86 to double
%88 = fcmp olt double %87, 0.000000e+00
br i1 %88, label %89, label %92
89: ; preds = %75
store float -1.000000e+00, float* %16
%90 = load float, float* %13
%91 = fneg float %90
store float %91, float* %15
br label %92
92: ; preds = %89, %75
%93 = load float*, float** %5
%94 = load i32, i32* %10
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds float, float* %93, i64 %95
%97 = load float, float* %96
store float %97, float* %12
%98 = load float, float* %12
store float %98, float* %14
%99 = load float, float* %14
%100 = fpext float %99 to double
%101 = fcmp olt double %100, 0.000000e+00
br i1 %101, label %102, label %105
102: ; preds = %92
%103 = load float, float* %12
%104 = fneg float %103
store float %104, float* %14
br label %105
105: ; preds = %102, %92
%106 = load float, float* %15
%107 = load float, float* %14
%108 = fdiv float %106, %107
store float %108, float* %17
%109 = load float, float* %17
%110 = load float, float* %18
%111 = fcmp ogt float %109, %110
br i1 %111, label %112, label %118
112: ; preds = %105
%113 = load float, float* %18
%114 = load float, float* %14
%115 = fmul float %113, %114
%116 = load float, float* %16
%117 = fmul float %115, %116
store float %117, float* %13
br label %130
118: ; preds = %105
%119 = load float, float* %17
%120 = load float, float* %19
%121 = fcmp olt float %119, %120
br i1 %121, label %122, label %128
122: ; preds = %118
%123 = load float, float* %19
%124 = load float, float* %14
%125 = fmul float %123, %124
%126 = load float, float* %16
%127 = fmul float %125, %126
store float %127, float* %13
br label %129
128: ; preds = %118
store i32 0, i32* %11
br label %129
129: ; preds = %128, %122
br label %130
130: ; preds = %129, %112
%131 = load float, float* %13
%132 = load float*, float** %5
%133 = load i32, i32* %10
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds float, float* %132, i64 %134
store float %131, float* %135
%136 = load i32, i32* %10
%137 = add nsw i32 %136, 1
store i32 %137, i32* %10
br label %30
138: ; preds = %30
%139 = load i32, i32* %10
%140 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([16 x i8], [16 x i8]* @.str, i64 0, i64 0), i32 %139)
ret void
}

LLVM IR (-O2):
@.str = [16 x i8] c"demean(): j=%d\0A\00"
define void @demean(float* nocapture %0, i32 %1, float* nocapture readnone %2, float* nocapture %3) {
%5 = zext i32 %1 to i64
br label %6
6: ; preds = %104, %4
%7 = phi i64 [ 0, %4 ], [ %70, %104 ]
%8 = phi i64 [ 1, %4 ], [ %109, %104 ]
%9 = icmp eq i64 %7, 0
br i1 %9, label %68, label %10
10: ; preds = %6
%11 = add i64 %7, -1
%12 = and i64 %7, 3
%13 = icmp ult i64 %11, 3
br i1 %13, label %50, label %14
14: ; preds = %10
%15 = and i64 %7, -4
br label %16
16: ; preds = %16, %14
%17 = phi i64 [ 0, %14 ], [ %47, %16 ]
%18 = phi float [ 0.000000e+00, %14 ], [ %46, %16 ]
%19 = phi i64 [ 0, %14 ], [ %48, %16 ]
%20 = trunc i64 %17 to i32
%21 = sub nsw i32 %1, %20
%22 = sitofp i32 %21 to float
%23 = getelementptr inbounds float, float* %0, i64 %17
%24 = load float, float* %23
%25 = tail call float @llvm.fmuladd.f32(float %22, float %24, float %18)
%26 = or i64 %17, 1
%27 = trunc i64 %26 to i32
%28 = sub nsw i32 %1, %27
%29 = sitofp i32 %28 to float
%30 = getelementptr inbounds float, float* %0, i64 %26
%31 = load float, float* %30
%32 = tail call float @llvm.fmuladd.f32(float %29, float %31, float %25)
%33 = or i64 %17, 2
%34 = trunc i64 %33 to i32
%35 = sub nsw i32 %1, %34
%36 = sitofp i32 %35 to float
%37 = getelementptr inbounds float, float* %0, i64 %33
%38 = load float, float* %37
%39 = tail call float @llvm.fmuladd.f32(float %36, float %38, float %32)
%40 = or i64 %17, 3
%41 = trunc i64 %40 to i32
%42 = sub nsw i32 %1, %41
%43 = sitofp i32 %42 to float
%44 = getelementptr inbounds float, float* %0, i64 %40
%45 = load float, float* %44
%46 = tail call float @llvm.fmuladd.f32(float %43, float %45, float %39)
%47 = add nuw nsw i64 %17, 4
%48 = add i64 %19, 4
%49 = icmp eq i64 %48, %15
br i1 %49, label %50, label %16
50: ; preds = %16, %10
%51 = phi float [ undef, %10 ], [ %46, %16 ]
%52 = phi i64 [ 0, %10 ], [ %47, %16 ]
%53 = phi float [ 0.000000e+00, %10 ], [ %46, %16 ]
%54 = icmp eq i64 %12, 0
br i1 %54, label %68, label %55
55: ; preds = %50, %55
%56 = phi i64 [ %65, %55 ], [ %52, %50 ]
%57 = phi float [ %64, %55 ], [ %53, %50 ]
%58 = phi i64 [ %66, %55 ], [ 0, %50 ]
%59 = trunc i64 %56 to i32
%60 = sub nsw i32 %1, %59
%61 = sitofp i32 %60 to float
%62 = getelementptr inbounds float, float* %0, i64 %56
%63 = load float, float* %62
%64 = tail call float @llvm.fmuladd.f32(float %61, float %63, float %57)
%65 = add nuw nsw i64 %56, 1
%66 = add i64 %58, 1
%67 = icmp eq i64 %66, %12
br i1 %67, label %68, label %55
68: ; preds = %50, %55, %6
%69 = phi float [ 0.000000e+00, %6 ], [ %51, %50 ], [ %64, %55 ]
%70 = add nuw i64 %7, 1
%71 = trunc i64 %70 to i32
%72 = icmp slt i32 %71, %1
br i1 %72, label %73, label %84
73: ; preds = %68, %73
%74 = phi i64 [ %82, %73 ], [ %8, %68 ]
%75 = phi float [ %81, %73 ], [ %69, %68 ]
%76 = trunc i64 %74 to i32
%77 = sub nsw i32 %1, %76
%78 = sitofp i32 %77 to float
%79 = getelementptr inbounds float, float* %0, i64 %74
%80 = load float, float* %79
%81 = tail call float @llvm.fmuladd.f32(float %78, float %80, float %75)
%82 = add i64 %74, 1
%83 = icmp eq i64 %82, %5
br i1 %83, label %84, label %73
84: ; preds = %73, %68
%85 = phi float [ %69, %68 ], [ %81, %73 ]
%86 = load float, float* %3
%87 = trunc i64 %7 to i32
%88 = sub nsw i32 %1, %87
%89 = sitofp i32 %88 to float
%90 = fdiv float %85, %89
%91 = fsub float %86, %90
%92 = fcmp olt float %91, 0.000000e+00
%93 = fneg float %91
%94 = select i1 %92, float %93, float %91
%95 = getelementptr inbounds float, float* %0, i64 %7
%96 = load float, float* %95
%97 = fcmp olt float %96, 0.000000e+00
%98 = fneg float %96
%99 = select i1 %97, float %98, float %96
%100 = fdiv float %94, %99
%101 = fcmp ogt float %100, 0x3FF19999A0000000
br i1 %101, label %104, label %102
102: ; preds = %84
%103 = fcmp olt float %100, 0x3FECCCCCC0000000
br i1 %103, label %104, label %110
104: ; preds = %102, %84
%105 = phi float [ 0x3FF19999A0000000, %84 ], [ 0x3FECCCCCC0000000, %102 ]
%106 = fmul float %99, %105
%107 = fneg float %106
%108 = select i1 %92, float %107, float %106
store float %108, float* %95
%109 = add nuw i64 %8, 1
br label %6
110: ; preds = %102
%111 = trunc i64 %70 to i32
store float %91, float* %95
%112 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([16 x i8], [16 x i8]* @.str, i64 0, i64 0), i32 %111)
ret void
}