Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

lk -> llk in *sytrd_sy2sb #812

Merged
merged 1 commit into from
May 13, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
30 changes: 15 additions & 15 deletions src/stdlib_linalg_lapack_c.fypp
Original file line number Diff line number Diff line change
Expand Up @@ -36507,7 +36507,7 @@ module stdlib_linalg_lapack_c

! Local Scalars
logical(lk) :: lquery, upper
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
tpos, wpos, s2pos, s1pos
! Intrinsic Functions
intrinsic :: min,max
Expand Down Expand Up @@ -36543,13 +36543,13 @@ module stdlib_linalg_lapack_c
if( n<=kd+1 ) then
if( upper ) then
do i = 1, n
lk = min( kd+1, i )
call stdlib_ccopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
llk = min( kd+1, i )
call stdlib_ccopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
end do
else
do i = 1, n
lk = min( kd+1, n-i+1 )
call stdlib_ccopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
llk = min( kd+1, n-i+1 )
call stdlib_ccopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
end do
endif
work( 1 ) = 1
Expand Down Expand Up @@ -36586,8 +36586,8 @@ module stdlib_linalg_lapack_c
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_ccopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min( kd, n-j ) + 1
call stdlib_ccopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
call stdlib_claset( 'LOWER', pk, pk, czero, cone,a( i, i+kd ), lda )
! form the matrix t
Expand All @@ -36609,8 +36609,8 @@ module stdlib_linalg_lapack_c
end do
! copy the upper band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_ccopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min(kd, n-j) + 1
call stdlib_ccopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
else
! reduce the lower triangle of a to lower band matrix
Expand All @@ -36622,8 +36622,8 @@ module stdlib_linalg_lapack_c
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_ccopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min( kd, n-j ) + 1
call stdlib_ccopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
call stdlib_claset( 'UPPER', pk, pk, czero, cone,a( i+kd, i ), lda )
! form the matrix t
Expand All @@ -36645,15 +36645,15 @@ module stdlib_linalg_lapack_c
! ==================================================================
! restore a for comparison and checking to be removed
! do 45 j = i, i+pk-1
! lk = min( kd, n-j ) + 1
! call stdlib_ccopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
! llk = min( kd, n-j ) + 1
! call stdlib_ccopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
45 continue
! ==================================================================
end do loop_40
! copy the lower band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_ccopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min(kd, n-j) + 1
call stdlib_ccopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
end if
work( 1 ) = lwmin
Expand Down
30 changes: 15 additions & 15 deletions src/stdlib_linalg_lapack_d.fypp
Original file line number Diff line number Diff line change
Expand Up @@ -57092,7 +57092,7 @@ module stdlib_linalg_lapack_d

! Local Scalars
logical(lk) :: lquery, upper
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
tpos, wpos, s2pos, s1pos
! Intrinsic Functions
intrinsic :: min,max
Expand Down Expand Up @@ -57128,13 +57128,13 @@ module stdlib_linalg_lapack_d
if( n<=kd+1 ) then
if( upper ) then
do i = 1, n
lk = min( kd+1, i )
call stdlib_dcopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
llk = min( kd+1, i )
call stdlib_dcopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
end do
else
do i = 1, n
lk = min( kd+1, n-i+1 )
call stdlib_dcopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
llk = min( kd+1, n-i+1 )
call stdlib_dcopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
end do
endif
work( 1 ) = 1
Expand Down Expand Up @@ -57171,8 +57171,8 @@ module stdlib_linalg_lapack_d
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_dcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min( kd, n-j ) + 1
call stdlib_dcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
call stdlib_dlaset( 'LOWER', pk, pk, zero, one,a( i, i+kd ), lda )
! form the matrix t
Expand All @@ -57194,8 +57194,8 @@ module stdlib_linalg_lapack_d
end do
! copy the upper band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_dcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min(kd, n-j) + 1
call stdlib_dcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
else
! reduce the lower triangle of a to lower band matrix
Expand All @@ -57207,8 +57207,8 @@ module stdlib_linalg_lapack_d
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_dcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min( kd, n-j ) + 1
call stdlib_dcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
call stdlib_dlaset( 'UPPER', pk, pk, zero, one,a( i+kd, i ), lda )
! form the matrix t
Expand All @@ -57230,15 +57230,15 @@ module stdlib_linalg_lapack_d
! ==================================================================
! restore a for comparison and checking to be removed
! do 45 j = i, i+pk-1
! lk = min( kd, n-j ) + 1
! call stdlib_dcopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
! llk = min( kd, n-j ) + 1
! call stdlib_dcopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
45 continue
! ==================================================================
end do loop_40
! copy the lower band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_dcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min(kd, n-j) + 1
call stdlib_dcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
end if
work( 1 ) = lwmin
Expand Down
30 changes: 15 additions & 15 deletions src/stdlib_linalg_lapack_q.fypp
Original file line number Diff line number Diff line change
Expand Up @@ -72917,7 +72917,7 @@ module stdlib_linalg_lapack_q

