user: 提供的汇编代码用C语言表示会是什么样子:

```asm
<compute_gv>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x230,%rsp
  mov    %rdi,-0x1e8(%rbp)
  mov    %rsi,-0x1f0(%rbp)
  movsd  %xmm0,-0x1f8(%rbp)
  movsd  %xmm1,-0x200(%rbp)
  movsd  %xmm2,-0x208(%rbp)
  movsd  %xmm3,-0x210(%rbp)
  mov    %rdx,-0x218(%rbp)
  mov    %rcx,-0x220(%rbp)
  mov    %r8d,-0x224(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  -0x208(%rbp),%xmm1
  movsd  D29(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D90(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1d0(%rbp)
  movsd  -0x208(%rbp),%xmm1
  movsd  D29(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D90(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1c8(%rbp)
  movsd  -0x1c8(%rbp),%xmm0
  movsd  %xmm0,-0xf0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xe8(%rbp)
  movsd  -0x1d0(%rbp),%xmm0
  movq   D34(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0xe0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xd8(%rbp)
  movsd  D56(%rip),%xmm0
  movsd  %xmm0,-0xd0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xc8(%rbp)
  movsd  -0x1d0(%rbp),%xmm0
  movsd  %xmm0,-0xc0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0xb8(%rbp)
  movsd  -0x1c8(%rbp),%xmm0
  movsd  %xmm0,-0xb0(%rbp)
  movsd  -0x210(%rbp),%xmm1
  movsd  D29(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D90(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1c0(%rbp)
  movsd  -0x210(%rbp),%xmm1
  movsd  D29(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D90(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1b8(%rbp)
  movsd  D56(%rip),%xmm0
  movsd  %xmm0,-0xa0(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x98(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  -0x1b8(%rbp),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x1c0(%rbp),%xmm0
  movq   D34(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0x1c0(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x1b8(%rbp),%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0xa0(%rbp),%xmm1
  movsd  -0xf0(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x88(%rbp),%xmm2
  movsd  -0xe8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x70(%rbp),%xmm2
  movsd  -0xe0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  -0xf0(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x80(%rbp),%xmm2
  movsd  -0xe8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x68(%rbp),%xmm2
  movsd  -0xe0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x90(%rbp),%xmm1
  movsd  -0xf0(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x78(%rbp),%xmm2
  movsd  -0xe8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x60(%rbp),%xmm2
  movsd  -0xe0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0xa0(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x88(%rbp),%xmm2
  movsd  -0xd0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x70(%rbp),%xmm2
  movsd  -0xc8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x80(%rbp),%xmm2
  movsd  -0xd0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x68(%rbp),%xmm2
  movsd  -0xc8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x90(%rbp),%xmm1
  movsd  -0xd8(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x78(%rbp),%xmm2
  movsd  -0xd0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x60(%rbp),%xmm2
  movsd  -0xc8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0xa0(%rbp),%xmm1
  movsd  -0xc0(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x88(%rbp),%xmm2
  movsd  -0xb8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x70(%rbp),%xmm2
  movsd  -0xb0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x98(%rbp),%xmm1
  movsd  -0xc0(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x80(%rbp),%xmm2
  movsd  -0xb8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x68(%rbp),%xmm2
  movsd  -0xb0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x90(%rbp),%xmm1
  movsd  -0xc0(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x78(%rbp),%xmm2
  movsd  -0xb8(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x60(%rbp),%xmm2
  movsd  -0xb0(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movl   $0x0,-0x1d4(%rbp)
  jmp    L18
L21:
  mov    -0x1d4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x1f0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  D175(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x1f8(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1b0(%rbp)
  mov    -0x1d4(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x1f0(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  D175(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  -0x1f8(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x1a8(%rbp)
  mov    -0x218(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  -0x1a8(%rbp),%xmm0
  mov    -0x218(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x1b0(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x190(%rbp)
  mov    -0x218(%rbp),%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x1b0(%rbp),%xmm1
  mov    -0x218(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mulsd  -0x1a8(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x188(%rbp)
  mov    -0x218(%rbp),%rax
  movsd  0x10(%rax),%xmm0
  movsd  %xmm0,-0x180(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x190(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x48(%rbp),%xmm2
  movsd  -0x188(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x40(%rbp),%xmm2
  movsd  -0x180(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x130(%rbp)
  movsd  -0x38(%rbp),%xmm1
  movsd  -0x190(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x30(%rbp),%xmm2
  movsd  -0x188(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x28(%rbp),%xmm2
  movsd  -0x180(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x128(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  -0x190(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x18(%rbp),%xmm2
  movsd  -0x188(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x10(%rbp),%xmm2
  movsd  -0x180(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x120(%rbp)
  mov    -0x1d4(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x1e8(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  -0x130(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x170(%rbp)
  mov    -0x1d4(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x1e8(%rbp),%rax
  add    %rdx,%rax
  movsd  0x8(%rax),%xmm0
  movsd  -0x128(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x168(%rbp)
  mov    -0x1d4(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x1e8(%rbp),%rax
  add    %rdx,%rax
  movsd  0x10(%rax),%xmm0
  movsd  -0x120(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x160(%rbp)
  movsd  -0x170(%rbp),%xmm1
  movsd  -0x170(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x168(%rbp),%xmm2
  movsd  -0x168(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x160(%rbp),%xmm2
  movsd  -0x160(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  movsd  D56(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x1a0(%rbp)
  movsd  D56(%rip),%xmm0
  divsd  -0x200(%rbp),%xmm0
  movsd  %xmm0,-0x198(%rbp)
  movsd  -0x170(%rbp),%xmm0
  mulsd  -0x1a0(%rbp),%xmm0
  movsd  D56(%rip),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  -0x198(%rbp),%xmm0
  movsd  %xmm0,-0x110(%rbp)
  movsd  -0x168(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  mulsd  -0x1a0(%rbp),%xmm0
  movsd  %xmm0,-0x108(%rbp)
  movsd  -0x160(%rbp),%xmm0
  mulsd  -0x198(%rbp),%xmm0
  mulsd  -0x1a0(%rbp),%xmm0
  movsd  %xmm0,-0x100(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  -0x110(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm2
  movsd  -0x108(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x20(%rbp),%xmm2
  movsd  -0x100(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x150(%rbp)
  movsd  -0x48(%rbp),%xmm1
  movsd  -0x110(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x30(%rbp),%xmm2
  movsd  -0x108(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x18(%rbp),%xmm2
  movsd  -0x100(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x148(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  -0x110(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x28(%rbp),%xmm2
  movsd  -0x108(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x10(%rbp),%xmm2
  movsd  -0x100(%rbp),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x140(%rbp)
  movsd  -0x150(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x1a8(%rbp),%xmm1
  movsd  -0x148(%rbp),%xmm0
  mulsd  -0x1b0(%rbp),%xmm0
  mov    -0x1d4(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x220(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  -0x1b0(%rbp),%xmm0
  movq   D34(%rip),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  -0x150(%rbp),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x148(%rbp),%xmm0
  mulsd  -0x1a8(%rbp),%xmm0
  mov    -0x1d4(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x220(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rax)
  mov    -0x1d4(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x220(%rbp),%rax
  add    %rdx,%rax
  movsd  -0x140(%rbp),%xmm0
  movsd  %xmm0,0x10(%rax)
  addl   $0x1,-0x1d4(%rbp)
L18:
  mov    -0x1d4(%rbp),%eax
  cmp    -0x224(%rbp),%eax
  jl     L21
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L20
  call   <__stack_chk_fail@plt>
L20:
  leave
  ret
```

