@@ -36398,13 +36398,13 @@ module stdlib_linalg_lapack_c
3639836398 call stdlib_clacpy( "A", kd+1, n, ab, ldab, work( apos ), lda )
3639936399 call stdlib_claset( "A", kd, n, czero, czero, work( awpos ), lda )
3640036400 ! 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
3640836408 ! main bulge chasing loop
3640936409 loop_100: do thgrid = 1, thgrnb
3641036410 stt = (thgrid-1)*thgrsiz+1
@@ -36438,24 +36438,24 @@ module stdlib_linalg_lapack_c
3643836438 endif
3643936439 endif
3644036440 ! 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()
3645436454 call stdlib_chb2st_kernels( uplo, wantq, ttype,stind, edind, &
3645536455 sweepid, n, kd, ib,work ( inda ), lda,hous( indv ), hous( &
3645636456 indtau ), ldv,work( indw + tid*kd ) )
36457- !$OMP END TASK
36458- !$ endif
36457+ !$OMP END TASK
36458+ !$ endif
3645936459 if ( blklastind>=(n-1) ) then
3646036460 stt = stt + 1
3646136461 exit
@@ -36465,8 +36465,8 @@ module stdlib_linalg_lapack_c
3646536465 end do loop_120
3646636466 end do loop_110
3646736467 end do loop_100
36468- !$OMP END MASTER
36469- !$OMP END PARALLEL
36468+ !$OMP END MASTER
36469+ !$OMP END PARALLEL
3647036470 ! copy the diagonal from a to d. note that d is real thus only
3647136471 ! the real part is needed, the imaginary part should be czero.
3647236472 do i = 1, n
0 commit comments