flip signs to be consistent with papaers
@@ -237,9 +237,9 @@ for i=0:max_iters-1
%}
%{%
epp = 0.;
- dF = F(:,1:end-1) - F(:,2:end);
+ dF = F(:,2:end) - F(:,1:end-1);
gg = (dF + epp * eye(size(dF))) \ F(:,end);
- dG = G(:,1:end-1) - G(:,2:end);
+ dG = G(:,2:end) - G(:,1:end-1);
uv = G(:,end) - dG * gg;
delta_accel = [delta_accel; norm([gg; zeros(anderson_lookback - length(gg),1)] - [a_last; zeros(anderson_lookback - length(a_last),1)])];
a_last = gg;
