// File: /workspace/github_c_makefile_star2_O0/PoulinV[P]AxiCLASS/build/../source/perturbations.c
// Function: perturbations_rsa_delta_and_theta
int perturbations_rsa_delta_and_theta(
                                      struct precision * ppr,
                                      struct background * pba,
                                      struct thermodynamics * pth,
                                      struct perturbations * ppt,
                                      double k,
                                      double * y,
                                      double a_prime_over_a,
                                      double * pvecthermo,
                                      struct perturbations_workspace * ppw,
                                      ErrorMsg error_message
                                      ) {
  /* - define local variables */

  double k2;

  k2 = k*k;

  class_test(ppw->approx[ppw->index_ap_rsa] == (int)rsa_off,
             "this function should not have been called now, bug was introduced",
             ppt->error_message,
             ppt->error_message);

  // formulas below TBC for curvaturema

  /* newtonian gauge */
  if (ppt->gauge == newtonian) {

    if (ppr->radiation_streaming_approximation == rsa_null) {
      ppw->rsa_delta_g = 0.;
      ppw->rsa_theta_g = 0.;
    }
    else {
      ppw->rsa_delta_g = -4.*y[ppw->pv->index_pt_phi];
      ppw->rsa_theta_g = 6.*ppw->pvecmetric[ppw->index_mt_phi_prime];
    }

    if (ppr->radiation_streaming_approximation == rsa_MD_with_reio) {

      ppw->rsa_delta_g +=
        -4./k2*ppw->pvecthermo[pth->index_th_dkappa]*y[ppw->pv->index_pt_theta_b];

      ppw->rsa_theta_g +=
        3./k2*(ppw->pvecthermo[pth->index_th_ddkappa]*y[ppw->pv->index_pt_theta_b]
               +ppw->pvecthermo[pth->index_th_dkappa]*
               (-a_prime_over_a*y[ppw->pv->index_pt_theta_b]
                +ppw->pvecthermo[pth->index_th_cb2]*k2*y[ppw->pv->index_pt_delta_b]
                +k2*y[ppw->pv->index_pt_phi]));
    }

    if (pba->has_ur == _TRUE_) {

      if (ppr->radiation_streaming_approximation == rsa_null) {
        ppw->rsa_delta_ur = 0.;
        ppw->rsa_theta_ur = 0.;
      }
      else {
        ppw->rsa_delta_ur = -4.*y[ppw->pv->index_pt_phi];
        ppw->rsa_theta_ur = 6.*ppw->pvecmetric[ppw->index_mt_phi_prime];
      }
    }
  }

  /* synchronous gauge */
  if (ppt->gauge == synchronous) {

    if (ppr->radiation_streaming_approximation == rsa_null) {
      ppw->rsa_delta_g = 0.;
      ppw->rsa_theta_g = 0.;
    }
    else {

      ppw->rsa_delta_g = 4./k2*(a_prime_over_a*ppw->pvecmetric[ppw->index_mt_h_prime]
                                -k2*y[ppw->pv->index_pt_eta]);
      ppw->rsa_theta_g = -0.5*ppw->pvecmetric[ppw->index_mt_h_prime];
    }

    if (ppr->radiation_streaming_approximation == rsa_MD_with_reio) {

      ppw->rsa_delta_g +=
        -4./k2*ppw->pvecthermo[pth->index_th_dkappa]*(y[ppw->pv->index_pt_theta_b]+0.5*ppw->pvecmetric[ppw->index_mt_h_prime]);

      ppw->rsa_theta_g +=
        3./k2*(ppw->pvecthermo[pth->index_th_ddkappa]*
               (y[ppw->pv->index_pt_theta_b]
                +0.5*ppw->pvecmetric[ppw->index_mt_h_prime])
               +ppw->pvecthermo[pth->index_th_dkappa]*
               (-a_prime_over_a*y[ppw->pv->index_pt_theta_b]
                + ppw->pvecthermo[pth->index_th_cb2]*k2*y[ppw->pv->index_pt_delta_b]
                -a_prime_over_a*ppw->pvecmetric[ppw->index_mt_h_prime]
                +k2*y[ppw->pv->index_pt_eta]));
    }

    if (pba->has_ur == _TRUE_) {

      if (ppr->radiation_streaming_approximation == rsa_null) {
        ppw->rsa_delta_ur = 0.;
        ppw->rsa_theta_ur = 0.;
      }
      else {
        ppw->rsa_delta_ur = 4./k2*(a_prime_over_a*ppw->pvecmetric[ppw->index_mt_h_prime]
                                   -k2*y[ppw->pv->index_pt_eta]);
        ppw->rsa_theta_ur = -0.5*ppw->pvecmetric[ppw->index_mt_h_prime];
      }
    }
  }

  /* update total delta and theta given rsa approximation results */

  ppw->delta_rho += ppw->pvecback[pba->index_bg_rho_g]*ppw->rsa_delta_g;
  ppw->delta_p += 1./3.*ppw->pvecback[pba->index_bg_rho_g]*ppw->rsa_delta_g;
  ppw->rho_plus_p_theta += 4./3.*ppw->pvecback[pba->index_bg_rho_g]*ppw->rsa_theta_g;

  if (pba->has_ur == _TRUE_) {
    ppw->delta_rho += ppw->pvecback[pba->index_bg_rho_ur]*ppw->rsa_delta_ur;
    ppw->delta_p += 1./3.*ppw->pvecback[pba->index_bg_rho_ur]*ppw->rsa_delta_ur;
    ppw->rho_plus_p_theta += 4./3.*ppw->pvecback[pba->index_bg_rho_ur]*ppw->rsa_theta_ur;
  }

  return _SUCCESS_;

}