! Local Scalars
logical(lk) :: lquery, upper
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
tpos, wpos, s2pos, s1pos
! Intrinsic Functions
intrinsic :: min,max
Expand Down Expand Up @@ -72953,13 +72953,13 @@ module stdlib_linalg_lapack_q
if( n<=kd+1 ) then
if( upper ) then
do i = 1, n
lk = min( kd+1, i )
call stdlib_qcopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
llk = min( kd+1, i )
call stdlib_qcopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
end do
else
do i = 1, n
lk = min( kd+1, n-i+1 )
call stdlib_qcopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
llk = min( kd+1, n-i+1 )
call stdlib_qcopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
end do
endif
work( 1 ) = 1
Expand Down Expand Up @@ -72996,8 +72996,8 @@ module stdlib_linalg_lapack_q
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_qcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min( kd, n-j ) + 1
call stdlib_qcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
call stdlib_qlaset( 'LOWER', pk, pk, zero, one,a( i, i+kd ), lda )
! form the matrix t
Expand All @@ -73019,8 +73019,8 @@ module stdlib_linalg_lapack_q
end do
! copy the upper band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_qcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min(kd, n-j) + 1
call stdlib_qcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
else
! reduce the lower triangle of a to lower band matrix
Expand All @@ -73032,8 +73032,8 @@ module stdlib_linalg_lapack_q
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_qcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min( kd, n-j ) + 1
call stdlib_qcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
call stdlib_qlaset( 'UPPER', pk, pk, zero, one,a( i+kd, i ), lda )
! form the matrix t
Expand All @@ -73055,15 +73055,15 @@ module stdlib_linalg_lapack_q
! ==================================================================
! restore a for comparison and checking to be removed
! do 45 j = i, i+pk-1
! lk = min( kd, n-j ) + 1
! call stdlib_qcopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
! llk = min( kd, n-j ) + 1
! call stdlib_qcopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
45 continue
! ==================================================================
end do loop_40
! copy the lower band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_qcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min(kd, n-j) + 1
call stdlib_qcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
end if
work( 1 ) = lwmin
Expand Down
30 changes: 15 additions & 15 deletions src/stdlib_linalg_lapack_s.fypp
Original file line number Diff line number Diff line change
Expand Up @@ -55567,7 +55567,7 @@ module stdlib_linalg_lapack_s

