// 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>::do_step_impl<RocketStage, std::array<double, 14ul>, std::array<double, 14ul>, std::array<double, 14ul>, std::array<double, 14ul>>(RocketStage, std::array<double, 14ul> const&, std::array<double, 14ul> const&, double, std::array<double, 14ul>&, std::array<double, 14ul>&, double)
void do_step_impl( System system , const StateIn &in , const DerivIn &dxdt_in , time_type t ,
            StateOut &out , DerivOut &dxdt_out , time_type dt )
    {
        const value_type a2 = static_cast<value_type> ( 1 ) / static_cast<value_type>( 5 );
        const value_type a3 = static_cast<value_type> ( 3 ) / static_cast<value_type> ( 10 );
        const value_type a4 = static_cast<value_type> ( 4 ) / static_cast<value_type> ( 5 );
        const value_type a5 = static_cast<value_type> ( 8 )/static_cast<value_type> ( 9 );

        const value_type b21 = static_cast<value_type> ( 1 ) / static_cast<value_type> ( 5 );

        const value_type b31 = static_cast<value_type> ( 3 ) / static_cast<value_type>( 40 );
        const value_type b32 = static_cast<value_type> ( 9 ) / static_cast<value_type>( 40 );

        const value_type b41 = static_cast<value_type> ( 44 ) / static_cast<value_type> ( 45 );
        const value_type b42 = static_cast<value_type> ( -56 ) / static_cast<value_type> ( 15 );
        const value_type b43 = static_cast<value_type> ( 32 ) / static_cast<value_type> ( 9 );

        const value_type b51 = static_cast<value_type> ( 19372 ) / static_cast<value_type>( 6561 );
        const value_type b52 = static_cast<value_type> ( -25360 ) / static_cast<value_type> ( 2187 );
        const value_type b53 = static_cast<value_type> ( 64448 ) / static_cast<value_type>( 6561 );
        const value_type b54 = static_cast<value_type> ( -212 ) / static_cast<value_type>( 729 );

        const value_type b61 = static_cast<value_type> ( 9017 ) / static_cast<value_type>( 3168 );
        const value_type b62 = static_cast<value_type> ( -355 ) / static_cast<value_type>( 33 );
        const value_type b63 = static_cast<value_type> ( 46732 ) / static_cast<value_type>( 5247 );
        const value_type b64 = static_cast<value_type> ( 49 ) / static_cast<value_type>( 176 );
        const value_type b65 = static_cast<value_type> ( -5103 ) / static_cast<value_type>( 18656 );

        const value_type c1 = static_cast<value_type> ( 35 ) / static_cast<value_type>( 384 );
        const value_type c3 = static_cast<value_type> ( 500 ) / static_cast<value_type>( 1113 );
        const value_type c4 = static_cast<value_type> ( 125 ) / static_cast<value_type>( 192 );
        const value_type c5 = static_cast<value_type> ( -2187 ) / static_cast<value_type>( 6784 );
        const value_type c6 = static_cast<value_type> ( 11 ) / static_cast<value_type>( 84 );

        typename odeint::unwrap_reference< System >::type &sys = system;

        m_k_x_tmp_resizer.adjust_size( in , detail::bind( &stepper_type::template resize_k_x_tmp_impl<StateIn> , detail::ref( *this ) , detail::_1 ) );

        //m_x_tmp = x + dt*b21*dxdt
        stepper_base_type::m_algebra.for_each3( m_x_tmp.m_v , in , dxdt_in ,
                typename operations_type::template scale_sum2< value_type , time_type >( 1.0 , dt*b21 ) );

        sys( m_x_tmp.m_v , m_k2.m_v , t + dt*a2 );
        // m_x_tmp = x + dt*b31*dxdt + dt*b32*m_k2
        stepper_base_type::m_algebra.for_each4( m_x_tmp.m_v , in , dxdt_in , m_k2.m_v ,
                typename operations_type::template scale_sum3< value_type , time_type , time_type >( 1.0 , dt*b31 , dt*b32 ));

        sys( m_x_tmp.m_v , m_k3.m_v , t + dt*a3 );
        // m_x_tmp = x + dt * (b41*dxdt + b42*m_k2 + b43*m_k3)
        stepper_base_type::m_algebra.for_each5( m_x_tmp.m_v , in , dxdt_in , m_k2.m_v , m_k3.m_v ,
                typename operations_type::template scale_sum4< value_type , time_type , time_type , time_type >( 1.0 , dt*b41 , dt*b42 , dt*b43 ));

        sys( m_x_tmp.m_v, m_k4.m_v , t + dt*a4 );
        stepper_base_type::m_algebra.for_each6( m_x_tmp.m_v , in , dxdt_in , m_k2.m_v , m_k3.m_v , m_k4.m_v ,
                typename operations_type::template scale_sum5< value_type , time_type , time_type , time_type , time_type >( 1.0 , dt*b51 , dt*b52 , dt*b53 , dt*b54 ));

        sys( m_x_tmp.m_v , m_k5.m_v , t + dt*a5 );
        stepper_base_type::m_algebra.for_each7( m_x_tmp.m_v , in , dxdt_in , m_k2.m_v , m_k3.m_v , m_k4.m_v , m_k5.m_v ,
                typename operations_type::template scale_sum6< value_type , time_type , time_type , time_type , time_type , time_type >( 1.0 , dt*b61 , dt*b62 , dt*b63 , dt*b64 , dt*b65 ));

        sys( m_x_tmp.m_v , m_k6.m_v , t + dt );
        stepper_base_type::m_algebra.for_each7( out , in , dxdt_in , m_k3.m_v , m_k4.m_v , m_k5.m_v , m_k6.m_v ,
                typename operations_type::template scale_sum6< value_type , time_type , time_type , time_type , time_type , time_type >( 1.0 , dt*c1 , dt*c3 , dt*c4 , dt*c5 , dt*c6 ));

        // the new derivative
        sys( out , dxdt_out , t + dt );
    }

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x98, %rsp
movaps %xmm1, 0x30(%rsp)
movq %r9, 0x58(%rsp)
movq %r8, %r14
movsd %xmm0, 0x20(%rsp)
movq %rcx, %r12
movq %rdx, %r13
movq %rsi, %rbx
movq %rdi, %rbp
addq $0x388, %rdi            # imm = 0x388
leaq 0x412(%rip), %rax       # 0x14c40
leaq 0x80(%rsp), %rdx
movq %rax, (%rdx)
andq $0x0, 0x8(%rdx)
movq %rbp, 0x10(%rdx)
movq %r13, %rsi
callq 0x147d8
leaq 0x78(%rbp), %r15
movsd 0xb6ea(%rip), %xmm1     # 0x1ff40
mulsd 0x30(%rsp), %xmm1
movsd %xmm1, 0x48(%rsp)
movsd 0xa7f6(%rip), %xmm0     # 0x1f060
movq %r15, %rdi
movq %r13, %rsi
movq %r12, %rdx
callq 0x14c44
movq %rbp, 0x50(%rsp)
leaq 0xe8(%rbp), %rdx
movsd 0x48(%rsp), %xmm0
addsd 0x20(%rsp), %xmm0
movq (%rbx), %rax
movq %rbx, 0x28(%rsp)
movq %rbx, %rdi
movq %r15, %rsi
movq %rdx, 0x48(%rsp)
callq *(%rax)
movapd 0x30(%rsp), %xmm4
movsd 0xb6cd(%rip), %xmm0     # 0x1ff80
mulsd %xmm4, %xmm0
movsd 0xb6c9(%rip), %xmm1     # 0x1ff88
mulsd %xmm4, %xmm1
xorl %eax, %eax
cmpq $0xe, %rax
je 0x148f6
movsd (%r12,%rax,8), %xmm2
mulsd %xmm0, %xmm2
addsd (%r13,%rax,8), %xmm2
movsd 0x70(%r15,%rax,8), %xmm3
mulsd %xmm1, %xmm3
addsd %xmm2, %xmm3
movsd %xmm3, (%r15,%rax,8)
incq %rax
jmp 0x148c5
movq 0x50(%rsp), %rbp
leaq 0x158(%rbp), %rbx
movsd 0xb686(%rip), %xmm0     # 0x1ff90
mulsd %xmm4, %xmm0
addsd 0x20(%rsp), %xmm0
movq 0x28(%rsp), %rdi
movq (%rdi), %rax
movq %r15, %rsi
movq %rbx, %rdx
callq *(%rax)
movapd 0x30(%rsp), %xmm0
unpcklpd %xmm0, %xmm0            # xmm0 = xmm0[0,0]
mulpd 0xb72a(%rip), %xmm0     # 0x20060
movsd 0xb65a(%rip), %xmm1     # 0x1ff98
mulsd 0x30(%rsp), %xmm1
movabsq $0x3ff0000000000000, %rax # imm = 0x3FF0000000000000
movq %rax, 0x60(%rsp)
movupd %xmm0, 0x68(%rsp)
movsd %xmm1, 0x78(%rsp)
movups 0x60(%rsp), %xmm0
movups 0x70(%rsp), %xmm1
movups %xmm1, 0x10(%rsp)
movups %xmm0, (%rsp)
movq %r15, %rdi
movq %r13, %rsi
movq %r12, %rdx
movq 0x48(%rsp), %rcx
movq %rbx, %r8
callq 0x14c6e
addq $0x1c8, %rbp            # imm = 0x1C8
movsd 0xb609(%rip), %xmm0     # 0x1ffa0
mulsd 0x30(%rsp), %xmm0
addsd 0x20(%rsp), %xmm0
movq 0x28(%rsp), %rdi
movq (%rdi), %rax
movq %r15, %rsi
movq %rbp, %rdx
callq *(%rax)
movapd 0x30(%rsp), %xmm6
movsd 0xb5e7(%rip), %xmm0     # 0x1ffa8
mulsd %xmm6, %xmm0
movsd 0xb5e3(%rip), %xmm1     # 0x1ffb0
mulsd %xmm6, %xmm1
movsd 0xb5df(%rip), %xmm2     # 0x1ffb8
mulsd %xmm6, %xmm2
movsd 0xb5db(%rip), %xmm3     # 0x1ffc0
mulsd %xmm6, %xmm3
xorl %eax, %eax
cmpq $0xe, %rax
je 0x14a3d
movsd (%r12,%rax,8), %xmm4
mulsd %xmm0, %xmm4
addsd (%r13,%rax,8), %xmm4
movsd -0xe0(%rbp,%rax,8), %xmm5
mulsd %xmm1, %xmm5
addsd %xmm4, %xmm5
movsd -0x70(%rbp,%rax,8), %xmm4
mulsd %xmm2, %xmm4
addsd %xmm5, %xmm4
movsd (%rbp,%rax,8), %xmm5
mulsd %xmm3, %xmm5
addsd %xmm4, %xmm5
movsd %xmm5, -0x150(%rbp,%rax,8)
incq %rax
jmp 0x149eb
movq 0x50(%rsp), %rax
leaq 0x238(%rax), %rbp
movsd 0xb577(%rip), %xmm0     # 0x1ffc8
mulsd %xmm6, %xmm0
addsd 0x20(%rsp), %xmm0
movq 0x28(%rsp), %rdi
movq (%rdi), %rax
movq %r15, %rsi
movq %rbp, %rdx
callq *(%rax)
movapd 0x30(%rsp), %xmm7
movsd 0xb557(%rip), %xmm0     # 0x1ffd0
mulsd %xmm7, %xmm0
movsd 0xb553(%rip), %xmm1     # 0x1ffd8
mulsd %xmm7, %xmm1
movsd 0xb54f(%rip), %xmm2     # 0x1ffe0
mulsd %xmm7, %xmm2
movsd 0xb54b(%rip), %xmm3     # 0x1ffe8
mulsd %xmm7, %xmm3
movsd 0xb547(%rip), %xmm4     # 0x1fff0
mulsd %xmm7, %xmm4
xorl %eax, %eax
cmpq $0xe, %rax
je 0x14b12
movsd (%r12,%rax,8), %xmm5
mulsd %xmm0, %xmm5
addsd (%r13,%rax,8), %xmm5
movsd -0x150(%rbp,%rax,8), %xmm6
mulsd %xmm1, %xmm6
addsd %xmm5, %xmm6
movsd -0xe0(%rbp,%rax,8), %xmm5
mulsd %xmm2, %xmm5
addsd %xmm6, %xmm5
movsd -0x70(%rbp,%rax,8), %xmm6
mulsd %xmm3, %xmm6
addsd %xmm5, %xmm6
movsd (%rbp,%rax,8), %xmm5
mulsd %xmm4, %xmm5
addsd %xmm6, %xmm5
movsd %xmm5, -0x1c0(%rbp,%rax,8)
incq %rax
jmp 0x14aaf
movq 0x50(%rsp), %rdx
addq $0x2a8, %rdx            # imm = 0x2A8
movsd 0x20(%rsp), %xmm0
addsd %xmm7, %xmm0
movq 0x28(%rsp), %rdi
movq (%rdi), %rax
movq %r15, %rsi
movsd %xmm0, 0x20(%rsp)
callq *(%rax)
movapd 0x30(%rsp), %xmm6
movsd 0xb3c7(%rip), %xmm0     # 0x1ff10
mulsd %xmm6, %xmm0
movsd 0xb3c3(%rip), %xmm1     # 0x1ff18
mulsd %xmm6, %xmm1
movsd 0xb3bf(%rip), %xmm2     # 0x1ff20
mulsd %xmm6, %xmm2
movsd 0xb3bb(%rip), %xmm3     # 0x1ff28
mulsd %xmm6, %xmm3
mulsd 0xb3b7(%rip), %xmm6     # 0x1ff30
xorl %eax, %eax
cmpq $0xe, %rax
je 0x14bda
movsd (%r12,%rax,8), %xmm4
mulsd %xmm0, %xmm4
addsd (%r13,%rax,8), %xmm4
movsd (%rbx,%rax,8), %xmm5
mulsd %xmm1, %xmm5
addsd %xmm4, %xmm5
movsd 0x70(%rbx,%rax,8), %xmm4
mulsd %xmm2, %xmm4
addsd %xmm5, %xmm4
movsd 0xe0(%rbx,%rax,8), %xmm5
mulsd %xmm3, %xmm5
addsd %xmm4, %xmm5
movsd 0x150(%rbx,%rax,8), %xmm4
mulsd %xmm6, %xmm4
addsd %xmm5, %xmm4
movsd %xmm4, (%r14,%rax,8)
incq %rax
jmp 0x14b7b
movq 0x28(%rsp), %rdi
movq (%rdi), %rax
movq %r14, %rsi
movq 0x58(%rsp), %rdx
movsd 0x20(%rsp), %xmm0
callq *(%rax)
addq $0x98, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq