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