// File: /workspace/github_cpp_makefile_star2_O1/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 0xc734(%rip), %xmm3     # 0x27d40
mulsd %xmm0, %xmm3
movsd 0xc72f(%rip), %xmm8     # 0x27d48
addsd %xmm3, %xmm8
mulsd 0xc729(%rip), %xmm8     # 0x27d50
divsd 0xc728(%rip), %xmm8     # 0x27d58
movsd 0xc728(%rip), %xmm6     # 0x27d60
mulsd %xmm0, %xmm6
addsd 0xc724(%rip), %xmm6     # 0x27d68
mulsd 0xb9cc(%rip), %xmm6     # 0x27018
divsd 0xc71c(%rip), %xmm6     # 0x27d70
addsd 0xc71c(%rip), %xmm3     # 0x27d78
mulsd 0xc71c(%rip), %xmm3     # 0x27d80
divsd 0xc71c(%rip), %xmm3     # 0x27d88
movsd 0xc71c(%rip), %xmm7     # 0x27d90
mulsd %xmm0, %xmm7
addsd 0xc718(%rip), %xmm7     # 0x27d98
mulsd 0xc718(%rip), %xmm7     # 0x27da0
divsd 0xc718(%rip), %xmm7     # 0x27da8
movsd 0xc717(%rip), %xmm9     # 0x27db0
mulsd %xmm0, %xmm9
addsd 0xc711(%rip), %xmm9     # 0x27db8
mulsd 0xc710(%rip), %xmm9     # 0x27dc0
divsd 0xc70f(%rip), %xmm9     # 0x27dc8
movsd 0xc70e(%rip), %xmm10    # 0x27dd0
mulsd %xmm0, %xmm10
addsd 0xc708(%rip), %xmm10    # 0x27dd8
mulsd 0xc657(%rip), %xmm10    # 0x27d30
divsd 0xc6fe(%rip), %xmm10    # 0x27de0
movsd 0xc6fe(%rip), %xmm4     # 0x27de8
addsd %xmm0, %xmm4
movapd %xmm0, %xmm11
mulsd %xmm0, %xmm11
movapd %xmm0, %xmm5
addsd %xmm0, %xmm5
movsd 0xc6e8(%rip), %xmm1     # 0x27df0
subsd %xmm5, %xmm1
mulsd %xmm11, %xmm1
mulsd %xmm4, %xmm11
mulsd %xmm4, %xmm0
mulsd %xmm4, %xmm0
mulsd %xmm11, %xmm4
mulsd %xmm4, %xmm8
movsd 0xc6c8(%rip), %xmm5     # 0x27df8
mulsd %xmm1, %xmm5
addsd %xmm8, %xmm5
addsd %xmm0, %xmm5
mulsd %xmm4, %xmm6
movsd 0xc6b7(%rip), %xmm0     # 0x27e00
mulsd %xmm1, %xmm0
addsd %xmm6, %xmm0
mulsd %xmm4, %xmm3
movsd 0xc6ab(%rip), %xmm6     # 0x27e08
mulsd %xmm1, %xmm6
addsd %xmm3, %xmm6
mulsd %xmm4, %xmm7
movsd 0xc69f(%rip), %xmm3     # 0x27e10
mulsd %xmm1, %xmm3
addsd %xmm7, %xmm3
mulsd %xmm4, %xmm9
mulsd 0xc692(%rip), %xmm1     # 0x27e18
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 0x1b7af
retq
nop