There appears to be a minor problem with the routine dgelss.f. If I run with array bound
checking switched on I get the following error message from a call to dgelss
Subscript 1 of A (value 2) is out of range (1:1)
Program terminated by fatal error
In DLASET, line 99 of dlaset.f
Called by DGELSS, line 284 of dgelss.f
This appears to be caused by an error in the dimensioning information being passed to DLASET
from DGELSS. The following edit appears to fix things -- S has MINMN rows and a single
column
284c284
< CALL DLASET( 'F', MINMN, 1, ZERO, ZERO, S, 1 )
---
> CALL DLASET( 'F', MINMN, 1, ZERO, ZERO, S, MINMN )
Best wishes,
Tim

