user: What would the C source look like for the assembly code provided here:

```asm
<make_double_block>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xb8,%rsp
  mov    %rdi,0x78(%rsp)
  mov    %rsi,0x80(%rsp)
  mov    %rdx,0x30(%rsp)
  mov    %rcx,0x38(%rsp)
  mov    %r8,0xa8(%rsp)
  mov    %r9,0x70(%rsp)
  mov    0x108(%rsp),%r15d
  mov    0x110(%rsp),%r13d
  mov    0x118(%rsp),%r14d
  mov    0x120(%rsp),%r12d
  mov    %r15d,%eax
  imul   0x128(%rsp),%eax
  imul   %r13d,%eax
  imul   %r15d,%eax
  imul   %r13d,%eax
  imul   %r12d,%eax
  imul   %r15d,%eax
  imul   %r13d,%eax
  imul   %r14d,%eax
  lea    0x0(,%rax,8),%ebx
  mov    %ebx,%edi
  call   <malloc@plt>
  mov    %rax,0x50(%rsp)
  mov    %ebx,%edi
  call   <malloc@plt>
  mov    %rax,0x58(%rsp)
  mov    %ebx,%edi
  call   <malloc@plt>
  mov    %rax,0x60(%rsp)
  mov    %ebx,%edi
  call   <malloc@plt>
  mov    %rax,0x68(%rsp)
  movl   $0x0,0x4(%rsp)
  cmpl   $0x0,0x128(%rsp)
  jle    L120
  mov    %r15d,%ebp
  mov    %r13d,%ebx
  mov    %r14d,0x118(%rsp)
  jmp    L498
L801:
  sub    $0x8,%rsp
  push   %r12
  mov    0x128(%rsp),%eax
  push   %rax
  push   %rbx
  push   %r15
  push   %r14
  push   %r13
  push   %rbp
  mov    0x50(%rsp),%r9d
  mov    0x54(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x90(%rsp),%r10
  movq   $0x0,(%r10,%rax,8)
  add    $0x38,%rsp
  push   %r12
  mov    0x128(%rsp),%eax
  push   %rax
  push   %rbx
  push   %r15
  push   %r14
  push   %r13
  push   %rbp
  mov    0x50(%rsp),%r9d
  mov    0x54(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x98(%rsp),%r11
  movq   $0x0,(%r11,%rax,8)
  add    $0x38,%rsp
  push   %r12
  mov    0x128(%rsp),%eax
  push   %rax
  push   %rbx
  push   %r15
  push   %r14
  push   %r13
  push   %rbp
  mov    0x50(%rsp),%r9d
  mov    0x54(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa0(%rsp),%r10
  movq   $0x0,(%r10,%rax,8)
  add    $0x38,%rsp
  push   %r12
  mov    0x128(%rsp),%eax
  push   %rax
  push   %rbx
  push   %r15
  push   %r14
  push   %r13
  push   %rbp
  mov    0x50(%rsp),%r9d
  mov    0x54(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa8(%rsp),%rcx
  movq   $0x0,(%rcx,%rax,8)
  add    $0x1,%r14d
  add    $0x40,%rsp
  cmp    %r14d,0x118(%rsp)
  jne    L801
L258:
  add    $0x1,%r13d
  cmp    %r13d,%ebx
  je     L543
L289:
  mov    $0x0,%r14d
  cmpl   $0x0,0x118(%rsp)
  jg     L801
  jmp    L258
L543:
  add    $0x1,%ebp
  cmp    %ebp,%r15d
  je     L85
L522:
  mov    $0x0,%r13d
  jmp    L289
L85:
  addl   $0x1,0x10(%rsp)
  mov    0x10(%rsp),%eax
  cmp    %eax,%r12d
  je     L309
L471:
  movl   $0x0,0x1c(%rsp)
  cmpl   $0x0,0x140(%rsp)
  jle    L619
L119:
  mov    $0x0,%ebp
L677:
  mov    $0x0,%r13d
L685:
  mov    $0x0,%r14d
L202:
  mov    0x140(%rsp),%eax
  push   %rax
  push   %r12
  push   %rbx
  push   %r15
  mov    0x3c(%rsp),%eax
  push   %rax
  push   %r14
  push   %r13
  push   %rbp
  mov    0x50(%rsp),%r9d
  mov    0x54(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Bdouble_index@plt>
  mov    0xb8(%rsp),%r11
  movq   $0x0,(%r11,%rax,8)
  add    $0x40,%rsp
  mov    0x140(%rsp),%r11d
  push   %r11
  push   %r12
  push   %rbx
  push   %r15
  mov    0x3c(%rsp),%eax
  push   %rax
  push   %r14
  push   %r13
  push   %rbp
  mov    0x50(%rsp),%r9d
  mov    0x54(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Bdouble_index@plt>
  mov    0xc0(%rsp),%rcx
  movq   $0x0,(%rcx,%rax,8)
  add    $0x1,%r14d
  add    $0x40,%rsp
  cmp    %r14d,%r12d
  jne    L202
  add    $0x1,%r13d
  cmp    %r13d,%ebx
  jne    L685
  add    $0x1,%ebp
  cmp    %ebp,%r15d
  jne    L677
  addl   $0x1,0x1c(%rsp)
  mov    0x1c(%rsp),%eax
  cmp    %eax,0x140(%rsp)
  jne    L119
L619:
  mov    $0x0,%ebp
  jmp    L522
L309:
  mov    0x14(%rsp),%eax
  add    $0x1,%eax
  cmp    %eax,%ebx
  je     L609
  mov    %eax,0x14(%rsp)
L98:
  movl   $0x0,0x10(%rsp)
  test   %r12d,%r12d
  jg     L471
  jmp    L309
L609:
  addl   $0x1,0xc(%rsp)
  mov    0xc(%rsp),%eax
  cmp    %eax,%r15d
  je     L582
L149:
  movl   $0x0,0x14(%rsp)
  jmp    L98
L582:
  mov    0x18(%rsp),%ecx
  lea    0x1(%rcx),%eax
  cmp    0x14(%rsp),%ecx
  je     L477
  mov    %eax,0x18(%rsp)
L550:
  movl   $0x0,0xc(%rsp)
  jmp    L149
L477:
  mov    %r15d,%ebp
L696:
  addl   $0x1,0x8(%rsp)
  mov    0x8(%rsp),%eax
  cmp    %eax,%ebp
  je     L692
L1:
  movl   $0x0,0x18(%rsp)
  mov    %ebp,%r15d
  test   %ebx,%ebx
  jg     L550
  jmp    L696
L692:
  addl   $0x1,0x4(%rsp)
  mov    0x4(%rsp),%eax
  cmp    %eax,0x128(%rsp)
  je     L817
L498:
  movl   $0x0,0x8(%rsp)
  test   %ebp,%ebp
  jg     L1
  jmp    L692
L817:
  mov    %ebp,%r15d
  mov    %ebx,%r13d
  mov    0x118(%rsp),%r14d
  cmpl   $0x0,0x130(%rsp)
  jle    L392
L592:
  mov    0x130(%rsp),%eax
  mov    %rax,0xa0(%rsp)
  movq   $0x0,0x20(%rsp)
  mov    %r12d,0x120(%rsp)
  mov    %r13d,%r12d
  mov    %r14d,%r13d
  mov    %r15d,%r14d
  mov    0x120(%rsp),%r15d
  jmp    L33
L722:
  mov    0xf0(%rsp),%rcx
  mov    0x20(%rsp),%rax
  movsd  (%rcx,%rax,8),%xmm3
  movsd  %xmm3,0x88(%rsp)
  sub    $0x8,%rsp
  mov    0x130(%rsp),%ecx
  push   %rcx
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  (%rdi,%rax,8),%xmm5
  movsd  %xmm5,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r11d
  push   %r11
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x88(%rsp),%ecx
  mov    0x80(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0xd0(%rsp),%xmm5
  mulsd  (%rdi,%rax,8),%xmm5
  movsd  %xmm5,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r10d
  push   %r10
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  (%r10,%rax,8),%xmm1
  movsd  %xmm1,0xd8(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r11d
  push   %r11
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x88(%rsp),%ecx
  mov    0x80(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0xd8(%rsp),%xmm1
  mulsd  (%r10,%rax,8),%xmm1
  movsd  0xd0(%rsp),%xmm5
  subsd  %xmm1,%xmm5
  movapd %xmm5,%xmm0
  mulsd  0xc8(%rsp),%xmm0
  movsd  %xmm0,0xc8(%rsp)
  add    $0x38,%rsp
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %r12
  mov    0x5c(%rsp),%eax
  push   %rax
  mov    0x60(%rsp),%r11d
  push   %r11
  mov    %ebp,%r9d
  mov    0x50(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x54(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x90(%rsp),%r8
  lea    (%r8,%rax,8),%rax
  movsd  0xc8(%rsp),%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  mov    0x130(%rsp),%r8
  mov    0x60(%rsp),%r9
  movsd  (%r8,%r9,8),%xmm7
  movsd  %xmm7,0xc8(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r8d
  push   %r8
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  (%rdi,%rax,8),%xmm2
  movsd  %xmm2,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%esi
  push   %rsi
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x88(%rsp),%ecx
  mov    0x80(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0xd0(%rsp),%xmm2
  mulsd  (%r10,%rax,8),%xmm2
  movsd  %xmm2,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%edx
  push   %rdx
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  (%r10,%rax,8),%xmm5
  movsd  %xmm5,0xd8(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%esi
  push   %rsi
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x88(%rsp),%ecx
  mov    0x80(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0xd8(%rsp),%xmm5
  mulsd  (%rdi,%rax,8),%xmm5
  movapd %xmm5,%xmm0
  addsd  0xd0(%rsp),%xmm0
  mulsd  0xc8(%rsp),%xmm0
  movsd  %xmm0,0xc8(%rsp)
  add    $0x38,%rsp
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %r12
  mov    0x5c(%rsp),%eax
  push   %rax
  mov    0x60(%rsp),%r11d
  push   %r11
  mov    %ebp,%r9d
  mov    0x50(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x54(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x98(%rsp),%r8
  lea    (%r8,%rax,8),%rax
  movsd  0xc8(%rsp),%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  mov    0x130(%rsp),%r8
  mov    0x60(%rsp),%r9
  movsd  (%r8,%r9,8),%xmm4
  movsd  %xmm4,0xc8(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r8d
  push   %r8
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x6c(%rsp),%ecx
  mov    0x68(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  (%rdi,%rax,8),%xmm6
  movsd  %xmm6,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r8d
  push   %r8
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0xd0(%rsp),%xmm6
  mulsd  (%rdi,%rax,8),%xmm6
  movsd  %xmm6,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%esi
  push   %rsi
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x6c(%rsp),%ecx
  mov    0x68(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  (%r10,%rax,8),%xmm7
  movsd  %xmm7,0xd8(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%esi
  push   %rsi
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0xd8(%rsp),%xmm1
  mulsd  (%r10,%rax,8),%xmm1
  movsd  0xd0(%rsp),%xmm6
  subsd  %xmm1,%xmm6
  movapd %xmm6,%xmm0
  mulsd  0xc8(%rsp),%xmm0
  movsd  %xmm0,0xc8(%rsp)
  add    $0x38,%rsp
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %r12
  mov    0x78(%rsp),%esi
  push   %rsi
  mov    0x78(%rsp),%ecx
  push   %rcx
  mov    %ebp,%r9d
  mov    0x50(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x54(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa0(%rsp),%r8
  lea    (%r8,%rax,8),%rax
  movsd  0xc8(%rsp),%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  mov    0x130(%rsp),%r8
  mov    0x60(%rsp),%r9
  movsd  (%r8,%r9,8),%xmm5
  movsd  %xmm5,0xc8(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r9d
  push   %r9
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x6c(%rsp),%ecx
  mov    0x68(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  (%rdi,%rax,8),%xmm6
  movsd  %xmm6,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r9d
  push   %r9
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0xd0(%rsp),%xmm6
  mulsd  (%r10,%rax,8),%xmm6
  movsd  %xmm6,0xd0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%esi
  push   %rsi
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x50(%rsp),%r9d
  mov    0x4c(%rsp),%r8d
  mov    0x6c(%rsp),%ecx
  mov    0x68(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  (%r10,%rax,8),%xmm3
  movsd  %xmm3,0xd8(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r10d
  push   %r10
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %rbp
  push   %r12
  mov    0x54(%rsp),%r9d
  mov    0x48(%rsp),%r8d
  mov    0x5c(%rsp),%ecx
  mov    0x58(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x1,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0xd8(%rsp),%xmm3
  mulsd  (%rdi,%rax,8),%xmm3
  movapd %xmm3,%xmm0
  addsd  0xd0(%rsp),%xmm0
  mulsd  0xc8(%rsp),%xmm0
  movsd  %xmm0,0xc8(%rsp)
  add    $0x38,%rsp
  push   %r15
  push   %r14
  push   %r13
  push   %rbx
  push   %r12
  mov    0x78(%rsp),%eax
  push   %rax
  mov    0x78(%rsp),%eax
  push   %rax
  mov    %ebp,%r9d
  mov    0x50(%rsp),%r8d
  mov    0x4c(%rsp),%ecx
  mov    0x54(%rsp),%edx
  mov    0x48(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa8(%rsp),%rcx
  lea    (%rcx,%rax,8),%rax
  movsd  0xc8(%rsp),%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x1,%ebp
  add    $0x40,%rsp
  cmp    %ebp,%r15d
  jne    L722
L511:
  add    $0x1,%r12d
  cmp    %r12d,%r14d
  je     L89
L225:
  mov    $0x0,%ebp
  test   %r15d,%r15d
  jg     L722
  jmp    L511
L89:
  mov    %ebx,%ebp
L438:
  mov    0x10(%rsp),%eax
  add    $0x1,%eax
  cmp    %eax,%r13d
  je     L6
  mov    %eax,0x10(%rsp)
L125:
  mov    $0x0,%r12d
  mov    %ebp,%ebx
  test   %r14d,%r14d
  jg     L225
  jmp    L438
L6:
  mov    %ebp,%r12d
  addl   $0x1,0xc(%rsp)
  mov    0xc(%rsp),%eax
  cmp    %eax,%ebp
  je     L274
L216:
  movl   $0x0,0x10(%rsp)
  mov    %r12d,%ebp
  jmp    L125
L274:
  mov    0x14(%rsp),%edx
  lea    0x1(%rdx),%eax
  cmp    0x10(%rsp),%edx
  je     L454
  mov    %eax,0x14(%rsp)
L694:
  movl   $0x0,0xc(%rsp)
  jmp    L216
L454:
  mov    %r13d,%eax
  mov    %ebp,%r13d
  mov    %eax,%r12d
L88:
  addl   $0x1,0x8(%rsp)
  mov    0x8(%rsp),%eax
  cmp    %eax,%r13d
  je     L142
L171:
  movl   $0x0,0x14(%rsp)
  test   %r12d,%r12d
  jle    L88
  mov    %r12d,%eax
  mov    %r13d,%r12d
  mov    %eax,%r13d
  jmp    L694
L142:
  mov    %r15d,%eax
  mov    %r14d,%r15d
  mov    %r13d,%r14d
  mov    %r12d,%r13d
  mov    %eax,%r12d
L366:
  addl   $0x1,0x4(%rsp)
  mov    0x4(%rsp),%eax
  cmp    %eax,0x128(%rsp)
  je     L624
L4:
  movl   $0x0,0x8(%rsp)
  test   %r14d,%r14d
  jle    L366
  mov    %r12d,0x120(%rsp)
  mov    %r13d,%r12d
  mov    %r14d,%r13d
  mov    %r15d,%r14d
  mov    0x120(%rsp),%r15d
  jmp    L171
L624:
  mov    %r12d,0x120(%rsp)
  mov    %r13d,%r12d
  mov    %r15d,%r13d
  mov    0x120(%rsp),%r15d
L224:
  addq   $0x1,0x20(%rsp)
  mov    0x20(%rsp),%rax
  cmp    %rax,0xa0(%rsp)
  je     L355
L33:
  mov    0x20(%rsp),%eax
  mov    0x138(%rsp),%edx
  mov    $0x0,%esi
  mov    %eax,%ebx
  mov    %eax,%edi
  call   <index_2d@plt>
  mov    0xa8(%rsp),%rbp
  mov    0x0(%rbp,%rax,4),%eax
  mov    %eax,0x28(%rsp)
  mov    0x138(%rsp),%edx
  mov    $0x0,%esi
  mov    %ebx,%edi
  call   <index_2d@plt>
  mov    0x70(%rsp),%rcx
  mov    (%rcx,%rax,4),%eax
  mov    %eax,0x2c(%rsp)
  mov    0x138(%rsp),%edx
  mov    $0x1,%esi
  mov    %ebx,%edi
  call   <index_2d@plt>
  mov    0x0(%rbp,%rax,4),%eax
  mov    %eax,0x18(%rsp)
  mov    0x138(%rsp),%edx
  mov    $0x1,%esi
  mov    %ebx,%edi
  call   <index_2d@plt>
  mov    0x70(%rsp),%rcx
  mov    (%rcx,%rax,4),%eax
  mov    %eax,0x1c(%rsp)
  mov    0x138(%rsp),%edx
  mov    $0x2,%esi
  mov    %ebx,%edi
  call   <index_2d@plt>
  mov    0x0(%rbp,%rax,4),%eax
  mov    %eax,0x40(%rsp)
  mov    0x138(%rsp),%edx
  mov    $0x2,%esi
  mov    %ebx,%edi
  call   <index_2d@plt>
  mov    0x70(%rsp),%rcx
  mov    (%rcx,%rax,4),%eax
  mov    %eax,0x48(%rsp)
  movl   $0x0,0x4(%rsp)
  cmpl   $0x0,0x128(%rsp)
  jle    L224
  mov    %r15d,0x120(%rsp)
  mov    %r13d,%r15d
  mov    %r12d,%r13d
  mov    0x120(%rsp),%r12d
  jmp    L4
L355:
  mov    %r15d,0x120(%rsp)
  mov    %r14d,%r15d
  mov    %r13d,%r14d
  mov    %r12d,%r13d
  mov    0x120(%rsp),%r12d
L392:
  movl   $0x0,0x4(%rsp)
  cmpl   $0x0,0x128(%rsp)
  jle    L317
  mov    %r13d,0x110(%rsp)
  mov    %r15d,%r13d
  mov    %r12d,%ebp
  mov    0x110(%rsp),%r15d
  jmp    L616
L464:
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %r13d,%edi
  call   <index_2d@plt>
  mov    0xf8(%rsp),%rsi
  movsd  0x40(%rsp),%xmm2
  mulsd  (%rsi,%rax,8),%xmm2
  movsd  %xmm2,0x70(%rsp)
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %r13d,%edi
  call   <index_2d@plt>
  mov    0x100(%rsp),%rcx
  movsd  0x48(%rsp),%xmm0
  mulsd  (%rcx,%rax,8),%xmm0
  movsd  0x70(%rsp),%xmm3
  subsd  %xmm0,%xmm3
  movsd  %xmm3,0x70(%rsp)
  push   %r15
  push   %r14
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x5c(%rsp),%edx
  push   %rdx
  mov    0x44(%rsp),%ecx
  push   %rcx
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x60(%rsp),%edx
  mov    0x5c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Bdouble_index@plt>
  mov    0xb8(%rsp),%r11
  lea    (%r11,%rax,8),%rax
  movsd  0xb0(%rsp),%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x40,%rsp
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %r13d,%edi
  call   <index_2d@plt>
  mov    0xf8(%rsp),%r11
  movsd  0x48(%rsp),%xmm0
  mulsd  (%r11,%rax,8),%xmm0
  movsd  %xmm0,0x70(%rsp)
  mov    %r15d,%edx
  mov    %ebx,%esi
  mov    %r13d,%edi
  call   <index_2d@plt>
  mov    0x100(%rsp),%r10
  movsd  0x40(%rsp),%xmm2
  mulsd  (%r10,%rax,8),%xmm2
  movapd %xmm2,%xmm7
  addsd  0x70(%rsp),%xmm7
  movsd  %xmm7,0x70(%rsp)
  push   %r15
  push   %r14
  push   %r12
  push   %rbp
  push   %rbx
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x5c(%rsp),%eax
  push   %rax
  mov    0x44(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x60(%rsp),%edx
  mov    0x5c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Bdouble_index@plt>
  mov    0xc0(%rsp),%rsi
  lea    (%rsi,%rax,8),%rax
  movsd  0xb0(%rsp),%xmm7
  addsd  (%rax),%xmm7
  movsd  %xmm7,(%rax)
  add    $0x1,%ebx
  add    $0x40,%rsp
  cmp    %ebx,%r15d
  jne    L464
  mov    0x88(%rsp),%r15d
  mov    0x90(%rsp),%ebx
L121:
  addl   $0x1,0x18(%rsp)
  mov    0x18(%rsp),%eax
  cmp    %eax,%r14d
  je     L449
L154:
  sub    $0x8,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x90(%rsp),%r11
  movsd  (%r11,%rax,8),%xmm4
  movsd  %xmm4,0x80(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%r10d
  push   %r10
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0x80(%rsp),%xmm4
  mulsd  (%rdi,%rax,8),%xmm4
  movsd  %xmm4,0x80(%rsp)
  add    $0x38,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x98(%rsp),%r10
  movsd  (%r10,%rax,8),%xmm7
  movsd  %xmm7,0x88(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%edx
  push   %rdx
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0x88(%rsp),%xmm7
  mulsd  (%r10,%rax,8),%xmm7
  movsd  0x80(%rsp),%xmm4
  subsd  %xmm7,%xmm4
  movsd  %xmm4,0xb0(%rsp)
  add    $0x38,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x90(%rsp),%rsi
  movsd  (%rsi,%rax,8),%xmm2
  movsd  %xmm2,0x80(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%ecx
  push   %rcx
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0x80(%rsp),%xmm2
  mulsd  (%r10,%rax,8),%xmm2
  movsd  %xmm2,0x80(%rsp)
  add    $0x38,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0x98(%rsp),%rsi
  movsd  (%rsi,%rax,8),%xmm5
  movsd  %xmm5,0x88(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%edx
  push   %rdx
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x0,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0x88(%rsp),%xmm5
  mulsd  (%rdi,%rax,8),%xmm5
  movapd %xmm5,%xmm0
  addsd  0x80(%rsp),%xmm0
  movsd  %xmm0,0xc8(%rsp)
  add    $0x38,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa0(%rsp),%rcx
  movsd  (%rcx,%rax,8),%xmm3
  movsd  %xmm3,0x80(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%edx
  push   %rdx
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0x80(%rsp),%xmm3
  mulsd  (%rdi,%rax,8),%xmm3
  movsd  %xmm3,0x80(%rsp)
  add    $0x38,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa8(%rsp),%rsi
  movsd  (%rsi,%rax,8),%xmm1
  movsd  %xmm1,0x88(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%edx
  push   %rdx
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0x88(%rsp),%xmm1
  mulsd  (%r10,%rax,8),%xmm1
  movsd  0x80(%rsp),%xmm3
  subsd  %xmm1,%xmm3
  movsd  0xb0(%rsp),%xmm4
  subsd  %xmm3,%xmm4
  movsd  %xmm4,0x80(%rsp)
  add    $0x38,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa0(%rsp),%rcx
  movsd  (%rcx,%rax,8),%xmm4
  movsd  %xmm4,0x88(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%ecx
  push   %rcx
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x78(%rsp),%r10
  movsd  0x88(%rsp),%xmm4
  mulsd  (%r10,%rax,8),%xmm4
  movsd  %xmm4,0x88(%rsp)
  add    $0x38,%rsp
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  push   %r13
  push   %r15
  mov    0x40(%rsp),%eax
  push   %rax
  mov    0x54(%rsp),%r9d
  mov    0x68(%rsp),%r8d
  mov    0x50(%rsp),%ecx
  mov    0x6c(%rsp),%edx
  mov    0x4c(%rsp),%esi
  mov    0x44(%rsp),%edi
  call   <Q_index@plt>
  mov    0xa8(%rsp),%rsi
  movsd  (%rsi,%rax,8),%xmm6
  movsd  %xmm6,0xb0(%rsp)
  add    $0x38,%rsp
  mov    0x130(%rsp),%ecx
  push   %rcx
  push   %r14
  push   %rbx
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r11d
  push   %r11
  push   %r13
  mov    0x60(%rsp),%r9d
  mov    0x5c(%rsp),%r8d
  mov    %r15d,%ecx
  mov    0x48(%rsp),%edx
  mov    0x44(%rsp),%esi
  mov    $0x2,%edi
  call   <prop_index@plt>
  mov    0x70(%rsp),%rdi
  movsd  0xb0(%rsp),%xmm6
  mulsd  (%rdi,%rax,8),%xmm6
  movapd %xmm6,%xmm0
  addsd  0x88(%rsp),%xmm0
  movsd  0xc8(%rsp),%xmm7
  subsd  %xmm0,%xmm7
  movsd  %xmm7,0x88(%rsp)
  add    $0x40,%rsp
  cmpl   $0x0,0x140(%rsp)
  jle    L121
  mov    $0x0,%eax
  mov    %r15d,0x88(%rsp)
  mov    %ebx,0x90(%rsp)
  mov    %eax,%ebx
  mov    0x140(%rsp),%r15d
  jmp    L464
L449:
  add    $0x1,%r13d
  cmp    %r13d,%ebx
  je     L636
L75:
  movl   $0x0,0x18(%rsp)
  jmp    L154
L636:
  mov    0x20(%rsp),%eax
  add    $0x1,%eax
  cmp    %eax,%r12d
  je     L572
  mov    %eax,0x20(%rsp)
L630:
  mov    $0x0,%r13d
  test   %ebx,%ebx
  jg     L75
  jmp    L636
L572:
  mov    %ebx,%r13d
  mov    0x1c(%rsp),%eax
  add    $0x1,%eax
  cmp    %eax,%ebp
  je     L674
  mov    %eax,0x1c(%rsp)
L585:
  movl   $0x0,0x20(%rsp)
  mov    %r13d,%ebx
  jmp    L630
L674:
  lea    0x1(%r15),%eax
  cmp    %r15d,0x20(%rsp)
  je     L381
  mov    %eax,%r15d
L725:
  movl   $0x0,0x1c(%rsp)
  jmp    L585
L381:
  mov    %ebp,%r15d
  mov    0x8(%rsp),%edx
  lea    0x1(%rdx),%eax
  cmp    %edx,0x1c(%rsp)
  je     L563
  mov    %eax,0x8(%rsp)
L821:
  mov    $0x0,%ebp
  mov    %r15d,%eax
  mov    %ebp,%r15d
  mov    %eax,%ebp
  jmp    L725
L563:
  mov    %r12d,%r13d
  mov    %ebp,%r12d
  mov    %ebx,%r15d
  addl   $0x1,0x14(%rsp)
  mov    0x14(%rsp),%eax
  cmp    %eax,%r14d
  je     L199
L285:
  movl   $0x0,0x8(%rsp)
  mov    %r15d,%eax
  mov    %r12d,%r15d
  mov    %r13d,%r12d
  mov    %eax,%r13d
  jmp    L821
L199:
  mov    %r14d,%ebp
  mov    %ebx,%r14d
  mov    %r13d,%r15d
  mov    %r12d,%r13d
L25:
  mov    0x28(%rsp),%eax
  add    $0x1,%eax
  cmp    %eax,%r15d
  je     L759
  mov    %eax,0x28(%rsp)
L86:
  movl   $0x0,0x14(%rsp)
  test   %ebp,%ebp
  jle    L25
  mov    %r13d,%r12d
  mov    %r15d,%r13d
  mov    %r14d,%r15d
  mov    %ebp,%r14d
  jmp    L285
L759:
  addl   $0x1,0x10(%rsp)
  mov    0x10(%rsp),%eax
  cmp    %eax,%r13d
  je     L8
L620:
  movl   $0x0,0x28(%rsp)
  jmp    L86
L8:
  mov    0x2c(%rsp),%esi
  lea    0x1(%rsi),%eax
  cmp    0x28(%rsp),%esi
  je     L719
  mov    %eax,0x2c(%rsp)
L20:
  movl   $0x0,0x10(%rsp)
  jmp    L620
L719:
  addl   $0x1,0xc(%rsp)
  mov    0xc(%rsp),%eax
  cmp    %eax,%r13d
  je     L105
L267:
  movl   $0x0,0x2c(%rsp)
  test   %r15d,%r15d
  jg     L20
  jmp    L719
L105:
  addl   $0x1,0x4(%rsp)
  mov    0x4(%rsp),%eax
  cmp    %eax,0x128(%rsp)
  je     L317
L616:
  movl   $0x0,0xc(%rsp)
  test   %r13d,%r13d
  jg     L267
  jmp    L105
L120:
  cmpl   $0x0,0x130(%rsp)
  jg     L592
L317:
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  mov    0x58(%rsp),%rdi
  call   <free@plt>
  mov    0x60(%rsp),%rdi
  call   <free@plt>
  mov    0x68(%rsp),%rdi
  call   <free@plt>
  add    $0xb8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void make_double_block(double *Bdouble_re, double *Bdouble_im, const double *prop_re, const double *prop_im, const int *color_weights, const int *spin_weights, const double *weights, const double *psi_re, const double *psi_im, const int Nc, const int Ns, const int Vsrc, const int Vsnk, const int Nt, const int Nw, const int Nq, const int Nsrc) {
  int iCprime, iSprime, jCprime, jSprime, kCprime, kSprime, iC, iS, jC, jS, kC, kS, x, x1, x2, t, y, wnum, n;
  double prop_prod_re;
  double prop_prod_im;
  double *Q12_re = malloc(Nt * Nc * Ns * Nc * Ns * Vsnk * Nc * Ns * Vsrc * sizeof(double));
  double *Q12_im = malloc(Nt * Nc * Ns * Nc * Ns * Vsnk * Nc * Ns * Vsrc * sizeof(double));
  double *Q01_re = malloc(Nt * Nc * Ns * Nc * Ns * Vsnk * Nc * Ns * Vsrc * sizeof(double));
  double *Q01_im = malloc(Nt * Nc * Ns * Nc * Ns * Vsnk * Nc * Ns * Vsrc * sizeof(double));
  for (t = 0; t < Nt; t++) {
    for (iCprime = 0; iCprime < Nc; iCprime++) {
      for (iSprime = 0; iSprime < Ns; iSprime++) {
        for (kCprime = 0; kCprime < Nc; kCprime++) {
          for (kSprime = 0; kSprime < Ns; kSprime++) {
            for (x1 = 0; x1 < Vsnk; x1++) {
              for (n = 0; n < Nsrc; n++) {
                for (jCprime = 0; jCprime < Nc; jCprime++) {
                  for (jSprime = 0; jSprime < Ns; jSprime++) {
                    for (x2 = 0; x2 < Vsnk; x2++) {
                      Bdouble_re[Bdouble_index(t, iCprime, iSprime, kCprime, kSprime, x1, jCprime, jSprime, x2, n, Nc, Ns, Vsnk, Nsrc)] = 0.0;
                      Bdouble_im[Bdouble_index(t, iCprime, iSprime, kCprime, kSprime, x1, jCprime, jSprime, x2, n, Nc, Ns, Vsnk, Nsrc)] = 0.0;
                    }
                  }
                }
              }
              for (jC = 0; jC < Nc; jC++) {
                for (jS = 0; jS < Ns; jS++) {
                  for (y = 0; y < Vsrc; y++) {
                    Q12_re[Q_index(t, iCprime, iSprime, kCprime, kSprime, x1, jC, jS, y, Nc, Ns, Vsrc, Vsnk)] = 0.0;
                    Q12_im[Q_index(t, iCprime, iSprime, kCprime, kSprime, x1, jC, jS, y, Nc, Ns, Vsrc, Vsnk)] = 0.0;
                    Q01_re[Q_index(t, iCprime, iSprime, kCprime, kSprime, x1, jC, jS, y, Nc, Ns, Vsrc, Vsnk)] = 0.0;
                    Q01_im[Q_index(t, iCprime, iSprime, kCprime, kSprime, x1, jC, jS, y, Nc, Ns, Vsrc, Vsnk)] = 0.0;
                  }
                }
              }
            }
          }
        }
      }
    }
  }
  for (wnum = 0; wnum < Nw; wnum++) {
    iC = color_weights[index_2d(wnum, 0, Nq)];
    iS = spin_weights[index_2d(wnum, 0, Nq)];
    jC = color_weights[index_2d(wnum, 1, Nq)];
    jS = spin_weights[index_2d(wnum, 1, Nq)];
    kC = color_weights[index_2d(wnum, 2, Nq)];
    kS = spin_weights[index_2d(wnum, 2, Nq)];
    for (t = 0; t < Nt; t++) {
      for (jCprime = 0; jCprime < Nc; jCprime++) {
        for (jSprime = 0; jSprime < Ns; jSprime++) {
          for (kCprime = 0; kCprime < Nc; kCprime++) {
            for (kSprime = 0; kSprime < Ns; kSprime++) {
              for (y = 0; y < Vsrc; y++) {
                for (x = 0; x < Vsnk; x++) {
                  Q12_re[Q_index(t, jCprime, jSprime, kCprime, kSprime, x, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] += weights[wnum] * (prop_re[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_re[prop_index(2, t, kC, kS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] - prop_im[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_im[prop_index(2, t, kC, kS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)]);
                  Q12_im[Q_index(t, jCprime, jSprime, kCprime, kSprime, x, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] += weights[wnum] * (prop_re[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_im[prop_index(2, t, kC, kS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] + prop_im[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_re[prop_index(2, t, kC, kS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)]);
                  Q01_re[Q_index(t, jCprime, jSprime, kCprime, kSprime, x, kC, kS, y, Nc, Ns, Vsrc, Vsnk)] += weights[wnum] * (prop_re[prop_index(0, t, iC, iS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_re[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] - prop_im[prop_index(0, t, iC, iS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_im[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)]);
                  Q01_im[Q_index(t, jCprime, jSprime, kCprime, kSprime, x, kC, kS, y, Nc, Ns, Vsrc, Vsnk)] += weights[wnum] * (prop_re[prop_index(0, t, iC, iS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_im[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] + prop_im[prop_index(0, t, iC, iS, kCprime, kSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)] * prop_re[prop_index(1, t, jC, jS, jCprime, jSprime, y, x, Nc, Ns, Vsrc, Vsnk, Nt)]);
                }
              }
            }
          }
        }
      }
    }
  }
  for (t = 0; t < Nt; t++) {
    for (jCprime = 0; jCprime < Nc; jCprime++) {
      for (jSprime = 0; jSprime < Ns; jSprime++) {
        for (kCprime = 0; kCprime < Nc; kCprime++) {
          for (kSprime = 0; kSprime < Ns; kSprime++) {
            for (x1 = 0; x1 < Vsnk; x1++) {
              for (iC = 0; iC < Nc; iC++) {
                for (iS = 0; iS < Ns; iS++) {
                  for (iCprime = 0; iCprime < Nc; iCprime++) {
                    for (iSprime = 0; iSprime < Ns; iSprime++) {
                      for (y = 0; y < Vsrc; y++) {
                        for (x2 = 0; x2 < Vsnk; x2++) {
                          prop_prod_re = Q12_re[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_re[prop_index(0, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)] - Q12_im[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_im[prop_index(0, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)];
                          prop_prod_im = Q12_re[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_im[prop_index(0, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)] + Q12_im[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_re[prop_index(0, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)];
                          prop_prod_re -= Q01_re[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_re[prop_index(2, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)] - Q01_im[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_im[prop_index(2, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)];
                          prop_prod_im -= Q01_re[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_im[prop_index(2, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)] + Q01_im[Q_index(t, jCprime, jSprime, kCprime, kSprime, x1, iC, iS, y, Nc, Ns, Vsrc, Vsnk)] * prop_re[prop_index(2, t, iC, iS, iCprime, iSprime, y, x2, Nc, Ns, Vsrc, Vsnk, Nt)];
                          for (n = 0; n < Nsrc; n++) {
                            Bdouble_re[Bdouble_index(t, iCprime, iSprime, kCprime, kSprime, x1, jCprime, jSprime, x2, n, Nc, Ns, Vsnk, Nsrc)] += psi_re[index_2d(y, n, Nsrc)] * prop_prod_re - psi_im[index_2d(y, n, Nsrc)] * prop_prod_im;
                            Bdouble_im[Bdouble_index(t, iCprime, iSprime, kCprime, kSprime, x1, jCprime, jSprime, x2, n, Nc, Ns, Vsnk, Nsrc)] += psi_re[index_2d(y, n, Nsrc)] * prop_prod_im + psi_im[index_2d(y, n, Nsrc)] * prop_prod_re;
                          }
                        }
                      }
                    }
                  }
                }
              }
            }
          }
        }
      }
    }
  }
  free(Q12_re);
  free(Q12_im);
  free(Q01_re);
  free(Q01_im);
}
```