suanPan
WilsonPenzienNewmark.h
Go to the documentation of this file.
1/*******************************************************************************
2 * Copyright (C) 2017-2024 Theodore Chang
3 *
4 * This program is free software: you can redistribute it and/or modify
5 * it under the terms of the GNU General Public License as published by
6 * the Free Software Foundation, either version 3 of the License, or
7 * (at your option) any later version.
8 *
9 * This program is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU General Public License for more details.
13 *
14 * You should have received a copy of the GNU General Public License
15 * along with this program. If not, see <http://www.gnu.org/licenses/>.
16 ******************************************************************************/
28#ifndef WILSONPENZIENNEWMARK_H
29#define WILSONPENZIENNEWMARK_H
30
32
33class WilsonPenzienNewmark final : public Newmark {
34 bool first_iteration = true;
35
36 const vec damping_ratio;
37
38 mat theta;
39 vec beta;
40
41public:
42 explicit WilsonPenzienNewmark(unsigned, vec&&, double = .25, double = .5);
43
44 int initialize() override;
45
46 [[nodiscard]] int process_constraint() override;
47
48 int solve(mat&, const mat&) override;
49 int solve(mat&, const sp_mat&) override;
50 int solve(mat&, mat&&) override;
51 int solve(mat&, sp_mat&&) override;
52
53 void commit_status() override;
54 void clear_status() override;
55 void reset_status() override;
56
57 void assemble_resistance() override;
58 void assemble_matrix() override;
59
60 void print() override;
61};
62
63#endif
64
A Newmark class defines a solver using Newmark algorithm.
Definition: Newmark.h:45
A WilsonPenzienNewmark class defines a solver using Newmark algorithm with Wilson-Penzien damping mod...
Definition: WilsonPenzienNewmark.h:33
void print() override
Definition: WilsonPenzienNewmark.cpp:181
void reset_status() override
Definition: WilsonPenzienNewmark.cpp:137
int solve(mat &, const mat &) override
Definition: WilsonPenzienNewmark.cpp:77
WilsonPenzienNewmark(unsigned, vec &&, double=.25, double=.5)
Definition: WilsonPenzienNewmark.cpp:23
int process_constraint() override
Definition: WilsonPenzienNewmark.cpp:39
void clear_status() override
Definition: WilsonPenzienNewmark.cpp:131
void commit_status() override
Definition: WilsonPenzienNewmark.cpp:125
void assemble_resistance() override
Definition: WilsonPenzienNewmark.cpp:143
int initialize() override
Definition: WilsonPenzienNewmark.cpp:27
void assemble_matrix() override
Definition: WilsonPenzienNewmark.cpp:162