user: What is the c source code of the assembly code below:

```asm
<sumGAMMAPROT>:
  endbr64
  sub    $0x8,%rsp
  mov    %edi,%r10d
  mov    %rsi,%rax
  mov    %rdx,%rsi
  mov    0x10(%rsp),%rdi
  mov    0x18(%rsp),%edx
  mov    0x20(%rsp),%r11d
  cmp    $0x1,%r10d
  je     L87
  cmp    $0x2,%r10d
  je     L10
  test   %r10d,%r10d
  je     L132
  mov    $0x0,%edi
  mov    $0x0,%eax
  call   <assert@plt>
L46:
  add    $0x8,%rsp
  ret
L10:
  cmp    %r11d,%edx
  jge    L46
  lea    (%rdx,%rdx,4),%r9d
  shl    $0x4,%r9d
  movslq %r9d,%rdx
  shl    $0x3,%rdx
  lea    (%rsi,%rdx,1),%r8
  lea    (%rcx,%rdx,1),%rsi
  lea    0x280(%rax,%rdx,1),%rdi
  lea    (%r11,%r11,4),%r10d
  shl    $0x4,%r10d
  jmp    L118
L87:
  cmp    %r11d,%edx
  jge    L46
  movslq %edx,%rsi
  add    %rsi,%r9
  lea    (%rdx,%rdx,4),%r10d
  shl    $0x4,%r10d
  movslq %r10d,%rdx
  shl    $0x3,%rdx
  lea    (%rcx,%rdx,1),%rsi
  lea    0x280(%rax,%rdx,1),%rdi
  lea    (%r11,%r11,4),%r11d
  shl    $0x4,%r11d
  jmp    L71
L132:
  cmp    %r11d,%edx
  jge    L46
  movslq %edx,%rsi
  lea    (%rdx,%rdx,4),%edx
  shl    $0x4,%edx
  movslq %edx,%rdx
  lea    0x280(%rax,%rdx,8),%r10
  jmp    L107
L45:
  add    $0x1,%rsi
  add    $0x280,%r10
  cmp    %esi,%r11d
  jle    L46
L107:
  movsbl (%r9,%rsi,1),%eax
  lea    (%rax,%rax,4),%eax
  shl    $0x2,%eax
  cltq
  lea    (%r8,%rax,8),%rcx
  movsbl (%rdi,%rsi,1),%eax
  lea    (%rax,%rax,4),%eax
  shl    $0x2,%eax
  cltq
  lea    (%r8,%rax,8),%rdx
  lea    -0x280(%r10),%rax
L173:
  movsd  (%rcx),%xmm0
  mulsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x8(%rcx),%xmm0
  mulsd  0x8(%rdx),%xmm0
  movsd  %xmm0,0x8(%rax)
  movsd  0x10(%rcx),%xmm0
  mulsd  0x10(%rdx),%xmm0
  movsd  %xmm0,0x10(%rax)
  movsd  0x18(%rcx),%xmm0
  mulsd  0x18(%rdx),%xmm0
  movsd  %xmm0,0x18(%rax)
  movsd  0x20(%rcx),%xmm0
  mulsd  0x20(%rdx),%xmm0
  movsd  %xmm0,0x20(%rax)
  movsd  0x28(%rcx),%xmm0
  mulsd  0x28(%rdx),%xmm0
  movsd  %xmm0,0x28(%rax)
  movsd  0x30(%rcx),%xmm0
  mulsd  0x30(%rdx),%xmm0
  movsd  %xmm0,0x30(%rax)
  movsd  0x38(%rcx),%xmm0
  mulsd  0x38(%rdx),%xmm0
  movsd  %xmm0,0x38(%rax)
  movsd  0x40(%rcx),%xmm0
  mulsd  0x40(%rdx),%xmm0
  movsd  %xmm0,0x40(%rax)
  movsd  0x48(%rcx),%xmm0
  mulsd  0x48(%rdx),%xmm0
  movsd  %xmm0,0x48(%rax)
  movsd  0x50(%rcx),%xmm0
  mulsd  0x50(%rdx),%xmm0
  movsd  %xmm0,0x50(%rax)
  movsd  0x58(%rcx),%xmm0
  mulsd  0x58(%rdx),%xmm0
  movsd  %xmm0,0x58(%rax)
  movsd  0x60(%rcx),%xmm0
  mulsd  0x60(%rdx),%xmm0
  movsd  %xmm0,0x60(%rax)
  movsd  0x68(%rcx),%xmm0
  mulsd  0x68(%rdx),%xmm0
  movsd  %xmm0,0x68(%rax)
  movsd  0x70(%rcx),%xmm0
  mulsd  0x70(%rdx),%xmm0
  movsd  %xmm0,0x70(%rax)
  movsd  0x78(%rcx),%xmm0
  mulsd  0x78(%rdx),%xmm0
  movsd  %xmm0,0x78(%rax)
  movsd  0x80(%rcx),%xmm0
  mulsd  0x80(%rdx),%xmm0
  movsd  %xmm0,0x80(%rax)
  movsd  0x88(%rcx),%xmm0
  mulsd  0x88(%rdx),%xmm0
  movsd  %xmm0,0x88(%rax)
  movsd  0x90(%rcx),%xmm0
  mulsd  0x90(%rdx),%xmm0
  movsd  %xmm0,0x90(%rax)
  add    $0xa0,%rax
  movsd  0x98(%rcx),%xmm0
  mulsd  0x98(%rdx),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %r10,%rax
  jne    L173
  jmp    L45
L44:
  add    $0x1,%r9
  add    $0x280,%rsi
  add    $0x50,%r10d
  add    $0x280,%rdi
  cmp    %r11d,%r10d
  je     L46
L71:
  movsbl (%r9),%eax
  lea    (%rax,%rax,4),%eax
  shl    $0x2,%eax
  cltq
  lea    (%r8,%rax,8),%rcx
  lea    -0x280(%rdi),%rax
  mov    %rsi,%rdx
L67:
  movsd  (%rcx),%xmm0
  mulsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x8(%rcx),%xmm0
  mulsd  0x8(%rdx),%xmm0
  movsd  %xmm0,0x8(%rax)
  movsd  0x10(%rcx),%xmm0
  mulsd  0x10(%rdx),%xmm0
  movsd  %xmm0,0x10(%rax)
  movsd  0x18(%rcx),%xmm0
  mulsd  0x18(%rdx),%xmm0
  movsd  %xmm0,0x18(%rax)
  movsd  0x20(%rcx),%xmm0
  mulsd  0x20(%rdx),%xmm0
  movsd  %xmm0,0x20(%rax)
  movsd  0x28(%rcx),%xmm0
  mulsd  0x28(%rdx),%xmm0
  movsd  %xmm0,0x28(%rax)
  movsd  0x30(%rcx),%xmm0
  mulsd  0x30(%rdx),%xmm0
  movsd  %xmm0,0x30(%rax)
  movsd  0x38(%rcx),%xmm0
  mulsd  0x38(%rdx),%xmm0
  movsd  %xmm0,0x38(%rax)
  movsd  0x40(%rcx),%xmm0
  mulsd  0x40(%rdx),%xmm0
  movsd  %xmm0,0x40(%rax)
  movsd  0x48(%rcx),%xmm0
  mulsd  0x48(%rdx),%xmm0
  movsd  %xmm0,0x48(%rax)
  movsd  0x50(%rcx),%xmm0
  mulsd  0x50(%rdx),%xmm0
  movsd  %xmm0,0x50(%rax)
  movsd  0x58(%rcx),%xmm0
  mulsd  0x58(%rdx),%xmm0
  movsd  %xmm0,0x58(%rax)
  movsd  0x60(%rcx),%xmm0
  mulsd  0x60(%rdx),%xmm0
  movsd  %xmm0,0x60(%rax)
  movsd  0x68(%rcx),%xmm0
  mulsd  0x68(%rdx),%xmm0
  movsd  %xmm0,0x68(%rax)
  movsd  0x70(%rcx),%xmm0
  mulsd  0x70(%rdx),%xmm0
  movsd  %xmm0,0x70(%rax)
  movsd  0x78(%rcx),%xmm0
  mulsd  0x78(%rdx),%xmm0
  movsd  %xmm0,0x78(%rax)
  movsd  0x80(%rcx),%xmm0
  mulsd  0x80(%rdx),%xmm0
  movsd  %xmm0,0x80(%rax)
  movsd  0x88(%rcx),%xmm0
  mulsd  0x88(%rdx),%xmm0
  movsd  %xmm0,0x88(%rax)
  movsd  0x90(%rcx),%xmm0
  mulsd  0x90(%rdx),%xmm0
  movsd  %xmm0,0x90(%rax)
  add    $0xa0,%rax
  movsd  0x98(%rcx),%xmm0
  mulsd  0x98(%rdx),%xmm0
  movsd  %xmm0,-0x8(%rax)
  add    $0xa0,%rdx
  cmp    %rdi,%rax
  jne    L67
  jmp    L44
L178:
  add    $0x280,%r8
  add    $0x280,%rsi
  add    $0x50,%r9d
  add    $0x280,%rdi
  cmp    %r10d,%r9d
  je     L46
L118:
  lea    -0x280(%rdi),%rax
  mov    %rsi,%rcx
  mov    %r8,%rdx
L68:
  movsd  (%rdx),%xmm0
  mulsd  (%rcx),%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x8(%rdx),%xmm0
  mulsd  0x8(%rcx),%xmm0
  movsd  %xmm0,0x8(%rax)
  movsd  0x10(%rdx),%xmm0
  mulsd  0x10(%rcx),%xmm0
  movsd  %xmm0,0x10(%rax)
  movsd  0x18(%rdx),%xmm0
  mulsd  0x18(%rcx),%xmm0
  movsd  %xmm0,0x18(%rax)
  movsd  0x20(%rdx),%xmm0
  mulsd  0x20(%rcx),%xmm0
  movsd  %xmm0,0x20(%rax)
  movsd  0x28(%rdx),%xmm0
  mulsd  0x28(%rcx),%xmm0
  movsd  %xmm0,0x28(%rax)
  movsd  0x30(%rdx),%xmm0
  mulsd  0x30(%rcx),%xmm0
  movsd  %xmm0,0x30(%rax)
  movsd  0x38(%rdx),%xmm0
  mulsd  0x38(%rcx),%xmm0
  movsd  %xmm0,0x38(%rax)
  movsd  0x40(%rdx),%xmm0
  mulsd  0x40(%rcx),%xmm0
  movsd  %xmm0,0x40(%rax)
  movsd  0x48(%rdx),%xmm0
  mulsd  0x48(%rcx),%xmm0
  movsd  %xmm0,0x48(%rax)
  movsd  0x50(%rdx),%xmm0
  mulsd  0x50(%rcx),%xmm0
  movsd  %xmm0,0x50(%rax)
  movsd  0x58(%rdx),%xmm0
  mulsd  0x58(%rcx),%xmm0
  movsd  %xmm0,0x58(%rax)
  movsd  0x60(%rdx),%xmm0
  mulsd  0x60(%rcx),%xmm0
  movsd  %xmm0,0x60(%rax)
  movsd  0x68(%rdx),%xmm0
  mulsd  0x68(%rcx),%xmm0
  movsd  %xmm0,0x68(%rax)
  movsd  0x70(%rdx),%xmm0
  mulsd  0x70(%rcx),%xmm0
  movsd  %xmm0,0x70(%rax)
  movsd  0x78(%rdx),%xmm0
  mulsd  0x78(%rcx),%xmm0
  movsd  %xmm0,0x78(%rax)
  movsd  0x80(%rdx),%xmm0
  mulsd  0x80(%rcx),%xmm0
  movsd  %xmm0,0x80(%rax)
  movsd  0x88(%rdx),%xmm0
  mulsd  0x88(%rcx),%xmm0
  movsd  %xmm0,0x88(%rax)
  movsd  0x90(%rdx),%xmm0
  mulsd  0x90(%rcx),%xmm0
  movsd  %xmm0,0x90(%rax)
  add    $0xa0,%rax
  movsd  0x98(%rdx),%xmm0
  mulsd  0x98(%rcx),%xmm0
  movsd  %xmm0,-0x8(%rax)
  add    $0xa0,%rdx
  add    $0xa0,%rcx
  cmp    %rdi,%rax
  jne    L68
  jmp    L178
```

