velocity_function.h

//
//  Deformabel Simplicial Complex (DSC) method
//  Copyright (C) 2013  Technical University of Denmark
//
//  This program is free software: you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or
//  (at your option) any later version.
//
//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  See licence.txt for a copy of the GNU General Public License.

#pragma once

#include <chrono>
#include <ctime>
#include "DSC.h"

namespace DSC {
    
    class VelocityFunc
    {
        double compute_time = 0.;
        double deform_time = 0.;
        
        double total_compute_time = 0.;
        double total_deform_time = 0.;
        
    protected:
        int time_step = 0;
        int MAX_TIME_STEPS;
        int deform_time_steps = 10;

        double VELOCITY; // Determines the distance each interface vertex moves at each iteration.
        double ACCURACY; // Determines the accuracy of the final result.
        
        std::vector<vec3> pos_old;
        
    public:
        VelocityFunc(double velocity, double accuracy, int max_time_steps);

        virtual ~VelocityFunc();

        virtual std::string get_name() const;

        int get_time_step() const;

        virtual void set_max_time_steps(int max_time_steps);

        double get_velocity() const;

        virtual void set_velocity(double vel);

        double get_accuracy() const;

        virtual void set_accuracy(double acc);

        double get_deform_time() const;

        double get_compute_time() const;

        double get_total_deform_time() const;

        double get_total_compute_time() const;
        
    protected:
        void update_compute_time(const std::chrono::time_point<std::chrono::system_clock>& start_time);

        void update_deform_time(const std::chrono::time_point<std::chrono::system_clock>& start_time);
        
        virtual void deform(DeformableSimplicialComplex& dsc);
        
    public:
        virtual bool is_motion_finished(DeformableSimplicialComplex& dsc);
        
        void take_time_step(DeformableSimplicialComplex& dsc);
        
        virtual void test(DeformableSimplicialComplex& dsc);
        
    };
    
}