From 19acd35a3d3c884efb5fbfea6fa9f7e24cb1ec44 Mon Sep 17 00:00:00 2001
From: Armando Luiz Nicolini Delgado <nicolui@inf.ufpr.br>
Date: Tue, 28 Jun 2022 10:45:05 -0300
Subject: [PATCH] =?UTF-8?q?/*=20Alterada=20a=20forma=20de=20gerar=20matriz?=
=?UTF-8?q?es=20diagonalmente=20dominantes=20ns=20fun=C3=A7=C3=A3o=20'iniS?=
=?UTF-8?q?isLin()'=20*/?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
utils/sislin.c | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/utils/sislin.c b/utils/sislin.c
index 875f6da..610b118 100644
--- a/utils/sislin.c
+++ b/utils/sislin.c
@@ -128,7 +128,10 @@ void iniSisLin (SistLinear_t *SL, tipoSistLinear_t tipo, real_t coef_max)
else if (tipo == diagDominante) {
// aumenta o expoente dos termos da diagonal principal
for (unsigned int i=0; i<n; ++i) {
- SL->A[i][i] *= (real_t)n;
+ real_t soma = 0.0;
+ for (unsigned int j=0; j < i; ++j) soma += SL->A[i][j];
+ for (unsigned int j=i+1; j < n; ++j) soma += SL->A[i][j];
+ SL->A[i][i] += soma;
}
}
}
--
GitLab