$title Test correctness of logit intrinsic (FNLOGIT,SEQ=880) $onText This test ensures the correct evaluation of logit Contributor: Lutz Westermann, October 2021 $offText * define some sets, etc. $include fnset_x.inc $gdxIn fnlogit $load aeps reps reps1 reps2 relToInput T data $gdxIn loop {T, data(T, 'f') = logit.value (data(T,'x')); data(T, 'fx') = logit.grad (data(T,'x')); data(T,'fxx') = logit.hess (data(T,'x')); data(T, 'rc') = mathlastrc; data(T, 'ec') = mathlastec; }; * the data set should contain some singular values abort$[execerror < 8] 'Should get execution errors evaluating logit.XXXX(x) for exceptional x'; execerror = 0; $include fntest_x.inc