Skip to content
Snippets Groups Projects
Commit 9fd8bf95 authored by Marco Ratto's avatar Marco Ratto
Browse files

bug fix in formula for analytic Hessian.

parent f7aae77f
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 1 addition and 1 deletion
matlab/getH.m
+ 1
− 1
View file @ 9fd8bf95
...@@ -161,7 +161,7 @@ if nargout>5, ...@@ -161,7 +161,7 @@ if nargout>5,
gp22(j,i,:,:)=gp22(j,i,:,:)+reshape(tmpu+tmpu',[1 1 M_.param_nbr M_.param_nbr]); gp22(j,i,:,:)=gp22(j,i,:,:)+reshape(tmpu+tmpu',[1 1 M_.param_nbr M_.param_nbr]);
end end
tmp0=reshape(g2(j,:),[nc nc]); tmp0=reshape(g2(j,:),[nc nc]);
gp22(j,:,:,:)=gp22(j,:,:,:)+reshape(tmp0(:,1:ns)*d2yssdtheta(I,:,:),[1 nc M_.param_nbr M_.param_nbr]); gp22(j,:,:,:)=gp22(j,:,:,:)+reshape(tmp0(:,1:ns)*d2yssdtheta(I,:),[1 nc M_.param_nbr M_.param_nbr]);
end end
g22 = g22+gp22; g22 = g22+gp22;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment