void bn_sqr_words(unsigned long *r, const unsigned long *a, int n) {
  if (n <= 0)
    return;
  while (n & ~3) {
    asm("mulq %2" : "=a"(r[0]), "=d"(r[1]) : "a"(a[0]) : "cc");
    ;
    asm("mulq %2" : "=a"(r[2]), "=d"(r[3]) : "a"(a[1]) : "cc");
    ;
    asm("mulq %2" : "=a"(r[4]), "=d"(r[5]) : "a"(a[2]) : "cc");
    ;
    asm("mulq %2" : "=a"(r[6]), "=d"(r[7]) : "a"(a[3]) : "cc");
    ;
    a += 4;
    r += 8;
    n -= 4;
  }
  if (n) {
    asm("mulq %2" : "=a"(r[0]), "=d"(r[1]) : "a"(a[0]) : "cc");
    ;
    if (--n == 0)
      return;
    asm("mulq %2" : "=a"(r[2]), "=d"(r[3]) : "a"(a[1]) : "cc");
    ;
    if (--n == 0)
      return;
    asm("mulq %2" : "=a"(r[4]), "=d"(r[5]) : "a"(a[2]) : "cc");
    ;
  }
}


x86 assembly (clang, -O0):
<bn_sqr_words>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jg     111e <bn_sqr_words+0x1e>
jmp    1296 <bn_sqr_words+0x196>
jmp    1123 <bn_sqr_words+0x23>
mov    -0x14(%rbp),%eax
and    $0xfffffffc,%eax
cmp    $0x0,%eax
je     11ee <bn_sqr_words+0xee>
mov    -0x8(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mul    %rax
mov    %rax,%rsi
mov    -0x38(%rbp),%rax
mov    %rsi,(%rcx)
mov    %rdx,0x8(%rax)
mov    -0x8(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mul    %rax
mov    %rax,%rsi
mov    -0x30(%rbp),%rax
mov    %rsi,0x10(%rcx)
mov    %rdx,0x18(%rax)
mov    -0x8(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    0x10(%rax),%rax
mul    %rax
mov    %rax,%rsi
mov    -0x28(%rbp),%rax
mov    %rsi,0x20(%rcx)
mov    %rdx,0x28(%rax)
mov    -0x8(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    0x18(%rax),%rax
mul    %rax
mov    %rax,%rsi
mov    -0x20(%rbp),%rax
mov    %rsi,0x30(%rcx)
mov    %rdx,0x38(%rax)
mov    -0x10(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x10(%rbp)
mov    -0x8(%rbp),%rax
add    $0x40,%rax
mov    %rax,-0x8(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x4,%eax
mov    %eax,-0x14(%rbp)
jmp    1123 <bn_sqr_words+0x23>
cmpl   $0x0,-0x14(%rbp)
je     1296 <bn_sqr_words+0x196>
mov    -0x8(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mul    %rax
mov    %rax,%rsi
mov    -0x40(%rbp),%rax
mov    %rsi,(%rcx)
mov    %rdx,0x8(%rax)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
cmp    $0x0,%eax
jne    1233 <bn_sqr_words+0x133>
jmp    1296 <bn_sqr_words+0x196>
mov    -0x8(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mul    %rax
mov    %rax,%rsi
mov    -0x48(%rbp),%rax
mov    %rsi,0x10(%rcx)
mov    %rdx,0x18(%rax)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
cmp    $0x0,%eax
jne    1270 <bn_sqr_words+0x170>
jmp    1296 <bn_sqr_words+0x196>
mov    -0x8(%rbp),%rcx
mov    -0x8(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    0x10(%rax),%rax
mul    %rax
mov    %rax,%rsi
mov    -0x50(%rbp),%rax
mov    %rsi,0x20(%rcx)
mov    %rdx,0x28(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<bn_sqr_words>:
test   %edx,%edx
jle    1193 <bn_sqr_words+0x93>
mov    %edx,%ecx
cmp    $0x4,%edx
jb     115e <bn_sqr_words+0x5e>
nop
mov    (%rsi),%rax
mul    %rax
mov    %rax,(%rdi)
mov    %rdx,0x8(%rdi)
mov    0x8(%rsi),%rax
mul    %rax
mov    %rax,0x10(%rdi)
mov    %rdx,0x18(%rdi)
mov    0x10(%rsi),%rax
mul    %rax
mov    %rax,0x20(%rdi)
mov    %rdx,0x28(%rdi)
mov    0x18(%rsi),%rax
mul    %rax
mov    %rax,0x30(%rdi)
mov    %rdx,0x38(%rdi)
add    $0x20,%rsi
add    $0x40,%rdi
add    $0xfffffffc,%ecx
cmp    $0x3,%ecx
ja     1110 <bn_sqr_words+0x10>
test   %ecx,%ecx
je     1193 <bn_sqr_words+0x93>
mov    (%rsi),%rax
mul    %rax
mov    %rax,(%rdi)
mov    %rdx,0x8(%rdi)
cmp    $0x1,%ecx
je     1193 <bn_sqr_words+0x93>
mov    0x8(%rsi),%rax
mul    %rax
mov    %rax,0x10(%rdi)
mov    %rdx,0x18(%rdi)
cmp    $0x2,%ecx
je     1193 <bn_sqr_words+0x93>
mov    0x10(%rsi),%rax
mul    %rax
mov    %rax,0x20(%rdi)
mov    %rdx,0x28(%rdi)
ret

LLVM IR (-O0):
define void @bn_sqr_words(i64* %0, i64* %1, i32 %2) {
%4 = alloca i64*
%5 = alloca i64*
%6 = alloca i32
store i64* %0, i64** %4
store i64* %1, i64** %5
store i32 %2, i32* %6
%7 = load i32, i32* %6
%8 = icmp sle i32 %7, 0
br i1 %8, label %9, label %10
9: ; preds = %3
br label %106
10: ; preds = %3
br label %11
11: ; preds = %15, %10
%12 = load i32, i32* %6
%13 = and i32 %12, -4
%14 = icmp ne i32 %13, 0
br i1 %14, label %15, label %62
15: ; preds = %11
%16 = load i64*, i64** %4
%17 = getelementptr inbounds i64, i64* %16, i64 0
%18 = load i64*, i64** %4
%19 = getelementptr inbounds i64, i64* %18, i64 1
%20 = load i64*, i64** %5
%21 = getelementptr inbounds i64, i64* %20, i64 0
%22 = load i64, i64* %21
%23 = call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %22)
%24 = extractvalue { i64, i64 } %23, 0
%25 = extractvalue { i64, i64 } %23, 1
store i64 %24, i64* %17
store i64 %25, i64* %19
%26 = load i64*, i64** %4
%27 = getelementptr inbounds i64, i64* %26, i64 2
%28 = load i64*, i64** %4
%29 = getelementptr inbounds i64, i64* %28, i64 3
%30 = load i64*, i64** %5
%31 = getelementptr inbounds i64, i64* %30, i64 1
%32 = load i64, i64* %31
%33 = call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %32)
%34 = extractvalue { i64, i64 } %33, 0
%35 = extractvalue { i64, i64 } %33, 1
store i64 %34, i64* %27
store i64 %35, i64* %29
%36 = load i64*, i64** %4
%37 = getelementptr inbounds i64, i64* %36, i64 4
%38 = load i64*, i64** %4
%39 = getelementptr inbounds i64, i64* %38, i64 5
%40 = load i64*, i64** %5
%41 = getelementptr inbounds i64, i64* %40, i64 2
%42 = load i64, i64* %41
%43 = call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %42)
%44 = extractvalue { i64, i64 } %43, 0
%45 = extractvalue { i64, i64 } %43, 1
store i64 %44, i64* %37
store i64 %45, i64* %39
%46 = load i64*, i64** %4
%47 = getelementptr inbounds i64, i64* %46, i64 6
%48 = load i64*, i64** %4
%49 = getelementptr inbounds i64, i64* %48, i64 7
%50 = load i64*, i64** %5
%51 = getelementptr inbounds i64, i64* %50, i64 3
%52 = load i64, i64* %51
%53 = call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %52)
%54 = extractvalue { i64, i64 } %53, 0
%55 = extractvalue { i64, i64 } %53, 1
store i64 %54, i64* %47
store i64 %55, i64* %49
%56 = load i64*, i64** %5
%57 = getelementptr inbounds i64, i64* %56, i64 4
store i64* %57, i64** %5
%58 = load i64*, i64** %4
%59 = getelementptr inbounds i64, i64* %58, i64 8
store i64* %59, i64** %4
%60 = load i32, i32* %6
%61 = sub nsw i32 %60, 4
store i32 %61, i32* %6
br label %11
62: ; preds = %11
%63 = load i32, i32* %6
%64 = icmp ne i32 %63, 0
br i1 %64, label %65, label %106
65: ; preds = %62
%66 = load i64*, i64** %4
%67 = getelementptr inbounds i64, i64* %66, i64 0
%68 = load i64*, i64** %4
%69 = getelementptr inbounds i64, i64* %68, i64 1
%70 = load i64*, i64** %5
%71 = getelementptr inbounds i64, i64* %70, i64 0
%72 = load i64, i64* %71
%73 = call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %72)
%74 = extractvalue { i64, i64 } %73, 0
%75 = extractvalue { i64, i64 } %73, 1
store i64 %74, i64* %67
store i64 %75, i64* %69
%76 = load i32, i32* %6
%77 = add nsw i32 %76, -1
store i32 %77, i32* %6
%78 = icmp eq i32 %77, 0
br i1 %78, label %79, label %80
79: ; preds = %65
br label %106
80: ; preds = %65
%81 = load i64*, i64** %4
%82 = getelementptr inbounds i64, i64* %81, i64 2
%83 = load i64*, i64** %4
%84 = getelementptr inbounds i64, i64* %83, i64 3
%85 = load i64*, i64** %5
%86 = getelementptr inbounds i64, i64* %85, i64 1
%87 = load i64, i64* %86
%88 = call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %87)
%89 = extractvalue { i64, i64 } %88, 0
%90 = extractvalue { i64, i64 } %88, 1
store i64 %89, i64* %82
store i64 %90, i64* %84
%91 = load i32, i32* %6
%92 = add nsw i32 %91, -1
store i32 %92, i32* %6
%93 = icmp eq i32 %92, 0
br i1 %93, label %94, label %95
94: ; preds = %80
br label %106
95: ; preds = %80
%96 = load i64*, i64** %4
%97 = getelementptr inbounds i64, i64* %96, i64 4
%98 = load i64*, i64** %4
%99 = getelementptr inbounds i64, i64* %98, i64 5
%100 = load i64*, i64** %5
%101 = getelementptr inbounds i64, i64* %100, i64 2
%102 = load i64, i64* %101
%103 = call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %102)
%104 = extractvalue { i64, i64 } %103, 0
%105 = extractvalue { i64, i64 } %103, 1
store i64 %104, i64* %97
store i64 %105, i64* %99
br label %106
106: ; preds = %9, %79, %94, %95, %62
ret void
}

LLVM IR (-O2):
define void @bn_sqr_words(i64* nocapture writeonly %0, i64* nocapture %1, i32 %2) {
%4 = icmp slt i32 %2, 1
br i1 %4, label %70, label %5
5: ; preds = %3
%6 = icmp ult i32 %2, 4
br i1 %6, label %43, label %7
7: ; preds = %5, %7
%8 = phi i32 [ %39, %7 ], [ %2, %5 ]
%9 = phi i64* [ %37, %7 ], [ %1, %5 ]
%10 = phi i64* [ %38, %7 ], [ %0, %5 ]
%11 = getelementptr inbounds i64, i64* %10, i64 1
%12 = load i64, i64* %9
%13 = tail call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %12)
%14 = extractvalue { i64, i64 } %13, 0
%15 = extractvalue { i64, i64 } %13, 1
store i64 %14, i64* %10
store i64 %15, i64* %11
%16 = getelementptr inbounds i64, i64* %10, i64 2
%17 = getelementptr inbounds i64, i64* %10, i64 3
%18 = getelementptr inbounds i64, i64* %9, i64 1
%19 = load i64, i64* %18
%20 = tail call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %19)
%21 = extractvalue { i64, i64 } %20, 0
%22 = extractvalue { i64, i64 } %20, 1
store i64 %21, i64* %16
store i64 %22, i64* %17
%23 = getelementptr inbounds i64, i64* %10, i64 4
%24 = getelementptr inbounds i64, i64* %10, i64 5
%25 = getelementptr inbounds i64, i64* %9, i64 2
%26 = load i64, i64* %25
%27 = tail call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %26)
%28 = extractvalue { i64, i64 } %27, 0
%29 = extractvalue { i64, i64 } %27, 1
store i64 %28, i64* %23
store i64 %29, i64* %24
%30 = getelementptr inbounds i64, i64* %10, i64 6
%31 = getelementptr inbounds i64, i64* %10, i64 7
%32 = getelementptr inbounds i64, i64* %9, i64 3
%33 = load i64, i64* %32
%34 = tail call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %33)
%35 = extractvalue { i64, i64 } %34, 0
%36 = extractvalue { i64, i64 } %34, 1
store i64 %35, i64* %30
store i64 %36, i64* %31
%37 = getelementptr inbounds i64, i64* %9, i64 4
%38 = getelementptr inbounds i64, i64* %10, i64 8
%39 = add nsw i32 %8, -4
%40 = icmp ult i32 %39, 4
br i1 %40, label %41, label %7
41: ; preds = %7
%42 = icmp eq i32 %39, 0
br i1 %42, label %70, label %43
43: ; preds = %5, %41
%44 = phi i32 [ %39, %41 ], [ %2, %5 ]
%45 = phi i64* [ %37, %41 ], [ %1, %5 ]
%46 = phi i64* [ %38, %41 ], [ %0, %5 ]
%47 = getelementptr inbounds i64, i64* %46, i64 1
%48 = load i64, i64* %45
%49 = tail call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %48)
%50 = extractvalue { i64, i64 } %49, 0
%51 = extractvalue { i64, i64 } %49, 1
store i64 %50, i64* %46
store i64 %51, i64* %47
%52 = icmp eq i32 %44, 1
br i1 %52, label %70, label %53
53: ; preds = %43
%54 = getelementptr inbounds i64, i64* %46, i64 2
%55 = getelementptr inbounds i64, i64* %46, i64 3
%56 = getelementptr inbounds i64, i64* %45, i64 1
%57 = load i64, i64* %56
%58 = tail call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %57)
%59 = extractvalue { i64, i64 } %58, 0
%60 = extractvalue { i64, i64 } %58, 1
store i64 %59, i64* %54
store i64 %60, i64* %55
%61 = icmp eq i32 %44, 2
br i1 %61, label %70, label %62
62: ; preds = %53
%63 = getelementptr inbounds i64, i64* %46, i64 4
%64 = getelementptr inbounds i64, i64* %46, i64 5
%65 = getelementptr inbounds i64, i64* %45, i64 2
%66 = load i64, i64* %65
%67 = tail call { i64, i64 } asm "mulq $2", "={ax},={dx},{ax},~{cc},~{dirflag},~{fpsr},~{flags}"(i64 %66)
%68 = extractvalue { i64, i64 } %67, 0
%69 = extractvalue { i64, i64 } %67, 1
store i64 %68, i64* %63
store i64 %69, i64* %64
br label %70
70: ; preds = %53, %43, %3, %62, %41
ret void
}