user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<vdeFun>:
  endbr64
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x88,%rsp
  mov    (%rdi),%rax
  test   %rax,%rax
  je     L76
  mov    (%rsi),%rdx
  movsd  0x10(%rax),%xmm0
  test   %rdx,%rdx
  je     L158
  movsd  %xmm0,(%rdx)
  movsd  0x18(%rax),%xmm6
  movsd  %xmm6,0x38(%rsp)
  movsd  %xmm6,0x8(%rdx)
L8:
  movsd  0x8(%rax),%xmm0
  lea    0x78(%rsp),%rdi
  lea    0x70(%rsp),%rsi
  call   <sincos@plt>
  movsd  D107(%rip),%xmm7
  movsd  0x78(%rsp),%xmm1
  movsd  0x70(%rsp),%xmm2
  movsd  %xmm7,(%rsp)
  mulsd  %xmm1,%xmm7
  movsd  %xmm7,0x18(%rsp)
L142:
  movsd  0x38(%rsp),%xmm0
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <vde_forw_pendulum_sq@plt>
  mov    0x18(%rbp),%rax
  pxor   %xmm6,%xmm6
  movsd  0x10(%rsp),%xmm1
  movsd  %xmm0,0x20(%rsp)
  movsd  0x28(%rsp),%xmm2
  mulsd  0x18(%rsp),%xmm0
  test   %rax,%rax
  movsd  %xmm6,0x8(%rsp)
  je     L177
  movsd  (%rax),%xmm4
  movsd  %xmm4,0x8(%rsp)
L177:
  movsd  D119(%rip),%xmm3
  addsd  0x8(%rsp),%xmm0
  movsd  %xmm1,0x40(%rsp)
  movsd  %xmm2,0x10(%rsp)
  mulsd  %xmm2,%xmm3
  movsd  %xmm3,0x28(%rsp)
  mulsd  %xmm1,%xmm3
  addsd  %xmm0,%xmm3
  movapd %xmm2,%xmm0
  movsd  %xmm3,0x30(%rsp)
  call   <vde_forw_pendulum_sq@plt>
  mulsd  D324(%rip),%xmm0
  mov    (%rbx),%rax
  movsd  D236(%rip),%xmm10
  movsd  0x30(%rsp),%xmm3
  test   %rax,%rax
  movsd  0x10(%rsp),%xmm2
  movsd  0x40(%rsp),%xmm1
  subsd  %xmm0,%xmm10
  divsd  %xmm10,%xmm3
  je     L198
  movsd  %xmm3,0x10(%rax)
