From 9abb2e9ff8d8190047b70fd0251fc8eedc54ed88 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?St=C3=A9phane=20Adjemian=20=28Charybdis=29?=
<stephane.adjemian@univ-lemans.fr>
Date: Mon, 23 Jan 2012 13:57:30 +0100
Subject: [PATCH] Adapted homotopy routine to SEP.
---
matlab/ep/homotopic_steps.m | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/matlab/ep/homotopic_steps.m b/matlab/ep/homotopic_steps.m
index 814f01ffa5..21d8d31adc 100644
--- a/matlab/ep/homotopic_steps.m
+++ b/matlab/ep/homotopic_steps.m
@@ -30,12 +30,11 @@ disp(' ')
if homotopy_1
while weight<1
iter = iter+1;
- oo_.exo_simul(2,:) = weight*exxo_simul(2,:);
+ oo_.exo_simul = weight*exxo_simul;
t0 = tic;
[flag,tmp] = bytecode('dynamic');
TeaTime = toc(t0);
ctime = ctime+TeaTime;
- %old_weight = weight;
info.convergence = ~flag;
if verbose
if ~info.convergence
@@ -80,7 +79,6 @@ if homotopy_1
break
end
weight = weight+step_length;
- %step_length = initial_step_length;
end
if iter>max_iter
info = NaN;
@@ -88,7 +86,7 @@ if homotopy_1
end
end
if weight<1 && homotopy_1
- oo_.exo_simul(2,:) = exxo_simul(2,:);
+ oo_.exo_simul = exxo_simul;
t0 = tic;
[flag,tmp] = bytecode('dynamic');
TeaTime = toc(t0);
--
GitLab