main.cpp
6973e704faffmaster

main.cpp
View Options

			#include <iostream>
			#include <cassert>
			#include <cmath>
			#include "Functions.hpp"
			#include "LU.hpp"


			int main(int argc, char *argv[]){
			int size=3;
			double** A=AllocateMemory(size,size);
			double* b=AllocateMemory(size);

			A[0][0]=1;A[0][1]=1;A[0][2]=1;
			A[1][0]=4;A[1][1]=3;A[1][2]=-1;
			A[2][0]=3;A[2][1]=5;A[2][2]=3;

			b[0]=1;b[1]=6;b[2]=4;
			std::cout<<"Matrix A:"<<std::endl;
			TerminalPrint(A,size,size);
			LU linsys(A,b,size);
			std::cout<<"Matrix U:"<<std::endl;
			TerminalPrint(linsys.get_U(),size,size);
			std::cout<<"Matrix L:"<<std::endl;
			TerminalPrint(linsys.get_L(),size,size);
			double** A_back=Multiply(linsys.get_L(),linsys.get_U(),size,size,size,size);
			std::cout<<"Matrix L * Matrix U:"<<std::endl;
			TerminalPrint(A_back,size,size);
			std::cout<<"Vector b:"<<std::endl;
			TerminalPrint(b,size);
			std::cout<<"Vector z:"<<std::endl;
			TerminalPrint(linsys.get_b_star(),size);
			std::cout<<"Solution of system Ax=b:"<<std::endl;
			TerminalPrint(linsys.Solve(),size);
			return 0;
			}