L198:
  movsd  (%rsp),%xmm8
  movsd  0x38(%rsp),%xmm0
  movsd  %xmm3,0x68(%rsp)
  movsd  %xmm10,0x60(%rsp)
  mulsd  %xmm2,%xmm8
  movsd  %xmm2,0x48(%rsp)
  movsd  %xmm1,0x40(%rsp)
  movapd %xmm8,%xmm6
  movsd  %xmm8,0x58(%rsp)
  mulsd  %xmm1,%xmm6
  movsd  %xmm6,0x10(%rsp)
  call   <vde_forw_pendulum_sq@plt>
  movsd  0x48(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm6
  movsd  0x8(%rsp),%xmm4
  movsd  0x40(%rsp),%xmm1
  movsd  %xmm0,0x30(%rsp)
  mulsd  %xmm0,%xmm6
  movsd  D119(%rip),%xmm0
  movsd  D301(%rip),%xmm5
  movsd  %xmm2,0x40(%rsp)
  mulsd  %xmm2,%xmm4
  movsd  %xmm1,0x50(%rsp)
  mulsd  %xmm1,%xmm0
  mulsd  %xmm1,%xmm5
  addsd  %xmm6,%xmm4
  addsd  %xmm4,%xmm0
  addsd  %xmm0,%xmm5
  movapd %xmm2,%xmm0
  movsd  %xmm5,0x48(%rsp)
  call   <vde_forw_pendulum_sq@plt>
  mulsd  D324(%rip),%xmm0
  mov    (%rbx),%rax
  movsd  D236(%rip),%xmm4
  movsd  0x48(%rsp),%xmm5
  test   %rax,%rax
  movsd  0x40(%rsp),%xmm2
  movsd  0x50(%rsp),%xmm1
  movsd  0x58(%rsp),%xmm8
  movsd  0x60(%rsp),%xmm10
  movsd  0x68(%rsp),%xmm3
  subsd  %xmm0,%xmm4
  mulsd  D308(%rip),%xmm4
  divsd  %xmm4,%xmm5
  je     L183
  movsd  %xmm5,0x18(%rax)
L183:
  mov    0x8(%rbp),%rdx
  mov    0x8(%rbx),%rax
  test   %rdx,%rdx
  je     L39
  movsd  0x10(%rdx),%xmm0
  test   %rax,%rax
  je     L227
  movsd  %xmm0,(%rax)
  movsd  0x38(%rsp),%xmm7
  movsd  0x18(%rdx),%xmm9
  addsd  %xmm7,%xmm7
  movsd  %xmm9,0x8(%rax)
  movsd  0x8(%rdx),%xmm11
  movapd %xmm11,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  %xmm7,%xmm9
L167:
  mulsd  %xmm1,%xmm11
  movsd  (%rsp),%xmm14
  movapd %xmm2,%xmm6
  movsd  D119(%rip),%xmm15
  movsd  0x28(%rsp),%xmm12
  addsd  %xmm2,%xmm6
  mulsd  %xmm0,%xmm14
  mulsd  %xmm0,%xmm12
  mulsd  0x20(%rsp),%xmm14
  movapd %xmm6,%xmm13
  mulsd  %xmm11,%xmm15
  mulsd  %xmm11,%xmm13
  divsd  %xmm10,%xmm3
  mulsd  D324(%rip),%xmm13
  mulsd  %xmm1,%xmm15
  subsd  %xmm15,%xmm12
  divsd  %xmm4,%xmm5
  movsd  %xmm3,0x38(%rsp)
  movsd  0x18(%rsp),%xmm3
  mulsd  %xmm9,%xmm3
  mulsd  0x10(%rsp),%xmm9
  addsd  %xmm3,%xmm14
  movsd  0x38(%rsp),%xmm3
  mulsd  %xmm13,%xmm3
  mulsd  D308(%rip),%xmm13
  addsd  %xmm14,%xmm12
  movsd  (%rsp),%xmm14
  mulsd  %xmm11,%xmm14
  divsd  %xmm10,%xmm12
  mulsd  0x8(%rsp),%xmm11
  mulsd  %xmm1,%xmm14
  mulsd  %xmm5,%xmm13
  subsd  %xmm3,%xmm12
  movsd  %xmm12,0x10(%rax)
  movapd %xmm8,%xmm12
  mulsd  %xmm0,%xmm12
  subsd  %xmm14,%xmm12
  mulsd  0x30(%rsp),%xmm12
  addsd  %xmm12,%xmm9
  subsd  %xmm11,%xmm9
  movsd  D119(%rip),%xmm11
  mulsd  %xmm0,%xmm11
  mulsd  D301(%rip),%xmm0
  addsd  %xmm9,%xmm11
  addsd  %xmm11,%xmm0
  divsd  %xmm4,%xmm0
  subsd  %xmm13,%xmm0
  movsd  %xmm0,0x18(%rax)
  test   %rdx,%rdx
  je     L60
  movsd  0x30(%rdx),%xmm0
  movsd  %xmm0,0x20(%rax)
  movsd  0x38(%rdx),%xmm0
  movsd  %xmm0,0x28(%rax)
  movsd  0x28(%rdx),%xmm11
  mulsd  %xmm7,%xmm0
  movapd %xmm11,%xmm9
  mulsd  %xmm1,%xmm11
  mulsd  %xmm2,%xmm9
  movapd %xmm11,%xmm13
  mulsd  %xmm6,%xmm13
  mulsd  D324(%rip),%xmm13
L34:
  movsd  D119(%rip),%xmm15
  movsd  (%rsp),%xmm14
  movsd  0x28(%rsp),%xmm12
  movsd  0x18(%rsp),%xmm3
  mulsd  %xmm11,%xmm15
  mulsd  %xmm9,%xmm14
  mulsd  %xmm0,%xmm3
  mulsd  %xmm9,%xmm12
  mulsd  0x20(%rsp),%xmm14
  mulsd  %xmm1,%xmm15
  mulsd  0x10(%rsp),%xmm0
  addsd  %xmm3,%xmm14
  movsd  0x38(%rsp),%xmm3
  subsd  %xmm15,%xmm12
  movsd  (%rsp),%xmm15
  mulsd  %xmm13,%xmm3
  mulsd  %xmm11,%xmm15
  addsd  %xmm14,%xmm12
  mulsd  0x8(%rsp),%xmm11
  mulsd  D308(%rip),%xmm13
  divsd  %xmm10,%xmm12
  mulsd  %xmm1,%xmm15
  mulsd  %xmm5,%xmm13
  subsd  %xmm3,%xmm12
  movsd  D119(%rip),%xmm3
  mulsd  %xmm9,%xmm3
  movsd  %xmm12,0x30(%rax)
  movsd  D301(%rip),%xmm12
  mulsd  %xmm9,%xmm12
  mulsd  %xmm8,%xmm9
  subsd  %xmm15,%xmm9
  mulsd  0x30(%rsp),%xmm9
  addsd  %xmm9,%xmm0
  subsd  %xmm11,%xmm0
  addsd  %xmm3,%xmm0
  addsd  %xmm12,%xmm0
  divsd  %xmm4,%xmm0
  subsd  %xmm13,%xmm0
  movsd  %xmm0,0x38(%rax)
  test   %rdx,%rdx
  je     L185
  movsd  0x50(%rdx),%xmm0
  movsd  %xmm0,0x40(%rax)
  movsd  0x58(%rdx),%xmm0
  movsd  %xmm0,0x48(%rax)
  movsd  0x48(%rdx),%xmm11
  mulsd  %xmm7,%xmm0
  movapd %xmm11,%xmm9
  mulsd  %xmm1,%xmm11
  mulsd  %xmm2,%xmm9
  movapd %xmm11,%xmm13
  mulsd  %xmm6,%xmm13
  mulsd  D324(%rip),%xmm13
L68:
  movsd  D119(%rip),%xmm15
  movsd  (%rsp),%xmm14
  movsd  0x28(%rsp),%xmm12
  movsd  0x18(%rsp),%xmm3
  mulsd  %xmm11,%xmm15
  mulsd  %xmm9,%xmm14
  mulsd  %xmm0,%xmm3
  mulsd  %xmm9,%xmm12
  mulsd  0x20(%rsp),%xmm14
  mulsd  %xmm1,%xmm15
  mulsd  0x10(%rsp),%xmm0
  addsd  %xmm3,%xmm14
  movsd  0x38(%rsp),%xmm3
  subsd  %xmm15,%xmm12
  movsd  (%rsp),%xmm15
  mulsd  %xmm13,%xmm3
  mulsd  %xmm11,%xmm15
  addsd  %xmm14,%xmm12
  mulsd  0x8(%rsp),%xmm11
  mulsd  D308(%rip),%xmm13
  divsd  %xmm10,%xmm12
  mulsd  %xmm1,%xmm15
  mulsd  %xmm5,%xmm13
  subsd  %xmm3,%xmm12
  movsd  D119(%rip),%xmm3
  mulsd  %xmm9,%xmm3
  movsd  %xmm12,0x50(%rax)
  movsd  D301(%rip),%xmm12
  mulsd  %xmm9,%xmm12
  mulsd  %xmm8,%xmm9
  subsd  %xmm15,%xmm9
  mulsd  0x30(%rsp),%xmm9
  addsd  %xmm9,%xmm0
  subsd  %xmm11,%xmm0
  addsd  %xmm3,%xmm0
  addsd  %xmm12,%xmm0
  divsd  %xmm4,%xmm0
  subsd  %xmm13,%xmm0
  movsd  %xmm0,0x58(%rax)
  test   %rdx,%rdx
  je     L63
  movsd  0x70(%rdx),%xmm0
  movsd  %xmm0,0x60(%rax)
  movsd  0x78(%rdx),%xmm9
  movsd  %xmm9,0x68(%rax)
  movsd  0x68(%rdx),%xmm11
  mulsd  %xmm7,%xmm9
  movapd %xmm11,%xmm0
  mulsd  %xmm1,%xmm11
  mulsd  %xmm2,%xmm0
  movapd %xmm11,%xmm13
  mulsd  %xmm6,%xmm13
  mulsd  D324(%rip),%xmm13
L224:
  movsd  D119(%rip),%xmm15
  movsd  (%rsp),%xmm14
  movsd  0x28(%rsp),%xmm12
  movsd  0x18(%rsp),%xmm3
  mulsd  %xmm11,%xmm15
  mov    0x10(%rbp),%rdx
  mulsd  %xmm0,%xmm14
  mulsd  %xmm9,%xmm3
  mulsd  %xmm0,%xmm12
  mulsd  0x20(%rsp),%xmm14
  mulsd  %xmm1,%xmm15
  mulsd  0x10(%rsp),%xmm9
  addsd  %xmm3,%xmm14
  movsd  0x38(%rsp),%xmm3
  subsd  %xmm15,%xmm12
  mulsd  %xmm13,%xmm3
  mulsd  D308(%rip),%xmm13
  addsd  %xmm14,%xmm12
  movsd  (%rsp),%xmm14
  mulsd  %xmm11,%xmm14
  divsd  %xmm10,%xmm12
  mulsd  0x8(%rsp),%xmm11
  mulsd  %xmm5,%xmm13
  mulsd  %xmm1,%xmm14
  subsd  %xmm3,%xmm12
  movsd  %xmm12,0x70(%rax)
  movapd %xmm8,%xmm12
  mulsd  %xmm0,%xmm12
  subsd  %xmm14,%xmm12
  mulsd  0x30(%rsp),%xmm12
  addsd  %xmm12,%xmm9
  subsd  %xmm11,%xmm9
  movsd  D119(%rip),%xmm11
  mulsd  %xmm0,%xmm11
  mulsd  D301(%rip),%xmm0
  addsd  %xmm9,%xmm11
  addsd  %xmm11,%xmm0
  divsd  %xmm4,%xmm0
  subsd  %xmm13,%xmm0
  movsd  %xmm0,0x78(%rax)
  mov    0x10(%rbx),%rax
  test   %rdx,%rdx
  je     L217
L175:
  movsd  0x10(%rdx),%xmm0
  test   %rax,%rax
  je     L203
  movsd  %xmm0,(%rax)
  movsd  0x18(%rdx),%xmm11
  movsd  %xmm11,0x8(%rax)
  movsd  0x8(%rdx),%xmm12
  mulsd  %xmm7,%xmm11
  movapd %xmm12,%xmm0
  mulsd  %xmm1,%xmm12
  mulsd  %xmm2,%xmm0
  mulsd  %xmm12,%xmm6
  mulsd  D324(%rip),%xmm6
L215:
  movsd  (%rsp),%xmm15
  movsd  0x28(%rsp),%xmm7
  mulsd  %xmm0,%xmm8
  movsd  0x18(%rsp),%xmm9
  movsd  0x38(%rsp),%xmm3
  divsd  %xmm4,%xmm2
  mulsd  %xmm0,%xmm7
  movapd %xmm15,%xmm14
  movsd  D190(%rip),%xmm13
  mulsd  %xmm0,%xmm14
  mulsd  0x20(%rsp),%xmm14
  mulsd  %xmm11,%xmm9
  mulsd  %xmm6,%xmm3
  mulsd  0x10(%rsp),%xmm11
  mulsd  D308(%rip),%xmm6
  addsd  %xmm14,%xmm9
  movsd  D119(%rip),%xmm14
  mulsd  %xmm12,%xmm14
  mulsd  %xmm5,%xmm6
  mulsd  %xmm1,%xmm14
  divsd  %xmm10,%xmm13
  subsd  %xmm14,%xmm7
  addsd  %xmm9,%xmm7
  divsd  %xmm10,%xmm7
  subsd  %xmm3,%xmm7
  movapd %xmm15,%xmm3
  mulsd  %xmm12,%xmm3
  mulsd  0x8(%rsp),%xmm12
  addsd  %xmm13,%xmm7
  mulsd  %xmm3,%xmm1
  movsd  %xmm7,0x10(%rax)
  subsd  %xmm1,%xmm8
  mulsd  0x30(%rsp),%xmm8
  movsd  D119(%rip),%xmm1
  mulsd  %xmm0,%xmm1
  mulsd  D301(%rip),%xmm0
  addsd  %xmm8,%xmm11
  subsd  %xmm12,%xmm11
  addsd  %xmm11,%xmm1
  addsd  %xmm1,%xmm0
  divsd  %xmm4,%xmm0
  subsd  %xmm6,%xmm0
  addsd  %xmm2,%xmm0
  movsd  %xmm0,0x18(%rax)
L203:
  add    $0x88,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  ret
  nopl   0x0(%rax)
L158:
  movsd  0x18(%rax),%xmm5
  movsd  %xmm5,0x38(%rsp)
  jmp    L8
L76:
  mov    (%rsi),%rax
  test   %rax,%rax
  je     L59
  movsd  D291(%rip),%xmm5
  movsd  D107(%rip),%xmm4
  pxor   %xmm1,%xmm1
  movq   $0x0,(%rax)
  movq   $0x0,0x8(%rax)
  movsd  D190(%rip),%xmm2
  movsd  %xmm5,0x18(%rsp)
  movsd  %xmm1,0x38(%rsp)
  movsd  %xmm4,(%rsp)
  jmp    L142
  nopl   (%rax)
L227:
  mov    0x10(%rbp),%rdx
  divsd  %xmm4,%xmm5
  movapd %xmm2,%xmm6
  mov    0x10(%rbx),%rax
  movsd  0x38(%rsp),%xmm7
  addsd  %xmm2,%xmm6
  addsd  %xmm7,%xmm7
  divsd  %xmm10,%xmm3
  movsd  %xmm3,0x38(%rsp)
  test   %rdx,%rdx
  jne    L175
L217:
  test   %rax,%rax
  je     L203
  pxor   %xmm12,%xmm12
  movapd %xmm2,%xmm0
  movq   $0x0,(%rax)
  mulsd  %xmm12,%xmm7
  movq   $0x0,0x8(%rax)
  mulsd  %xmm12,%xmm0
  mulsd  %xmm1,%xmm12
  movapd %xmm7,%xmm11
  mulsd  %xmm12,%xmm6
  mulsd  D324(%rip),%xmm6
  jmp    L215
  nopl   0x0(%rax)
L39:
  test   %rax,%rax
  je     L227
  movsd  0x38(%rsp),%xmm7
  pxor   %xmm11,%xmm11
  movapd %xmm2,%xmm0
  movq   $0x0,(%rax)
  mulsd  %xmm11,%xmm0
  movq   $0x0,0x8(%rax)
  addsd  %xmm7,%xmm7
  movapd %xmm7,%xmm9
  mulsd  %xmm11,%xmm9
  jmp    L167
  nop
L63:
  pxor   %xmm11,%xmm11
  movapd %xmm2,%xmm0
  movapd %xmm7,%xmm9
  movq   $0x0,0x60(%rax)
  mulsd  %xmm11,%xmm0
  movapd %xmm6,%xmm13
  movq   $0x0,0x68(%rax)
  mulsd  %xmm11,%xmm9
  mulsd  %xmm1,%xmm11
  mulsd  %xmm11,%xmm13
  mulsd  D324(%rip),%xmm13
  jmp    L224
  nopl   (%rax)
L185:
  pxor   %xmm11,%xmm11
  movapd %xmm2,%xmm9
  movapd %xmm7,%xmm0
  movq   $0x0,0x40(%rax)
  mulsd  %xmm11,%xmm9
  movapd %xmm6,%xmm13
  movq   $0x0,0x48(%rax)
  mulsd  %xmm11,%xmm0
  mulsd  %xmm1,%xmm11
  mulsd  %xmm11,%xmm13
  mulsd  D324(%rip),%xmm13
  jmp    L68
  nopl   (%rax)
L60:
  pxor   %xmm11,%xmm11
  movapd %xmm2,%xmm9
  movapd %xmm7,%xmm0
  movq   $0x0,0x20(%rax)
  mulsd  %xmm11,%xmm9
  movapd %xmm6,%xmm13
  movq   $0x0,0x28(%rax)
  mulsd  %xmm11,%xmm0
  mulsd  %xmm1,%xmm11
  mulsd  %xmm11,%xmm13
  mulsd  D324(%rip),%xmm13
  jmp    L34
  nopl   (%rax)
L59:
  movsd  D291(%rip),%xmm4
  movsd  D107(%rip),%xmm5
  pxor   %xmm1,%xmm1
  movsd  D190(%rip),%xmm2
  movsd  %xmm1,0x38(%rsp)
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm5,(%rsp)
  jmp    L142
```

tool_calls: name: read_data
arguments: label: D107
datatype: byte[8]
name: read_data
arguments: label: D119
datatype: byte[8]
name: read_data
arguments: label: D324
datatype: byte[8]
name: read_data
arguments: label: D236
datatype: byte[8]
name: read_data
arguments: label: D301
datatype: byte[8]
name: read_data
arguments: label: D308
datatype: byte[8]
name: read_data
arguments: label: D190
datatype: f64
name: read_data
arguments: label: D291
datatype: f32[2]

tool: {"D107": [0x7c,0x14,0xae,0x47,0xe1,0x7a,0xb4,0xbf]}

tool: {"D119": [0x32,0x08,0xac,0x1c,0x5a,0x64,0xef,0x3f]}

tool: {"D324": [0x9a,0x99,0x99,0x99,0x99,0x99,0xb9,0x3f]}

tool: {"D236": [0x9a,0x99,0x99,0x99,0x99,0x99,0xf1,0x3f]}

tool: {"D301": [0x1f,0x85,0xeb,0x51,0xb8,0x9e,0x23,0x40]}

tool: {"D308": [0x9a,0x99,0x99,0x99,0x99,0x99,0xe9,0x3f]}

tool: D190: 1

tool: {"D291": [0.0,-0.0]}

assistant: ```c
int vdeFun(const double **arg, double **res, int *iw, double *w, int mem) {
  double a0 = arg[0] ? arg[0][2] : 0;
  if (res[0] != 0)
    res[0][0] = a0;
  a0 = arg[0] ? arg[0][3] : 0;
  if (res[0] != 0)
    res[0][1] = a0;
  double a1 = arg[0] ? arg[0][1] : 0;
  double a2 = sin(a1);
  double a3 = -8.0000000000000016e-02;
  a2 = (a3 * a2);
  double a4 = vde_forw_pendulum_sq(a0);
  double a5 = (a2 * a4);
  double a6 = arg[3] ? arg[3][0] : 0;
  a5 = (a5 + a6);
  double a7 = cos(a1);
  double a8 = 9.8100000000000009e-01;
  a7 = (a8 * a7);
  double a9 = sin(a1);
  double a10 = (a7 * a9);
  a5 = (a5 + a10);
  a10 = cos(a1);
  double a11 = vde_forw_pendulum_sq(a10);
  double a12 = 1.0000000000000001e-01;
  a11 = (a12 * a11);
  double a13 = 1.1000000000000001e+00;
  a11 = (a13 - a11);
  a5 = (a5 / a11);
  if (res[0] != 0)
    res[0][2] = a5;
  double a14 = cos(a1);
  a14 = (a3 * a14);
  double a15 = sin(a1);
  double a16 = (a14 * a15);
  double a17 = vde_forw_pendulum_sq(a0);
  double a18 = (a16 * a17);
  double a19 = cos(a1);
  double a20 = (a6 * a19);
  a18 = (a18 + a20);
  a20 = sin(a1);
  a20 = (a8 * a20);
  a18 = (a18 + a20);
  a20 = sin(a1);
  double a21 = 9.8100000000000005e+00;
  a20 = (a21 * a20);
  a18 = (a18 + a20);
  a20 = cos(a1);
  double a22 = vde_forw_pendulum_sq(a20);
  a22 = (a12 * a22);
  a13 = (a13 - a22);
  a22 = 8.0000000000000004e-01;
  a13 = (a22 * a13);
  a18 = (a18 / a13);
  if (res[0] != 0)
    res[0][3] = a18;
  double a23 = arg[1] ? arg[1][2] : 0;
  if (res[1] != 0)
    res[1][0] = a23;
  a23 = arg[1] ? arg[1][3] : 0;
  if (res[1] != 0)
    res[1][1] = a23;
  double a24 = cos(a1);
  double a25 = arg[1] ? arg[1][1] : 0;
  double a26 = (a24 * a25);
  a26 = (a3 * a26);
  a26 = (a4 * a26);
  double a27 = (a0 + a0);
  double a28 = (a27 * a23);
  a28 = (a2 * a28);
  a26 = (a26 + a28);
  a28 = cos(a1);
  double a29 = (a28 * a25);
  a29 = (a7 * a29);
  double a30 = sin(a1);
  double a31 = (a30 * a25);
  a31 = (a8 * a31);
  a31 = (a9 * a31);
  a29 = (a29 - a31);
  a26 = (a26 + a29);
  a26 = (a26 / a11);
  a29 = (a5 / a11);
  a31 = (a10 + a10);
  double a32 = sin(a1);
  double a33 = (a32 * a25);
  a33 = (a31 * a33);
  a33 = (a12 * a33);
  a33 = (a29 * a33);
  a26 = (a26 - a33);
  if (res[1] != 0)
    res[1][2] = a26;
  a26 = cos(a1);
  a33 = (a26 * a25);
  a33 = (a14 * a33);
  double a34 = sin(a1);
  double a35 = (a34 * a25);
  a35 = (a3 * a35);
  a35 = (a15 * a35);
  a33 = (a33 - a35);
  a33 = (a17 * a33);
  a35 = (a0 + a0);
  a23 = (a35 * a23);
  a23 = (a16 * a23);
  a33 = (a33 + a23);
  a23 = sin(a1);
  double a36 = (a23 * a25);
  a36 = (a6 * a36);
  a33 = (a33 - a36);
  a36 = cos(a1);
  double a37 = (a36 * a25);
  a37 = (a8 * a37);
  a33 = (a33 + a37);
  a37 = cos(a1);
  double a38 = (a37 * a25);
  a38 = (a21 * a38);
  a33 = (a33 + a38);
  a33 = (a33 / a13);
  a38 = (a18 / a13);
  double a39 = (a20 + a20);
  double a40 = sin(a1);
  a25 = (a40 * a25);
  a25 = (a39 * a25);
  a25 = (a12 * a25);
  a25 = (a22 * a25);
  a25 = (a38 * a25);
  a33 = (a33 - a25);
  if (res[1] != 0)
    res[1][3] = a33;
  a33 = arg[1] ? arg[1][6] : 0;
  if (res[1] != 0)
    res[1][4] = a33;
  a33 = arg[1] ? arg[1][7] : 0;
  if (res[1] != 0)
    res[1][5] = a33;
  a25 = arg[1] ? arg[1][5] : 0;
  double a41 = (a24 * a25);
  a41 = (a3 * a41);
  a41 = (a4 * a41);
  double a42 = (a27 * a33);
  a42 = (a2 * a42);
  a41 = (a41 + a42);
  a42 = (a28 * a25);
  a42 = (a7 * a42);
  double a43 = (a30 * a25);
  a43 = (a8 * a43);
  a43 = (a9 * a43);
  a42 = (a42 - a43);
  a41 = (a41 + a42);
  a41 = (a41 / a11);
  a42 = (a32 * a25);
  a42 = (a31 * a42);
  a42 = (a12 * a42);
  a42 = (a29 * a42);
  a41 = (a41 - a42);
  if (res[1] != 0)
    res[1][6] = a41;
  a41 = (a26 * a25);
  a41 = (a14 * a41);
  a42 = (a34 * a25);
  a42 = (a3 * a42);
  a42 = (a15 * a42);
  a41 = (a41 - a42);
  a41 = (a17 * a41);
  a33 = (a35 * a33);
  a33 = (a16 * a33);
  a41 = (a41 + a33);
  a33 = (a23 * a25);
  a33 = (a6 * a33);
  a41 = (a41 - a33);
  a33 = (a36 * a25);
  a33 = (a8 * a33);
  a41 = (a41 + a33);
  a33 = (a37 * a25);
  a33 = (a21 * a33);
  a41 = (a41 + a33);
  a41 = (a41 / a13);
  a25 = (a40 * a25);
  a25 = (a39 * a25);
  a25 = (a12 * a25);
  a25 = (a22 * a25);
  a25 = (a38 * a25);
  a41 = (a41 - a25);
  if (res[1] != 0)
    res[1][7] = a41;
  a41 = arg[1] ? arg[1][10] : 0;
  if (res[1] != 0)
    res[1][8] = a41;
  a41 = arg[1] ? arg[1][11] : 0;
  if (res[1] != 0)
    res[1][9] = a41;
  a25 = arg[1] ? arg[1][9] : 0;
  a33 = (a24 * a25);
  a33 = (a3 * a33);
  a33 = (a4 * a33);
  a42 = (a27 * a41);
  a42 = (a2 * a42);
  a33 = (a33 + a42);
  a42 = (a28 * a25);
  a42 = (a7 * a42);
  a43 = (a30 * a25);
  a43 = (a8 * a43);
  a43 = (a9 * a43);
  a42 = (a42 - a43);
  a33 = (a33 + a42);
  a33 = (a33 / a11);
  a42 = (a32 * a25);
  a42 = (a31 * a42);
  a42 = (a12 * a42);
  a42 = (a29 * a42);
  a33 = (a33 - a42);
  if (res[1] != 0)
    res[1][10] = a33;
  a33 = (a26 * a25);
  a33 = (a14 * a33);
  a42 = (a34 * a25);
  a42 = (a3 * a42);
  a42 = (a15 * a42);
  a33 = (a33 - a42);
  a33 = (a17 * a33);
  a41 = (a35 * a41);
  a41 = (a16 * a41);
  a33 = (a33 + a41);
  a41 = (a23 * a25);
  a41 = (a6 * a41);
  a33 = (a33 - a41);
  a41 = (a36 * a25);
  a41 = (a8 * a41);
  a33 = (a33 + a41);
  a41 = (a37 * a25);
  a41 = (a21 * a41);
  a33 = (a33 + a41);
  a33 = (a33 / a13);
  a25 = (a40 * a25);
  a25 = (a39 * a25);
  a25 = (a12 * a25);
  a25 = (a22 * a25);
  a25 = (a38 * a25);
  a33 = (a33 - a25);
  if (res[1] != 0)
    res[1][11] = a33;
  a33 = arg[1] ? arg[1][14] : 0;
  if (res[1] != 0)
    res[1][12] = a33;
  a33 = arg[1] ? arg[1][15] : 0;
  if (res[1] != 0)
    res[1][13] = a33;
  a25 = arg[1] ? arg[1][13] : 0;
  a24 = (a24 * a25);
  a24 = (a3 * a24);
  a24 = (a4 * a24);
  a27 = (a27 * a33);
  a27 = (a2 * a27);
  a24 = (a24 + a27);
  a28 = (a28 * a25);
  a28 = (a7 * a28);
  a30 = (a30 * a25);
  a30 = (a8 * a30);
  a30 = (a9 * a30);
  a28 = (a28 - a30);
  a24 = (a24 + a28);
  a24 = (a24 / a11);
  a32 = (a32 * a25);
  a31 = (a31 * a32);
  a31 = (a12 * a31);
  a29 = (a29 * a31);
  a24 = (a24 - a29);
  if (res[1] != 0)
    res[1][14] = a24;
  a26 = (a26 * a25);
  a26 = (a14 * a26);
  a34 = (a34 * a25);
  a34 = (a3 * a34);
  a34 = (a15 * a34);
  a26 = (a26 - a34);
  a26 = (a17 * a26);
  a35 = (a35 * a33);
  a35 = (a16 * a35);
  a26 = (a26 + a35);
  a23 = (a23 * a25);
  a23 = (a6 * a23);
  a26 = (a26 - a23);
  a36 = (a36 * a25);
  a36 = (a8 * a36);
  a26 = (a26 + a36);
  a37 = (a37 * a25);
  a37 = (a21 * a37);
  a26 = (a26 + a37);
  a26 = (a26 / a13);
  a40 = (a40 * a25);
  a39 = (a39 * a40);
  a39 = (a12 * a39);
  a39 = (a22 * a39);
  a38 = (a38 * a39);
  a26 = (a26 - a38);
  if (res[1] != 0)
    res[1][15] = a26;
  a26 = arg[2] ? arg[2][2] : 0;
  if (res[2] != 0)
    res[2][0] = a26;
  a26 = arg[2] ? arg[2][3] : 0;
  if (res[2] != 0)
    res[2][1] = a26;
  a38 = cos(a1);
  a39 = arg[2] ? arg[2][1] : 0;
  a38 = (a38 * a39);
  a38 = (a3 * a38);
  a4 = (a4 * a38);
  a38 = (a0 + a0);
  a38 = (a38 * a26);
  a2 = (a2 * a38);
  a4 = (a4 + a2);
  a2 = cos(a1);
  a2 = (a2 * a39);
  a7 = (a7 * a2);
  a2 = sin(a1);
  a2 = (a2 * a39);
  a2 = (a8 * a2);
  a9 = (a9 * a2);
  a7 = (a7 - a9);
  a4 = (a4 + a7);
  a4 = (a4 / a11);
  a5 = (a5 / a11);
  a10 = (a10 + a10);
  a7 = sin(a1);
  a7 = (a7 * a39);
  a10 = (a10 * a7);
  a10 = (a12 * a10);
  a5 = (a5 * a10);
  a4 = (a4 - a5);
  a11 = (1. / a11);
  a11 = (a11 + a4);
  if (res[2] != 0)
    res[2][2] = a11;
  a19 = (a19 / a13);
  a11 = cos(a1);
  a11 = (a11 * a39);
  a14 = (a14 * a11);
  a11 = sin(a1);
  a11 = (a11 * a39);
  a3 = (a3 * a11);
  a15 = (a15 * a3);
  a14 = (a14 - a15);
  a17 = (a17 * a14);
  a0 = (a0 + a0);
  a0 = (a0 * a26);
  a16 = (a16 * a0);
  a17 = (a17 + a16);
  a16 = sin(a1);
  a16 = (a16 * a39);
  a6 = (a6 * a16);
  a17 = (a17 - a6);
  a6 = cos(a1);
  a6 = (a6 * a39);
  a8 = (a8 * a6);
  a17 = (a17 + a8);
  a8 = cos(a1);
  a8 = (a8 * a39);
  a21 = (a21 * a8);
  a17 = (a17 + a21);
  a17 = (a17 / a13);
  a18 = (a18 / a13);
  a20 = (a20 + a20);
  a1 = sin(a1);
  a1 = (a1 * a39);
  a20 = (a20 * a1);
  a12 = (a12 * a20);
  a22 = (a22 * a12);
  a18 = (a18 * a22);
  a17 = (a17 - a18);
  a19 = (a19 + a17);
  if (res[2] != 0)
    res[2][3] = a19;
  return 0;
}
```