void detrend(s, nt, dt, iv) float *s, *dt, *iv;
int nt;
{
  int it, np;
  float *s1, b0;
  np = nt / 4;
  s1 = (float *)check_malloc(nt * sizeof(float));
  for (it = 0; it < nt; it++)
    s1[it] = s[it];
  integrate(s1, nt, dt, iv);
  get_trend(s1 + (nt - np), np, dt, &b0);
  s[0] = s[0] - b0 / (*dt);
  printf("detrend(): b0=%f\n", b0);
  free(s1);
}


x86 assembly (clang, -O0):
<detrend>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,-0x28(%rbp)
movslq -0xc(%rbp),%rdi
shl    $0x2,%rdi
mov    $0x0,%al
call   1070 <check_malloc@plt>
cltq
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11c5 <detrend+0x75>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1191 <detrend+0x41>
mov    -0x30(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x18(%rbp),%rdx
mov    -0x20(%rbp),%rcx
mov    $0x0,%al
call   1040 <integrate@plt>
mov    -0x30(%rbp),%rdi
mov    -0xc(%rbp),%eax
sub    -0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x28(%rbp),%esi
mov    -0x18(%rbp),%rdx
lea    -0x34(%rbp),%rcx
mov    $0x0,%al
call   1060 <get_trend@plt>
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  -0x34(%rbp),%xmm1
mov    -0x18(%rbp),%rax
divss  (%rax),%xmm1
subss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xdcf(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<detrend>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %esi,%ebp
lea    0x3(%rbp),%r12d
test   %esi,%esi
cmovns %esi,%r12d
mov    %rcx,0x10(%rsp)
mov    %rdx,%r14
mov    %rdi,%r15
sar    $0x2,%r12d
movslq %esi,%r13
lea    0x0(,%r13,4),%rdi
xor    %eax,%eax
call   1070 <check_malloc@plt>
movslq %eax,%rbx
test   %r13d,%r13d
jle    1304 <detrend+0x1b4>
mov    %ebp,%eax
cmp    $0x8,%ebp
jb     11b8 <detrend+0x68>
lea    (%r15,%rax,4),%rcx
cmp    %rcx,%rbx
jae    1230 <detrend+0xe0>
lea    (%rbx,%rax,4),%rcx
cmp    %r15,%rcx
jbe    1230 <detrend+0xe0>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rax,%rsi
and    $0x3,%rsi
je     11e5 <detrend+0x95>
nopl   0x0(%rax)
movss  (%r15,%rcx,4),%xmm0
movss  %xmm0,(%rbx,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    11d0 <detrend+0x80>
cmp    $0x3,%rdx
jb     1304 <detrend+0x1b4>
nop
movss  (%r15,%rcx,4),%xmm0
movss  %xmm0,(%rbx,%rcx,4)
movss  0x4(%r15,%rcx,4),%xmm0
movss  %xmm0,0x4(%rbx,%rcx,4)
movss  0x8(%r15,%rcx,4),%xmm0
movss  %xmm0,0x8(%rbx,%rcx,4)
movss  0xc(%r15,%rcx,4),%xmm0
movss  %xmm0,0xc(%rbx,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    11f0 <detrend+0xa0>
jmp    1304 <detrend+0x1b4>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rdx
mov    %rdx,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %esi,%r8d
and    $0x3,%r8d
cmp    $0x18,%rdx
jae    1255 <detrend+0x105>
xor    %edi,%edi
jmp    12c0 <detrend+0x170>
and    $0xfffffffffffffffc,%rsi
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movups (%r15,%rdi,4),%xmm0
movups 0x10(%r15,%rdi,4),%xmm1
movups %xmm0,(%rbx,%rdi,4)
movups %xmm1,0x10(%rbx,%rdi,4)
movups 0x20(%r15,%rdi,4),%xmm0
movups 0x30(%r15,%rdi,4),%xmm1
movups %xmm0,0x20(%rbx,%rdi,4)
movups %xmm1,0x30(%rbx,%rdi,4)
movups 0x40(%r15,%rdi,4),%xmm0
movups 0x50(%r15,%rdi,4),%xmm1
movups %xmm0,0x40(%rbx,%rdi,4)
movups %xmm1,0x50(%rbx,%rdi,4)
movups 0x60(%r15,%rdi,4),%xmm0
movups 0x70(%r15,%rdi,4),%xmm1
movups %xmm0,0x60(%rbx,%rdi,4)
movups %xmm1,0x70(%rbx,%rdi,4)
add    $0x20,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1260 <detrend+0x110>
test   %r8,%r8
je     12fb <detrend+0x1ab>
lea    (%rbx,%rdi,4),%rsi
add    $0x10,%rsi
lea    (%r15,%rdi,4),%rdi
add    $0x10,%rdi
shl    $0x5,%r8
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movups -0x10(%rdi,%rdx,1),%xmm0
movups (%rdi,%rdx,1),%xmm1
movups %xmm0,-0x10(%rsi,%rdx,1)
movups %xmm1,(%rsi,%rdx,1)
add    $0x20,%rdx
cmp    %rdx,%r8
jne    12e0 <detrend+0x190>
cmp    %rax,%rcx
jne    11ba <detrend+0x6a>
mov    %rbx,%rdi
mov    %ebp,%esi
mov    %r14,%rdx
mov    0x10(%rsp),%rcx
xor    %eax,%eax
call   1040 <integrate@plt>
sub    %r12d,%ebp
movslq %ebp,%rax
lea    (%rbx,%rax,4),%rdi
lea    0xc(%rsp),%rcx
mov    %r12d,%esi
mov    %r14,%rdx
xor    %eax,%eax
call   1060 <get_trend@plt>
movss  (%r15),%xmm1
movss  0xc(%rsp),%xmm2
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
divss  (%r14),%xmm2
subss  %xmm2,%xmm1
movss  %xmm1,(%r15)
lea    0xca5(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [18 x i8] c"detrend(): b0=%f\0A\00"
define void @detrend(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 float*
%12 = alloca float
store float* %0, float** %5
store i32 %1, i32* %6
store float* %2, float** %7
store float* %3, float** %8
%13 = load i32, i32* %6
%14 = sdiv i32 %13, 4
store i32 %14, i32* %10
%15 = load i32, i32* %6
%16 = sext i32 %15 to i64
%17 = mul i64 %16, 4
%18 = call i32 (i64, ...) bitcast (i32 (...)* @check_malloc to i32 (i64, ...)*)(i64 %17)
%19 = sext i32 %18 to i64
%20 = inttoptr i64 %19 to float*
store float* %20, float** %11
store i32 0, i32* %9
br label %21
21: ; preds = %35, %4
%22 = load i32, i32* %9
%23 = load i32, i32* %6
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %38
25: ; preds = %21
%26 = load float*, float** %5
%27 = load i32, i32* %9
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %26, i64 %28
%30 = load float, float* %29
%31 = load float*, float** %11
%32 = load i32, i32* %9
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %31, i64 %33
store float %30, float* %34
br label %35
35: ; preds = %25
%36 = load i32, i32* %9
%37 = add nsw i32 %36, 1
store i32 %37, i32* %9
br label %21
38: ; preds = %21
%39 = load float*, float** %11
%40 = load i32, i32* %6
%41 = load float*, float** %7
%42 = load float*, float** %8
%43 = call i32 (float*, i32, float*, float*, ...) bitcast (i32 (...)* @integrate to i32 (float*, i32, float*, float*, ...)*)(float* %39, i32 %40, float* %41, float* %42)
%44 = load float*, float** %11
%45 = load i32, i32* %6
%46 = load i32, i32* %10
%47 = sub nsw i32 %45, %46
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %44, i64 %48
%50 = load i32, i32* %10
%51 = load float*, float** %7
%52 = call i32 (float*, i32, float*, float*, ...) bitcast (i32 (...)* @get_trend to i32 (float*, i32, float*, float*, ...)*)(float* %49, i32 %50, float* %51, float* %12)
%53 = load float*, float** %5
%54 = getelementptr inbounds float, float* %53, i64 0
%55 = load float, float* %54
%56 = load float, float* %12
%57 = load float*, float** %7
%58 = load float, float* %57
%59 = fdiv float %56, %58
%60 = fsub float %55, %59
%61 = load float*, float** %5
%62 = getelementptr inbounds float, float* %61, i64 0
store float %60, float* %62
%63 = load float, float* %12
%64 = fpext float %63 to double
%65 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str, i64 0, i64 0), double %64)
%66 = load float*, float** %11
%67 = bitcast float* %66 to i8*
ret void
}

LLVM IR (-O2):
@.str = [18 x i8] c"detrend(): b0=%f\0A\00"
define void @detrend(float* nocapture %0, i32 %1, float* %2, float* %3) {
%5 = alloca float
%6 = bitcast float* %5 to i8*
%7 = sdiv i32 %1, 4
%8 = sext i32 %1 to i64
%9 = shl nsw i64 %8, 2
%10 = tail call i32 (i64, ...) bitcast (i32 (...)* @check_malloc to i32 (i64, ...)*)(i64 %9)
%11 = sext i32 %10 to i64
%12 = inttoptr i64 %11 to float*
%13 = inttoptr i64 %11 to i8*
%14 = icmp sgt i32 %1, 0
br i1 %14, label %15, label %141
15: ; preds = %4
%16 = zext i32 %1 to i64
%17 = icmp ult i32 %1, 8
br i1 %17, label %104, label %18
18: ; preds = %15
%19 = getelementptr float, float* %12, i64 %16
%20 = getelementptr float, float* %0, i64 %16
%21 = bitcast float* %20 to i8*
%22 = icmp ult i8* %13, %21
%23 = icmp ugt float* %19, %0
%24 = and i1 %22, %23
br i1 %24, label %104, label %25
25: ; preds = %18
%26 = and i64 %16, 4294967288
%27 = add nsw i64 %26, -8
%28 = lshr exact i64 %27, 3
%29 = add nuw nsw i64 %28, 1
%30 = and i64 %29, 3
%31 = icmp ult i64 %27, 24
br i1 %31, label %83, label %32
32: ; preds = %25
%33 = and i64 %29, 4611686018427387900
br label %34
34: ; preds = %34, %32
%35 = phi i64 [ 0, %32 ], [ %80, %34 ]
%36 = phi i64 [ 0, %32 ], [ %81, %34 ]
%37 = getelementptr inbounds float, float* %0, i64 %35
%38 = bitcast float* %37 to <4 x float>*
%39 = load <4 x float>, <4 x float>* %38
%40 = getelementptr inbounds float, float* %37, i64 4
%41 = bitcast float* %40 to <4 x float>*
%42 = load <4 x float>, <4 x float>* %41
%43 = getelementptr inbounds float, float* %12, i64 %35
%44 = bitcast float* %43 to <4 x float>*
store <4 x float> %39, <4 x float>* %44
%45 = getelementptr inbounds float, float* %43, i64 4
%46 = bitcast float* %45 to <4 x float>*
store <4 x float> %42, <4 x float>* %46
%47 = or i64 %35, 8
%48 = getelementptr inbounds float, float* %0, i64 %47
%49 = bitcast float* %48 to <4 x float>*
%50 = load <4 x float>, <4 x float>* %49
%51 = getelementptr inbounds float, float* %48, i64 4
%52 = bitcast float* %51 to <4 x float>*
%53 = load <4 x float>, <4 x float>* %52
%54 = getelementptr inbounds float, float* %12, i64 %47
%55 = bitcast float* %54 to <4 x float>*
store <4 x float> %50, <4 x float>* %55
%56 = getelementptr inbounds float, float* %54, i64 4
%57 = bitcast float* %56 to <4 x float>*
store <4 x float> %53, <4 x float>* %57
%58 = or i64 %35, 16
%59 = getelementptr inbounds float, float* %0, i64 %58
%60 = bitcast float* %59 to <4 x float>*
%61 = load <4 x float>, <4 x float>* %60
%62 = getelementptr inbounds float, float* %59, i64 4
%63 = bitcast float* %62 to <4 x float>*
%64 = load <4 x float>, <4 x float>* %63
%65 = getelementptr inbounds float, float* %12, i64 %58
%66 = bitcast float* %65 to <4 x float>*
store <4 x float> %61, <4 x float>* %66
%67 = getelementptr inbounds float, float* %65, i64 4
%68 = bitcast float* %67 to <4 x float>*
store <4 x float> %64, <4 x float>* %68
%69 = or i64 %35, 24
%70 = getelementptr inbounds float, float* %0, i64 %69
%71 = bitcast float* %70 to <4 x float>*
%72 = load <4 x float>, <4 x float>* %71
%73 = getelementptr inbounds float, float* %70, i64 4
%74 = bitcast float* %73 to <4 x float>*
%75 = load <4 x float>, <4 x float>* %74
%76 = getelementptr inbounds float, float* %12, i64 %69
%77 = bitcast float* %76 to <4 x float>*
store <4 x float> %72, <4 x float>* %77
%78 = getelementptr inbounds float, float* %76, i64 4
%79 = bitcast float* %78 to <4 x float>*
store <4 x float> %75, <4 x float>* %79
%80 = add nuw i64 %35, 32
%81 = add i64 %36, 4
%82 = icmp eq i64 %81, %33
br i1 %82, label %83, label %34
83: ; preds = %34, %25
%84 = phi i64 [ 0, %25 ], [ %80, %34 ]
%85 = icmp eq i64 %30, 0
br i1 %85, label %102, label %86
86: ; preds = %83, %86
%87 = phi i64 [ %99, %86 ], [ %84, %83 ]
%88 = phi i64 [ %100, %86 ], [ 0, %83 ]
%89 = getelementptr inbounds float, float* %0, i64 %87
%90 = bitcast float* %89 to <4 x float>*
%91 = load <4 x float>, <4 x float>* %90
%92 = getelementptr inbounds float, float* %89, i64 4
%93 = bitcast float* %92 to <4 x float>*
%94 = load <4 x float>, <4 x float>* %93
%95 = getelementptr inbounds float, float* %12, i64 %87
%96 = bitcast float* %95 to <4 x float>*
store <4 x float> %91, <4 x float>* %96
%97 = getelementptr inbounds float, float* %95, i64 4
%98 = bitcast float* %97 to <4 x float>*
store <4 x float> %94, <4 x float>* %98
%99 = add nuw i64 %87, 8
%100 = add i64 %88, 1
%101 = icmp eq i64 %100, %30
br i1 %101, label %102, label %86
102: ; preds = %86, %83
%103 = icmp eq i64 %26, %16
br i1 %103, label %141, label %104
104: ; preds = %18, %15, %102
%105 = phi i64 [ 0, %18 ], [ 0, %15 ], [ %26, %102 ]
%106 = xor i64 %105, -1
%107 = add nsw i64 %106, %16
%108 = and i64 %16, 3
%109 = icmp eq i64 %108, 0
br i1 %109, label %119, label %110
110: ; preds = %104, %110
%111 = phi i64 [ %116, %110 ], [ %105, %104 ]
%112 = phi i64 [ %117, %110 ], [ 0, %104 ]
%113 = getelementptr inbounds float, float* %0, i64 %111
%114 = load float, float* %113
%115 = getelementptr inbounds float, float* %12, i64 %111
store float %114, float* %115
%116 = add nuw nsw i64 %111, 1
%117 = add i64 %112, 1
%118 = icmp eq i64 %117, %108
br i1 %118, label %119, label %110
119: ; preds = %110, %104
%120 = phi i64 [ %105, %104 ], [ %116, %110 ]
%121 = icmp ult i64 %107, 3
br i1 %121, label %141, label %122
122: ; preds = %119, %122
%123 = phi i64 [ %139, %122 ], [ %120, %119 ]
%124 = getelementptr inbounds float, float* %0, i64 %123
%125 = load float, float* %124
%126 = getelementptr inbounds float, float* %12, i64 %123
store float %125, float* %126
%127 = add nuw nsw i64 %123, 1
%128 = getelementptr inbounds float, float* %0, i64 %127
%129 = load float, float* %128
%130 = getelementptr inbounds float, float* %12, i64 %127
store float %129, float* %130
%131 = add nuw nsw i64 %123, 2
%132 = getelementptr inbounds float, float* %0, i64 %131
%133 = load float, float* %132
%134 = getelementptr inbounds float, float* %12, i64 %131
store float %133, float* %134
%135 = add nuw nsw i64 %123, 3
%136 = getelementptr inbounds float, float* %0, i64 %135
%137 = load float, float* %136
%138 = getelementptr inbounds float, float* %12, i64 %135
store float %137, float* %138
%139 = add nuw nsw i64 %123, 4
%140 = icmp eq i64 %139, %16
br i1 %140, label %141, label %122
141: ; preds = %119, %122, %102, %4
%142 = tail call i32 (float*, i32, float*, float*, ...) bitcast (i32 (...)* @integrate to i32 (float*, i32, float*, float*, ...)*)(float* %12, i32 %1, float* %2, float* %3)
%143 = sub nsw i32 %1, %7
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds float, float* %12, i64 %144
%146 = call i32 (float*, i32, float*, float*, ...) bitcast (i32 (...)* @get_trend to i32 (float*, i32, float*, float*, ...)*)(float* %145, i32 %7, float* %2, float* %5)
%147 = load float, float* %0
%148 = load float, float* %5
%149 = load float, float* %2
%150 = fdiv float %148, %149
%151 = fsub float %147, %150
store float %151, float* %0
%152 = fpext float %148 to double
%153 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str, i64 0, i64 0), double %152)
%154 = inttoptr i64 %11 to i8*
ret void
}