! Local Scalars
logical(lk) :: lquery, upper
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
tpos, wpos, s2pos, s1pos
! Intrinsic Functions
intrinsic :: min,max
Expand Down Expand Up @@ -55603,13 +55603,13 @@ module stdlib_linalg_lapack_s
if( n<=kd+1 ) then
if( upper ) then
do i = 1, n
lk = min( kd+1, i )
call stdlib_scopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
llk = min( kd+1, i )
call stdlib_scopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
end do
else
do i = 1, n
lk = min( kd+1, n-i+1 )
call stdlib_scopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
llk = min( kd+1, n-i+1 )
call stdlib_scopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
end do
endif
work( 1 ) = 1
Expand Down Expand Up @@ -55646,8 +55646,8 @@ module stdlib_linalg_lapack_s
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_scopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min( kd, n-j ) + 1
call stdlib_scopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
call stdlib_slaset( 'LOWER', pk, pk, zero, one,a( i, i+kd ), lda )
! form the matrix t
Expand All @@ -55669,8 +55669,8 @@ module stdlib_linalg_lapack_s
end do
! copy the upper band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_scopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min(kd, n-j) + 1
call stdlib_scopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
else
! reduce the lower triangle of a to lower band matrix
Expand All @@ -55682,8 +55682,8 @@ module stdlib_linalg_lapack_s
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_scopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min( kd, n-j ) + 1
call stdlib_scopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
call stdlib_slaset( 'UPPER', pk, pk, zero, one,a( i+kd, i ), lda )
! form the matrix t
Expand All @@ -55705,15 +55705,15 @@ module stdlib_linalg_lapack_s
! ==================================================================
! restore a for comparison and checking to be removed
! do 45 j = i, i+pk-1
! lk = min( kd, n-j ) + 1
! call stdlib_scopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
! llk = min( kd, n-j ) + 1
! call stdlib_scopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
45 continue
! ==================================================================
end do loop_40
! copy the lower band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_scopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min(kd, n-j) + 1
call stdlib_scopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
end if
work( 1 ) = lwmin
Expand Down
30 changes: 15 additions & 15 deletions src/stdlib_linalg_lapack_w.fypp
Original file line number Diff line number Diff line change
Expand Up @@ -28868,7 +28868,7 @@ module stdlib_linalg_lapack_w

! Local Scalars
logical(lk) :: lquery, upper
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
tpos, wpos, s2pos, s1pos
! Intrinsic Functions
intrinsic :: min,max
Expand Down Expand Up @@ -28904,13 +28904,13 @@ module stdlib_linalg_lapack_w
if( n<=kd+1 ) then
if( upper ) then
do i = 1, n
lk = min( kd+1, i )
call stdlib_wcopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
llk = min( kd+1, i )
call stdlib_wcopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
end do
else
do i = 1, n
lk = min( kd+1, n-i+1 )
call stdlib_wcopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
llk = min( kd+1, n-i+1 )
call stdlib_wcopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
end do
endif
work( 1 ) = 1
Expand Down Expand Up @@ -28947,8 +28947,8 @@ module stdlib_linalg_lapack_w
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_wcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min( kd, n-j ) + 1
call stdlib_wcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
call stdlib_wlaset( 'LOWER', pk, pk, czero, cone,a( i, i+kd ), lda )
! form the matrix t
Expand All @@ -28970,8 +28970,8 @@ module stdlib_linalg_lapack_w
end do
! copy the upper band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_wcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
llk = min(kd, n-j) + 1
call stdlib_wcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
end do
else
! reduce the lower triangle of a to lower band matrix
Expand All @@ -28983,8 +28983,8 @@ module stdlib_linalg_lapack_w
iinfo )
! copy the upper portion of a into ab
do j = i, i+pk-1
lk = min( kd, n-j ) + 1
call stdlib_wcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min( kd, n-j ) + 1
call stdlib_wcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
call stdlib_wlaset( 'UPPER', pk, pk, czero, cone,a( i+kd, i ), lda )
! form the matrix t
Expand All @@ -29006,15 +29006,15 @@ module stdlib_linalg_lapack_w
! ==================================================================
! restore a for comparison and checking to be removed
! do 45 j = i, i+pk-1
! lk = min( kd, n-j ) + 1
! call stdlib_wcopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
! llk = min( kd, n-j ) + 1
! call stdlib_wcopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
45 continue
! ==================================================================
end do loop_40
! copy the lower band to ab which is the band storage matrix
do j = n-kd+1, n
lk = min(kd, n-j) + 1
call stdlib_wcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
llk = min(kd, n-j) + 1
call stdlib_wcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
end do
end if
work( 1 ) = lwmin
Expand Down
Loading
Loading