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

```asm
<make_gaussian>:
  endbr64
  push   %r15
  pxor   %xmm7,%xmm7
  push   %r14
  cvtss2sd %xmm0,%xmm7
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  movsd  D15(%rip),%xmm1
  mov    %edi,0x2c(%rsp)
  comisd %xmm7,%xmm1
  mov    %rsi,0x38(%rsp)
  ja     L131
L33:
  mov    0x2c(%rsp),%esi
  mov    %esi,%edx
  shr    $0x1f,%edx
  lea    (%rsi,%rdx,1),%eax
  and    $0x1,%eax
  sub    %edx,%eax
  cmp    $0x1,%eax
  jne    L101
  lea    -0x3(%rsi),%eax
  cmp    $0x7c,%eax
  ja     L101
  mov    %esi,%edx
  mov    %esi,%eax
  mov    0x38(%rsp),%rsi
  pxor   %xmm7,%xmm7
  sar    %edx
  movl   $0x0,0x24(%rsp)
  xor    %ebp,%ebp
  imul   %edx,%eax
  movslq %edx,%rcx
  movss  %xmm7,0xc(%rsp)
  pxor   %xmm7,%xmm7
  cvtss2sd %xmm0,%xmm7
  movsd  %xmm7,0x18(%rsp)
  cltq
  add    %rcx,%rax
  lea    (%rsi,%rax,4),%rax
  mov    %rax,0x30(%rsp)
  lea    0x1(%rdx),%eax
  mov    %eax,0x28(%rsp)
  nopl   0x0(%rax,%rax,1)
L93:
  mov    %ebp,%eax
  mov    0x30(%rsp),%rcx
  xor    %r14d,%r14d
  xor    %r15d,%r15d
  imul   %ebp,%eax
  mov    %eax,0x20(%rsp)
  movslq 0x28(%rsp),%rax
  mov    %rax,0x10(%rsp)
  movslq 0x24(%rsp),%rax
  shl    $0x2,%rax
  lea    (%rcx,%rax,1),%r13
  sub    %rax,%rcx
  mov    %rcx,%r12
  jmp    L38
  nopl   0x0(%rax,%rax,1)
L150:
  movss  0xc(%rsp),%xmm4
  addss  %xmm0,%xmm4
  movss  %xmm4,0xc(%rsp)
L108:
  movss  %xmm0,0x0(%r13,%r15,4)
  movss  %xmm0,(%r12,%r15,4)
  add    $0x1,%r15
  movss  %xmm0,(%r12,%r14,1)
  movss  %xmm0,0x0(%r13,%r14,1)
  sub    $0x4,%r14
  cmp    0x10(%rsp),%r15
  je     L128
L38:
  mov    %r15d,%edi
  pxor   %xmm0,%xmm0
  pxor   %xmm2,%xmm2
  mov    %r15d,%ebx
  imul   %r15d,%edi
  add    0x20(%rsp),%edi
  cvtsi2sd %edi,%xmm0
  ucomisd %xmm0,%xmm2
  ja     L73
  sqrtsd %xmm0,%xmm0
L15:
  xorpd  D9(%rip),%xmm0
  divsd  0x18(%rsp),%xmm0
  call   <exp@plt>
  or     %ebp,%ebx
  cvtsd2ss %xmm0,%xmm0
  je     L150
  test   %r15d,%r15d
  je     L42
  test   %ebp,%ebp
  je     L42
  movss  D24(%rip),%xmm1
  mulss  %xmm0,%xmm1
  addss  0xc(%rsp),%xmm1
  movss  %xmm1,0xc(%rsp)
  jmp    L108
  cs nopw 0x0(%rax,%rax,1)
L42:
  movaps %xmm0,%xmm1
  addss  %xmm0,%xmm1
  addss  0xc(%rsp),%xmm1
  movss  %xmm1,0xc(%rsp)
  jmp    L108
  nopl   0x0(%rax,%rax,1)
L128:
  mov    0x2c(%rsp),%esi
  add    $0x1,%ebp
  add    %esi,0x24(%rsp)
  cmp    0x28(%rsp),%ebp
  jne    L93
  mov    %esi,%ecx
  movss  0xc(%rsp),%xmm1
  imul   %esi,%ecx
  mov    0x38(%rsp),%rsi
  shufps $0x0,%xmm1,%xmm1
  mov    %rsi,%rax
  mov    %ecx,%edx
  shr    $0x2,%edx
  sub    $0x1,%edx
  shl    $0x4,%rdx
  lea    0x10(%rsi,%rdx,1),%rdx
  nopl   (%rax)
L47:
  movups (%rax),%xmm0
  add    $0x10,%rax
  divps  %xmm1,%xmm0
  movups %xmm0,-0x10(%rax)
  cmp    %rdx,%rax
  jne    L47
  mov    %ecx,%eax
  and    $0xfffffffc,%eax
  test   $0x3,%cl
  je     L36
  mov    0x38(%rsp),%rbx
  mov    %eax,%edx
  movss  0xc(%rsp),%xmm7
  shl    $0x2,%rdx
  lea    (%rbx,%rdx,1),%rsi
  movss  (%rsi),%xmm0
  divss  %xmm7,%xmm0
  movss  %xmm0,(%rsi)
  lea    0x1(%rax),%esi
  cmp    %esi,%ecx
  jle    L36
  lea    0x4(%rbx,%rdx,1),%rsi
  add    $0x2,%eax
  movss  (%rsi),%xmm0
  divss  %xmm7,%xmm0
  movss  %xmm0,(%rsi)
  cmp    %eax,%ecx
  jle    L36
  lea    0x8(%rbx,%rdx,1),%rax
  movss  (%rax),%xmm0
  divss  %xmm7,%xmm0
  movss  %xmm0,(%rax)
L36:
  xor    %eax,%eax
L187:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L131:
  lea    D27(%rip),%rdi
  xor    %eax,%eax
  call   <err_printf@plt>
  movss  D16(%rip),%xmm0
  jmp    L33
L73:
  call   <sqrt@plt>
  jmp    L15
L101:
  lea    D36(%rip),%rdi
  xor    %eax,%eax
  call   <err_printf@plt>
  mov    $0xffffffff,%eax
  jmp    L187
```