tool_calls: name: read_data
arguments: label: D29
type_infered: f64
name: read_data
arguments: label: D90
type_infered: f64
name: read_data
arguments: label: D34
type_infered: f32[2]
name: read_data
arguments: label: D56
type_infered: f64
name: read_data
arguments: label: D175
type_infered: byte[8]

tool: D29: 3.141592653589793

tool: D90: 180

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

tool: D56: 1

tool: {"D175": [0x39,0x9d,0x52,0xa2,0x46,0xdf,0x91,0x3f]}

assistant: ```c
void compute_gv(double xlylzl[][3], double omega[], double omegasign, double wvln, double wedge, double chi, double t[3], double gv[][3], int n) {
  double sc, cc, sw, cw, wmat[9], cmat[9], mat[9], u[3], d[3], v[3];
  double modyz, o[3], co, so, ds, k[3];
  int i;
  sw = sin(wedge * 3.14159265358979323846 / 180.0);
  cw = cos(wedge * 3.14159265358979323846 / 180.0);
  wmat[0] = cw;
  wmat[1] = 0.0;
  wmat[2] = -sw;
  wmat[3] = 0.;
  wmat[4] = 1.0;
  wmat[5] = 0.;
  wmat[6] = sw;
  wmat[7] = 0.0;
  wmat[8] = cw;
  sc = sin(chi * 3.14159265358979323846 / 180.0);
  cc = cos(chi * 3.14159265358979323846 / 180.0);
  cmat[0] = 1.;
  cmat[1] = 0.0;
  cmat[2] = 0.;
  cmat[3] = 0.;
  cmat[4] = cc;
  cmat[5] = -sc;
  cmat[6] = 0.;
  cmat[7] = sc;
  cmat[8] = cc;
  (mat)[0] = (cmat)[0] * (wmat)[0] + (cmat)[3] * (wmat)[1] + (cmat)[6] * (wmat)[2];
  (mat)[1] = (cmat)[1] * (wmat)[0] + (cmat)[4] * (wmat)[1] + (cmat)[7] * (wmat)[2];
  (mat)[2] = (cmat)[2] * (wmat)[0] + (cmat)[5] * (wmat)[1] + (cmat)[8] * (wmat)[2];
  (mat)[3] = (cmat)[0] * (wmat)[3] + (cmat)[3] * (wmat)[4] + (cmat)[6] * (wmat)[5];
  (mat)[4] = (cmat)[1] * (wmat)[3] + (cmat)[4] * (wmat)[4] + (cmat)[7] * (wmat)[5];
  (mat)[5] = (cmat)[2] * (wmat)[3] + (cmat)[5] * (wmat)[4] + (cmat)[8] * (wmat)[5];
  (mat)[6] = (cmat)[0] * (wmat)[6] + (cmat)[3] * (wmat)[7] + (cmat)[6] * (wmat)[8];
  (mat)[7] = (cmat)[1] * (wmat)[6] + (cmat)[4] * (wmat)[7] + (cmat)[7] * (wmat)[8];
  (mat)[8] = (cmat)[2] * (wmat)[6] + (cmat)[5] * (wmat)[7] + (cmat)[8] * (wmat)[8];
  ;
#pragma omp parallel for private(so, co, u, o, d, modyz, ds, v, k)
  for (i = 0; i < n; i++) {
    so = sin(3.14159265358979323846 / 180.0 * omega[i] * omegasign);
    co = cos(3.14159265358979323846 / 180.0 * omega[i] * omegasign);
    u[0] = co * t[0] - so * t[1];
    u[1] = so * t[0] + co * t[1];
    u[2] = t[2];
    (o)[0] = (mat)[0] * (u)[0] + (mat)[1] * (u)[1] + (mat)[2] * (u)[2];
    (o)[1] = (mat)[3] * (u)[0] + (mat)[4] * (u)[1] + (mat)[5] * (u)[2];
    (o)[2] = (mat)[6] * (u)[0] + (mat)[7] * (u)[1] + (mat)[8] * (u)[2];
    ;
    (d)[0] = (xlylzl[i])[0] - (o)[0];
    (d)[1] = (xlylzl[i])[1] - (o)[1];
    (d)[2] = (xlylzl[i])[2] - (o)[2];
    ;
    modyz = 1. / sqrt(d[0] * d[0] + d[1] * d[1] + d[2] * d[2]);
    ds = 1. / wvln;
    k[0] = ds * (d[0] * modyz - 1.);
    k[1] = ds * d[1] * modyz;
    k[2] = ds * d[2] * modyz;
    (v)[0] = (mat)[0] * (k)[0] + (mat)[3] * (k)[1] + (mat)[6] * (k)[2];
    (v)[1] = (mat)[1] * (k)[0] + (mat)[4] * (k)[1] + (mat)[7] * (k)[2];
    (v)[2] = (mat)[2] * (k)[0] + (mat)[5] * (k)[1] + (mat)[8] * (k)[2];
    ;
    gv[i][0] = co * v[0] + so * v[1];
    gv[i][1] = -so * v[0] + co * v[1];
    gv[i][2] = v[2];
  }
}
```