assistant: ```c
void sumGAMMAPROT(int tipCase, double *sumtable, double *x1, double *x2, double *tipVector, char *tipX1, char *tipX2, int lower, int n) {
  int i, l;
  double *left, *right, *sum;
  switch (tipCase) {
  case 0:
    for (i = lower; i < n; i++) {
      left = &(tipVector[20 * tipX1[i]]);
      right = &(tipVector[20 * tipX2[i]]);
      sum = &(sumtable[i * 80]);
      for (l = 0; l < 4; l++) {
        *sum++ = left[0] * right[0];
        *sum++ = left[1] * right[1];
        *sum++ = left[2] * right[2];
        *sum++ = left[3] * right[3];
        *sum++ = left[4] * right[4];
        *sum++ = left[5] * right[5];
        *sum++ = left[6] * right[6];
        *sum++ = left[7] * right[7];
        *sum++ = left[8] * right[8];
        *sum++ = left[9] * right[9];
        *sum++ = left[10] * right[10];
        *sum++ = left[11] * right[11];
        *sum++ = left[12] * right[12];
        *sum++ = left[13] * right[13];
        *sum++ = left[14] * right[14];
        *sum++ = left[15] * right[15];
        *sum++ = left[16] * right[16];
        *sum++ = left[17] * right[17];
        *sum++ = left[18] * right[18];
        *sum++ = left[19] * right[19];
      }
    }
    break;
  case 1:
    for (i = lower; i < n; i++) {
      left = &(tipVector[20 * tipX1[i]]);
      sum = &(sumtable[i * 80]);
      for (l = 0; l < 4; l++) {
        right = &(x2[80 * i + l * 20]);
        *sum++ = left[0] * right[0];
        *sum++ = left[1] * right[1];
        *sum++ = left[2] * right[2];
        *sum++ = left[3] * right[3];
        *sum++ = left[4] * right[4];
        *sum++ = left[5] * right[5];
        *sum++ = left[6] * right[6];
        *sum++ = left[7] * right[7];
        *sum++ = left[8] * right[8];
        *sum++ = left[9] * right[9];
        *sum++ = left[10] * right[10];
        *sum++ = left[11] * right[11];
        *sum++ = left[12] * right[12];
        *sum++ = left[13] * right[13];
        *sum++ = left[14] * right[14];
        *sum++ = left[15] * right[15];
        *sum++ = left[16] * right[16];
        *sum++ = left[17] * right[17];
        *sum++ = left[18] * right[18];
        *sum++ = left[19] * right[19];
      }
    }
    break;
  case 2:
    for (i = lower; i < n; i++) {
      sum = &sumtable[i * 80];
      for (l = 0; l < 4; l++) {
        left = &(x1[80 * i + l * 20]);
        right = &(x2[80 * i + l * 20]);
        *sum++ = left[0] * right[0];
        *sum++ = left[1] * right[1];
        *sum++ = left[2] * right[2];
        *sum++ = left[3] * right[3];
        *sum++ = left[4] * right[4];
        *sum++ = left[5] * right[5];
        *sum++ = left[6] * right[6];
        *sum++ = left[7] * right[7];
        *sum++ = left[8] * right[8];
        *sum++ = left[9] * right[9];
        *sum++ = left[10] * right[10];
        *sum++ = left[11] * right[11];
        *sum++ = left[12] * right[12];
        *sum++ = left[13] * right[13];
        *sum++ = left[14] * right[14];
        *sum++ = left[15] * right[15];
        *sum++ = left[16] * right[16];
        *sum++ = left[17] * right[17];
        *sum++ = left[18] * right[18];
        *sum++ = left[19] * right[19];
      }
    }
    break;
  default:
    assert(0);
  }
}
```