Integrator.h

Go to the documentation of this file.

/*******************************************************************************
 * Copyright (C) 2017-2022 Theodore Chang
 *
 * 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.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 ******************************************************************************/
#ifndef INTEGRATOR_H
#define INTEGRATOR_H

#include <Domain/Tag.h>

class DomainBase;

class Integrator : public Tag {
    bool time_step_switch = true;

    weak_ptr<DomainBase> database;

public:
    explicit Integrator(unsigned = 0);
    Integrator(const Integrator&) = delete;            // copy forbidden
    Integrator(Integrator&&) = delete;                 // move forbidden
    Integrator& operator=(const Integrator&) = delete; // assign forbidden
    Integrator& operator=(Integrator&&) = delete;      // assign forbidden
    ~Integrator() override;

    void set_domain(const weak_ptr<DomainBase>&);
    [[nodiscard]] const weak_ptr<DomainBase>& get_domain() const;

    virtual int initialize();

    // ! some multistep integrators may require fixed time step for some consecutive sub-steps
    void set_time_step_switch(bool);
    [[nodiscard]] bool allow_to_change_time_step() const;

    [[nodiscard]] virtual int process_load();
    [[nodiscard]] virtual int process_constraint();
    [[nodiscard]] virtual int process_criterion();
    [[nodiscard]] virtual int process_modifier();
    [[nodiscard]] virtual int process_load_resistance();
    [[nodiscard]] virtual int process_constraint_resistance();

    void record() const;

    virtual void assemble_resistance();
    virtual void assemble_matrix();

    virtual vec get_force_residual();
    virtual vec get_displacement_residual();
    virtual vec get_auxiliary_residual();
    virtual sp_mat get_reference_load();
    virtual sp_mat get_auxiliary_stiffness();

    virtual void update_load();
    virtual void update_constraint();

    virtual void update_trial_load_factor(double);
    virtual void update_trial_load_factor(const vec&);
    virtual void update_trial_displacement(const vec&);

    void update_trial_time(double);
    void update_incre_time(double);
    virtual int update_trial_status();
    virtual int update_incre_status();

    virtual int update_internal(const mat&);

    mat solve(const mat&);
    mat solve(const sp_mat&);
    mat solve(mat&&);
    mat solve(sp_mat&&);
    virtual int solve(mat&, const mat&);
    virtual int solve(mat&, const sp_mat&);
    virtual int solve(mat&, mat&&);
    virtual int solve(mat&, sp_mat&&);

    virtual void erase_machine_error() const;

    virtual void stage_and_commit_status();

    virtual void stage_status();
    virtual void commit_status();
    virtual void clear_status();
    virtual void reset_status();

    virtual void update_parameter(double);
    virtual void update_compatibility() const;

    virtual vec from_incre_velocity(const vec&, const uvec&);     // obtain target displacement from increment of velocity
    virtual vec from_incre_acceleration(const vec&, const uvec&); // obtain target displacement from increment of acceleration
    vec from_total_velocity(const vec&, const uvec&);
    vec from_total_acceleration(const vec&, const uvec&);
    vec from_incre_velocity(double, const uvec&);
    vec from_incre_acceleration(double, const uvec&);
    vec from_total_velocity(double, const uvec&);
    vec from_total_acceleration(double, const uvec&);
};

#endif