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

```asm
<mirror_ripples>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  movdqa D9(%rip),%xmm0
  mov    %esi,0x24(%rsp)
  lea    0x30(%rsp),%rcx
  mov    %fs:0x28,%rax
  mov    %rax,0x48(%rsp)
  xor    %eax,%eax
  lea    0x0(,%rdi,4),%eax
  movl   $0x0,0x1c(%rsp)
  mov    %eax,0x20(%rsp)
  mov    %rcx,0x28(%rsp)
  movaps %xmm0,0x30(%rsp)
  test   %eax,%eax
  jle    L23
  nopl   0x0(%rax,%rax,1)
L24:
  xor    %edi,%edi
  xor    %ebx,%ebx
  call   <fill@plt>
  movzbl 0x1c(%rsp),%eax
  mov    %al,0x1b(%rsp)
  mov    0x28(%rsp),%rax
  mov    %rax,0x10(%rsp)
L36:
  mov    $0x7,%r13d
  xor    %r14d,%r14d
  sub    %ebx,%r13d
L5:
  mov    0x10(%rsp),%rax
  movzbl 0x1b(%rsp),%edi
  mov    $0x7,%r15d
  sub    %r14d,%r15d
  add    (%rax,%r14,4),%dil
  movzbl %dil,%edi
  call   <totty_sin@plt>
  mov    %r14d,%esi
  mov    %ebx,%edi
  add    $0xc4,%eax
  movslq %eax,%rbp
  sar    $0x1f,%eax
  imul   $0x5397829d,%rbp,%rbp
  sar    $0x24,%rbp
  sub    %eax,%ebp
  movzbl %bpl,%ebp
  mov    %ebp,%edx
  call   <setvoxel@plt>
  mov    %r14d,%esi
  mov    %ebp,%edx
  mov    %r13d,%edi
  call   <setvoxel@plt>
  mov    %ebp,%edx
  mov    %r15d,%esi
  mov    %ebx,%edi
  call   <setvoxel@plt>
  mov    %ebp,%edx
  mov    %r15d,%esi
  mov    %r13d,%edi
  call   <setvoxel@plt>
  mov    $0x7,%edx
  mov    %r14d,%esi
  mov    %ebx,%edi
  sub    %ebp,%edx
  mov    %edx,0xc(%rsp)
  call   <setvoxel@plt>
  mov    0xc(%rsp),%edx
  mov    %r14d,%esi
  mov    %r13d,%edi
  add    $0x1,%r14
  call   <setvoxel@plt>
  mov    0xc(%rsp),%edx
  mov    %r15d,%esi
  mov    %ebx,%edi
  call   <setvoxel@plt>
  mov    0xc(%rsp),%edx
  mov    %r15d,%esi
  mov    %r13d,%edi
  call   <setvoxel@plt>
  cmp    $0x4,%r14
  jne    L5
  add    $0x1,%ebx
  addq   $0x1,0x10(%rsp)
  cmp    $0x4,%ebx
  jne    L36
  mov    0x24(%rsp),%edi
  call   <delay_ms@plt>
  addl   $0x4,0x1c(%rsp)
  mov    0x1c(%rsp),%eax
  cmp    %eax,0x20(%rsp)
  jg     L24
L23:
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L31
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L31:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_tag: D9
data_type: byte[16]

tool: {"D9": [0x31,0x2b,0x26,0x23,0x2b,0x23,0x1d,0x1a,0x26,0x1d,0x15,0x10,0x23,0x19,0x10,0x07]}

assistant: ```c
void mirror_ripples(int iterations, int delay) {
  unsigned char sqrt_LUT[] = {49, 43, 38, 35, 43, 35, 29, 26, 38, 29, 21, 16, 35, 25, 16, 7};
  unsigned char x, y, height, distance;
  int i;
  for (i = 0; i < iterations * 4; i += 4) {
    fill(0x00);
    for (x = 0; x < 4; x++)
      for (y = 0; y < 4; y++) {
        distance = sqrt_LUT[x + y * 4];
        height = (196 + totty_sin(distance + i)) / 49;
        setvoxel(x, y, height);
        setvoxel(7 - x, y, height);
        setvoxel(x, 7 - y, height);
        setvoxel(7 - x, 7 - y, height);
        setvoxel(x, y, 7 - height);
        setvoxel(7 - x, y, 7 - height);
        setvoxel(x, 7 - y, 7 - height);
        setvoxel(7 - x, 7 - y, 7 - height);
      }
    delay_ms(delay);
  }
}
```