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

```asm
<color_feature_standard_deviation>:
  endbr64
  push   %r15
  mov    %edx,%r15d
  push   %r14
  mov    %rdi,%r14
  push   %r13
  mov    %esi,%r13d
  lea    0x0(,%rcx,8),%esi
  push   %r12
  mov    %esi,%edi
  push   %rbp
  push   %rbx
  mov    %r9,%rbx
  sub    $0x18,%rsp
  mov    %esi,0x4(%rsp)
  mov    0x50(%rsp),%rbp
  mov    %ecx,(%rsp)
  mov    %r8,0x8(%rsp)
  call   <malloc@plt>
  mov    0x4(%rsp),%esi
  test   %rax,%rax
  mov    %rax,%r12
  je     L95
L12:
  xor    %esi,%esi
  mov    $0x8,%edx
  mov    %r12,%rdi
  call   <memset@plt>
  test   %r15d,%r15d
  je     L109
  test   %r13d,%r13d
  je     L109
  lea    0x0(%r13,%r13,2),%r11d
  lea    -0x1(%r15),%eax
  movsd  D15(%rip),%xmm3
  movsd  D21(%rip),%xmm2
  movsd  D0(%rip),%xmm1
  lea    0x8(%r14,%rax,8),%r15
  mov    %r11d,%r9d
  xor    %r10d,%r10d
  nopw   0x0(%rax,%rax,1)
L106:
  mov    (%r14),%rcx
  lea    (%r10,%r10,2),%edx
  nopw   0x0(%rax,%rax,1)
L28:
  movslq (%rcx),%rsi
  test   %esi,%esi
  je     L41
  lea    0x1(%rdx),%eax
  pxor   %xmm0,%xmm0
  movzbl (%rbx,%rax,1),%eax
  cvtsi2sd %eax,%xmm0
  mov    %edx,%eax
  movzbl (%rbx,%rax,1),%eax
  mulsd  %xmm3,%xmm0
  cvttsd2si %xmm0,%edi
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  %xmm2,%xmm0
  cvttsd2si %xmm0,%eax
  pxor   %xmm0,%xmm0
  add    %edi,%eax
  lea    0x2(%rdx),%edi
  movzbl (%rbx,%rdi,1),%edi
  cvtsi2sd %edi,%xmm0
  mulsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%edi
  pxor   %xmm0,%xmm0
  add    %edi,%eax
  movzbl %al,%eax
  cvtsi2sd %eax,%xmm0
  subsd  -0x8(%rbp,%rsi,8),%xmm0
  lea    -0x8(%r12,%rsi,8),%rax
  mulsd  %xmm0,%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
L41:
  add    $0x3,%edx
  add    $0x4,%rcx
  cmp    %edx,%r9d
  jne    L28
  add    $0x8,%r14
  add    %r11d,%r9d
  add    %r13d,%r10d
  cmp    %r14,%r15
  jne    L106
L109:
  mov    (%rsp),%eax
  test   %eax,%eax
  je     L40
  mov    %eax,%ebp
  mov    0x8(%rsp),%r13
  mov    %r12,%rbx
  pxor   %xmm2,%xmm2
  lea    (%r12,%rbp,8),%r14
  mov    0x58(%rsp),%rbp
  sub    %r12,%rbp
  nopl   0x0(%rax)
L54:
  pxor   %xmm1,%xmm1
  movsd  (%rbx),%xmm0
  cvtsi2sdl 0x0(%r13),%xmm1
  divsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm2
  ja     L116
  sqrtsd %xmm0,%xmm0
  add    $0x4,%r13
  movsd  %xmm0,0x0(%rbp,%rbx,1)
  add    $0x8,%rbx
  cmp    %rbx,%r14
  jne    L54
L40:
  mov    %r12,%rdi
  call   <free@plt>
  add    $0x18,%rsp
  mov    $0x1,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L95:
  lea    D16(%rip),%rdi
  call   <printf@plt>
  mov    $0xffffffff,%edi
  call   <exit@plt>
  jmp    L12
L116:
  call   <sqrt@plt>
  add    $0x4,%r13
  pxor   %xmm2,%xmm2
  movsd  %xmm0,0x0(%rbp,%rbx,1)
  add    $0x8,%rbx
  cmp    %r14,%rbx
  jne    L54
  jmp    L40
```

tool_calls: name: read_data
arguments: label: D15
type_guessed: f64
name: read_data
arguments: label: D21
type_guessed: f64
name: read_data
arguments: label: D0
type_guessed: f64
name: read_data
arguments: label: D16
type_guessed: string

tool: D15: 0.7152

tool: D21: 0.2126

tool: D0: 0.0722

tool: D16: color_feature_standard_deviation:Cannot allocate %d bytes for memory.


assistant: ```c
int color_feature_standard_deviation(int **obj_id, int width, int height, int n_object, int *area, unsigned char *image, double *mean, double *color_std_dev) {
  unsigned int i = 0, j = 0;
  double *sum = NULL;
  sum = (double *)malloc(n_object * sizeof(double));
  if (sum == NULL) {
    printf("color_feature_standard_deviation:Cannot allocate %d bytes for memory.\n", (n_object * sizeof(double)));
    exit(-1);
  }
  memset(sum, 0, sizeof(n_object * sizeof(double)));
  for (i = 0; i < height; i++) {
    for (j = 0; j < width; j++) {
      if (!obj_id[i][j])
        continue;
      register unsigned int idx = (j + i * width) * 3;
      unsigned char val = ((unsigned char)(0.2126 * (double)(image[idx])) + (unsigned char)(0.7152 * (double)(image[idx + 1])) + (unsigned char)(0.0722 * (double)(image[idx + 2])));
      double tmp = ((double)val - mean[obj_id[i][j] - 1]);
      sum[obj_id[i][j] - 1] += (tmp * tmp);
    }
  }
  for (i = 0; i < n_object; i++) {
    color_std_dev[i] = sqrt(sum[i] / area[i]);
  }
  free(sum);
  return 1;
}
```