Page 1 of 1

Minor bug in DGELSS

PostPosted: Fri Jul 31, 2009 8:26 am
by t.r.hopkins
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

Re: Minor bug in DGELSS

PostPosted: Fri Jul 31, 2009 3:57 pm
by Julien Langou
Thanks Tim. We will correct this. I let you know when this is done. (I fear this is a busy Summer.) Julien.