From c7ff86c752d7d8521d739baf07113a8e9cf7cf3c Mon Sep 17 00:00:00 2001 From: ParkSuMin Date: Mon, 8 Sep 2025 22:22:57 +0300 Subject: [PATCH] =?UTF-8?q?=D0=94=D0=BE=D0=B1=D0=B0=D0=B2=D0=BB=D0=B5?= =?UTF-8?q?=D0=BD=D0=BE=20=D1=80=D0=B5=D1=88=D0=B5=D0=BD=D0=B8=D0=B5=20?= =?UTF-8?q?=D0=A1=D0=9B=D0=90=D0=A3?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit + вывод локальной матрицы жёсткости --- MemMAPR-MKE.cpp | 15 ++++++++++----- 1 file changed, 10 insertions(+), 5 deletions(-) diff --git a/MemMAPR-MKE.cpp b/MemMAPR-MKE.cpp index 571ce7f..7ae95fa 100644 --- a/MemMAPR-MKE.cpp +++ b/MemMAPR-MKE.cpp @@ -4,9 +4,9 @@ using namespace Eigen; int main() { std::ifstream file("matrix.txt"); - int N = 20; - int upper = 6, lower = 1; - double A = 1, B = 1, C = 1, L = upper - lower, du_dx = 1.; + int N = 3; + int upper = 1, lower = 0; + double A = 1, B = 0, C = 1, L = upper - lower; double dx = L / N; @@ -15,10 +15,11 @@ int main() { local(0, 1) = -A / dx + B / 2; local(1, 0) = -A / dx - B / 2; local(1, 1) = A / dx + B / 2; + std::cout << "Local matrix:\n" << local << std::endl; VectorXd local_load(2); - local_load(0) = C * dx / 2 - A * du_dx; - local_load(1) = C * dx / 2 + A * du_dx; + local_load(0) = C * dx / 2; + local_load(1) = C * dx / 2; // Ансамблированные матрица и вектор MatrixXd ansamb = MatrixXd::Zero(N + 1, N + 1); @@ -72,5 +73,9 @@ int main() { std::cout << "Modified matrix:\n" << ansamb << std::endl; std::cout << "Modified load vector:\n" << global_load << std::endl; + VectorXd solution = ansamb.fullPivLu().solve(global_load); + std::cout << "\nSolution:" << std::endl; + std::cout << solution << std::endl; + return 0; } \ No newline at end of file