Skip to content
Snippets Groups Projects
Commit 0c15ba09 authored by Sylvain CALINON's avatar Sylvain CALINON
Browse files

asym matrix corrected

parent aab6c5f5
Branches
No related tags found
No related merge requests found
......@@ -228,7 +228,7 @@ function [f, J] = f_curvature(x, param)
for t=1:param.nbData
a = x(:,t)' * Sa * x(:,t);
b = x(:,t)' * Sb * x(:,t) + 1E-8;
Jtmp = 2 * b^(-3/2) * Sa * x(:,t) - 3 * a * b^(-5/2) * Sb * x(:,t);
Jtmp = b^(-3/2) * (Sa + Sa') * x(:,t) - 3 * a * b^(-5/2) * Sb * x(:,t);
J = blkdiag(J, Jtmp');
end
end
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment