void psdes(unsigned long *lword, unsigned long *irword) {
  unsigned long i, ia, ib, iswap, itmph = 0, itmpl = 0;
  static unsigned long c1[4] = {0xbaa96887L, 0x1e17d32cL, 0x03bcdc3cL,
                                0x0f33d1b2L};
  static unsigned long c2[4] = {0x4b0f3b58L, 0xe874f0c3L, 0x6955c5a6L,
                                0x55a7ca46L};
  for (i = 0; i < 4; i++) {
    ia = (iswap = (*irword)) ^ c1[i];
    itmpl = ia & 0xffff;
    itmph = ia >> 16;
    ib = itmpl * itmpl + ~(itmph * itmph);
    *irword = (*lword) ^ (((ia = (ib >> 16) | ((ib & 0xffff) << 16)) ^ c2[i]) +
                          itmpl * itmph);
    *lword = iswap;
  }
}


x86 assembly (clang, -O0):
<psdes>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movq   $0x0,-0x38(%rbp)
movq   $0x0,-0x40(%rbp)
movq   $0x0,-0x18(%rbp)
cmpq   $0x4,-0x18(%rbp)
jae    11ea <psdes+0xea>
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rdx
lea    0x2eeb(%rip),%rcx
xor    (%rcx,%rdx,8),%rax
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rax
and    $0xffff,%rax
mov    %rax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
shr    $0x10,%rax
mov    %rax,-0x38(%rbp)
mov    -0x40(%rbp),%rax
imul   -0x40(%rbp),%rax
mov    -0x38(%rbp),%rcx
imul   -0x38(%rbp),%rcx
xor    $0xffffffffffffffff,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%rcx
mov    -0x28(%rbp),%rax
shr    $0x10,%rax
mov    -0x28(%rbp),%rdx
and    $0xffff,%rdx
shl    $0x10,%rdx
or     %rdx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x18(%rbp),%rsi
lea    0x2e9c(%rip),%rdx
xor    (%rdx,%rsi,8),%rax
mov    -0x40(%rbp),%rdx
imul   -0x38(%rbp),%rdx
add    %rdx,%rax
xor    %rax,%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x30(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x18(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x18(%rbp)
jmp    1124 <psdes+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<psdes>:
mov    (%rsi),%r8
mov    $0xbaa96887,%ecx
mov    %r8,%rdx
xor    %rcx,%rdx
movzwl %dx,%r9d
shr    $0x10,%rdx
mov    %r9,%rax
imul   %r9,%rax
imul   %rdx,%r9
imul   %rdx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rdx,%rax
shr    $0x10,%rax
shl    $0x10,%edx
or     %rax,%rdx
xor    $0x4b0f3b58,%rdx
add    %r9,%rdx
xor    (%rdi),%rdx
mov    %rdx,(%rsi)
mov    %r8,(%rdi)
mov    (%rsi),%r9
mov    %r9,%rax
xor    $0x1e17d32c,%rax
movzwl %ax,%r10d
shr    $0x10,%rax
mov    %r10,%rdx
imul   %r10,%rdx
imul   %rax,%r10
imul   %rax,%rax
not    %rax
add    %rdx,%rax
mov    %rax,%rdx
shr    $0x10,%rdx
shl    $0x10,%eax
or     %rdx,%rax
add    $0x2dcb883c,%rcx
xor    %rax,%rcx
add    %r10,%rcx
xor    %r8,%rcx
mov    %rcx,(%rsi)
mov    %r9,(%rdi)
mov    (%rsi),%r8
mov    %r8,%rcx
xor    $0x3bcdc3c,%rcx
movzwl %cx,%eax
shr    $0x10,%rcx
mov    %rax,%rdx
imul   %rax,%rdx
imul   %rcx,%rax
imul   %rcx,%rcx
not    %rcx
add    %rdx,%rcx
mov    %rcx,%rdx
shr    $0x10,%rdx
shl    $0x10,%ecx
or     %rdx,%rcx
xor    $0x6955c5a6,%rcx
add    %rax,%rcx
xor    %r9,%rcx
mov    %rcx,(%rsi)
mov    %r8,(%rdi)
mov    (%rsi),%r9
mov    %r9,%rdx
xor    $0xf33d1b2,%rdx
movzwl %dx,%ecx
shr    $0x10,%rdx
mov    %rcx,%rax
imul   %rcx,%rax
imul   %rdx,%rcx
imul   %rdx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rdx,%rax
shr    $0x10,%rax
shl    $0x10,%edx
or     %rax,%rdx
xor    $0x55a7ca46,%rdx
add    %rcx,%rdx
xor    %r8,%rdx
mov    %rdx,(%rsi)
mov    %r9,(%rdi)
ret

LLVM IR (-O0):
@psdes.c1 = internal global [4 x i64] [i64 3131664519, i64 504877868, i64 62708796, i64 255054258]
@psdes.c2 = internal global [4 x i64] [i64 1259289432, i64 3899977923, i64 1767228838, i64 1437059654]
define void @psdes(i64* %0, i64* %1) {
%3 = alloca i64*
%4 = alloca i64*
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
store i64* %0, i64** %3
store i64* %1, i64** %4
store i64 0, i64* %9
store i64 0, i64* %10
store i64 0, i64* %5
br label %11
11: ; preds = %53, %2
%12 = load i64, i64* %5
%13 = icmp ult i64 %12, 4
br i1 %13, label %14, label %56
14: ; preds = %11
%15 = load i64*, i64** %4
%16 = load i64, i64* %15
store i64 %16, i64* %8
%17 = load i64, i64* %5
%18 = getelementptr inbounds [4 x i64], [4 x i64]* @psdes.c1, i64 0, i64 %17
%19 = load i64, i64* %18
%20 = xor i64 %16, %19
store i64 %20, i64* %6
%21 = load i64, i64* %6
%22 = and i64 %21, 65535
store i64 %22, i64* %10
%23 = load i64, i64* %6
%24 = lshr i64 %23, 16
store i64 %24, i64* %9
%25 = load i64, i64* %10
%26 = load i64, i64* %10
%27 = mul i64 %25, %26
%28 = load i64, i64* %9
%29 = load i64, i64* %9
%30 = mul i64 %28, %29
%31 = xor i64 %30, -1
%32 = add i64 %27, %31
store i64 %32, i64* %7
%33 = load i64*, i64** %3
%34 = load i64, i64* %33
%35 = load i64, i64* %7
%36 = lshr i64 %35, 16
%37 = load i64, i64* %7
%38 = and i64 %37, 65535
%39 = shl i64 %38, 16
%40 = or i64 %36, %39
store i64 %40, i64* %6
%41 = load i64, i64* %5
%42 = getelementptr inbounds [4 x i64], [4 x i64]* @psdes.c2, i64 0, i64 %41
%43 = load i64, i64* %42
%44 = xor i64 %40, %43
%45 = load i64, i64* %10
%46 = load i64, i64* %9
%47 = mul i64 %45, %46
%48 = add i64 %44, %47
%49 = xor i64 %34, %48
%50 = load i64*, i64** %4
store i64 %49, i64* %50
%51 = load i64, i64* %8
%52 = load i64*, i64** %3
store i64 %51, i64* %52
br label %53
53: ; preds = %14
%54 = load i64, i64* %5
%55 = add i64 %54, 1
store i64 %55, i64* %5
br label %11
56: ; preds = %11
ret void
}

LLVM IR (-O2):
define void @psdes(i64* nocapture %0, i64* nocapture %1) {
%3 = load i64, i64* %1
%4 = xor i64 %3, 3131664519
%5 = and i64 %4, 65535
%6 = lshr i64 %4, 16
%7 = mul nuw nsw i64 %5, %5
%8 = mul i64 %6, %6
%9 = xor i64 %8, -1
%10 = add i64 %7, %9
%11 = load i64, i64* %0
%12 = lshr i64 %10, 16
%13 = shl i64 %10, 16
%14 = and i64 %13, 4294901760
%15 = or i64 %14, %12
%16 = xor i64 %15, 1259289432
%17 = mul nuw i64 %5, %6
%18 = add i64 %16, %17
%19 = xor i64 %18, %11
store i64 %19, i64* %1
store i64 %3, i64* %0
%20 = load i64, i64* %1
%21 = xor i64 %20, 504877868
%22 = and i64 %21, 65535
%23 = lshr i64 %21, 16
%24 = mul nuw nsw i64 %22, %22
%25 = mul i64 %23, %23
%26 = xor i64 %25, -1
%27 = add i64 %24, %26
%28 = lshr i64 %27, 16
%29 = shl i64 %27, 16
%30 = and i64 %29, 4294901760
%31 = or i64 %30, %28
%32 = xor i64 %31, 3899977923
%33 = mul nuw i64 %22, %23
%34 = add i64 %32, %33
%35 = xor i64 %34, %3
store i64 %35, i64* %1
store i64 %20, i64* %0
%36 = load i64, i64* %1
%37 = xor i64 %36, 62708796
%38 = and i64 %37, 65535
%39 = lshr i64 %37, 16
%40 = mul nuw nsw i64 %38, %38
%41 = mul i64 %39, %39
%42 = xor i64 %41, -1
%43 = add i64 %40, %42
%44 = lshr i64 %43, 16
%45 = shl i64 %43, 16
%46 = and i64 %45, 4294901760
%47 = or i64 %46, %44
%48 = xor i64 %47, 1767228838
%49 = mul nuw i64 %38, %39
%50 = add i64 %48, %49
%51 = xor i64 %50, %20
store i64 %51, i64* %1
store i64 %36, i64* %0
%52 = load i64, i64* %1
%53 = xor i64 %52, 255054258
%54 = and i64 %53, 65535
%55 = lshr i64 %53, 16
%56 = mul nuw nsw i64 %54, %54
%57 = mul i64 %55, %55
%58 = xor i64 %57, -1
%59 = add i64 %56, %58
%60 = lshr i64 %59, 16
%61 = shl i64 %59, 16
%62 = and i64 %61, 4294901760
%63 = or i64 %62, %60
%64 = xor i64 %63, 1437059654
%65 = mul nuw i64 %54, %55
%66 = add i64 %64, %65
%67 = xor i64 %66, %36
store i64 %67, i64* %1
store i64 %52, i64* %0
ret void
}