void CSPL_Array_clip(double *inval, const double min, const double max,
                     const long n) {
  long i;
  for (i = 0; i < n; i++) {
    inval[i] = (inval[i] < min) ? min : inval[i];
    inval[i] = (inval[i] > max) ? max : inval[i];
  }
}


x86 assembly (clang, -O0):
<CSPL_Array_clip>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    11d7 <CSPL_Array_clip+0xd7>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1157 <CSPL_Array_clip+0x57>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    1169 <CSPL_Array_clip+0x69>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    11a2 <CSPL_Array_clip+0xa2>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    11b4 <CSPL_Array_clip+0xb4>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    111e <CSPL_Array_clip+0x1e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<CSPL_Array_clip>:
test   %rsi,%rsi
jle    122f <CSPL_Array_clip+0x12f>
cmp    $0x4,%rsi
jae    1116 <CSPL_Array_clip+0x16>
xor    %eax,%eax
jmp    1210 <CSPL_Array_clip+0x110>
mov    %rsi,%rax
and    $0xfffffffffffffffc,%rax
movapd %xmm0,%xmm8
unpcklpd %xmm0,%xmm8
movapd %xmm1,%xmm3
unpcklpd %xmm1,%xmm3
lea    -0x4(%rax),%rcx
mov    %rcx,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rcx,%rcx
je     1230 <CSPL_Array_clip+0x130>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %ecx,%ecx
movupd (%rdi,%rcx,8),%xmm4
movupd 0x10(%rdi,%rcx,8),%xmm5
movupd 0x20(%rdi,%rcx,8),%xmm6
movupd 0x30(%rdi,%rcx,8),%xmm7
movapd %xmm8,%xmm2
maxpd  %xmm4,%xmm2
movapd %xmm8,%xmm4
maxpd  %xmm5,%xmm4
movapd %xmm3,%xmm5
minpd  %xmm2,%xmm5
movapd %xmm3,%xmm2
minpd  %xmm4,%xmm2
movupd %xmm5,(%rdi,%rcx,8)
movupd %xmm2,0x10(%rdi,%rcx,8)
movapd %xmm8,%xmm2
maxpd  %xmm6,%xmm2
movapd %xmm8,%xmm4
maxpd  %xmm7,%xmm4
movapd %xmm3,%xmm5
minpd  %xmm2,%xmm5
movapd %xmm3,%xmm2
minpd  %xmm4,%xmm2
movupd %xmm5,0x20(%rdi,%rcx,8)
movupd %xmm2,0x30(%rdi,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdx
jne    1150 <CSPL_Array_clip+0x50>
test   $0x1,%r8b
je     1203 <CSPL_Array_clip+0x103>
movupd (%rdi,%rcx,8),%xmm2
movupd 0x10(%rdi,%rcx,8),%xmm4
movapd %xmm8,%xmm5
maxpd  %xmm2,%xmm5
maxpd  %xmm4,%xmm8
movapd %xmm3,%xmm2
minpd  %xmm5,%xmm2
minpd  %xmm8,%xmm3
movupd %xmm2,(%rdi,%rcx,8)
movupd %xmm3,0x10(%rdi,%rcx,8)
cmp    %rsi,%rax
je     122f <CSPL_Array_clip+0x12f>
nopl   0x0(%rax,%rax,1)
movapd %xmm0,%xmm2
maxsd  (%rdi,%rax,8),%xmm2
movapd %xmm1,%xmm3
minsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rax,8)
add    $0x1,%rax
cmp    %rax,%rsi
jne    1210 <CSPL_Array_clip+0x110>
ret
xor    %ecx,%ecx
test   $0x1,%r8b
jne    11d2 <CSPL_Array_clip+0xd2>
jmp    1203 <CSPL_Array_clip+0x103>

LLVM IR (-O0):
define void @CSPL_Array_clip(double* %0, double %1, double %2, i64 %3) {
%5 = alloca double*
%6 = alloca double
%7 = alloca double
%8 = alloca i64
%9 = alloca i64
store double* %0, double** %5
store double %1, double* %6
store double %2, double* %7
store i64 %3, i64* %8
store i64 0, i64* %9
br label %10
10: ; preds = %51, %4
%11 = load i64, i64* %9
%12 = load i64, i64* %8
%13 = icmp slt i64 %11, %12
br i1 %13, label %14, label %54
14: ; preds = %10
%15 = load double*, double** %5
%16 = load i64, i64* %9
%17 = getelementptr inbounds double, double* %15, i64 %16
%18 = load double, double* %17
%19 = load double, double* %6
%20 = fcmp olt double %18, %19
br i1 %20, label %21, label %23
21: ; preds = %14
%22 = load double, double* %6
br label %28
23: ; preds = %14
%24 = load double*, double** %5
%25 = load i64, i64* %9
%26 = getelementptr inbounds double, double* %24, i64 %25
%27 = load double, double* %26
br label %28
28: ; preds = %23, %21
%29 = phi double [ %22, %21 ], [ %27, %23 ]
%30 = load double*, double** %5
%31 = load i64, i64* %9
%32 = getelementptr inbounds double, double* %30, i64 %31
store double %29, double* %32
%33 = load double*, double** %5
%34 = load i64, i64* %9
%35 = getelementptr inbounds double, double* %33, i64 %34
%36 = load double, double* %35
%37 = load double, double* %7
%38 = fcmp ogt double %36, %37
br i1 %38, label %39, label %41
39: ; preds = %28
%40 = load double, double* %7
br label %46
41: ; preds = %28
%42 = load double*, double** %5
%43 = load i64, i64* %9
%44 = getelementptr inbounds double, double* %42, i64 %43
%45 = load double, double* %44
br label %46
46: ; preds = %41, %39
%47 = phi double [ %40, %39 ], [ %45, %41 ]
%48 = load double*, double** %5
%49 = load i64, i64* %9
%50 = getelementptr inbounds double, double* %48, i64 %49
store double %47, double* %50
br label %51
51: ; preds = %46
%52 = load i64, i64* %9
%53 = add nsw i64 %52, 1
store i64 %53, i64* %9
br label %10
54: ; preds = %10
ret void
}

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