// Assembly:
subq $0x1058, %rsp           # imm = 0x1058
movq 0x1068(%rsp), %rax
movq 0x1060(%rsp), %rax
movq %rdi, 0x1048(%rsp)
movq %rsi, 0x1040(%rsp)
movq %rdx, 0x1038(%rsp)
movq %rcx, 0x1030(%rsp)
movsd %xmm0, 0x1028(%rsp)
movq %r8, 0x1020(%rsp)
movsd %xmm1, 0x1018(%rsp)
movq %r9, 0x1010(%rsp)
movsd 0x1028(%rsp), %xmm0
mulsd 0x1028(%rsp), %xmm0
movsd %xmm0, 0x1008(%rsp)
movq 0x1060(%rsp), %rax
movq 0x180(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x164(%rcx), %rcx
cmpl $0x0, (%rax,%rcx,4)
jne 0x10731d
leaq 0x800(%rsp), %rdi
movq 0x1030(%rsp), %rsi
addq $0x6c8a, %rsi           # imm = 0x6C8A
movq 0x1030(%rsp), %rdx
addq $0x6c8a, %rdx           # imm = 0x6C8A
movb $0x0, %al
callq 0x5b9c0
movq %rsp, %rdi
leaq 0x800(%rsp), %rcx
leaq 0xcada7(%rip), %rsi     # 0x1d2082
leaq 0xf4e30(%rip), %rdx     # 0x1fc112
movb $0x0, %al
callq 0x5b9c0
movq %rsp, %r8
leaq 0xf4e4e(%rip), %rdi     # 0x1fc141
leaq 0xcad3a(%rip), %rsi     # 0x1d2034
leaq 0xf4e82(%rip), %rdx     # 0x1fc183
movl $0x2cff, %ecx           # imm = 0x2CFF
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x1054(%rsp)
jmp 0x107d23
movq 0x1030(%rsp), %rax
cmpl $0x0, 0x83c(%rax)
jne 0x107673
movq 0x1048(%rsp), %rax
cmpl $0x0, 0x620(%rax)
jne 0x10736b
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xa8(%rax)
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xb0(%rax)
jmp 0x1073d8
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0xd4(%rcx), %rcx
movsd 0xce132(%rip), %xmm0    # 0x1d54c0
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xa8(%rax)
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x4(%rcx), %rcx
movsd 0xcc2c5(%rip), %xmm0    # 0x1d3688
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xb0(%rax)
movq 0x1048(%rsp), %rax
cmpl $0x2, 0x620(%rax)
jne 0x1075b1
movsd 0xce0cb(%rip), %xmm0    # 0x1d54c0
divsd 0x1008(%rsp), %xmm0
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xe80(%rcx), %rcx
mulsd (%rax,%rcx,8), %xmm0
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x2c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x1060(%rsp), %rax
movsd 0xa8(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0xa8(%rax)
movsd 0xcc235(%rip), %xmm0    # 0x1d3698
divsd 0x1008(%rsp), %xmm0
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xe88(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x2c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm3
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xe80(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm1
movsd 0x1018(%rsp), %xmm4
movq %xmm4, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm5
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x2c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm6
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xef0(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm4
mulsd 0x1008(%rsp), %xmm4
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x28(%rcx), %rcx
mulsd (%rax,%rcx,8), %xmm4
mulsd %xmm6, %xmm5
addsd %xmm4, %xmm5
movsd 0x1008(%rsp), %xmm4
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0xd4(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm6
mulsd %xmm6, %xmm4
addsd %xmm5, %xmm4
mulsd %xmm4, %xmm1
mulsd %xmm3, %xmm2
addsd %xmm1, %xmm2
movq 0x1060(%rsp), %rax
movsd 0xb0(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0xb0(%rax)
movq 0x1040(%rsp), %rax
movswl 0x4c4(%rax), %eax
cmpl $0x1, %eax
jne 0x107671
movq 0x1048(%rsp), %rax
cmpl $0x0, 0x620(%rax)
jne 0x107602
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xb8(%rax)
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xc0(%rax)
jmp 0x10766f
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0xd4(%rcx), %rcx
movsd 0xcde9b(%rip), %xmm0    # 0x1d54c0
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xb8(%rax)
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x4(%rcx), %rcx
movsd 0xcc02e(%rip), %xmm0    # 0x1d3688
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xc0(%rax)
jmp 0x107671
jmp 0x107673
movq 0x1030(%rsp), %rax
cmpl $0x1, 0x83c(%rax)
jne 0x107b24
movq 0x1048(%rsp), %rax
cmpl $0x0, 0x620(%rax)
jne 0x1076c4
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xa8(%rax)
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xb0(%rax)
jmp 0x10778c
movsd 0xcddfc(%rip), %xmm0    # 0x1d54c8
divsd 0x1008(%rsp), %xmm0
movsd 0x1018(%rsp), %xmm1
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm3
movsd 0x1008(%rsp), %xmm2
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0xd0(%rcx), %rcx
mulsd (%rax,%rcx,8), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xa8(%rax)
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x8(%rcx), %rcx
movsd 0xcbf29(%rip), %xmm0    # 0x1d36a0
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xb0(%rax)
movq 0x1048(%rsp), %rax
cmpl $0x2, 0x620(%rax)
jne 0x107a04
movsd 0xcdd17(%rip), %xmm0    # 0x1d54c0
divsd 0x1008(%rsp), %xmm0
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xe80(%rcx), %rcx
mulsd (%rax,%rcx,8), %xmm0
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x2c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm2
movsd 0xcb474(%rip), %xmm3    # 0x1d2c88
mulsd %xmm3, %xmm2
addsd %xmm1, %xmm2
movq 0x1060(%rsp), %rax
movsd 0xa8(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0xa8(%rax)
movsd 0xcbe54(%rip), %xmm0    # 0x1d3698
divsd 0x1008(%rsp), %xmm0
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xe88(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm2
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x2c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm3
movsd 0xcb3d9(%rip), %xmm4    # 0x1d2c88
mulsd %xmm4, %xmm3
addsd %xmm1, %xmm3
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xe80(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm1
movsd 0x1018(%rsp), %xmm4
movq %xmm4, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm4
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x2c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm6
movq 0x1060(%rsp), %rax
movq 0x38(%rax), %rax
movq 0x1038(%rsp), %rcx
movslq 0xef0(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm5
mulsd 0x1008(%rsp), %xmm5
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0x28(%rcx), %rcx
mulsd (%rax,%rcx,8), %xmm5
mulsd %xmm6, %xmm4
addsd %xmm5, %xmm4
movsd 0x1018(%rsp), %xmm5
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm6
movq %xmm5, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm5
mulsd %xmm6, %xmm5
addsd %xmm4, %xmm5
movsd 0x1008(%rsp), %xmm4
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0xd0(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm6
mulsd %xmm6, %xmm4
addsd %xmm5, %xmm4
mulsd %xmm4, %xmm1
mulsd %xmm3, %xmm2
addsd %xmm1, %xmm2
movq 0x1060(%rsp), %rax
movsd 0xb0(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0xb0(%rax)
movq 0x1040(%rsp), %rax
movswl 0x4c4(%rax), %eax
cmpl $0x1, %eax
jne 0x107b22
movq 0x1048(%rsp), %rax
cmpl $0x0, 0x620(%rax)
jne 0x107a58
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xb8(%rax)
movq 0x1060(%rsp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, 0xc0(%rax)
jmp 0x107b20
movsd 0xcda68(%rip), %xmm0    # 0x1d54c8
divsd 0x1008(%rsp), %xmm0
movsd 0x1018(%rsp), %xmm1
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm3
movsd 0x1008(%rsp), %xmm2
movq 0x1020(%rsp), %rax
movq 0x1060(%rsp), %rcx
movq 0x48(%rcx), %rcx
movslq 0xd0(%rcx), %rcx
mulsd (%rax,%rcx,8), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xb8(%rax)
movq 0x1060(%rsp), %rax
movq 0x40(%rax), %rax
movq 0x1060(%rsp), %rcx
movslq 0x8(%rcx), %rcx
movsd 0xcbb95(%rip), %xmm0    # 0x1d36a0
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd %xmm0, 0xc0(%rax)
jmp 0x107b22
jmp 0x107b24
movq 0x1060(%rsp), %rax
movq 0x30(%rax), %rax
movq 0x1040(%rsp), %rcx
movslq 0x38c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd 0xa8(%rax), %xmm2
movq 0x1060(%rsp), %rax
movsd 0x50(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x50(%rax)
movq 0x1060(%rsp), %rax
movq 0x30(%rax), %rax
movq 0x1040(%rsp), %rcx
movslq 0x38c(%rcx), %rcx
movsd 0xcd907(%rip), %xmm0    # 0x1d5498
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd 0xa8(%rax), %xmm2
movq 0x1060(%rsp), %rax
movsd 0x68(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x68(%rax)
movq 0x1060(%rsp), %rax
movq 0x30(%rax), %rax
movq 0x1040(%rsp), %rcx
movslq 0x38c(%rcx), %rcx
movsd 0xcd92d(%rip), %xmm0    # 0x1d5510
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd 0xb0(%rax), %xmm2
movq 0x1060(%rsp), %rax
movsd 0x58(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x58(%rax)
movq 0x1040(%rsp), %rax
movswl 0x4c4(%rax), %eax
cmpl $0x1, %eax
jne 0x107d18
movq 0x1060(%rsp), %rax
movq 0x30(%rax), %rax
movq 0x1040(%rsp), %rcx
movslq 0x3ac(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd 0xb8(%rax), %xmm2
movq 0x1060(%rsp), %rax
movsd 0x50(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x50(%rax)
movq 0x1060(%rsp), %rax
movq 0x30(%rax), %rax
movq 0x1040(%rsp), %rcx
movslq 0x3ac(%rcx), %rcx
movsd 0xcd801(%rip), %xmm0    # 0x1d5498
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd 0xb8(%rax), %xmm2
movq 0x1060(%rsp), %rax
movsd 0x68(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x68(%rax)
movq 0x1060(%rsp), %rax
movq 0x30(%rax), %rax
movq 0x1040(%rsp), %rcx
movslq 0x3ac(%rcx), %rcx
movsd 0xcd827(%rip), %xmm0    # 0x1d5510
mulsd (%rax,%rcx,8), %xmm0
movq 0x1060(%rsp), %rax
movsd 0xc0(%rax), %xmm2
movq 0x1060(%rsp), %rax
movsd 0x58(%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x58(%rax)
movl $0x0, 0x1054(%rsp)
movl 0x1054(%rsp), %eax
addq $0x1058, %rsp           # imm = 0x1058
retq
nopw %cs:(%rax,%rax)