// File: /workspace/github_c_makefile_star2_O2/YunyiQian[P]multitel3/radiats.c
// Function: fk2mtg
void fk2mtg(float az, int npt, float *fk[12], float *mt[18]) {
  int i;
  float sa, ca, sa2, ca2;
  az = az*DEG2RAD;
  sa = sin(az); ca = cos(az);
  sa2 = sin(2*az); ca2 = cos(2*az);
  // mxx.zr = EX/3-DD/6-SS*cos2Az/2
  for(i=0;i<npt;i++) mt[0][i] =  fk[9][i]/3. - fk[0][i]/6. - fk[6][i]*ca2/2;
  for(i=0;i<npt;i++) mt[1][i] = fk[10][i]/3. - fk[1][i]/6. - fk[7][i]*ca2/2;
  // mxx.t = -SS*sin2Az/2
  for(i=0;i<npt;i++) mt[2][i] = -fk[8][i]*sa2/2;
  // mxy.zr = -SS*sin2Az
  for(i=0;i<npt;i++) mt[3][i] = -fk[6][i]*sa2;
  for(i=0;i<npt;i++) mt[4][i] = -fk[7][i]*sa2;
  // mxy.t =  SS*cos2Az
  for(i=0;i<npt;i++) mt[5][i] = fk[8][i]*ca2;
  // mxz.zr = -DS*cosAz
  for(i=0;i<npt;i++) mt[6][i] = -fk[3][i]*ca;
  for(i=0;i<npt;i++) mt[7][i] = -fk[4][i]*ca;
  // mxz.t = -DS*sinAz
  for(i=0;i<npt;i++) mt[8][i] = -fk[5][i]*sa;
  // myy.zr = Ex/3-DD/6+SS*cos2Az/2
  for(i=0;i<npt;i++) mt[9][i]  =  fk[9][i]/3. - fk[0][i]/6. + fk[6][i]*ca2/2;
  for(i=0;i<npt;i++) mt[10][i] = fk[10][i]/3. - fk[1][i]/6. + fk[7][i]*ca2/2;
  // myy.t = -mxx.t
  for(i=0;i<npt;i++) mt[11][i] = -mt[2][i];
  // myz.zr = -DS*sinAz
  for(i=0;i<npt;i++) mt[12][i] = -fk[3][i]*sa;
  for(i=0;i<npt;i++) mt[13][i] = -fk[4][i]*sa;
  // myz.t =  DS*cosAz
  for(i=0;i<npt;i++) mt[14][i] = fk[5][i]*ca;
  // mzz.zr = Ex/3+DD/3
  for(i=0;i<npt;i++) mt[15][i] =  (fk[9][i] + fk[0][i])/3.;
  for(i=0;i<npt;i++) mt[16][i] = (fk[10][i] + fk[1][i])/3.;
  // mzz.t = 0
  for(i=0;i<npt;i++) mt[17][i] = 0.;
  return;
}

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
subq $0x20, %rsp
movq %rdx, %rbx
movq %rsi, %r15
cvtss2sd %xmm0, %xmm0
mulsd 0x220a(%rip), %xmm0     # 0x6e90
movl %edi, %r14d
cvtsd2ss %xmm0, %xmm0
movss %xmm0, 0x8(%rsp)
cvtss2sd %xmm0, %xmm0
movsd %xmm0, (%rsp)
callq 0x11d0
movsd %xmm0, 0x18(%rsp)
movsd (%rsp), %xmm0
callq 0x10d0
movsd %xmm0, (%rsp)
movss 0x8(%rsp), %xmm0
addss %xmm0, %xmm0
cvtss2sd %xmm0, %xmm0
movsd %xmm0, 0x10(%rsp)
callq 0x11d0
movsd %xmm0, 0x8(%rsp)
movsd 0x10(%rsp), %xmm0
callq 0x10d0
cvtsd2ss %xmm0, %xmm1
xorl %eax, %eax
testl %r14d, %r14d
cmovlel %eax, %r14d
movapd 0x21db(%rip), %xmm2     # 0x6ed0
movss 0x21eb(%rip), %xmm3     # 0x6ee8
cmpq %rax, %r14
je 0x4d53
movq 0x48(%r15), %rcx
movq (%r15), %rdx
movss (%rdx,%rax,4), %xmm0
movss (%rcx,%rax,4), %xmm4
unpcklps %xmm0, %xmm4            # xmm4 = xmm4[0],xmm0[0],xmm4[1],xmm0[1]
cvtps2pd %xmm4, %xmm0
movq 0x30(%r15), %rcx
divpd %xmm2, %xmm0
movapd %xmm0, %xmm4
unpckhpd %xmm0, %xmm4            # xmm4 = xmm4[1],xmm0[1]
subsd %xmm4, %xmm0
movss (%rcx,%rax,4), %xmm4
mulss %xmm1, %xmm4
mulss %xmm3, %xmm4
cvtss2sd %xmm4, %xmm4
addsd %xmm0, %xmm4
cvtsd2ss %xmm4, %xmm0
movq (%rbx), %rcx
movss %xmm0, (%rcx,%rax,4)
incq %rax
jmp 0x4cfd
cvtsd2ss (%rsp), %xmm0
xorl %eax, %eax
cmpq %rax, %r14
je 0x4db2
movq 0x50(%r15), %rcx
movq 0x8(%r15), %rdx
movss (%rdx,%rax,4), %xmm4
movss (%rcx,%rax,4), %xmm5
unpcklps %xmm4, %xmm5            # xmm5 = xmm5[0],xmm4[0],xmm5[1],xmm4[1]
cvtps2pd %xmm5, %xmm4
movq 0x38(%r15), %rcx
divpd %xmm2, %xmm4
movapd %xmm4, %xmm5
unpckhpd %xmm4, %xmm5            # xmm5 = xmm5[1],xmm4[1]
subsd %xmm5, %xmm4
movss (%rcx,%rax,4), %xmm5
mulss %xmm1, %xmm5
mulss %xmm3, %xmm5
cvtss2sd %xmm5, %xmm5
addsd %xmm4, %xmm5
cvtsd2ss %xmm5, %xmm4
movq 0x8(%rbx), %rcx
movss %xmm4, (%rcx,%rax,4)
incq %rax
jmp 0x4d5a
cvtsd2ss 0x8(%rsp), %xmm4
xorl %eax, %eax
cmpq %rax, %r14
je 0x4dde
movq 0x40(%r15), %rcx
movss (%rcx,%rax,4), %xmm5
mulss %xmm4, %xmm5
mulss %xmm3, %xmm5
movq 0x10(%rbx), %rcx
movss %xmm5, (%rcx,%rax,4)
incq %rax
jmp 0x4dba
xorl %eax, %eax
movaps 0x1289(%rip), %xmm3     # 0x6070
cmpq %rax, %r14
je 0x4e0a
movq 0x30(%r15), %rcx
movss (%rcx,%rax,4), %xmm5
xorps %xmm3, %xmm5
mulss %xmm4, %xmm5
movq 0x18(%rbx), %rcx
movss %xmm5, (%rcx,%rax,4)
incq %rax
jmp 0x4de7
xorl %eax, %eax
cmpq %rax, %r14
je 0x4e2f
movq 0x38(%r15), %rcx
movss (%rcx,%rax,4), %xmm5
xorps %xmm3, %xmm5
mulss %xmm4, %xmm5
movq 0x20(%rbx), %rcx
movss %xmm5, (%rcx,%rax,4)
incq %rax
jmp 0x4e0c
xorl %eax, %eax
cmpq %rax, %r14
je 0x4e51
movq 0x40(%r15), %rcx
movss (%rcx,%rax,4), %xmm4
mulss %xmm1, %xmm4
movq 0x28(%rbx), %rcx
movss %xmm4, (%rcx,%rax,4)
incq %rax
jmp 0x4e31
xorl %eax, %eax
cmpq %rax, %r14
je 0x4e76
movq 0x18(%r15), %rcx
movss (%rcx,%rax,4), %xmm4
xorps %xmm3, %xmm4
mulss %xmm0, %xmm4
movq 0x30(%rbx), %rcx
movss %xmm4, (%rcx,%rax,4)
incq %rax
jmp 0x4e53
xorl %eax, %eax
cmpq %rax, %r14
je 0x4e9b
movq 0x20(%r15), %rcx
movss (%rcx,%rax,4), %xmm4
xorps %xmm3, %xmm4
mulss %xmm0, %xmm4
movq 0x38(%rbx), %rcx
movss %xmm4, (%rcx,%rax,4)
incq %rax
jmp 0x4e78
cvtsd2ss 0x18(%rsp), %xmm4
xorl %eax, %eax
cmpq %rax, %r14
je 0x4ec6
movq 0x28(%r15), %rcx
movss (%rcx,%rax,4), %xmm5
xorps %xmm3, %xmm5
mulss %xmm4, %xmm5
movq 0x40(%rbx), %rcx
movss %xmm5, (%rcx,%rax,4)
incq %rax
jmp 0x4ea3
xorl %eax, %eax
movss 0x1b18(%rip), %xmm5     # 0x69e8
cmpq %rax, %r14
je 0x4f27
movq 0x48(%r15), %rcx
movq (%r15), %rdx
movss (%rdx,%rax,4), %xmm6
movss (%rcx,%rax,4), %xmm7
unpcklps %xmm6, %xmm7            # xmm7 = xmm7[0],xmm6[0],xmm7[1],xmm6[1]
cvtps2pd %xmm7, %xmm6
movq 0x30(%r15), %rcx
divpd %xmm2, %xmm6
movapd %xmm6, %xmm7
unpckhpd %xmm6, %xmm7            # xmm7 = xmm7[1],xmm6[1]
subsd %xmm7, %xmm6
movss (%rcx,%rax,4), %xmm7
mulss %xmm1, %xmm7
mulss %xmm5, %xmm7
cvtss2sd %xmm7, %xmm7
addsd %xmm6, %xmm7
cvtsd2ss %xmm7, %xmm6
movq 0x48(%rbx), %rcx
movss %xmm6, (%rcx,%rax,4)
incq %rax
jmp 0x4ed0
xorl %eax, %eax
cmpq %rax, %r14
je 0x4f81
movq 0x50(%r15), %rcx
movq 0x8(%r15), %rdx
movss (%rdx,%rax,4), %xmm6
movss (%rcx,%rax,4), %xmm7
unpcklps %xmm6, %xmm7            # xmm7 = xmm7[0],xmm6[0],xmm7[1],xmm6[1]
cvtps2pd %xmm7, %xmm6
movq 0x38(%r15), %rcx
divpd %xmm2, %xmm6
movapd %xmm6, %xmm7
unpckhpd %xmm6, %xmm7            # xmm7 = xmm7[1],xmm6[1]
subsd %xmm7, %xmm6
movss (%rcx,%rax,4), %xmm7
mulss %xmm1, %xmm7
mulss %xmm5, %xmm7
cvtss2sd %xmm7, %xmm7
addsd %xmm6, %xmm7
cvtsd2ss %xmm7, %xmm6
movq 0x50(%rbx), %rcx
movss %xmm6, (%rcx,%rax,4)
incq %rax
jmp 0x4f29
xorl %eax, %eax
cmpq %rax, %r14
je 0x4fa2
movq 0x10(%rbx), %rcx
movq 0x58(%rbx), %rdx
movss (%rcx,%rax,4), %xmm1
xorps %xmm3, %xmm1
movss %xmm1, (%rdx,%rax,4)
incq %rax
jmp 0x4f83
xorl %eax, %eax
cmpq %rax, %r14
je 0x4fc7
movq 0x18(%r15), %rcx
movss (%rcx,%rax,4), %xmm1
xorps %xmm3, %xmm1
mulss %xmm4, %xmm1
movq 0x60(%rbx), %rcx
movss %xmm1, (%rcx,%rax,4)
incq %rax
jmp 0x4fa4
xorl %eax, %eax
cmpq %rax, %r14
je 0x4fec
movq 0x20(%r15), %rcx
movss (%rcx,%rax,4), %xmm1
xorps %xmm3, %xmm1
mulss %xmm4, %xmm1
movq 0x68(%rbx), %rcx
movss %xmm1, (%rcx,%rax,4)
incq %rax
jmp 0x4fc9
xorl %eax, %eax
cmpq %rax, %r14
je 0x500e
movq 0x28(%r15), %rcx
movss (%rcx,%rax,4), %xmm1
mulss %xmm0, %xmm1
movq 0x70(%rbx), %rcx
movss %xmm1, (%rcx,%rax,4)
incq %rax
jmp 0x4fee
xorl %eax, %eax
movss 0x1ecc(%rip), %xmm0     # 0x6ee4
cmpq %rax, %r14
je 0x5040
movq (%r15), %rcx
movq 0x48(%r15), %rdx
movss (%rdx,%rax,4), %xmm1
addss (%rcx,%rax,4), %xmm1
divss %xmm0, %xmm1
movq 0x78(%rbx), %rcx
movss %xmm1, (%rcx,%rax,4)
incq %rax
jmp 0x5018
xorl %eax, %eax
cmpq %rax, %r14
je 0x506e
movq 0x8(%r15), %rcx
movq 0x50(%r15), %rdx
movss (%rdx,%rax,4), %xmm1
addss (%rcx,%rax,4), %xmm1
divss %xmm0, %xmm1
movq 0x80(%rbx), %rcx
movss %xmm1, (%rcx,%rax,4)
incq %rax
jmp 0x5042
xorl %eax, %eax
cmpq %rax, %r14
je 0x5085
movq 0x88(%rbx), %rcx
andl $0x0, (%rcx,%rax,4)
incq %rax
jmp 0x5070
addq $0x20, %rsp
popq %rbx
popq %r14
popq %r15
retq