// File: /workspace/github_cpp_makefile_star2_O3/istellartech[P]OpenTsiolkovsky/./boost/numeric/odeint/stepper/runge_kutta_dopri5.hpp
// Function: void boost::numeric::odeint::runge_kutta_dopri5<std::array<double, 14ul>, double, std::array<double, 14ul>, double, boost::numeric::odeint::array_algebra, boost::numeric::odeint::default_operations, boost::numeric::odeint::initially_resizer>::calc_state<std::array<double, 14ul>, std::array<double, 14ul>, std::array<double, 14ul>, std::array<double, 14ul>, std::array<double, 14ul>>(double, std::array<double, 14ul>&, std::array<double, 14ul> const&, std::array<double, 14ul> const&, double, std::array<double, 14ul> const&, std::array<double, 14ul> const&, double) const
void calc_state( time_type t , StateOut &x ,
                     const StateIn1 &x_old , const DerivIn1 &deriv_old , time_type t_old ,
                     const StateIn2 & /* x_new */ , const DerivIn2 &deriv_new , time_type t_new ) const
    {
        const value_type b1 = static_cast<value_type> ( 35 ) / static_cast<value_type>( 384 );
        const value_type b3 = static_cast<value_type> ( 500 ) / static_cast<value_type>( 1113 );
        const value_type b4 = static_cast<value_type> ( 125 ) / static_cast<value_type>( 192 );
        const value_type b5 = static_cast<value_type> ( -2187 ) / static_cast<value_type>( 6784 );
        const value_type b6 = static_cast<value_type> ( 11 ) / static_cast<value_type>( 84 );

        const time_type dt = ( t_new - t_old );
        const value_type theta = ( t - t_old ) / dt;
        const value_type X1 = static_cast< value_type >( 5 ) * ( static_cast< value_type >( 2558722523LL ) - static_cast< value_type >( 31403016 ) * theta ) / static_cast< value_type >( 11282082432LL );
        const value_type X3 = static_cast< value_type >( 100 ) * ( static_cast< value_type >( 882725551 ) - static_cast< value_type >( 15701508 ) * theta ) / static_cast< value_type >( 32700410799LL );
        const value_type X4 = static_cast< value_type >( 25 ) * ( static_cast< value_type >( 443332067 ) - static_cast< value_type >( 31403016 ) * theta ) / static_cast< value_type >( 1880347072LL ) ;
        const value_type X5 = static_cast< value_type >( 32805 ) * ( static_cast< value_type >( 23143187 ) - static_cast< value_type >( 3489224 ) * theta ) / static_cast< value_type >( 199316789632LL );
        const value_type X6 = static_cast< value_type >( 55 ) * ( static_cast< value_type >( 29972135 ) - static_cast< value_type >( 7076736 ) * theta ) / static_cast< value_type >( 822651844 );
        const value_type X7 = static_cast< value_type >( 10 ) * ( static_cast< value_type >( 7414447 ) - static_cast< value_type >( 829305 ) * theta ) / static_cast< value_type >( 29380423 );

        const value_type theta_m_1 = theta - static_cast< value_type >( 1 );
        const value_type theta_sq = theta * theta;
        const value_type A = theta_sq * ( static_cast< value_type >( 3 ) - static_cast< value_type >( 2 ) * theta );
        const value_type B = theta_sq * theta_m_1;
        const value_type C = theta_sq * theta_m_1 * theta_m_1;
        const value_type D = theta * theta_m_1 * theta_m_1;

        const value_type b1_theta = A * b1 - C * X1 + D;
        const value_type b3_theta = A * b3 + C * X3;
        const value_type b4_theta = A * b4 - C * X4;
        const value_type b5_theta = A * b5 + C * X5;
        const value_type b6_theta = A * b6 - C * X6;
        const value_type b7_theta = B + C * X7;

        // const state_type &k1 = *m_old_deriv;
        // const state_type &k3 = dopri5().m_k3;
        // const state_type &k4 = dopri5().m_k4;
        // const state_type &k5 = dopri5().m_k5;
        // const state_type &k6 = dopri5().m_k6;
        // const state_type &k7 = *m_current_deriv;

        stepper_base_type::m_algebra.for_each8( x , x_old , deriv_old , m_k3.m_v , m_k4.m_v , m_k5.m_v , m_k6.m_v , deriv_new ,
                typename operations_type::template scale_sum7< value_type , time_type , time_type , time_type , time_type , time_type , time_type >( 1.0 , dt * b1_theta , dt * b3_theta , dt * b4_theta , dt * b5_theta , dt * b6_theta , dt * b7_theta ) );
    }

// Assembly:
subsd %xmm1, %xmm2
subsd %xmm1, %xmm0
divsd %xmm2, %xmm0
movsd 0xc518(%rip), %xmm3     # 0x25d78
mulsd %xmm0, %xmm3
movsd 0xc513(%rip), %xmm8     # 0x25d80
addsd %xmm3, %xmm8
mulsd 0xc50d(%rip), %xmm8     # 0x25d88
divsd 0xc50c(%rip), %xmm8     # 0x25d90
movsd 0xc50c(%rip), %xmm6     # 0x25d98
mulsd %xmm0, %xmm6
addsd 0xc508(%rip), %xmm6     # 0x25da0
mulsd 0xb7d8(%rip), %xmm6     # 0x25078
divsd 0xc500(%rip), %xmm6     # 0x25da8
addsd 0xc500(%rip), %xmm3     # 0x25db0
mulsd 0xc500(%rip), %xmm3     # 0x25db8
divsd 0xc500(%rip), %xmm3     # 0x25dc0
movsd 0xc500(%rip), %xmm7     # 0x25dc8
mulsd %xmm0, %xmm7
addsd 0xc4fc(%rip), %xmm7     # 0x25dd0
mulsd 0xc4fc(%rip), %xmm7     # 0x25dd8
divsd 0xc4fc(%rip), %xmm7     # 0x25de0
movsd 0xc4fb(%rip), %xmm9     # 0x25de8
mulsd %xmm0, %xmm9
addsd 0xc4f5(%rip), %xmm9     # 0x25df0
mulsd 0xc4f4(%rip), %xmm9     # 0x25df8
divsd 0xc4f3(%rip), %xmm9     # 0x25e00
movsd 0xc4f2(%rip), %xmm10    # 0x25e08
mulsd %xmm0, %xmm10
addsd 0xc4ec(%rip), %xmm10    # 0x25e10
mulsd 0xc43b(%rip), %xmm10    # 0x25d68
divsd 0xc4e2(%rip), %xmm10    # 0x25e18
movsd 0xc4e2(%rip), %xmm4     # 0x25e20
addsd %xmm0, %xmm4
movapd %xmm0, %xmm11
mulsd %xmm0, %xmm11
movapd %xmm0, %xmm5
addsd %xmm0, %xmm5
movsd 0xc4cc(%rip), %xmm1     # 0x25e28
subsd %xmm5, %xmm1
mulsd %xmm11, %xmm1
mulsd %xmm4, %xmm11
mulsd %xmm4, %xmm0
mulsd %xmm4, %xmm0
mulsd %xmm11, %xmm4
mulsd %xmm4, %xmm8
movsd 0xc4ac(%rip), %xmm5     # 0x25e30
mulsd %xmm1, %xmm5
addsd %xmm8, %xmm5
addsd %xmm0, %xmm5
mulsd %xmm4, %xmm6
movsd 0xc49b(%rip), %xmm0     # 0x25e38
mulsd %xmm1, %xmm0
addsd %xmm6, %xmm0
mulsd %xmm4, %xmm3
movsd 0xc48f(%rip), %xmm6     # 0x25e40
mulsd %xmm1, %xmm6
addsd %xmm3, %xmm6
mulsd %xmm4, %xmm7
movsd 0xc483(%rip), %xmm3     # 0x25e48
mulsd %xmm1, %xmm3
addsd %xmm7, %xmm3
mulsd %xmm4, %xmm9
mulsd 0xc476(%rip), %xmm1     # 0x25e50
addsd %xmm9, %xmm1
mulsd %xmm10, %xmm4
addsd %xmm11, %xmm4
mulsd %xmm2, %xmm5
mulsd %xmm2, %xmm0
mulsd %xmm2, %xmm6
mulsd %xmm2, %xmm3
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm4
xorl %eax, %eax
movsd (%rcx,%rax), %xmm2
mulsd %xmm5, %xmm2
addsd (%rdx,%rax), %xmm2
movsd 0x158(%rdi,%rax), %xmm7
mulsd %xmm0, %xmm7
addsd %xmm2, %xmm7
movsd 0x1c8(%rdi,%rax), %xmm2
mulsd %xmm6, %xmm2
addsd %xmm7, %xmm2
movsd 0x238(%rdi,%rax), %xmm7
mulsd %xmm3, %xmm7
addsd %xmm2, %xmm7
movsd 0x2a8(%rdi,%rax), %xmm2
mulsd %xmm1, %xmm2
addsd %xmm7, %xmm2
movsd (%r9,%rax), %xmm7
mulsd %xmm4, %xmm7
addsd %xmm2, %xmm7
movsd %xmm7, (%rsi,%rax)
addq $0x8, %rax
cmpq $0x70, %rax
jne 0x19a03
retq
nop