Skip to content

Commit 1dc58c9

Browse files
jvdp1jvdp1
authored
Merge pull request #812 from perazz/fix_lk
`lk` -> `llk` in `*sytrd_sy2sb`
2 parents f63585d + 4d76130 commit 1dc58c9

File tree

6 files changed

+90
-90
lines changed

6 files changed

+90
-90
lines changed

‎src/stdlib_linalg_lapack_c.fypp‎

Lines changed: 15 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -36507,7 +36507,7 @@ module stdlib_linalg_lapack_c
3650736507

3650836508
! Local Scalars
3650936509
logical(lk) :: lquery, upper
36510-
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
36510+
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
3651136511
tpos, wpos, s2pos, s1pos
3651236512
! Intrinsic Functions
3651336513
intrinsic :: min,max
@@ -36543,13 +36543,13 @@ module stdlib_linalg_lapack_c
3654336543
if( n<=kd+1 ) then
3654436544
if( upper ) then
3654536545
do i = 1, n
36546-
lk = min( kd+1, i )
36547-
call stdlib_ccopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
36546+
llk = min( kd+1, i )
36547+
call stdlib_ccopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
3654836548
end do
3654936549
else
3655036550
do i = 1, n
36551-
lk = min( kd+1, n-i+1 )
36552-
call stdlib_ccopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
36551+
llk = min( kd+1, n-i+1 )
36552+
call stdlib_ccopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
3655336553
end do
3655436554
endif
3655536555
work( 1 ) = 1
@@ -36586,8 +36586,8 @@ module stdlib_linalg_lapack_c
3658636586
iinfo )
3658736587
! copy the upper portion of a into ab
3658836588
do j = i, i+pk-1
36589-
lk = min( kd, n-j ) + 1
36590-
call stdlib_ccopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
36589+
llk = min( kd, n-j ) + 1
36590+
call stdlib_ccopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
3659136591
end do
3659236592
call stdlib_claset( 'LOWER', pk, pk, czero, cone,a( i, i+kd ), lda )
3659336593
! form the matrix t
@@ -36609,8 +36609,8 @@ module stdlib_linalg_lapack_c
3660936609
end do
3661036610
! copy the upper band to ab which is the band storage matrix
3661136611
do j = n-kd+1, n
36612-
lk = min(kd, n-j) + 1
36613-
call stdlib_ccopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
36612+
llk = min(kd, n-j) + 1
36613+
call stdlib_ccopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
3661436614
end do
3661536615
else
3661636616
! reduce the lower triangle of a to lower band matrix
@@ -36622,8 +36622,8 @@ module stdlib_linalg_lapack_c
3662236622
iinfo )
3662336623
! copy the upper portion of a into ab
3662436624
do j = i, i+pk-1
36625-
lk = min( kd, n-j ) + 1
36626-
call stdlib_ccopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
36625+
llk = min( kd, n-j ) + 1
36626+
call stdlib_ccopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
3662736627
end do
3662836628
call stdlib_claset( 'UPPER', pk, pk, czero, cone,a( i+kd, i ), lda )
3662936629
! form the matrix t
@@ -36645,15 +36645,15 @@ module stdlib_linalg_lapack_c
3664536645
! ==================================================================
3664636646
! restore a for comparison and checking to be removed
3664736647
! do 45 j = i, i+pk-1
36648-
! lk = min( kd, n-j ) + 1
36649-
! call stdlib_ccopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
36648+
! llk = min( kd, n-j ) + 1
36649+
! call stdlib_ccopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
3665036650
45 continue
3665136651
! ==================================================================
3665236652
end do loop_40
3665336653
! copy the lower band to ab which is the band storage matrix
3665436654
do j = n-kd+1, n
36655-
lk = min(kd, n-j) + 1
36656-
call stdlib_ccopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
36655+
llk = min(kd, n-j) + 1
36656+
call stdlib_ccopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
3665736657
end do
3665836658
end if
3665936659
work( 1 ) = lwmin

‎src/stdlib_linalg_lapack_d.fypp‎

Lines changed: 15 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -57092,7 +57092,7 @@ module stdlib_linalg_lapack_d
5709257092

