Diffusion PCSC_project (master) Browse

R7934/Richardson/Richardson.cpp
68db15c70096master

Richardson.cpp
View Options

			#include <iostream>
			#include "Richardson.h"

			using namespace std;

			int main() {

			double x[5];
			x[1]=1;
			x[2]=1;
			x[3]=1;
			x[4]=1;
			x[5]=1;
			double A[5][5];

			cout << "La matrice A vaut : " << endl;
			for(int i=0;i<5;++i) {
			for (int j = 0; j < 5; ++j) {
			if(i==j) {
			A[i][j]=5;
			cout << A[i][j] << " ";
			}else{
			A[i][j] = 0;
			cout << A[i][j] << " ";
			}
			}
			cout << endl;
			}
			cout << endl;

			double b[5];

			for(int j=0;j<5;++j)
			{
			b[j]=2;
			}
			double *solution;

			solution = Richardson(A,b);
			cout << endl;
			cout << "La valeur de x est: " << endl;
			cout << solution[0] << endl;
			cout << solution[1] << endl;
			cout << solution[2] << endl;
			cout << solution[3] << endl;
			cout << solution[4] << endl;

			return 0;


			return 0;
			}