MatricesC-/Matriz.cpp at master · pieroMC98/MatricesC- · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
#include "Matriz.h"

#include <iostream>

Matriz::Matriz(/* args */) {}
Matriz &Matriz::setI(int N) {
	this->i = N;
	return *this;
}

int Matriz::getI() { return i; }

Matriz &Matriz::setJ(int N) {
	this->j = N;
	return *this;
}

int Matriz::getJ() { return j; }

void Matriz::setMatriz(long double A[][SIZE]) {
	this->matriz = new long double *[this->i];
	for (int k = 0; k < this->i; k++) (this->matriz)[k] = new long double[this->j];

	for (int i = 0; i < this->i; i++)
		for (int j = 0; j < this->j; j++) this->matriz[i][j] = A[i][j];
}

void Matriz::setMatriz() {
	this->matriz = new long double *[this->i];
	for (int k = 0; k < this->i; k++) (this->matriz)[k] = new long double[this->j];

	for (int i = 0; i < this->i; i++)
		for (int j = 0; j < this->j; j++) {
			cout << "Introduzca [" << i + 1 << "][" << j + 1 << "]: " << endl;
			cin >> (this->matriz)[i][j];
		}
}

void Matriz::setMatrizFraccion(Fraccion **A) {
	this->matriz_gauss = new Fraccion *[this->i];
	for (int k = 0; k < this->i; k++) this->matriz_gauss[k] = new Fraccion[this->j];

	for (int i = 0; i < this->i; i++)
		for (int j = 0; j < this->j; j++) this->matriz_gauss[i][j] = A[i][j];
}

void Matriz::getMatriz() {
	for (int i = 0; i < this->i; i++) {
		for (int j = 0; j < this->j; j++) cout << setprecision(4) << "\t " << (this->matriz)[i][j];
		cout << endl << endl << endl;
	}
}

void Matriz::getMatrizGauss() {
	for (int i = 0; i < this->i; i++) {
		for (int j = 0; j < this->j; j++) cout << this->matriz_gauss[i][j].imprimir();
		cout << endl << endl << endl;
	}
}

double *resultado(int i, int j, double (*)[SIZE]) { /*
											   double* aux =
						       new double[7]; if( i < 7 ){

												   return
						       aux[i] + resultado(); }else return aux;
											*/
	return 0;
}

void Matriz::gauss(int i, int j) {
	if (j < this->i) {
		if (i < (this->i)) {
			long double m = this->matriz[i][j] / this->matriz[j][j];

			if (!m) gauss(j + 2, j + 1);

			for (int k = 0; k < this->j; k++) this->matriz[i][k] -= m * this->matriz[j][k];

			gauss(i + 1, j);
		} else
			gauss(j + 2, j + 1);
	}
}

void Matriz::gaussFraccion(int i, int j) {
	if (j < this->i) {
		if (i < (this->i)) {
			Fraccion *aux = new Fraccion;
			Fraccion diferencia, multiplicar, cociente;

			// cociente = *cociente.cociente(this->matriz_gauss[i][j], this->matriz_gauss[j][j]);
			aux = aux->cociente(this->matriz_gauss[i][j], this->matriz_gauss[j][j]);
			string c = this->matriz_gauss[i][j].imprimir(), c1 = this->matriz_gauss[j][j].imprimir(),
			       c3 = aux->imprimir(), M;

			if (!aux->getN()) gaussFraccion(j + 2, j + 1);
			// if (!cociente.getN()) gaussFraccion(j + 2, j + 1);

			for (int k = 0; k < this->j; k++) {
				this->matriz_gauss[i][k] = *(
					(new Fraccion())
						->diferencia(this->matriz_gauss[i][k],
							     *(new Fraccion())
								      ->multiplicar(*(aux), this->matriz_gauss[j][k])));
				M = this->matriz_gauss[i][k].imprimir();
			}

			// for (int k = 0; k < this->j; k++) {
			// 	this->matriz_gauss[i][k] = *(diferencia.diferencia(
			// 		this->matriz_gauss[i][k],
			// 		*(multiplicar).multiplicar(*(aux), this->matriz_gauss[j][k])));
			// }
			this->getMatrizGauss();
			Error(i);
			Error(j);
			jump;
			jump;
			gaussFraccion(i + 1, j);

		} else
			gaussFraccion(j + 2, j + 1);
	}
}

Matriz::~Matriz() { std::cout << "Objeto eliminado" << std::endl; }