tool_calls: name: unpack_data
arguments: data_label: D15
type: f64
name: unpack_data
arguments: data_label: D9
type: byte[16]
name: unpack_data
arguments: data_label: D24
type: byte[4]
name: unpack_data
arguments: data_label: D27
type: string
name: unpack_data
arguments: data_label: D16
type: byte[4]
name: unpack_data
arguments: data_label: D36
type: string

tool: D15: 0.01

tool: {"D9": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: {"D24": [0x00,0x00,0x80,0x40]}

tool: D27: make_gaussian: clipping sigma to 0.01


tool: {"D16": [0x0a,0xd7,0x23,0x3c]}

tool: D36: make_gaussian: bad size %d


assistant: ```c
int make_gaussian(float sigma, int size, float *kern) {
  int mid, all;
  float *mp;
  float sum, v;
  int x, y;
  if (sigma < 0.01) {
    err_printf("make_gaussian: clipping sigma to 0.01\n");
    sigma = 0.01;
  }
  if (size % 2 != 1 || size < 3 || size > 127) {
    err_printf("make_gaussian: bad size %d\n");
    return -1;
  }
  mid = size / 2;
  mp = kern + mid * size + mid;
  sum = 0;
  for (y = 0; y < mid + 1; y++)
    for (x = 0; x < mid + 1; x++) {
      v = exp(-sqrt(1.0 * (((x) * (x)) + ((y) * (y)))) / sigma);
      if (x == 0 && y == 0)
        sum += v;
      else if (x == 0 || y == 0)
        sum += 2 * v;
      else
        sum += 4 * v;
      *(mp + x + size * y) = v;
      *(mp + x - size * y) = v;
      *(mp - x - size * y) = v;
      *(mp - x + size * y) = v;
    }
  all = size * size;
  for (y = 0; y < all; y++)
    kern[y] /= sum;
  return 0;
}
```