// File: occivink[P]pyct/monocypher.cpp
// Function: poly_block(crypto_poly1305_ctx*, unsigned char const*, unsigned int)
static void poly_block(crypto_poly1305_ctx *ctx, const u8 in[16], unsigned end)
{
    u32 s[4];
    load32_le_buf(s, in, 4);

    // s = h + c, without carry propagation
    const u64 s0 = ctx->h[0] + (u64)s[0]; // s0 <= 1_fffffffe
    const u64 s1 = ctx->h[1] + (u64)s[1]; // s1 <= 1_fffffffe
    const u64 s2 = ctx->h[2] + (u64)s[2]; // s2 <= 1_fffffffe
    const u64 s3 = ctx->h[3] + (u64)s[3]; // s3 <= 1_fffffffe
    const u32 s4 = ctx->h[4] + end;       // s4 <=          5

    // Local all the things!
    const u32 r0 = ctx->r[0];       // r0  <= 0fffffff
    const u32 r1 = ctx->r[1];       // r1  <= 0ffffffc
    const u32 r2 = ctx->r[2];       // r2  <= 0ffffffc
    const u32 r3 = ctx->r[3];       // r3  <= 0ffffffc
    const u32 rr0 = (r0 >> 2) * 5;  // rr0 <= 13fffffb // lose 2 bits...
    const u32 rr1 = (r1 >> 2) + r1; // rr1 <= 13fffffb // rr1 == (r1 >> 2) * 5
    const u32 rr2 = (r2 >> 2) + r2; // rr2 <= 13fffffb // rr1 == (r2 >> 2) * 5
    const u32 rr3 = (r3 >> 2) + r3; // rr3 <= 13fffffb // rr1 == (r3 >> 2) * 5

    // (h + c) * r, without carry propagation
    const u64 x0 = s0*r0+ s1*rr3+ s2*rr2+ s3*rr1+ s4*rr0; // <= 97ffffe007fffff8
    const u64 x1 = s0*r1+ s1*r0 + s2*rr3+ s3*rr2+ s4*rr1; // <= 8fffffe20ffffff6
    const u64 x2 = s0*r2+ s1*r1 + s2*r0 + s3*rr3+ s4*rr2; // <= 87ffffe417fffff4
    const u64 x3 = s0*r3+ s1*r2 + s2*r1 + s3*r0 + s4*rr3; // <= 7fffffe61ffffff2
    const u32 x4 = s4 * (r0 & 3); // ...recover 2 bits    // <=                f

    // partial reduction modulo 2^130 - 5
    const u32 u5 = x4 + (x3 >> 32); // u5 <= 7ffffff5
    const u64 u0 = (u5 >>  2) * 5 + (x0 & 0xffffffff);
    const u64 u1 = (u0 >> 32)     + (x1 & 0xffffffff) + (x0 >> 32);
    const u64 u2 = (u1 >> 32)     + (x2 & 0xffffffff) + (x1 >> 32);
    const u64 u3 = (u2 >> 32)     + (x3 & 0xffffffff) + (x2 >> 32);
    const u64 u4 = (u3 >> 32)     + (u5 & 3);

    // Update the hash
    ctx->h[0] = (u32)u0; // u0 <= 1_9ffffff0
    ctx->h[1] = (u32)u1; // u1 <= 1_97ffffe0
    ctx->h[2] = (u32)u2; // u2 <= 1_8fffffe2
    ctx->h[3] = (u32)u3; // u3 <= 1_87ffffe4
    ctx->h[4] = (u32)u4; // u4 <=          4
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xd0, %rsp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movl %edx, -0x14(%rbp)
leaq -0x30(%rbp), %rdi
movq -0x10(%rbp), %rsi
movl $0x4, %edx
callq 0x3990
movq -0x8(%rbp), %rax
movl 0x10(%rax), %eax
movl -0x30(%rbp), %ecx
addq %rcx, %rax
movq %rax, -0x38(%rbp)
movq -0x8(%rbp), %rax
movl 0x14(%rax), %eax
movl -0x2c(%rbp), %ecx
addq %rcx, %rax
movq %rax, -0x40(%rbp)
movq -0x8(%rbp), %rax
movl 0x18(%rax), %eax
movl -0x28(%rbp), %ecx
addq %rcx, %rax
movq %rax, -0x48(%rbp)
movq -0x8(%rbp), %rax
movl 0x1c(%rax), %eax
movl -0x24(%rbp), %ecx
addq %rcx, %rax
movq %rax, -0x50(%rbp)
movq -0x8(%rbp), %rax
movl 0x20(%rax), %eax
addl -0x14(%rbp), %eax
movl %eax, -0x54(%rbp)
movq -0x8(%rbp), %rax
movl (%rax), %eax
movl %eax, -0x58(%rbp)
movq -0x8(%rbp), %rax
movl 0x4(%rax), %eax
movl %eax, -0x5c(%rbp)
movq -0x8(%rbp), %rax
movl 0x8(%rax), %eax
movl %eax, -0x60(%rbp)
movq -0x8(%rbp), %rax
movl 0xc(%rax), %eax
movl %eax, -0x64(%rbp)
movl -0x58(%rbp), %eax
shrl $0x2, %eax
imull $0x5, %eax, %eax
movl %eax, -0x68(%rbp)
movl -0x5c(%rbp), %eax
shrl $0x2, %eax
addl -0x5c(%rbp), %eax
movl %eax, -0x6c(%rbp)
movl -0x60(%rbp), %eax
shrl $0x2, %eax
addl -0x60(%rbp), %eax
movl %eax, -0x70(%rbp)
movl -0x64(%rbp), %eax
shrl $0x2, %eax
addl -0x64(%rbp), %eax
movl %eax, -0x74(%rbp)
movq -0x38(%rbp), %rax
movl -0x58(%rbp), %ecx
imulq %rcx, %rax
movq -0x40(%rbp), %rcx
movl -0x74(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x48(%rbp), %rcx
movl -0x70(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x50(%rbp), %rcx
movl -0x6c(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movl -0x54(%rbp), %ecx
imull -0x68(%rbp), %ecx
movl %ecx, %ecx
addq %rcx, %rax
movq %rax, -0x80(%rbp)
movq -0x38(%rbp), %rax
movl -0x5c(%rbp), %ecx
imulq %rcx, %rax
movq -0x40(%rbp), %rcx
movl -0x58(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x48(%rbp), %rcx
movl -0x74(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x50(%rbp), %rcx
movl -0x70(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movl -0x54(%rbp), %ecx
imull -0x6c(%rbp), %ecx
movl %ecx, %ecx
addq %rcx, %rax
movq %rax, -0x88(%rbp)
movq -0x38(%rbp), %rax
movl -0x60(%rbp), %ecx
imulq %rcx, %rax
movq -0x40(%rbp), %rcx
movl -0x5c(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x48(%rbp), %rcx
movl -0x58(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x50(%rbp), %rcx
movl -0x74(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movl -0x54(%rbp), %ecx
imull -0x70(%rbp), %ecx
movl %ecx, %ecx
addq %rcx, %rax
movq %rax, -0x90(%rbp)
movq -0x38(%rbp), %rax
movl -0x64(%rbp), %ecx
imulq %rcx, %rax
movq -0x40(%rbp), %rcx
movl -0x60(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x48(%rbp), %rcx
movl -0x5c(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movq -0x50(%rbp), %rcx
movl -0x58(%rbp), %edx
imulq %rdx, %rcx
addq %rcx, %rax
movl -0x54(%rbp), %ecx
imull -0x74(%rbp), %ecx
movl %ecx, %ecx
addq %rcx, %rax
movq %rax, -0x98(%rbp)
movl -0x54(%rbp), %eax
movl -0x58(%rbp), %ecx
andl $0x3, %ecx
imull %ecx, %eax
movl %eax, -0x9c(%rbp)
movl -0x9c(%rbp), %eax
movq -0x98(%rbp), %rcx
shrq $0x20, %rcx
addq %rcx, %rax
movl %eax, -0xa0(%rbp)
movl -0xa0(%rbp), %eax
shrl $0x2, %eax
imull $0x5, %eax, %eax
movl %eax, %eax
movabsq $0xffffffff, %rcx       # imm = 0xFFFFFFFF
andq -0x80(%rbp), %rcx
addq %rcx, %rax
movq %rax, -0xa8(%rbp)
movq -0xa8(%rbp), %rax
shrq $0x20, %rax
movabsq $0xffffffff, %rcx       # imm = 0xFFFFFFFF
andq -0x88(%rbp), %rcx
addq %rcx, %rax
movq -0x80(%rbp), %rcx
shrq $0x20, %rcx
addq %rcx, %rax
movq %rax, -0xb0(%rbp)
movq -0xb0(%rbp), %rax
shrq $0x20, %rax
movabsq $0xffffffff, %rcx       # imm = 0xFFFFFFFF
andq -0x90(%rbp), %rcx
addq %rcx, %rax
movq -0x88(%rbp), %rcx
shrq $0x20, %rcx
addq %rcx, %rax
movq %rax, -0xb8(%rbp)
movq -0xb8(%rbp), %rax
shrq $0x20, %rax
movabsq $0xffffffff, %rcx       # imm = 0xFFFFFFFF
andq -0x98(%rbp), %rcx
addq %rcx, %rax
movq -0x90(%rbp), %rcx
shrq $0x20, %rcx
addq %rcx, %rax
movq %rax, -0xc0(%rbp)
movq -0xc0(%rbp), %rax
shrq $0x20, %rax
movl -0xa0(%rbp), %ecx
andl $0x3, %ecx
movl %ecx, %ecx
addq %rcx, %rax
movq %rax, -0xc8(%rbp)
movq -0xa8(%rbp), %rax
movl %eax, %ecx
movq -0x8(%rbp), %rax
movl %ecx, 0x10(%rax)
movq -0xb0(%rbp), %rax
movl %eax, %ecx
movq -0x8(%rbp), %rax
movl %ecx, 0x14(%rax)
movq -0xb8(%rbp), %rax
movl %eax, %ecx
movq -0x8(%rbp), %rax
movl %ecx, 0x18(%rax)
movq -0xc0(%rbp), %rax
movl %eax, %ecx
movq -0x8(%rbp), %rax
movl %ecx, 0x1c(%rax)
movq -0xc8(%rbp), %rax
movl %eax, %ecx
movq -0x8(%rbp), %rax
movl %ecx, 0x20(%rax)
addq $0xd0, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)