Node.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 NODE_H
#define NODE_H

#include <Domain/Tag.h>

class DomainBase;
enum class OutputType;
enum class DOF : unsigned short;

struct NodeData {
    unsigned num_dof = 0; // number of DoFs

    vec coordinate; // coordinates of the node

    uvec original_dof;  // original indices
    uvec reordered_dof; // renumbered indices

    vec current_resistance;     // current resistance
    vec current_damping_force;  // current damping force
    vec current_inertial_force; // current inertial force
    vec current_displacement;   // current displacement
    vec current_velocity;       // current velocity
    vec current_acceleration;   // current acceleration

    vec incre_resistance;     // incremental resistance
    vec incre_damping_force;  // incremental damping force
    vec incre_inertial_force; // incremental inertial force
    vec incre_displacement;   // incremental displacement
    vec incre_velocity;       // incremental velocity
    vec incre_acceleration;   // incremental acceleration

    vec trial_resistance;     // trial resistance
    vec trial_damping_force;  // trial damping force
    vec trial_inertial_force; // trial inertial force
    vec trial_displacement;   // trial displacement
    vec trial_velocity;       // trial velocity
    vec trial_acceleration;   // trial acceleration
};

class Node final : protected NodeData, public Tag {
    bool initialized = false;

    std::mutex node_mutex;

    std::vector<DOF> dof_identifier;

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

    void initialize(const shared_ptr<DomainBase>&);

    void set_initialized(bool);
    [[nodiscard]] bool get_initialized() const;

    void set_dof_number(unsigned);
    [[nodiscard]] unsigned get_dof_number() const;

    void set_dof_identifier(const std::vector<DOF>&);
    [[nodiscard]] const std::vector<DOF>& get_dof_identifier() const;

    void set_original_dof(unsigned&);
    void set_original_dof(const uvec&);
    [[nodiscard]] const uvec& get_original_dof() const;

    void set_reordered_dof(const uvec&);
    [[nodiscard]] const uvec& get_reordered_dof() const;

    void set_coordinate(const vec&);
    [[nodiscard]] const vec& get_coordinate() const;

    void set_current_resistance(const vec&);
    void set_current_damping_force(const vec&);
    void set_current_inertial_force(const vec&);
    void set_current_displacement(const vec&);
    void set_current_velocity(const vec&);
    void set_current_acceleration(const vec&);

    void set_incre_resistance(const vec&);
    void set_incre_damping_force(const vec&);
    void set_incre_inertial_force(const vec&);
    void set_incre_displacement(const vec&);
    void set_incre_velocity(const vec&);
    void set_incre_acceleration(const vec&);

    void set_trial_resistance(const vec&);
    void set_trial_damping_force(const vec&);
    void set_trial_inertial_force(const vec&);
    void set_trial_displacement(const vec&);
    void set_trial_velocity(const vec&);
    void set_trial_acceleration(const vec&);

    [[nodiscard]] const vec& get_current_resistance() const;
    [[nodiscard]] const vec& get_current_damping_force() const;
    [[nodiscard]] const vec& get_current_inertial_force() const;
    [[nodiscard]] const vec& get_current_displacement() const;
    [[nodiscard]] const vec& get_current_velocity() const;
    [[nodiscard]] const vec& get_current_acceleration() const;

    [[nodiscard]] const vec& get_incre_resistance() const;
    [[nodiscard]] const vec& get_incre_damping_force() const;
    [[nodiscard]] const vec& get_incre_inertial_force() const;
    [[nodiscard]] const vec& get_incre_displacement() const;
    [[nodiscard]] const vec& get_incre_velocity() const;
    [[nodiscard]] const vec& get_incre_acceleration() const;

    [[nodiscard]] const vec& get_trial_resistance() const;
    [[nodiscard]] const vec& get_trial_damping_force() const;
    [[nodiscard]] const vec& get_trial_inertial_force() const;
    [[nodiscard]] const vec& get_trial_displacement() const;
    [[nodiscard]] const vec& get_trial_velocity() const;
    [[nodiscard]] const vec& get_trial_acceleration() const;

    void update_current_resistance(const vec&);
    void update_current_damping_force(const vec&);
    void update_current_inertial_force(const vec&);
    void update_current_displacement(const vec&);
    void update_current_velocity(const vec&);
    void update_current_acceleration(const vec&);

    void update_incre_resistance(const vec&);
    void update_incre_damping_force(const vec&);
    void update_incre_inertial_force(const vec&);
    void update_incre_displacement(const vec&);
    void update_incre_velocity(const vec&);
    void update_incre_acceleration(const vec&);

    void update_trial_resistance(const vec&);
    void update_trial_damping_force(const vec&);
    void update_trial_inertial_force(const vec&);
    void update_trial_displacement(const vec&);
    void update_trial_velocity(const vec&);
    void update_trial_acceleration(const vec&);

    void update_current_status(const vec&);
    void update_current_status(const vec&, const vec&);
    void update_current_status(const vec&, const vec&, const vec&);

    void update_incre_status(const vec&);
    void update_incre_status(const vec&, const vec&);
    void update_incre_status(const vec&, const vec&, const vec&);

    void update_trial_status(const vec&);
    void update_trial_status(const vec&, const vec&);
    void update_trial_status(const vec&, const vec&, const vec&);

    void commit_status();
    void reset_status();
    void clear_status();

    [[nodiscard]] std::vector<vec> record(OutputType) const;

    void print() override;
};

#endif