From f17dc240ef3d9690e8402fd30635c2f69992f1a5 Mon Sep 17 00:00:00 2001 From: Weslley S Pereira Date: Tue, 23 May 2023 18:25:27 -0600 Subject: [PATCH 1/2] Adds comment related to zero-sized matrix B in GELSY --- SRC/cgelsy.f | 1 + SRC/dgelsy.f | 1 + SRC/sgelsy.f | 1 + SRC/zgelsy.f | 1 + 4 files changed, 4 insertions(+) diff --git a/SRC/cgelsy.f b/SRC/cgelsy.f index 67140f1913..f32315cbbf 100644 --- a/SRC/cgelsy.f +++ b/SRC/cgelsy.f @@ -116,6 +116,7 @@ *> B is COMPLEX array, dimension (LDB,NRHS) *> On entry, the M-by-NRHS right hand side matrix B. *> On exit, the N-by-NRHS solution matrix X. +*> If M = 0 or N = 0, B is not referenced. *> \endverbatim *> *> \param[in] LDB diff --git a/SRC/dgelsy.f b/SRC/dgelsy.f index aebab92640..538e90c540 100644 --- a/SRC/dgelsy.f +++ b/SRC/dgelsy.f @@ -115,6 +115,7 @@ *> B is DOUBLE PRECISION array, dimension (LDB,NRHS) *> On entry, the M-by-NRHS right hand side matrix B. *> On exit, the N-by-NRHS solution matrix X. +*> If M = 0 or N = 0, B is not referenced. *> \endverbatim *> *> \param[in] LDB diff --git a/SRC/sgelsy.f b/SRC/sgelsy.f index 9c60f78a7a..92853176c3 100644 --- a/SRC/sgelsy.f +++ b/SRC/sgelsy.f @@ -115,6 +115,7 @@ *> B is REAL array, dimension (LDB,NRHS) *> On entry, the M-by-NRHS right hand side matrix B. *> On exit, the N-by-NRHS solution matrix X. +*> If M = 0 or N = 0, B is not referenced. *> \endverbatim *> *> \param[in] LDB diff --git a/SRC/zgelsy.f b/SRC/zgelsy.f index 65fa87ae98..bc76330efc 100644 --- a/SRC/zgelsy.f +++ b/SRC/zgelsy.f @@ -116,6 +116,7 @@ *> B is COMPLEX*16 array, dimension (LDB,NRHS) *> On entry, the M-by-NRHS right hand side matrix B. *> On exit, the N-by-NRHS solution matrix X. +*> If M = 0 or N = 0, B is not referenced. *> \endverbatim *> *> \param[in] LDB From b788b4ad205e5107f19f1bee8d108a774636b570 Mon Sep 17 00:00:00 2001 From: Weslley S Pereira Date: Tue, 23 May 2023 18:45:34 -0600 Subject: [PATCH 2/2] Fix comments to close issue #259 --- SRC/cgelsy.f | 1 + SRC/dgelsy.f | 1 + SRC/sgelsy.f | 1 + SRC/zgelsy.f | 1 + 4 files changed, 4 insertions(+) diff --git a/SRC/cgelsy.f b/SRC/cgelsy.f index f32315cbbf..64bb100231 100644 --- a/SRC/cgelsy.f +++ b/SRC/cgelsy.f @@ -149,6 +149,7 @@ *> The effective rank of A, i.e., the order of the submatrix *> R11. This is the same as the order of the submatrix T11 *> in the complete orthogonal factorization of A. +*> If NRHS = 0, RANK = 0 on output. *> \endverbatim *> *> \param[out] WORK diff --git a/SRC/dgelsy.f b/SRC/dgelsy.f index 538e90c540..e9fcd9682d 100644 --- a/SRC/dgelsy.f +++ b/SRC/dgelsy.f @@ -148,6 +148,7 @@ *> The effective rank of A, i.e., the order of the submatrix *> R11. This is the same as the order of the submatrix T11 *> in the complete orthogonal factorization of A. +*> If NRHS = 0, RANK = 0 on output. *> \endverbatim *> *> \param[out] WORK diff --git a/SRC/sgelsy.f b/SRC/sgelsy.f index 92853176c3..89dd39e80d 100644 --- a/SRC/sgelsy.f +++ b/SRC/sgelsy.f @@ -148,6 +148,7 @@ *> The effective rank of A, i.e., the order of the submatrix *> R11. This is the same as the order of the submatrix T11 *> in the complete orthogonal factorization of A. +*> If NRHS = 0, RANK = 0 on output. *> \endverbatim *> *> \param[out] WORK diff --git a/SRC/zgelsy.f b/SRC/zgelsy.f index bc76330efc..497becf8bd 100644 --- a/SRC/zgelsy.f +++ b/SRC/zgelsy.f @@ -149,6 +149,7 @@ *> The effective rank of A, i.e., the order of the submatrix *> R11. This is the same as the order of the submatrix T11 *> in the complete orthogonal factorization of A. +*> If NRHS = 0, RANK = 0 on output. *> \endverbatim *> *> \param[out] WORK