kitten/MemMAPR-MKE
1
0
Fork 0
forked from ParkSuMin/MemMAPR-MKE
Code Pull Requests Activity

Ready linear MKE solver

Browse Source
This commit is contained in:
ParkSuMin
2025-09-10 14:34:30 +03:00
parent b66a923d26
commit e5ed6001f8
4 changed files with 26 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
Solver.cpp
View File

@@ -9,7 +9,7 @@ Solver::Solver(double _A, double _B, double _C, int _N, int _l, int _u) {
dx = L / N;
}
void Solver::Execute() {
void Solver::Execute_Linear() {
MatrixXd local = MatrixXd::Zero(2, 2);
local(0, 0) = A / dx - B / 2;
local(0, 1) = -A / dx + B / 2;
@@ -76,4 +76,24 @@ void Solver::Execute() {
VectorXd solution = ansamb.fullPivLu().solve(global_load);
std::cout << "\nSolution:" << std::endl;
std::cout << solution << std::endl;
std::ofstream file("matrix_2.txt");
for (int i = 0; i < N + 1; i++) {
file << solution(i) << ' ';
}
file << std::endl;
}
void Solver::Execute_Cubic() {
MatrixXd local = MatrixXd::Zero(2, 2);
local(0, 0) = A * 37/(10 * dx) + B / 2;
local(0, 1) = -A / dx + B / 2;
local(0, 1) = A;
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;
local_load(1) = C * dx / 2;
}