int fp_precision_reduction(double *in, unsigned int len,
                           unsigned int precision) {
  unsigned int i, b;
  for (i = 0; i < len; i++) {
    for (b = 0; b < (52 - precision); b++) {
      *(unsigned long *)&in[i] &= ~((unsigned long)1 << b);
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<fp_precision_reduction>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    117c <fp_precision_reduction+0x7c>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x34,%ecx
sub    -0x10(%rbp),%ecx
cmp    %ecx,%eax
jae    1169 <fp_precision_reduction+0x69>
mov    -0x18(%rbp),%eax
mov    %eax,%ecx
mov    $0x1,%edx
shl    %cl,%rdx
xor    $0xffffffffffffffff,%rdx
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
and    (%rax,%rcx,8),%rdx
mov    %rdx,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1128 <fp_precision_reduction+0x28>
jmp    116e <fp_precision_reduction+0x6e>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1115 <fp_precision_reduction+0x15>
mov    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<fp_precision_reduction>:
push   %r14
push   %rbx
test   %esi,%esi
je     1337 <fp_precision_reduction+0x237>
mov    $0x34,%ecx
sub    %edx,%ecx
cmp    $0x2,%ecx
mov    $0x1,%eax
cmovae %ecx,%eax
mov    %esi,%r14d
mov    %eax,%r11d
and    $0xfffffffc,%r11d
lea    -0x4(%r11),%r8
mov    %r8,%r9
shr    $0x2,%r9
add    $0x1,%r9
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %esi,%esi
movapd 0xeb8(%rip),%xmm8
movdqa 0xebf(%rip),%xmm9
movdqa 0xec6(%rip),%xmm13
movdqa 0xecd(%rip),%xmm14
movdqa 0xed4(%rip),%xmm10
movdqa 0xedb(%rip),%xmm11
movdqa 0xee2(%rip),%xmm12
jmp    1191 <fp_precision_reduction+0x91>
mov    %rbx,(%rdi,%rsi,8)
add    $0x1,%rsi
cmp    %r14,%rsi
je     1337 <fp_precision_reduction+0x237>
cmp    $0x34,%edx
je     1184 <fp_precision_reduction+0x84>
mov    (%rdi,%rsi,8),%rbx
cmp    $0x4,%eax
jae    11b0 <fp_precision_reduction+0xb0>
xor    %ecx,%ecx
jmp    1310 <fp_precision_reduction+0x210>
cs nopw 0x0(%rax,%rax,1)
movq   %rbx,%xmm3
shufpd $0x2,%xmm8,%xmm3
test   %r8,%r8
je     1322 <fp_precision_reduction+0x222>
pcmpeqd %xmm15,%xmm15
mov    %r10,%rbx
movdqa %xmm9,%xmm2
movapd %xmm3,%xmm5
cs nopw 0x0(%rax,%rax,1)
nop
movdqa %xmm14,%xmm0
psllq  %xmm2,%xmm0
movdqa %xmm2,%xmm3
pshufd $0xee,%xmm2,%xmm6
movdqa %xmm14,%xmm1
psllq  %xmm6,%xmm1
paddq  %xmm13,%xmm3
movsd  %xmm0,%xmm1
movdqa %xmm14,%xmm0
psllq  %xmm3,%xmm0
pshufd $0xee,%xmm3,%xmm3
movdqa %xmm14,%xmm6
psllq  %xmm3,%xmm6
movsd  %xmm0,%xmm6
movdqa %xmm2,%xmm0
paddq  %xmm10,%xmm0
movdqa %xmm2,%xmm7
paddq  %xmm11,%xmm7
movdqa %xmm14,%xmm4
psllq  %xmm0,%xmm4
pshufd $0xee,%xmm0,%xmm0
movdqa %xmm14,%xmm3
psllq  %xmm0,%xmm3
movsd  %xmm4,%xmm3
orpd   %xmm1,%xmm3
andnpd %xmm5,%xmm3
movdqa %xmm14,%xmm1
psllq  %xmm7,%xmm1
pshufd $0xee,%xmm7,%xmm4
movdqa %xmm14,%xmm0
psllq  %xmm4,%xmm0
movsd  %xmm1,%xmm0
orpd   %xmm6,%xmm0
andnpd %xmm15,%xmm0
paddq  %xmm12,%xmm2
movapd %xmm3,%xmm5
movapd %xmm0,%xmm15
add    $0xfffffffffffffffe,%rbx
jne    11e0 <fp_precision_reduction+0xe0>
test   $0x1,%r9b
je     12e5 <fp_precision_reduction+0x1e5>
movdqa %xmm2,%xmm1
paddq  %xmm13,%xmm1
movdqa %xmm14,%xmm4
psllq  %xmm2,%xmm4
pshufd $0xee,%xmm2,%xmm2
movdqa %xmm14,%xmm5
psllq  %xmm2,%xmm5
movsd  %xmm4,%xmm5
andnpd %xmm3,%xmm5
movdqa %xmm14,%xmm2
psllq  %xmm1,%xmm2
pshufd $0xee,%xmm1,%xmm1
movdqa %xmm14,%xmm4
psllq  %xmm1,%xmm4
movsd  %xmm2,%xmm4
andnpd %xmm0,%xmm4
movapd %xmm5,%xmm3
movapd %xmm4,%xmm0
andpd  %xmm0,%xmm3
pshufd $0xee,%xmm3,%xmm0
pand   %xmm3,%xmm0
movq   %xmm0,%rbx
mov    %r11,%rcx
cmp    %rax,%r11
je     1180 <fp_precision_reduction+0x80>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
btr    %rcx,%rbx
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1310 <fp_precision_reduction+0x210>
jmp    1180 <fp_precision_reduction+0x80>
pcmpeqd %xmm0,%xmm0
movdqa %xmm9,%xmm2
test   $0x1,%r9b
jne    1296 <fp_precision_reduction+0x196>
jmp    12e5 <fp_precision_reduction+0x1e5>
mov    $0x1,%eax
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define i32 @fp_precision_reduction(double* %0, i32 %1, i32 %2) {
%4 = alloca double*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
store double* %0, double** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
store i32 0, i32* %7
br label %9
9: ; preds = %35, %3
%10 = load i32, i32* %7
%11 = load i32, i32* %5
%12 = icmp ult i32 %10, %11
br i1 %12, label %13, label %38
13: ; preds = %9
store i32 0, i32* %8
br label %14
14: ; preds = %31, %13
%15 = load i32, i32* %8
%16 = load i32, i32* %6
%17 = sub i32 52, %16
%18 = icmp ult i32 %15, %17
br i1 %18, label %19, label %34
19: ; preds = %14
%20 = load i32, i32* %8
%21 = zext i32 %20 to i64
%22 = shl i64 1, %21
%23 = xor i64 %22, -1
%24 = load double*, double** %4
%25 = load i32, i32* %7
%26 = zext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
%28 = bitcast double* %27 to i64*
%29 = load i64, i64* %28
%30 = and i64 %29, %23
store i64 %30, i64* %28
br label %31
31: ; preds = %19
%32 = load i32, i32* %8
%33 = add i32 %32, 1
store i32 %33, i32* %8
br label %14
34: ; preds = %14
br label %35
35: ; preds = %34
%36 = load i32, i32* %7
%37 = add i32 %36, 1
store i32 %37, i32* %7
br label %9
38: ; preds = %9
ret i32 1
}

LLVM IR (-O2):
define i32 @fp_precision_reduction(double* nocapture %0, i32 %1, i32 %2) {
%4 = icmp eq i32 %1, 0
br i1 %4, label %86, label %5
5: ; preds = %3
%6 = sub i32 52, %2
%7 = icmp eq i32 %2, 52
%8 = call i32 @llvm.umax.i32(i32 %6, i32 1)
%9 = zext i32 %1 to i64
%10 = zext i32 %8 to i64
%11 = and i64 %10, 4294967292
%12 = add nsw i64 %11, -4
%13 = lshr exact i64 %12, 2
%14 = add nuw nsw i64 %13, 1
%15 = icmp ult i32 %8, 4
%16 = and i64 %10, 4294967292
%17 = and i64 %14, 1
%18 = icmp eq i64 %12, 0
%19 = and i64 %14, 9223372036854775806
%20 = icmp eq i64 %17, 0
%21 = icmp eq i64 %16, %10
br label %22
22: ; preds = %5, %83
%23 = phi i64 [ 0, %5 ], [ %84, %83 ]
br i1 %7, label %83, label %24
24: ; preds = %22
%25 = getelementptr inbounds double, double* %0, i64 %23
%26 = bitcast double* %25 to i64*
%27 = load i64, i64* %26
br i1 %15, label %70, label %28
28: ; preds = %24
%29 = insertelement <2 x i64> <i64 poison, i64 -1>, i64 %27, i64 0
br i1 %18, label %51, label %30
30: ; preds = %28, %30
%31 = phi <2 x i64> [ %48, %30 ], [ <i64 0, i64 1>, %28 ]
%32 = phi <2 x i64> [ %44, %30 ], [ %29, %28 ]
%33 = phi <2 x i64> [ %47, %30 ], [ <i64 -1, i64 -1>, %28 ]
%34 = phi i64 [ %49, %30 ], [ 0, %28 ]
%35 = add <2 x i64> %31, <i64 2, i64 2>
%36 = shl nuw <2 x i64> <i64 1, i64 1>, %31
%37 = shl nuw <2 x i64> <i64 1, i64 1>, %35
%38 = add <2 x i64> %31, <i64 4, i64 4>
%39 = add <2 x i64> %31, <i64 6, i64 6>
%40 = shl nuw <2 x i64> <i64 1, i64 1>, %38
%41 = shl nuw <2 x i64> <i64 1, i64 1>, %39
%42 = or <2 x i64> %36, %40
%43 = xor <2 x i64> %42, <i64 -1, i64 -1>
%44 = and <2 x i64> %32, %43
%45 = or <2 x i64> %37, %41
%46 = xor <2 x i64> %45, <i64 -1, i64 -1>
%47 = and <2 x i64> %33, %46
%48 = add <2 x i64> %31, <i64 8, i64 8>
%49 = add i64 %34, 2
%50 = icmp eq i64 %49, %19
br i1 %50, label %51, label %30
51: ; preds = %30, %28
%52 = phi <2 x i64> [ undef, %28 ], [ %44, %30 ]
%53 = phi <2 x i64> [ undef, %28 ], [ %47, %30 ]
%54 = phi <2 x i64> [ <i64 0, i64 1>, %28 ], [ %48, %30 ]
%55 = phi <2 x i64> [ %29, %28 ], [ %44, %30 ]
%56 = phi <2 x i64> [ <i64 -1, i64 -1>, %28 ], [ %47, %30 ]
br i1 %20, label %65, label %57
57: ; preds = %51
%58 = add <2 x i64> %54, <i64 2, i64 2>
%59 = shl nuw <2 x i64> <i64 1, i64 1>, %54
%60 = shl nuw <2 x i64> <i64 1, i64 1>, %58
%61 = xor <2 x i64> %59, <i64 -1, i64 -1>
%62 = xor <2 x i64> %60, <i64 -1, i64 -1>
%63 = and <2 x i64> %55, %61
%64 = and <2 x i64> %56, %62
br label %65
65: ; preds = %51, %57
%66 = phi <2 x i64> [ %52, %51 ], [ %63, %57 ]
%67 = phi <2 x i64> [ %53, %51 ], [ %64, %57 ]
%68 = and <2 x i64> %67, %66
%69 = call i64 @llvm.vector.reduce.and.v2i64(<2 x i64> %68)
br i1 %21, label %81, label %70
70: ; preds = %24, %65
%71 = phi i64 [ 0, %24 ], [ %16, %65 ]
%72 = phi i64 [ %27, %24 ], [ %69, %65 ]
br label %73
73: ; preds = %70, %73
%74 = phi i64 [ %79, %73 ], [ %71, %70 ]
%75 = phi i64 [ %78, %73 ], [ %72, %70 ]
%76 = shl nuw i64 1, %74
%77 = xor i64 %76, -1
%78 = and i64 %75, %77
%79 = add nuw nsw i64 %74, 1
%80 = icmp eq i64 %79, %10
br i1 %80, label %81, label %73
81: ; preds = %73, %65
%82 = phi i64 [ %69, %65 ], [ %78, %73 ]
store i64 %82, i64* %26
br label %83
83: ; preds = %81, %22
%84 = add nuw nsw i64 %23, 1
%85 = icmp eq i64 %84, %9
br i1 %85, label %86, label %22
86: ; preds = %83, %3
ret i32 1
}