int psd_tonemasker(float *psd, int psd_len, int splnum1, int splnum2,
                   float *ptm, int *tone_flag) {
  int k;
  memset(ptm, 0, sizeof(float) * psd_len);
  memset(tone_flag, 0, sizeof(int) * psd_len);
  for (k = 0; k < psd_len; k++) {
    if (istone(psd, psd_len, k, splnum1, splnum2, tone_flag)) {
      ptm[k] = db_pos(db_inv(psd[k - 1]) + db_inv(psd[k]) + db_inv(psd[k + 1]));
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<psd_tonemasker>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x20(%rbp),%rdi
movslq -0xc(%rbp),%rdx
shl    $0x2,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
mov    -0x28(%rbp),%rdi
movslq -0xc(%rbp),%rdx
shl    $0x2,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    124a <psd_tonemasker+0x10a>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x2c(%rbp),%edx
mov    -0x10(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x28(%rbp),%r9
mov    $0x0,%al
call   1050 <istone@plt>
cmp    $0x0,%eax
je     1237 <psd_tonemasker+0xf7>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1060 <db_inv@plt>
mov    %eax,-0x34(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1060 <db_inv@plt>
mov    %eax,%ecx
mov    -0x34(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1060 <db_inv@plt>
mov    -0x30(%rbp),%edi
add    %eax,%edi
mov    $0x0,%al
call   1040 <db_pos@plt>
cvtsi2ss %eax,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    123c <psd_tonemasker+0xfc>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    118a <psd_tonemasker+0x4a>
xor    %eax,%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<psd_tonemasker>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r9,%r12
mov    %ecx,0xc(%rsp)
mov    %edx,0x8(%rsp)
mov    %rdi,%rbx
mov    %esi,0x4(%rsp)
movslq %esi,%r15
lea    0x0(,%r15,4),%r14
mov    %r8,0x10(%rsp)
mov    %r8,%rdi
xor    %esi,%esi
mov    %r14,%rdx
call   1030 <memset@plt>
mov    %r12,%rbp
mov    %r12,%rdi
xor    %esi,%esi
mov    %r14,%rdx
call   1030 <memset@plt>
test   %r15d,%r15d
jle    1234 <psd_tonemasker+0xf4>
mov    0x4(%rsp),%r12d
xor    %r14d,%r14d
jmp    1206 <psd_tonemasker+0xc6>
movss  -0x4(%rbx,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1060 <db_inv@plt>
mov    %eax,%r15d
movss  (%rbx,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1060 <db_inv@plt>
add    %eax,%r15d
lea    0x1(%r14),%r13
movss  0x4(%rbx,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1060 <db_inv@plt>
add    %eax,%r15d
mov    %r15d,%edi
xor    %eax,%eax
call   1040 <db_pos@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
mov    0x10(%rsp),%rax
movss  %xmm0,(%rax,%r14,4)
mov    %r13,%r14
cmp    %r12,%r14
je     1234 <psd_tonemasker+0xf4>
mov    %rbx,%rdi
mov    0x4(%rsp),%esi
mov    %r14d,%edx
mov    0x8(%rsp),%ecx
mov    0xc(%rsp),%r8d
mov    %rbp,%r9
xor    %eax,%eax
call   1050 <istone@plt>
test   %eax,%eax
jne    11a0 <psd_tonemasker+0x60>
add    $0x1,%r14
cmp    %r12,%r14
jne    1206 <psd_tonemasker+0xc6>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @psd_tonemasker(float* %0, i32 %1, i32 %2, i32 %3, float* %4, i32* %5) {
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca float*
%12 = alloca i32*
%13 = alloca i32
store float* %0, float** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store float* %4, float** %11
store i32* %5, i32** %12
%14 = load float*, float** %11
%15 = bitcast float* %14 to i8*
%16 = load i32, i32* %8
%17 = sext i32 %16 to i64
%18 = mul i64 4, %17
%19 = load i32*, i32** %12
%20 = bitcast i32* %19 to i8*
%21 = load i32, i32* %8
%22 = sext i32 %21 to i64
%23 = mul i64 4, %22
store i32 0, i32* %13
br label %24
24: ; preds = %70, %6
%25 = load i32, i32* %13
%26 = load i32, i32* %8
%27 = icmp slt i32 %25, %26
br i1 %27, label %28, label %73
28: ; preds = %24
%29 = load float*, float** %7
%30 = load i32, i32* %8
%31 = load i32, i32* %13
%32 = load i32, i32* %9
%33 = load i32, i32* %10
%34 = load i32*, i32** %12
%35 = call i32 (float*, i32, i32, i32, i32, i32*, ...) bitcast (i32 (...)* @istone to i32 (float*, i32, i32, i32, i32, i32*, ...)*)(float* %29, i32 %30, i32 %31, i32 %32, i32 %33, i32* %34)
%36 = icmp ne i32 %35, 0
br i1 %36, label %37, label %69
37: ; preds = %28
%38 = load float*, float** %7
%39 = load i32, i32* %13
%40 = sub nsw i32 %39, 1
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %38, i64 %41
%43 = load float, float* %42
%44 = fpext float %43 to double
%45 = call i32 (double, ...) bitcast (i32 (...)* @db_inv to i32 (double, ...)*)(double %44)
%46 = load float*, float** %7
%47 = load i32, i32* %13
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %46, i64 %48
%50 = load float, float* %49
%51 = fpext float %50 to double
%52 = call i32 (double, ...) bitcast (i32 (...)* @db_inv to i32 (double, ...)*)(double %51)
%53 = add nsw i32 %45, %52
%54 = load float*, float** %7
%55 = load i32, i32* %13
%56 = add nsw i32 %55, 1
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds float, float* %54, i64 %57
%59 = load float, float* %58
%60 = fpext float %59 to double
%61 = call i32 (double, ...) bitcast (i32 (...)* @db_inv to i32 (double, ...)*)(double %60)
%62 = add nsw i32 %53, %61
%63 = call i32 (i32, ...) bitcast (i32 (...)* @db_pos to i32 (i32, ...)*)(i32 %62)
%64 = sitofp i32 %63 to float
%65 = load float*, float** %11
%66 = load i32, i32* %13
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds float, float* %65, i64 %67
store float %64, float* %68
br label %69
69: ; preds = %37, %28
br label %70
70: ; preds = %69
%71 = load i32, i32* %13
%72 = add nsw i32 %71, 1
store i32 %72, i32* %13
br label %24
73: ; preds = %24
ret i32 0
}

LLVM IR (-O2):
define i32 @psd_tonemasker(float* %0, i32 %1, i32 %2, i32 %3, float* nocapture writeonly %4, i32* %5) {
%7 = bitcast float* %4 to i8*
%8 = sext i32 %1 to i64
%9 = shl nsw i64 %8, 2
tail call void @llvm.memset.p0i8.i64(i8* align 4 %7, i8 0, i64 %9, i1 false)
%10 = bitcast i32* %5 to i8*
tail call void @llvm.memset.p0i8.i64(i8* align 4 %10, i8 0, i64 %9, i1 false)
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %44
12: ; preds = %6
%13 = zext i32 %1 to i64
br label %14
14: ; preds = %12, %41
%15 = phi i64 [ 0, %12 ], [ %42, %41 ]
%16 = trunc i64 %15 to i32
%17 = tail call i32 (float*, i32, i32, i32, i32, i32*, ...) bitcast (i32 (...)* @istone to i32 (float*, i32, i32, i32, i32, i32*, ...)*)(float* %0, i32 %1, i32 %16, i32 %2, i32 %3, i32* %5)
%18 = icmp eq i32 %17, 0
br i1 %18, label %19, label %21
19: ; preds = %14
%20 = add nuw nsw i64 %15, 1
br label %41
21: ; preds = %14
%22 = add nsw i64 %15, -1
%23 = getelementptr inbounds float, float* %0, i64 %22
%24 = load float, float* %23
%25 = fpext float %24 to double
%26 = tail call i32 (double, ...) bitcast (i32 (...)* @db_inv to i32 (double, ...)*)(double %25)
%27 = getelementptr inbounds float, float* %0, i64 %15
%28 = load float, float* %27
%29 = fpext float %28 to double
%30 = tail call i32 (double, ...) bitcast (i32 (...)* @db_inv to i32 (double, ...)*)(double %29)
%31 = add nsw i32 %30, %26
%32 = add nuw nsw i64 %15, 1
%33 = getelementptr inbounds float, float* %0, i64 %32
%34 = load float, float* %33
%35 = fpext float %34 to double
%36 = tail call i32 (double, ...) bitcast (i32 (...)* @db_inv to i32 (double, ...)*)(double %35)
%37 = add nsw i32 %31, %36
%38 = tail call i32 (i32, ...) bitcast (i32 (...)* @db_pos to i32 (i32, ...)*)(i32 %37)
%39 = sitofp i32 %38 to float
%40 = getelementptr inbounds float, float* %4, i64 %15
store float %39, float* %40
br label %41
41: ; preds = %19, %21
%42 = phi i64 [ %20, %19 ], [ %32, %21 ]
%43 = icmp eq i64 %42, %13
br i1 %43, label %44, label %14
44: ; preds = %41, %6
ret i32 0
}