From 5ab3cf7986c255aed3dede55ac20d1489ce44eb9 Mon Sep 17 00:00:00 2001
From: ParkSuMin <kamalovap@paksurion.ru>
Date: Mon, 8 Sep 2025 20:10:23 +0300
Subject: [PATCH] Test solution

---
 Header.h                    |  3 ++
 MemMAPR-MKE.cpp             | 79 ++++++++++++++++++++++++++++++-------
 MemMAPR-MKE.vcxproj         |  3 ++
 MemMAPR-MKE.vcxproj.filters |  5 +++
 4 files changed, 76 insertions(+), 14 deletions(-)
 create mode 100644 Header.h

diff --git a/Header.h b/Header.h
new file mode 100644
index 0000000..7a1d8de
--- /dev/null
+++ b/Header.h
@@ -0,0 +1,3 @@
+#include <iostream>
+#include <fstream>
+#include <Eigen/Dense>
\ No newline at end of file
diff --git a/MemMAPR-MKE.cpp b/MemMAPR-MKE.cpp
index e2dbb99..571ce7f 100644
--- a/MemMAPR-MKE.cpp
+++ b/MemMAPR-MKE.cpp
@@ -1,25 +1,76 @@
-#include <iostream>
-#include <fstream>
-#include <Eigen/Dense>
+#include "Header.h"
 
 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.;
 
-    // Если размеры известны заранее (3x3 в вашем случае)
-    Matrix <double, 3, 3> matrix;
+    double dx = L / N;
 
-    for (int i = 0; i < 3; ++i) {
-        for (int j = 0; j < 3; ++j) {
-            if (!(file >> matrix(i, j))) {
-                std::cerr << "Error in matrix initialization";
-                return -1;
-            }
-            
-        }
+    MatrixXd local = MatrixXd::Zero(2, 2);
+    local(0, 0) = A / dx - B / 2;
+    local(0, 1) = -A / dx + B / 2;
+    local(1, 0) = -A / dx - B / 2;
+    local(1, 1) = A / dx + B / 2;
+
+    VectorXd local_load(2);
+    local_load(0) = C * dx / 2 - A * du_dx; 
+    local_load(1) = C * dx / 2 + A * du_dx;
+
+    // Ансамблированные матрица и вектор
+    MatrixXd ansamb = MatrixXd::Zero(N + 1, N + 1);
+    VectorXd global_load = VectorXd::Zero(N + 1);
+
+    // Ансамблирование для каждого конечного элемента
+    for (int elem = 0; elem < N; ++elem) {
+        int node_i = elem;
+        int node_j = elem + 1;
+
+        // Матрица жесткости
+        ansamb(node_i, node_i) += local(0, 0);
+        ansamb(node_i, node_j) += local(0, 1);
+        ansamb(node_j, node_i) += local(1, 0);
+        ansamb(node_j, node_j) += local(1, 1);
+
+        // Вектор нагрузки
+        global_load(node_i) += local_load(0);
+        global_load(node_j) += local_load(1);
     }
 
-    std::cout << "Matrix:\n" << matrix << std::endl;
+    std::cout << "Before:" << std::endl;
+    std::cout << "Ansamb matrix:\n" << ansamb << std::endl;
+    std::cout << "Ansamb load vector:\n" << global_load << std::endl;
+
+    // ГУ 1-го рода
+    double u_left = 5.0;   // u(1) = 5
+    double u_right = 15.0; // u(6) = 15
+
+    // 1. Учёт условия u(1)=5
+    for (int i = 1; i < N + 1; ++i) {
+        global_load(i) -= ansamb(i, 0) * u_left;
+    }
+    // Обнуляем первую строку и столбец матрицы, устанавливаем диагональ = 1
+    ansamb.row(0).setZero();
+    ansamb.col(0).setZero();
+    ansamb(0, 0) = 1.0;
+    global_load(0) = u_left;
+
+    // 2. Учёт условия u(6)=15
+    for (int i = 0; i < N; ++i) {
+        global_load(i) -= ansamb(i, N) * u_right;
+    }
+    // Обнуляем последнюю строку и столбец матрицы, устанавливаем диагональ = 1
+    ansamb.row(N).setZero();
+    ansamb.col(N).setZero();
+    ansamb(N, N) = 1.0;
+    global_load(N) = u_right;
+
+    std::cout << "\nAfter:" << std::endl;
+    std::cout << "Modified matrix:\n" << ansamb << std::endl;
+    std::cout << "Modified load vector:\n" << global_load << std::endl;
+
     return 0;
 }
\ No newline at end of file
diff --git a/MemMAPR-MKE.vcxproj b/MemMAPR-MKE.vcxproj
index 8e46e19..d5110a0 100644
--- a/MemMAPR-MKE.vcxproj
+++ b/MemMAPR-MKE.vcxproj
@@ -128,6 +128,9 @@
   <ItemGroup>
     <None Include="packages.config" />
   </ItemGroup>
+  <ItemGroup>
+    <ClInclude Include="Header.h" />
+  </ItemGroup>
   <Import Project="$(VCTargetsPath)\Microsoft.Cpp.targets" />
   <ImportGroup Label="ExtensionTargets">
     <Import Project="packages\Eigen.3.3.3\build\native\Eigen.targets" Condition="Exists('packages\Eigen.3.3.3\build\native\Eigen.targets')" />
diff --git a/MemMAPR-MKE.vcxproj.filters b/MemMAPR-MKE.vcxproj.filters
index 069fe56..09afe7f 100644
--- a/MemMAPR-MKE.vcxproj.filters
+++ b/MemMAPR-MKE.vcxproj.filters
@@ -22,4 +22,9 @@
   <ItemGroup>
     <None Include="packages.config" />
   </ItemGroup>
+  <ItemGroup>
+    <ClInclude Include="Header.h">
+      <Filter>Исходные файлы</Filter>
+    </ClInclude>
+  </ItemGroup>
 </Project>
\ No newline at end of file