5709357093
! Local Scalars
5709457094
logical(lk) :: lquery, upper
57095-
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
57095+
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
5709657096
tpos, wpos, s2pos, s1pos
5709757097
! Intrinsic Functions
5709857098
intrinsic :: min,max
@@ -57128,13 +57128,13 @@ module stdlib_linalg_lapack_d
5712857128
if( n<=kd+1 ) then
5712957129
if( upper ) then
5713057130
do i = 1, n
57131-
lk = min( kd+1, i )
57132-
call stdlib_dcopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
57131+
llk = min( kd+1, i )
57132+
call stdlib_dcopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
5713357133
end do
5713457134
else
5713557135
do i = 1, n
57136-
lk = min( kd+1, n-i+1 )
57137-
call stdlib_dcopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
57136+
llk = min( kd+1, n-i+1 )
57137+
call stdlib_dcopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
5713857138
end do
5713957139
endif
5714057140
work( 1 ) = 1
@@ -57171,8 +57171,8 @@ module stdlib_linalg_lapack_d
5717157171
iinfo )
5717257172
! copy the upper portion of a into ab
5717357173
do j = i, i+pk-1
57174-
lk = min( kd, n-j ) + 1
57175-
call stdlib_dcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
57174+
llk = min( kd, n-j ) + 1
57175+
call stdlib_dcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
5717657176
end do
5717757177
call stdlib_dlaset( 'LOWER', pk, pk, zero, one,a( i, i+kd ), lda )
5717857178
! form the matrix t
@@ -57194,8 +57194,8 @@ module stdlib_linalg_lapack_d
5719457194
end do
5719557195
! copy the upper band to ab which is the band storage matrix
5719657196
do j = n-kd+1, n
57197-
lk = min(kd, n-j) + 1
57198-
call stdlib_dcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
57197+
llk = min(kd, n-j) + 1
57198+
call stdlib_dcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
5719957199
end do
5720057200
else
5720157201
! reduce the lower triangle of a to lower band matrix
@@ -57207,8 +57207,8 @@ module stdlib_linalg_lapack_d
5720757207
iinfo )
5720857208
! copy the upper portion of a into ab
5720957209
do j = i, i+pk-1
57210-
lk = min( kd, n-j ) + 1
57211-
call stdlib_dcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
57210+
llk = min( kd, n-j ) + 1
57211+
call stdlib_dcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
5721257212
end do
5721357213
call stdlib_dlaset( 'UPPER', pk, pk, zero, one,a( i+kd, i ), lda )
5721457214
! form the matrix t
@@ -57230,15 +57230,15 @@ module stdlib_linalg_lapack_d
5723057230
! ==================================================================
5723157231
! restore a for comparison and checking to be removed
5723257232
! do 45 j = i, i+pk-1
57233-
! lk = min( kd, n-j ) + 1
57234-
! call stdlib_dcopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
57233+
! llk = min( kd, n-j ) + 1
57234+
! call stdlib_dcopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
5723557235
45 continue
5723657236
! ==================================================================
5723757237
end do loop_40
5723857238
! copy the lower band to ab which is the band storage matrix
5723957239
do j = n-kd+1, n
57240-
lk = min(kd, n-j) + 1
57241-
call stdlib_dcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
57240+
llk = min(kd, n-j) + 1
57241+
call stdlib_dcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
5724257242
end do
5724357243
end if
5724457244
work( 1 ) = lwmin

‎src/stdlib_linalg_lapack_q.fypp‎

Lines changed: 15 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -72917,7 +72917,7 @@ module stdlib_linalg_lapack_q
7291772917

7291872918
! Local Scalars
7291972919
logical(lk) :: lquery, upper
72920-
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
72920+
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
7292172921
tpos, wpos, s2pos, s1pos
7292272922
! Intrinsic Functions
7292372923
intrinsic :: min,max
@@ -72953,13 +72953,13 @@ module stdlib_linalg_lapack_q
7295372953
if( n<=kd+1 ) then
7295472954
if( upper ) then
7295572955
do i = 1, n
72956-
lk = min( kd+1, i )
72957-
call stdlib_qcopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
72956+
llk = min( kd+1, i )
72957+
call stdlib_qcopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
7295872958
end do
7295972959
else
7296072960
do i = 1, n
72961-
lk = min( kd+1, n-i+1 )
72962-
call stdlib_qcopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
72961+
llk = min( kd+1, n-i+1 )
72962+
call stdlib_qcopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
7296372963
end do
7296472964
endif
7296572965
work( 1 ) = 1
@@ -72996,8 +72996,8 @@ module stdlib_linalg_lapack_q
7299672996
iinfo )
7299772997
! copy the upper portion of a into ab
7299872998
do j = i, i+pk-1
72999-
lk = min( kd, n-j ) + 1
73000-
call stdlib_qcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
72999+
llk = min( kd, n-j ) + 1
73000+
call stdlib_qcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
7300173001
end do
7300273002
call stdlib_qlaset( 'LOWER', pk, pk, zero, one,a( i, i+kd ), lda )
7300373003
! form the matrix t
@@ -73019,8 +73019,8 @@ module stdlib_linalg_lapack_q
7301973019
end do
7302073020
! copy the upper band to ab which is the band storage matrix
7302173021
do j = n-kd+1, n
73022-
lk = min(kd, n-j) + 1
73023-
call stdlib_qcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
73022+
llk = min(kd, n-j) + 1
73023+
call stdlib_qcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
7302473024
end do
7302573025
else
7302673026
! reduce the lower triangle of a to lower band matrix
@@ -73032,8 +73032,8 @@ module stdlib_linalg_lapack_q
7303273032
iinfo )
7303373033
! copy the upper portion of a into ab
7303473034
do j = i, i+pk-1
73035-
lk = min( kd, n-j ) + 1
73036-
call stdlib_qcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
73035+
llk = min( kd, n-j ) + 1
73036+
call stdlib_qcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
7303773037
end do
7303873038
call stdlib_qlaset( 'UPPER', pk, pk, zero, one,a( i+kd, i ), lda )
7303973039
! form the matrix t
@@ -73055,15 +73055,15 @@ module stdlib_linalg_lapack_q
7305573055
! ==================================================================
7305673056
! restore a for comparison and checking to be removed
7305773057
! do 45 j = i, i+pk-1
73058-
! lk = min( kd, n-j ) + 1
73059-
! call stdlib_qcopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
73058+
! llk = min( kd, n-j ) + 1
73059+
! call stdlib_qcopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
7306073060
45 continue
7306173061
! ==================================================================
7306273062
end do loop_40
7306373063
! copy the lower band to ab which is the band storage matrix
7306473064
do j = n-kd+1, n
73065-
lk = min(kd, n-j) + 1
73066-
call stdlib_qcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
73065+
llk = min(kd, n-j) + 1
73066+
call stdlib_qcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
7306773067
end do
7306873068
end if
7306973069
work( 1 ) = lwmin

