@@ -36398,13 +36398,13 @@ module stdlib_linalg_lapack_c
36398
36398
call stdlib_clacpy( "A", kd+1, n, ab, ldab, work( apos ), lda )
36399
36399
call stdlib_claset( "A", kd, n, czero, czero, work( awpos ), lda )
36400
36400
! openmp parallelisation start here
36401
- !$OMP PARALLEL PRIVATE( TID, THGRID, BLKLASTIND ) &
36402
- !$OMP& PRIVATE( THED, I, M, K, ST, ED, STT, SWEEPID ) &
36403
- !$OMP& PRIVATE( MYID, TTYPE, COLPT, STIND, EDIND ) &
36404
- !$OMP& SHARED ( UPLO, WANTQ, INDV, INDTAU, HOUS, WORK) &
36405
- !$OMP& SHARED ( N, KD, IB, NBTILES, LDA, LDV, INDA ) &
36406
- !$OMP& SHARED ( STEPERCOL, THGRNB, THGRSIZ, GRSIZ, SHIFT )
36407
- !$OMP MASTER
36401
+ !$OMP PARALLEL PRIVATE( TID, THGRID, BLKLASTIND ) &
36402
+ !$OMP& PRIVATE( THED, I, M, K, ST, ED, STT, SWEEPID ) &
36403
+ !$OMP& PRIVATE( MYID, TTYPE, COLPT, STIND, EDIND ) &
36404
+ !$OMP& SHARED ( UPLO, WANTQ, INDV, INDTAU, HOUS, WORK) &
36405
+ !$OMP& SHARED ( N, KD, IB, NBTILES, LDA, LDV, INDA ) &
36406
+ !$OMP& SHARED ( STEPERCOL, THGRNB, THGRSIZ, GRSIZ, SHIFT )
36407
+ !$OMP MASTER
36408
36408
! main bulge chasing loop
36409
36409
loop_100: do thgrid = 1, thgrnb
36410
36410
stt = (thgrid-1)*thgrsiz+1
@@ -36438,24 +36438,24 @@ module stdlib_linalg_lapack_c
36438
36438
endif
36439
36439
endif
36440
36440
! call the kernel
36441
- !$ if( ttype/=1 ) then
36442
- !$OMP TASK DEPEND(in:WORK(MYID+SHIFT-1)) &
36443
- !$OMP& DEPEND(in:WORK(MYID-1)) &
36444
- !$OMP& DEPEND(out:WORK(MYID))
36445
- !$ tid = omp_get_thread_num()
36446
- !$ call stdlib_chb2st_kernels( uplo, wantq, ttype,stind, edind, &
36447
- !$ sweepid, n, kd, ib,work ( inda ), lda,hous( indv ), hous( &
36448
- !$ indtau ), ldv,work( indw + tid*kd ) )
36449
- !$OMP END TASK
36450
- !$ else
36451
- !$OMP TASK DEPEND(in:WORK(MYID+SHIFT-1)) &
36452
- !$OMP& DEPEND(out:WORK(MYID))
36453
- !$ tid = omp_get_thread_num()
36441
+ !$ if( ttype/=1 ) then
36442
+ !$OMP TASK DEPEND(in:WORK(MYID+SHIFT-1)) &
36443
+ !$OMP& DEPEND(in:WORK(MYID-1)) &
36444
+ !$OMP& DEPEND(out:WORK(MYID))
36445
+ !$ tid = omp_get_thread_num()
36446
+ !$ call stdlib_chb2st_kernels( uplo, wantq, ttype,stind, edind, &
36447
+ !$ sweepid, n, kd, ib,work ( inda ), lda,hous( indv ), hous( &
36448
+ !$ indtau ), ldv,work( indw + tid*kd ) )
36449
+ !$OMP END TASK
36450
+ !$ else
36451
+ !$OMP TASK DEPEND(in:WORK(MYID+SHIFT-1)) &
36452
+ !$OMP& DEPEND(out:WORK(MYID))
36453
+ !$ tid = omp_get_thread_num()
36454
36454
call stdlib_chb2st_kernels( uplo, wantq, ttype,stind, edind, &
36455
36455
sweepid, n, kd, ib,work ( inda ), lda,hous( indv ), hous( &
36456
36456
indtau ), ldv,work( indw + tid*kd ) )
36457
- !$OMP END TASK
36458
- !$ endif
36457
+ !$OMP END TASK
36458
+ !$ endif
36459
36459
if ( blklastind>=(n-1) ) then
36460
36460
stt = stt + 1
36461
36461
exit
@@ -36465,8 +36465,8 @@ module stdlib_linalg_lapack_c
36465
36465
end do loop_120
36466
36466
end do loop_110
36467
36467
end do loop_100
36468
- !$OMP END MASTER
36469
- !$OMP END PARALLEL
36468
+ !$OMP END MASTER
36469
+ !$OMP END PARALLEL
36470
36470
! copy the diagonal from a to d. note that d is real thus only
36471
36471
! the real part is needed, the imaginary part should be czero.
36472
36472
do i = 1, n
0 commit comments