LAPACK  3.6.1
LAPACK: Linear Algebra PACKage
subroutine dlacsg ( integer  M,
integer  P,
integer  Q,
double precision, dimension( * )  THETA,
integer, dimension( 4 )  ISEED,
double precision, dimension( ldx, * )  X,
integer  LDX,
double precision, dimension( * )  WORK 
)

Definition at line 353 of file dckcsd.f.

353  IMPLICIT NONE
354 *
355  INTEGER ldx, m, p, q
356  INTEGER iseed( 4 )
357  DOUBLE PRECISION theta( * )
358  DOUBLE PRECISION work( * ), x( ldx, * )
359 *
360  DOUBLE PRECISION one, zero
361  parameter ( one = 1.0d0, zero = 0.0d0 )
362 *
363  INTEGER i, info, r
364 *
365  r = min( p, m-p, q, m-q )
366 *
367  CALL dlaset( 'Full', m, m, zero, zero, x, ldx )
368 *
369  DO i = 1, min(p,q)-r
370  x(i,i) = one
371  END DO
372  DO i = 1, r
373  x(min(p,q)-r+i,min(p,q)-r+i) = cos(theta(i))
374  END DO
375  DO i = 1, min(p,m-q)-r
376  x(p-i+1,m-i+1) = -one
377  END DO
378  DO i = 1, r
379  x(p-(min(p,m-q)-r)+1-i,m-(min(p,m-q)-r)+1-i) =
380  $ -sin(theta(r-i+1))
381  END DO
382  DO i = 1, min(m-p,q)-r
383  x(m-i+1,q-i+1) = one
384  END DO
385  DO i = 1, r
386  x(m-(min(m-p,q)-r)+1-i,q-(min(m-p,q)-r)+1-i) =
387  $ sin(theta(r-i+1))
388  END DO
389  DO i = 1, min(m-p,m-q)-r
390  x(p+i,q+i) = one
391  END DO
392  DO i = 1, r
393  x(p+(min(m-p,m-q)-r)+i,q+(min(m-p,m-q)-r)+i) =
394  $ cos(theta(i))
395  END DO
396  CALL dlaror( 'Left', 'No init', p, m, x, ldx, iseed, work, info )
397  CALL dlaror( 'Left', 'No init', m-p, m, x(p+1,1), ldx,
398  $ iseed, work, info )
399  CALL dlaror( 'Right', 'No init', m, q, x, ldx, iseed,
400  $ work, info )
401  CALL dlaror( 'Right', 'No init', m, m-q,
402  $ x(1,q+1), ldx, iseed, work, info )
403 *
subroutine dlaset(UPLO, M, N, ALPHA, BETA, A, LDA)
DLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
Definition: dlaset.f:112
subroutine dlaror(SIDE, INIT, M, N, A, LDA, ISEED, X, INFO)
DLAROR
Definition: dlaror.f:148

Here is the call graph for this function:

Here is the caller graph for this function: