// File: /workspace/github_cpp_makefile_star2_O2/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 0xc018(%rip), %xmm4     # 0x1fe58
mulsd %xmm0, %xmm4
movsd 0xc013(%rip), %xmm8     # 0x1fe60
addsd %xmm4, %xmm8
mulsd 0xc00d(%rip), %xmm8     # 0x1fe68
divsd 0xc00c(%rip), %xmm8     # 0x1fe70
movsd 0xc00c(%rip), %xmm6     # 0x1fe78
mulsd %xmm0, %xmm6
addsd 0xc008(%rip), %xmm6     # 0x1fe80
mulsd 0xb1e8(%rip), %xmm6     # 0x1f068
divsd 0xc000(%rip), %xmm6     # 0x1fe88
addsd 0xc000(%rip), %xmm4     # 0x1fe90
mulsd 0xc000(%rip), %xmm4     # 0x1fe98
divsd 0xc000(%rip), %xmm4     # 0x1fea0
movsd 0xc000(%rip), %xmm7     # 0x1fea8
mulsd %xmm0, %xmm7
addsd 0xbffc(%rip), %xmm7     # 0x1feb0
mulsd 0xbffc(%rip), %xmm7     # 0x1feb8
divsd 0xbffc(%rip), %xmm7     # 0x1fec0
movsd 0xbffb(%rip), %xmm9     # 0x1fec8
mulsd %xmm0, %xmm9
addsd 0xbff5(%rip), %xmm9     # 0x1fed0
mulsd 0xbff4(%rip), %xmm9     # 0x1fed8
divsd 0xbff3(%rip), %xmm9     # 0x1fee0
movsd 0xbff2(%rip), %xmm10    # 0x1fee8
mulsd %xmm0, %xmm10
addsd 0xbfec(%rip), %xmm10    # 0x1fef0
mulsd 0xbf3b(%rip), %xmm10    # 0x1fe48
divsd 0xbfe2(%rip), %xmm10    # 0x1fef8
movsd 0xbfe2(%rip), %xmm3     # 0x1ff00
addsd %xmm0, %xmm3
movapd %xmm0, %xmm11
mulsd %xmm0, %xmm11
movapd %xmm0, %xmm5
addsd %xmm0, %xmm5
movsd 0xbfcc(%rip), %xmm1     # 0x1ff08
subsd %xmm5, %xmm1
mulsd %xmm11, %xmm1
mulsd %xmm3, %xmm11
mulsd %xmm3, %xmm0
mulsd %xmm3, %xmm0
mulsd %xmm11, %xmm3
mulsd %xmm3, %xmm8
movsd 0xbfac(%rip), %xmm5     # 0x1ff10
mulsd %xmm1, %xmm5
addsd %xmm8, %xmm5
addsd %xmm0, %xmm5
mulsd %xmm3, %xmm6
movsd 0xbf9b(%rip), %xmm0     # 0x1ff18
mulsd %xmm1, %xmm0
addsd %xmm6, %xmm0
mulsd %xmm3, %xmm4
movsd 0xbf8f(%rip), %xmm6     # 0x1ff20
mulsd %xmm1, %xmm6
addsd %xmm4, %xmm6
mulsd %xmm3, %xmm7
movsd 0xbf83(%rip), %xmm4     # 0x1ff28
mulsd %xmm1, %xmm4
addsd %xmm7, %xmm4
mulsd %xmm3, %xmm9
mulsd 0xbf76(%rip), %xmm1     # 0x1ff30
addsd %xmm9, %xmm1
mulsd %xmm10, %xmm3
addsd %xmm11, %xmm3
mulsd %xmm2, %xmm5
mulsd %xmm2, %xmm0
mulsd %xmm2, %xmm6
mulsd %xmm2, %xmm4
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm3
xorl %eax, %eax
cmpq $0x70, %rax
je 0x14054
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 %xmm4, %xmm7
addsd %xmm2, %xmm7
movsd 0x2a8(%rdi,%rax), %xmm2
mulsd %xmm1, %xmm2
addsd %xmm7, %xmm2
movsd (%r9,%rax), %xmm7
mulsd %xmm3, %xmm7
addsd %xmm2, %xmm7
movsd %xmm7, (%rsi,%rax)
addq $0x8, %rax
jmp 0x13fe3
retq
nop