Matrix notation: fix bug in numerical Hessian of external functions

2020-07-06 18:42:05 +02:00 · 2020-07-06 18:42:05 +02:00 · 2acacf61f1
parent 1975e467b5
commit 2acacf61f1
1 changed files with 1 additions and 1 deletions
--- a/matlab/hess_element.m
+++ b/matlab/hess_element.m
@ -59,7 +59,7 @@ m11{arg2}(elem2) = m11{arg2}(elem2) - h;

 % From Abramowitz and Stegun. Handbook of Mathematical Functions (1965)
 % formulas 25.3.24 and 25.3.27 p. 884
-if arg1==arg2
+if arg1 == arg2 && elem1 == elem2
    d = (16*func(p10{:})...
         +16*func(m10{:})...
         -30*func(args{:})...