GaussSeidelfunction.cpp
2a17905ce737master

GaussSeidelfunction.cpp
View Options

			//
			// Created by choitel on 12/3/18.
			//
			#include <iostream>
			#include <cmath>
			#include "GaussSeidel.h"

			using namespace std;

			double* GaussSeidel(double A[5][5],double b[5])
			{
			static double sol[5];

			double x[5];
			double xold[5];
			double e[5];
			double norm_e=1;
			for(int i=0;i<5;++i){
			x[i]=0;
			xold[i]=0;
			e[i]=1;
			}
			int count = 0;
			while(norm_e>1/100 && count<30) {

			for (int i = 0; i < 5; ++i) {
			x[i] = b[i] / A[i][i];
			for (int j = 0; j < 5; ++j) {
			if (i == j) {
			} else {
			x[i] = x[i] - (A[i][j] * x[j]) / A[i][i];
			}
			}
			if(xold[i]==0){
			e[i]=1;
			}else {
			e[i] = abs((x[i] - xold[i]) / xold[i]);
			}
			xold[i] = x[i];
			}
			norm_e=0;
			for(int i=0;i<5;++i){
			norm_e=norm_e + e[i]*e[i];
			}
			norm_e=sqrt(norm_e);

			count =count +1;
			}
			cout << "le nombre d'itérations du code est : " << endl;
			cout << count << endl;

			for(int i = 0;i<5;++i)
			{
			sol[i]=x[i];
			}

			return sol;
			}