The subroutine DSYCONV stores the supdiagonal of the symmetric 1-by-1 or 2-by-2
authorlangou <langou@users.noreply.github.com>
Thu, 29 Oct 2015 04:58:50 +0000 (04:58 +0000)
committerlangou <langou@users.noreply.github.com>
Thu, 29 Oct 2015 04:58:50 +0000 (04:58 +0000)
commit794766767994706d17cb1d0fe20e92acf8e267b2
treee24564a6d2fe7ef3ff3d7d572a555fb73173d30e
parentb2f92b663a60f844694fbc38fbe046e01f7bb7ac
The subroutine DSYCONV stores the supdiagonal of the symmetric 1-by-1 or 2-by-2
block matrix in the array WORK. (Or the subdiagonal, same array.) Before
commit, only the nonzero elements of WORK were initialized to their correct
nonzero values. The zero entries were not initialized and not referenced and
assumed to be zero in subsequent subroutines. This is fine, but I think a
better behavior for the subroutine is to initialize these zero entries to zero.
So the commit is initializing these entries to zero.

Also a (double precision) ZERO in the D routine should have been an (integer)
0.

Also fix some formatting to have C, D, and S consistent.

Note: the indentations, comments of the Z routine is not consistent with C, D and
S. I did not have time to fix this.
SRC/csyconv.f
SRC/dsyconv.f
SRC/ssyconv.f
SRC/zsyconv.f