From 5af1effd8ffb8887fa39beeddd99cde84a60a6e8 Mon Sep 17 00:00:00 2001
From: Michel Juillard <michel.juillard@mjui.fr>
Date: Wed, 30 May 2012 10:25:10 +0200
Subject: [PATCH] removed debugging code from nonlinear solver when Jacobian is
near singular. The treatment of this case still needs to be further improved.
---
matlab/solve1.m | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/matlab/solve1.m b/matlab/solve1.m
index dfe00e53d7..5705e2c793 100644
--- a/matlab/solve1.m
+++ b/matlab/solve1.m
@@ -116,8 +116,7 @@ for its = 1:maxit
p = e*[-r(1:end-n,1:end-n)\fvec(1:end-n);zeros(n,1)];
end
elseif bad_cond_flag && cond(fjac) > 1/sqrt(eps)
- fjac
- fjac2=fjac'*fjac
+ fjac2=fjac'*fjac;
p=-(fjac2+1e6*sqrt(nn*eps)*max(sum(abs(fjac2)))*eye(nn))\(fjac'*fvec);
else
p = -fjac\fvec ;
--
GitLab