Hi, I've noticed the following behavior of the invpd function: invpd(-1) yields -1; instead and as expected, invpd(-I(2)) throws an error about non-pd-ness of the input. Shouldn't this apply to scalar inputs as well?

On Thu, Dec 18, 2025 at 11:58 AM Riccardo (Jack) Lucchetti <p002264(a)staff.univpm.it&gt; wrote: > On 18/12/2025 17:15, Sven Schreiber wrote: >> Hi, >> >> I've noticed the following behavior of the invpd function: invpd(-1) >> yields -1; instead and as expected, invpd(-I(2)) throws an error about >> non-pd-ness of the input. Shouldn't this apply to scalar inputs as well? > Yes, IMO. I just pushed a fix to git. Allin, could you please review? Looks good, thanks.