void findpartshift(parts, explosion, partcenter, partshift) int parts;
double explosion;
double *partcenter;
double *partshift;
{
  int i;
  for (i = 0; i < parts; i++) {
    partshift[i << 1] =
        explosion * (partcenter[i << 1] - partcenter[parts << 1]);
    partshift[(i << 1) + 1] =
        explosion * (partcenter[(i << 1) + 1] - partcenter[(parts << 1) + 1]);
  }
}


x86 assembly (clang, -O0):
<findpartshift>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11bc <findpartshift+0xbc>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111b <findpartshift+0x1b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<findpartshift>:
push   %rbp
push   %r15
push   %r14
push   %rbx
test   %edi,%edi
jle    117e <findpartshift+0x7e>
mov    %rdx,%r10
lea    (%rdi,%rdi,1),%eax
movslq %eax,%r11
lea    (%rsi,%r11,8),%r8
lea    (%rdi,%rdi,1),%eax
add    $0x1,%eax
movslq %eax,%r9
mov    %edi,%ecx
cmp    $0x1,%edi
jne    1185 <findpartshift+0x85>
xor    %ebx,%ebx
lea    (%rbx,%rbx,1),%eax
add    $0x1,%eax
sub    %rbx,%rcx
shl    $0x21,%rbx
movabs $0x200000000,%rdi
mov    %rbx,%rdx
sar    $0x1d,%rdx
movsd  (%rsi,%rdx,1),%xmm1
subsd  (%r8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r10,%rdx,1)
cltq
movsd  (%rsi,%rax,8),%xmm1
subsd  (%rsi,%r9,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r10,%rax,8)
add    $0x2,%eax
add    %rdi,%rbx
add    $0xffffffffffffffff,%rcx
jne    1140 <findpartshift+0x40>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
lea    -0x1(%rcx),%rdi
mov    %edi,%ebp
add    %edi,%ebp
setb   %r15b
mov    $0x10,%edx
mov    %rdi,%rax
mul    %rdx
seto   %r14b
xor    %ebx,%ebx
test   %ebp,%ebp
js     1129 <findpartshift+0x29>
test   %r15b,%r15b
jne    1129 <findpartshift+0x29>
shr    $0x20,%rdi
jne    1129 <findpartshift+0x29>
lea    (%r10,%rax,1),%rdi
cmp    %r10,%rdi
jb     1129 <findpartshift+0x29>
test   %r14b,%r14b
jne    1129 <findpartshift+0x29>
lea    0x8(%r10),%rdi
add    %r10,%rax
add    $0x8,%rax
cmp    %rdi,%rax
jb     1129 <findpartshift+0x29>
test   %r14b,%r14b
jne    1129 <findpartshift+0x29>
lea    (%rsi,%r11,8),%rax
add    $0x10,%rax
mov    %rcx,%rdx
shl    $0x4,%rdx
lea    (%r10,%rdx,1),%rdi
add    %rsi,%rdx
cmp    %rdi,%r8
setb   %r11b
cmp    %r10,%rax
seta   %bpl
cmp    %rsi,%rdi
seta   %al
cmp    %r10,%rdx
seta   %dil
xor    %ebx,%ebx
test   %bpl,%r11b
jne    1129 <findpartshift+0x29>
and    %dil,%al
jne    1129 <findpartshift+0x29>
mov    %ecx,%ebx
and    $0xfffffffe,%ebx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    $0x1,%eax
xor    %ebp,%ebp
movabs $0x400000000,%r11
mov    %rbx,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %rbp,%rdi
sar    $0x1d,%rdi
movupd (%rsi,%rdi,1),%xmm2
movupd 0x10(%rsi,%rdi,1),%xmm3
movapd %xmm2,%xmm4
unpcklpd %xmm3,%xmm4
unpckhpd %xmm3,%xmm2
movsd  (%r8),%xmm3
unpcklpd %xmm3,%xmm3
subpd  %xmm3,%xmm4
mulpd  %xmm1,%xmm4
cltq
movsd  (%rsi,%r9,8),%xmm3
unpcklpd %xmm3,%xmm3
subpd  %xmm3,%xmm2
mulpd  %xmm1,%xmm2
movapd %xmm4,%xmm3
unpcklpd %xmm2,%xmm3
movhlps %xmm4,%xmm2
movups %xmm2,0x8(%r10,%rax,8)
movupd %xmm3,-0x8(%r10,%rax,8)
add    %r11,%rbp
add    $0x4,%eax
add    $0xfffffffffffffffe,%rdx
jne    1260 <findpartshift+0x160>
cmp    %rcx,%rbx
jne    1129 <findpartshift+0x29>
jmp    117e <findpartshift+0x7e>

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

LLVM IR (-O2):
define void @findpartshift(i32 %0, double %1, double* nocapture %2, double* nocapture writeonly %3) {
%5 = bitcast double* %3 to i8*
%6 = icmp sgt i32 %0, 0
br i1 %6, label %7, label %114
7: ; preds = %4
%8 = shl nuw i32 %0, 1
%9 = sext i32 %8 to i64
%10 = getelementptr double, double* %2, i64 %9
%11 = or i32 %8, 1
%12 = sext i32 %11 to i64
%13 = getelementptr inbounds double, double* %2, i64 %12
%14 = zext i32 %0 to i64
%15 = icmp eq i32 %0, 1
br i1 %15, label %90, label %16
16: ; preds = %7
%17 = add nsw i64 %14, -1
%18 = trunc i64 %17 to i32
%19 = call { i32, i1 } @llvm.umul.with.overflow.i32(i32 %18, i32 2)
%20 = extractvalue { i32, i1 } %19, 0
%21 = extractvalue { i32, i1 } %19, 1
%22 = icmp slt i32 %20, 0
%23 = or i1 %22, %21
%24 = icmp ugt i64 %17, 4294967295
%25 = or i1 %23, %24
%26 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %17, i64 16)
%27 = extractvalue { i64, i1 } %26, 0
%28 = extractvalue { i64, i1 } %26, 1
%29 = getelementptr i8, i8* %5, i64 %27
%30 = icmp ult i8* %29, %5
%31 = or i1 %30, %28
%32 = getelementptr double, double* %3, i64 1
%33 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %17, i64 16)
%34 = extractvalue { i64, i1 } %33, 0
%35 = extractvalue { i64, i1 } %33, 1
%36 = bitcast double* %32 to i8*
%37 = getelementptr i8, i8* %36, i64 %34
%38 = icmp ult i8* %37, %36
%39 = or i1 %38, %35
%40 = or i1 %25, %31
%41 = or i1 %40, %39
br i1 %41, label %90, label %42
42: ; preds = %16
%43 = add nsw i64 %9, 2
%44 = getelementptr double, double* %2, i64 %43
%45 = shl nuw nsw i64 %14, 1
%46 = getelementptr double, double* %3, i64 %45
%47 = getelementptr double, double* %2, i64 %45
%48 = icmp ult double* %10, %46
%49 = icmp ugt double* %44, %3
%50 = and i1 %48, %49
%51 = icmp ugt double* %46, %2
%52 = icmp ugt double* %47, %3
%53 = and i1 %51, %52
%54 = or i1 %50, %53
br i1 %54, label %90, label %55
55: ; preds = %42
%56 = and i64 %14, 4294967294
%57 = insertelement <2 x double> poison, double %1, i64 0
%58 = shufflevector <2 x double> %57, <2 x double> poison, <2 x i32> zeroinitializer
%59 = getelementptr inbounds double, double* %3, i64 -1
br label %60
60: ; preds = %60, %55
%61 = phi i64 [ 0, %55 ], [ %86, %60 ]
%62 = trunc i64 %61 to i32
%63 = shl i64 %61, 33
%64 = ashr exact i64 %63, 32
%65 = getelementptr inbounds double, double* %2, i64 %64
%66 = bitcast double* %65 to <4 x double>*
%67 = load <4 x double>, <4 x double>* %66
%68 = shufflevector <4 x double> %67, <4 x double> poison, <2 x i32> <i32 0, i32 2>
%69 = shufflevector <4 x double> %67, <4 x double> poison, <2 x i32> <i32 1, i32 3>
%70 = load double, double* %10
%71 = insertelement <2 x double> poison, double %70, i64 0
%72 = shufflevector <2 x double> %71, <2 x double> poison, <2 x i32> zeroinitializer
%73 = fsub <2 x double> %68, %72
%74 = fmul <2 x double> %73, %58
%75 = shl i32 %62, 1
%76 = or i32 %75, 1
%77 = sext i32 %76 to i64
%78 = load double, double* %13
%79 = insertelement <2 x double> poison, double %78, i64 0
%80 = shufflevector <2 x double> %79, <2 x double> poison, <2 x i32> zeroinitializer
%81 = fsub <2 x double> %69, %80
%82 = fmul <2 x double> %81, %58
%83 = getelementptr inbounds double, double* %59, i64 %77
%84 = bitcast double* %83 to <4 x double>*
%85 = shufflevector <2 x double> %74, <2 x double> %82, <4 x i32> <i32 0, i32 2, i32 1, i32 3>
store <4 x double> %85, <4 x double>* %84
%86 = add nuw i64 %61, 2
%87 = icmp eq i64 %86, %56
br i1 %87, label %88, label %60
88: ; preds = %60
%89 = icmp eq i64 %56, %14
br i1 %89, label %114, label %90
90: ; preds = %42, %16, %7, %88
%91 = phi i64 [ 0, %42 ], [ 0, %16 ], [ 0, %7 ], [ %56, %88 ]
br label %92
92: ; preds = %90, %92
%93 = phi i64 [ %112, %92 ], [ %91, %90 ]
%94 = shl i64 %93, 33
%95 = ashr exact i64 %94, 32
%96 = getelementptr inbounds double, double* %2, i64 %95
%97 = load double, double* %96
%98 = load double, double* %10
%99 = fsub double %97, %98
%100 = fmul double %99, %1
%101 = getelementptr inbounds double, double* %3, i64 %95
store double %100, double* %101
%102 = trunc i64 %93 to i32
%103 = shl i32 %102, 1
%104 = or i32 %103, 1
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds double, double* %2, i64 %105
%107 = load double, double* %106
%108 = load double, double* %13
%109 = fsub double %107, %108
%110 = fmul double %109, %1
%111 = getelementptr inbounds double, double* %3, i64 %105
store double %110, double* %111
%112 = add nuw nsw i64 %93, 1
%113 = icmp eq i64 %112, %14
br i1 %113, label %114, label %92
114: ; preds = %92, %88, %4
ret void
}