SLICOT / SLICOT-Reference

SLICOT - A Fortran subroutines library for systems and control
BSD 3-Clause "New" or "Revised" License
49 stars 23 forks source link

Fix: IB01AD: Incorrect calculation of minimal working space #18

Open vlasovdl opened 1 month ago

vlasovdl commented 1 month ago

When subroutine IB01AD call with parameters ALG='F' and BATCH='O' minimal needed workspace should be calculated by expression MINWRK = 2*NR*( M + L + 1 ) + NR (line 584) but it calculating in line 582 MINWRK = NR*( M + L + 1 ). It is because logical FIRST is always true when BATCH='O' .