‎src/stdlib_linalg_lapack_s.fypp‎

Lines changed: 15 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -55567,7 +55567,7 @@ module stdlib_linalg_lapack_s
5556755567

5556855568
! Local Scalars
5556955569
logical(lk) :: lquery, upper
55570-
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
55570+
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
5557155571
tpos, wpos, s2pos, s1pos
5557255572
! Intrinsic Functions
5557355573
intrinsic :: min,max
@@ -55603,13 +55603,13 @@ module stdlib_linalg_lapack_s
5560355603
if( n<=kd+1 ) then
5560455604
if( upper ) then
5560555605
do i = 1, n
55606-
lk = min( kd+1, i )
55607-
call stdlib_scopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
55606+
llk = min( kd+1, i )
55607+
call stdlib_scopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
5560855608
end do
5560955609
else
5561055610
do i = 1, n
55611-
lk = min( kd+1, n-i+1 )
55612-
call stdlib_scopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
55611+
llk = min( kd+1, n-i+1 )
55612+
call stdlib_scopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
5561355613
end do
5561455614
endif
5561555615
work( 1 ) = 1
@@ -55646,8 +55646,8 @@ module stdlib_linalg_lapack_s
5564655646
iinfo )
5564755647
! copy the upper portion of a into ab
5564855648
do j = i, i+pk-1
55649-
lk = min( kd, n-j ) + 1
55650-
call stdlib_scopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
55649+
llk = min( kd, n-j ) + 1
55650+
call stdlib_scopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
5565155651
end do
5565255652
call stdlib_slaset( 'LOWER', pk, pk, zero, one,a( i, i+kd ), lda )
5565355653
! form the matrix t
@@ -55669,8 +55669,8 @@ module stdlib_linalg_lapack_s
5566955669
end do
5567055670
! copy the upper band to ab which is the band storage matrix
5567155671
do j = n-kd+1, n
55672-
lk = min(kd, n-j) + 1
55673-
call stdlib_scopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
55672+
llk = min(kd, n-j) + 1
55673+
call stdlib_scopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
5567455674
end do
5567555675
else
5567655676
! reduce the lower triangle of a to lower band matrix
@@ -55682,8 +55682,8 @@ module stdlib_linalg_lapack_s
5568255682
iinfo )
5568355683
! copy the upper portion of a into ab
5568455684
do j = i, i+pk-1
55685-
lk = min( kd, n-j ) + 1
55686-
call stdlib_scopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
55685+
llk = min( kd, n-j ) + 1
55686+
call stdlib_scopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
5568755687
end do
5568855688
call stdlib_slaset( 'UPPER', pk, pk, zero, one,a( i+kd, i ), lda )
5568955689
! form the matrix t
@@ -55705,15 +55705,15 @@ module stdlib_linalg_lapack_s
5570555705
! ==================================================================
5570655706
! restore a for comparison and checking to be removed
5570755707
! do 45 j = i, i+pk-1
55708-
! lk = min( kd, n-j ) + 1
55709-
! call stdlib_scopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
55708+
! llk = min( kd, n-j ) + 1
55709+
! call stdlib_scopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
5571055710
45 continue
5571155711
! ==================================================================
5571255712
end do loop_40
5571355713
! copy the lower band to ab which is the band storage matrix
5571455714
do j = n-kd+1, n
55715-
lk = min(kd, n-j) + 1
55716-
call stdlib_scopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
55715+
llk = min(kd, n-j) + 1
55716+
call stdlib_scopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
5571755717
end do
5571855718
end if
5571955719
work( 1 ) = lwmin

‎src/stdlib_linalg_lapack_w.fypp‎

Lines changed: 15 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -28868,7 +28868,7 @@ module stdlib_linalg_lapack_w
2886828868

2886928869
! Local Scalars
2887028870
logical(lk) :: lquery, upper
28871-
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, lk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
28871+
integer(ilp) :: i, j, iinfo, lwmin, pn, pk, llk, ldt, ldw, lds2, lds1, ls2, ls1, lw, lt,&
2887228872
tpos, wpos, s2pos, s1pos
2887328873
! Intrinsic Functions
2887428874
intrinsic :: min,max
@@ -28904,13 +28904,13 @@ module stdlib_linalg_lapack_w
2890428904
if( n<=kd+1 ) then
2890528905
if( upper ) then
2890628906
do i = 1, n
28907-
lk = min( kd+1, i )
28908-
call stdlib_wcopy( lk, a( i-lk+1, i ), 1,ab( kd+1-lk+1, i ), 1 )
28907+
llk = min( kd+1, i )
28908+
call stdlib_wcopy( llk, a( i-llk+1, i ), 1,ab( kd+1-llk+1, i ), 1 )
2890928909
end do
2891028910
else
2891128911
do i = 1, n
28912-
lk = min( kd+1, n-i+1 )
28913-
call stdlib_wcopy( lk, a( i, i ), 1, ab( 1, i ), 1 )
28912+
llk = min( kd+1, n-i+1 )
28913+
call stdlib_wcopy( llk, a( i, i ), 1, ab( 1, i ), 1 )
2891428914
end do
2891528915
endif
2891628916
work( 1 ) = 1
@@ -28947,8 +28947,8 @@ module stdlib_linalg_lapack_w
2894728947
iinfo )
2894828948
! copy the upper portion of a into ab
2894928949
do j = i, i+pk-1
28950-
lk = min( kd, n-j ) + 1
28951-
call stdlib_wcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
28950+
llk = min( kd, n-j ) + 1
28951+
call stdlib_wcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
2895228952
end do
2895328953
call stdlib_wlaset( 'LOWER', pk, pk, czero, cone,a( i, i+kd ), lda )
2895428954
! form the matrix t
@@ -28970,8 +28970,8 @@ module stdlib_linalg_lapack_w
2897028970
end do
2897128971
! copy the upper band to ab which is the band storage matrix
2897228972
do j = n-kd+1, n
28973-
lk = min(kd, n-j) + 1
28974-
call stdlib_wcopy( lk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
28973+
llk = min(kd, n-j) + 1
28974+
call stdlib_wcopy( llk, a( j, j ), lda, ab( kd+1, j ), ldab-1 )
2897528975
end do
2897628976
else
2897728977
! reduce the lower triangle of a to lower band matrix
@@ -28983,8 +28983,8 @@ module stdlib_linalg_lapack_w
2898328983
iinfo )
2898428984
! copy the upper portion of a into ab
2898528985
do j = i, i+pk-1
28986-
lk = min( kd, n-j ) + 1
28987-
call stdlib_wcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
28986+
llk = min( kd, n-j ) + 1
28987+
call stdlib_wcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
2898828988
end do
2898928989
call stdlib_wlaset( 'UPPER', pk, pk, czero, cone,a( i+kd, i ), lda )
2899028990
! form the matrix t
@@ -29006,15 +29006,15 @@ module stdlib_linalg_lapack_w
2900629006
! ==================================================================
2900729007
! restore a for comparison and checking to be removed
2900829008
! do 45 j = i, i+pk-1
29009-
! lk = min( kd, n-j ) + 1
29010-
! call stdlib_wcopy( lk, ab( 1, j ), 1, a( j, j ), 1 )
29009+
! llk = min( kd, n-j ) + 1
29010+
! call stdlib_wcopy( llk, ab( 1, j ), 1, a( j, j ), 1 )
2901129011
45 continue
2901229012
! ==================================================================
2901329013
end do loop_40
2901429014
! copy the lower band to ab which is the band storage matrix
2901529015
do j = n-kd+1, n
29016-
lk = min(kd, n-j) + 1
29017-
call stdlib_wcopy( lk, a( j, j ), 1, ab( 1, j ), 1 )
29016+
llk = min(kd, n-j) + 1
29017+
call stdlib_wcopy( llk, a( j, j ), 1, ab( 1, j ), 1 )
2901829018
end do
2901929019
end if
2902029020
work( 1 ) = lwmin

0 commit comments

Comments
 (0)