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

```asm
<vdeFun>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x60,%rsp
  mov    %rdi,%rbp
  mov    %rsi,%rbx
  mov    (%rdi),%rax
  pxor   %xmm0,%xmm0
  test   %rax,%rax
  je     L83
  movsd  0x10(%rax),%xmm0
L83:
  mov    (%rbx),%rax
  test   %rax,%rax
  je     L624
  movsd  %xmm0,(%rax)
  mov    0x0(%rbp),%rax
  test   %rax,%rax
  je     L84
L513:
  mov    0x18(%rax),%r14
  mov    (%rbx),%rdx
  test   %rdx,%rdx
  je     L489
L540:
  mov    %r14,0x8(%rdx)
  mov    0x0(%rbp),%rax
  pxor   %xmm0,%xmm0
  test   %rax,%rax
  je     L319
L489:
  movsd  0x8(%rax),%xmm0
L319:
  lea    0x58(%rsp),%rdi
  lea    0x50(%rsp),%rsi
  call   <sincos@plt>
  movsd  0x50(%rsp),%xmm2
  movsd  %xmm2,(%rsp)
  movsd  0x58(%rsp),%xmm3
  movsd  %xmm3,0x8(%rsp)
  mulsd  D136(%rip),%xmm3
  movsd  %xmm3,0x10(%rsp)
  movq   %r14,%xmm0
  call   <vde_forw_pendulum_sq@plt>
  movsd  %xmm0,0x28(%rsp)
  movsd  0x10(%rsp),%xmm3
  mulsd  %xmm0,%xmm3
  movapd %xmm3,%xmm0
  mov    0x18(%rbp),%rax
  mov    D310(%rip),%r12
  test   %rax,%rax
  je     L380
  mov    (%rax),%r12
L380:
  movq   %r12,%xmm4
  addsd  %xmm4,%xmm0
  movsd  (%rsp),%xmm5
  movapd %xmm5,%xmm4
  mulsd  D248(%rip),%xmm4
  movsd  %xmm4,0x30(%rsp)
  movsd  0x8(%rsp),%xmm1
  mulsd  %xmm4,%xmm1
  addsd  %xmm0,%xmm1
  movq   %xmm1,%r13
  movapd %xmm5,%xmm0
  call   <vde_forw_pendulum_sq@plt>
  mulsd  D91(%rip),%xmm0
  movsd  D378(%rip),%xmm1
  movapd %xmm1,%xmm5
  subsd  %xmm0,%xmm5
  movsd  %xmm5,0x20(%rsp)
  movq   %r13,%xmm3
  divsd  %xmm5,%xmm3
  movq   %xmm3,%r13
  mov    (%rbx),%rax
  test   %rax,%rax
  je     L566
  movsd  %xmm3,0x10(%rax)
L566:
  movsd  (%rsp),%xmm5
  mulsd  D136(%rip),%xmm5
  movsd  %xmm5,0x38(%rsp)
  movsd  0x8(%rsp),%xmm3
  mulsd  %xmm5,%xmm3
  movsd  %xmm3,0x18(%rsp)
  movq   %r14,%xmm0
  call   <vde_forw_pendulum_sq@plt>
  movsd  %xmm0,0x40(%rsp)
  movsd  0x18(%rsp),%xmm4
  mulsd  %xmm0,%xmm4
  movapd %xmm4,%xmm0
  movq   %r12,%xmm1
  movsd  (%rsp),%xmm4
  mulsd  %xmm4,%xmm1
  addsd  %xmm0,%xmm1
  movsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm2
  mulsd  D248(%rip),%xmm2
  movapd %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mulsd  D43(%rip),%xmm5
  addsd  %xmm0,%xmm5
  movsd  %xmm5,0x48(%rsp)
  movapd %xmm4,%xmm0
  call   <vde_forw_pendulum_sq@plt>
  mulsd  D91(%rip),%xmm0
  movsd  D378(%rip),%xmm4
  subsd  %xmm0,%xmm4
  mulsd  D214(%rip),%xmm4
  movsd  0x48(%rsp),%xmm5
  divsd  %xmm4,%xmm5
  mov    (%rbx),%rax
  test   %rax,%rax
  je     L256
  movsd  %xmm5,0x18(%rax)
L256:
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L464
  movsd  0x10(%rax),%xmm0
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L509
L162:
  movsd  %xmm0,(%rax)
L509:
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L627
  movsd  0x18(%rax),%xmm6
  mov    0x8(%rbx),%rdx
  test   %rdx,%rdx
  je     L147
L221:
  movsd  %xmm6,0x8(%rdx)
  mov    0x8(%rbp),%rax
  pxor   %xmm9,%xmm9
  test   %rax,%rax
  je     L265
L147:
  movsd  0x8(%rax),%xmm9
L265:
  movsd  (%rsp),%xmm7
  movapd %xmm7,%xmm0
  mulsd  %xmm9,%xmm0
  movq   %r14,%xmm3
  addsd  %xmm3,%xmm3
  mulsd  %xmm3,%xmm6
  movsd  0x8(%rsp),%xmm13
  mulsd  %xmm13,%xmm9
  movq   %r13,%xmm1
  movsd  0x20(%rsp),%xmm14
  divsd  %xmm14,%xmm1
  movapd %xmm7,%xmm2
  addsd  %xmm7,%xmm2
  movapd %xmm9,%xmm8
  mulsd  %xmm2,%xmm8
  mulsd  D91(%rip),%xmm8
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L484
  movapd %xmm0,%xmm7
  mulsd  D136(%rip),%xmm7
  mulsd  0x28(%rsp),%xmm7
  movsd  0x10(%rsp),%xmm11
  mulsd  %xmm6,%xmm11
  addsd  %xmm7,%xmm11
  movapd %xmm9,%xmm10
  mulsd  D248(%rip),%xmm10
  mulsd  %xmm13,%xmm10
  movsd  0x30(%rsp),%xmm7
  mulsd  %xmm0,%xmm7
  subsd  %xmm10,%xmm7
  addsd  %xmm11,%xmm7
  divsd  %xmm14,%xmm7
  movapd %xmm1,%xmm10
  mulsd  %xmm8,%xmm10
  subsd  %xmm10,%xmm7
  movsd  %xmm7,0x10(%rax)
  divsd  %xmm4,%xmm5
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L403
  mulsd  D214(%rip),%xmm8
  movapd %xmm9,%xmm10
  mulsd  D136(%rip),%xmm10
  movsd  0x38(%rsp),%xmm7
  mulsd  %xmm0,%xmm7
  mulsd  %xmm13,%xmm10
  subsd  %xmm10,%xmm7
  mulsd  0x40(%rsp),%xmm7
  mulsd  0x18(%rsp),%xmm6
  addsd  %xmm7,%xmm6
  movq   %r12,%xmm7
  mulsd  %xmm7,%xmm9
  subsd  %xmm9,%xmm6
  movapd %xmm0,%xmm7
  mulsd  D248(%rip),%xmm7
  addsd  %xmm7,%xmm6
  mulsd  D43(%rip),%xmm0
  addsd  %xmm6,%xmm0
  divsd  %xmm4,%xmm0
  mulsd  %xmm5,%xmm8
  subsd  %xmm8,%xmm0
  movsd  %xmm0,0x18(%rax)
L403:
  mov    0x8(%rbp),%rax
  pxor   %xmm0,%xmm0
  test   %rax,%rax
  je     L215
  movsd  0x30(%rax),%xmm0
L215:
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L1
  movsd  %xmm0,0x20(%rax)
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L377
  movsd  0x38(%rax),%xmm0
  mov    0x8(%rbx),%rdx
  test   %rdx,%rdx
  je     L345
L581:
  movsd  %xmm0,0x28(%rdx)
  mov    0x8(%rbp),%rax
  pxor   %xmm7,%xmm7
  test   %rax,%rax
  je     L225
L345:
  movsd  0x28(%rax),%xmm7
L225:
  movsd  (%rsp),%xmm6
  mulsd  %xmm7,%xmm6
  mulsd  %xmm3,%xmm0
  movsd  0x8(%rsp),%xmm14
  mulsd  %xmm14,%xmm7
  movapd %xmm7,%xmm9
  mulsd  %xmm2,%xmm9
  mulsd  D91(%rip),%xmm9
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L284
  movapd %xmm6,%xmm8
  mulsd  D136(%rip),%xmm8
  mulsd  0x28(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm11
  mulsd  %xmm0,%xmm11
  addsd  %xmm8,%xmm11
  movapd %xmm7,%xmm10
  mulsd  D248(%rip),%xmm10
  mulsd  %xmm14,%xmm10
  movsd  0x30(%rsp),%xmm8
  mulsd  %xmm6,%xmm8
  subsd  %xmm10,%xmm8
  addsd  %xmm11,%xmm8
  divsd  0x20(%rsp),%xmm8
  movapd %xmm9,%xmm10
  mulsd  %xmm1,%xmm10
  subsd  %xmm10,%xmm8
  movsd  %xmm8,0x30(%rax)
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L284
  mulsd  D214(%rip),%xmm9
  movapd %xmm6,%xmm8
  mulsd  D43(%rip),%xmm8
  movapd %xmm6,%xmm10
  mulsd  D248(%rip),%xmm10
  movapd %xmm7,%xmm11
  mulsd  D136(%rip),%xmm11
  mulsd  0x38(%rsp),%xmm6
  mulsd  %xmm14,%xmm11
  subsd  %xmm11,%xmm6
  mulsd  0x40(%rsp),%xmm6
  mulsd  0x18(%rsp),%xmm0
  addsd  %xmm6,%xmm0
  movq   %r12,%xmm6
  mulsd  %xmm6,%xmm7
  subsd  %xmm7,%xmm0
  addsd  %xmm10,%xmm0
  addsd  %xmm8,%xmm0
  divsd  %xmm4,%xmm0
  mulsd  %xmm5,%xmm9
  subsd  %xmm9,%xmm0
  movsd  %xmm0,0x38(%rax)
L284:
  mov    0x8(%rbp),%rax
  pxor   %xmm0,%xmm0
  test   %rax,%rax
  je     L206
  movsd  0x50(%rax),%xmm0
L206:
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L281
  movsd  %xmm0,0x40(%rax)
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L243
  movsd  0x58(%rax),%xmm0
  mov    0x8(%rbx),%rdx
  test   %rdx,%rdx
  je     L126
L505:
  movsd  %xmm0,0x48(%rdx)
  mov    0x8(%rbp),%rax
  pxor   %xmm7,%xmm7
  test   %rax,%rax
  je     L361
L126:
  movsd  0x48(%rax),%xmm7
L361:
  movsd  (%rsp),%xmm6
  mulsd  %xmm7,%xmm6
  mulsd  %xmm3,%xmm0
  movsd  0x8(%rsp),%xmm15
  mulsd  %xmm15,%xmm7
  movapd %xmm7,%xmm9
  mulsd  %xmm2,%xmm9
  mulsd  D91(%rip),%xmm9
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L588
  movapd %xmm6,%xmm8
  mulsd  D136(%rip),%xmm8
  mulsd  0x28(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm11
  mulsd  %xmm0,%xmm11
  addsd  %xmm8,%xmm11
  movapd %xmm7,%xmm10
  mulsd  D248(%rip),%xmm10
  mulsd  %xmm15,%xmm10
  movsd  0x30(%rsp),%xmm8
  mulsd  %xmm6,%xmm8
  subsd  %xmm10,%xmm8
  addsd  %xmm11,%xmm8
  divsd  0x20(%rsp),%xmm8
  movapd %xmm9,%xmm10
  mulsd  %xmm1,%xmm10
  subsd  %xmm10,%xmm8
  movsd  %xmm8,0x50(%rax)
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L588
  mulsd  D214(%rip),%xmm9
  movapd %xmm6,%xmm8
  mulsd  D43(%rip),%xmm8
  movapd %xmm6,%xmm10
  mulsd  D248(%rip),%xmm10
  movapd %xmm7,%xmm11
  mulsd  D136(%rip),%xmm11
  mulsd  0x38(%rsp),%xmm6
  mulsd  %xmm15,%xmm11
  subsd  %xmm11,%xmm6
  mulsd  0x40(%rsp),%xmm6
  mulsd  0x18(%rsp),%xmm0
  addsd  %xmm6,%xmm0
  movq   %r12,%xmm6
  mulsd  %xmm6,%xmm7
  subsd  %xmm7,%xmm0
  addsd  %xmm10,%xmm0
  addsd  %xmm8,%xmm0
  divsd  %xmm4,%xmm0
  mulsd  %xmm5,%xmm9
  subsd  %xmm9,%xmm0
  movsd  %xmm0,0x58(%rax)
L588:
  mov    0x8(%rbp),%rax
  pxor   %xmm0,%xmm0
  test   %rax,%rax
  je     L34
  movsd  0x70(%rax),%xmm0
L34:
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L92
  movsd  %xmm0,0x60(%rax)
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L532
  movsd  0x78(%rax),%xmm6
  mov    0x8(%rbx),%rdx
  test   %rdx,%rdx
  je     L472
L12:
  movsd  %xmm6,0x68(%rdx)
  mov    0x8(%rbp),%rax
  pxor   %xmm7,%xmm7
  test   %rax,%rax
  je     L213
L472:
  movsd  0x68(%rax),%xmm7
L213:
  movsd  (%rsp),%xmm0
  mulsd  %xmm7,%xmm0
  mulsd  %xmm3,%xmm6
  movsd  0x8(%rsp),%xmm12
  mulsd  %xmm12,%xmm7
  movapd %xmm7,%xmm9
  mulsd  %xmm2,%xmm9
  mulsd  D91(%rip),%xmm9
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L111
  movapd %xmm0,%xmm8
  mulsd  D136(%rip),%xmm8
  mulsd  0x28(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm11
  mulsd  %xmm6,%xmm11
  addsd  %xmm8,%xmm11
  movapd %xmm7,%xmm10
  mulsd  D248(%rip),%xmm10
  mulsd  %xmm12,%xmm10
  movsd  0x30(%rsp),%xmm8
  mulsd  %xmm0,%xmm8
  subsd  %xmm10,%xmm8
  addsd  %xmm11,%xmm8
  divsd  0x20(%rsp),%xmm8
  movapd %xmm9,%xmm10
  mulsd  %xmm1,%xmm10
  subsd  %xmm10,%xmm8
  movsd  %xmm8,0x70(%rax)
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L111
  movapd %xmm7,%xmm10
  mulsd  D136(%rip),%xmm10
  movsd  0x38(%rsp),%xmm8
  mulsd  %xmm0,%xmm8
  mulsd  %xmm12,%xmm10
  subsd  %xmm10,%xmm8
  mulsd  0x40(%rsp),%xmm8
  mulsd  0x18(%rsp),%xmm6
  addsd  %xmm8,%xmm6
  movq   %r12,%xmm12
  mulsd  %xmm12,%xmm7
  subsd  %xmm7,%xmm6
  movapd %xmm0,%xmm7
  mulsd  D248(%rip),%xmm7
  addsd  %xmm7,%xmm6
  mulsd  D43(%rip),%xmm0
  addsd  %xmm6,%xmm0
  divsd  %xmm4,%xmm0
  mulsd  D214(%rip),%xmm9
  mulsd  %xmm5,%xmm9
  subsd  %xmm9,%xmm0
  movsd  %xmm0,0x78(%rax)
L111:
  mov    0x10(%rbp),%rax
  pxor   %xmm0,%xmm0
  test   %rax,%rax
  je     L5
  movsd  0x10(%rax),%xmm0
L5:
  mov    0x10(%rbx),%rax
  test   %rax,%rax
  je     L606
  movsd  %xmm0,(%rax)
  mov    0x10(%rbp),%rax
  test   %rax,%rax
  je     L46
  movsd  0x18(%rax),%xmm6
  mov    0x10(%rbx),%rdx
  test   %rdx,%rdx
  je     L119
L515:
  movsd  %xmm6,0x8(%rdx)
  mov    0x10(%rbp),%rax
  pxor   %xmm7,%xmm7
  test   %rax,%rax
  je     L317
L119:
  movsd  0x8(%rax),%xmm7
L317:
  movsd  (%rsp),%xmm12
  movapd %xmm12,%xmm0
  mulsd  %xmm7,%xmm0
  mulsd  %xmm3,%xmm6
  movapd %xmm7,%xmm3
  movsd  0x8(%rsp),%xmm13
  mulsd  %xmm13,%xmm3
  mulsd  %xmm3,%xmm2
  mulsd  D91(%rip),%xmm2
  mov    0x10(%rbx),%rax
  test   %rax,%rax
  je     L45
  movsd  D133(%rip),%xmm9
  movsd  0x20(%rsp),%xmm15
  divsd  %xmm15,%xmm9
  movsd  0x30(%rsp),%xmm8
  mulsd  %xmm0,%xmm8
  movapd %xmm0,%xmm10
  mulsd  D136(%rip),%xmm10
  mulsd  0x28(%rsp),%xmm10
  movsd  0x10(%rsp),%xmm7
  mulsd  %xmm6,%xmm7
  addsd  %xmm10,%xmm7
  movapd %xmm3,%xmm10
  mulsd  D248(%rip),%xmm10
  mulsd  %xmm13,%xmm10
  subsd  %xmm10,%xmm8
  addsd  %xmm8,%xmm7
  divsd  %xmm15,%xmm7
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm7
  addsd  %xmm9,%xmm7
  movsd  %xmm7,0x10(%rax)
  mov    0x10(%rbx),%rax
  test   %rax,%rax
  je     L45
  movsd  0x38(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  movapd %xmm3,%xmm7
  mulsd  D136(%rip),%xmm7
  mulsd  %xmm13,%xmm7
  subsd  %xmm7,%xmm1
  mulsd  0x40(%rsp),%xmm1
  mulsd  0x18(%rsp),%xmm6
  addsd  %xmm1,%xmm6
  movq   %r12,%xmm1
  mulsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm6
  movapd %xmm0,%xmm1
  mulsd  D248(%rip),%xmm1
  addsd  %xmm6,%xmm1
  mulsd  D43(%rip),%xmm0
  addsd  %xmm1,%xmm0
  divsd  %xmm4,%xmm0
  movapd %xmm12,%xmm1
  divsd  %xmm4,%xmm1
  mulsd  D214(%rip),%xmm2
  mulsd  %xmm5,%xmm2
  subsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%rax)
L45:
  mov    $0x0,%eax
  add    $0x60,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L46:
  mov    0x10(%rbx),%rdx
  pxor   %xmm6,%xmm6
  test   %rdx,%rdx
  jne    L515
  jmp    L45
L532:
  mov    0x8(%rbx),%rdx
  pxor   %xmm6,%xmm6
  test   %rdx,%rdx
  jne    L12
  jmp    L111
L243:
  mov    0x8(%rbx),%rdx
  pxor   %xmm0,%xmm0
  test   %rdx,%rdx
  jne    L505
  jmp    L588
L377:
  mov    0x8(%rbx),%rdx
  pxor   %xmm0,%xmm0
  test   %rdx,%rdx
  jne    L581
  jmp    L284
L84:
  mov    (%rbx),%rdx
  test   %rdx,%rdx
  je     L398
  mov    D310(%rip),%r14
  jmp    L540
L398:
  mov    D310(%rip),%r14
  movq   %r14,%xmm0
  jmp    L319
L606:
  mov    0x10(%rbp),%rax
  test   %rax,%rax
  je     L46
  movsd  0x18(%rax),%xmm6
  jmp    L119
L92:
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L532
  movsd  0x78(%rax),%xmm6
  jmp    L472
L281:
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L243
  movsd  0x58(%rax),%xmm0
  jmp    L126
L1:
  mov    0x8(%rbp),%rax
  test   %rax,%rax
  je     L377
  movsd  0x38(%rax),%xmm0
  jmp    L345
L464:
  mov    0x8(%rbx),%rax
  pxor   %xmm0,%xmm0
  test   %rax,%rax
  jne    L162
L627:
  mov    0x8(%rbx),%rdx
  pxor   %xmm6,%xmm6
  test   %rdx,%rdx
  jne    L221
  movq   %r14,%xmm3
  addsd  %xmm3,%xmm3
  movq   %r13,%xmm1
  divsd  0x20(%rsp),%xmm1
  movsd  (%rsp),%xmm6
  movapd %xmm6,%xmm2
  addsd  %xmm6,%xmm2
L484:
  divsd  %xmm4,%xmm5
  jmp    L403
L624:
  mov    0x0(%rbp),%rax
  mov    D310(%rip),%r14
  movq   %r14,%xmm0
  test   %rax,%rax
  jne    L513
  jmp    L319
```

tool_calls: name: unpack_data
arguments: data_label: D136
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D310
type_guessed: f64
name: unpack_data
arguments: data_label: D248
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D91
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D378
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D43
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D214
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D133
type_guessed: f64

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

tool: D310: 0

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

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

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

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

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

tool: D133: 1

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;
}
```