Ready linear MKE solver

.gitignore

@@ -636,3 +636,4 @@ FodyWeavers.xsd
 *.msix
 *.msm
 *.msp
+*.txt

MemMAPR-MKE.cpp

@@ -2,9 +2,8 @@
 #include "Solver.h"
 
 int main() {
-    std::ifstream file("matrix.txt");
-    Solver slv(1., 0., 1., 3, 0, 1);
-    slv.Execute();
+    Solver slv(1., 0., 1., 20, 0, 1);
+    slv.Execute_Linear();
 
     return 0;
 }

Solver.cpp

@@ -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;
 }

Solver.h

@@ -5,5 +5,6 @@ private:
 	int N, upper, lower;
 public:
 	Solver(double _A, double _B, double _C, int _N, int _l, int _u);
-	void Execute();
+	void Execute_Linear();
+	void Execute_Cubic();
 };

matrix.txt

@@ -1,3 +0,0 @@
-1 2 4
-4 5 6
-7 8 9