From 272c508d88b4f56ec3942303294afef83784c905 Mon Sep 17 00:00:00 2001
From: Emilie Marchandise <emilie.marchandise@uclouvain.be>
Date: Wed, 7 Apr 2010 15:16:16 +0000
Subject: [PATCH] Riemann solver bifurcation
---
Numeric/Numeric.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/Numeric/Numeric.cpp b/Numeric/Numeric.cpp
index 5988610d3f..9d5628ffa6 100644
--- a/Numeric/Numeric.cpp
+++ b/Numeric/Numeric.cpp
@@ -642,15 +642,15 @@ bool newton_fd(void (*func)(fullVector<double> &, fullVector<double> &, void *),
if (isZero == false) break;
}
if (isZero) break;
- //printf("**** fffffff0=%g %g %g %g %g %g %g %g %g\n", f(0), f(1), f(2), f(3), f(4), f(5), f(6), f(7), f(8));
for (int j = 0; j < N; j++){
double h = EPS * fabs(x(j));
if(h == 0.) h = EPS;
x(j) += h;
func(x, feps, data);
- for (int i = 0; i < N; i++)
+ for (int i = 0; i < N; i++){
J(i, j) = (feps(i) - f(i)) / h;
+ }
x(j) -= h;
}
--
GitLab