SpringSim.cpp

#include "SpringSim.h"

bool SpringSim::advance() {
    // perform time integration with different integrators

	// use p_cube, m_spring, m_dt, m_gravity

	Eigen::Vector3d v = p_cube->getLinearVelocity();
	Eigen::Vector3d p = p_cube->getPosition();

    // TODO

	// note that it is required to update both m_sptring.end and p_cube's position
    switch (m_method) {
        case 0: // analytical solution
        {
            break;
        }
        case 1: // explicit euler
			
            break;

        case 2: // symplectic euler
        
            break;

        case 3: // midpoint

            break;

        default:
            std::cerr << m_method << " is not a valid integrator method."
                        << std::endl;
    }

	// update spring end position
	m_spring.end = p_cube->getPosition();


    // advance m_time
    m_time += m_dt;
    m_step++;

    // log
    if ((m_step % m_log_frequency) == 0) {
        m_trajectories.back().push_back(p_cube->getPosition());
    }

    return false;
}