3486
3487
3488
3489
3490
3491
3492
3493
3494
3495
3496
3497
3498
3499
3500 INTEGER LDA, M, N
3501 CHARACTER*1 UPLO
3502 CHARACTER*2 TYPE
3503
3504 COMPLEX*16 AA( LDA, * ), AS( LDA, * )
3505
3506 INTEGER I, IBEG, IEND, J
3507 LOGICAL UPPER
3508
3509 upper = uplo.EQ.'U'
3510 IF( type.EQ.'GE' )THEN
3511 DO 20 j = 1, n
3512 DO 10 i = m + 1, lda
3513 IF( aa( i, j ).NE.as( i, j ) )
3514 $ GO TO 70
3515 10 CONTINUE
3516 20 CONTINUE
3517 ELSE IF( type.EQ.'HE'.OR.type.EQ.'SY' )THEN
3518 DO 50 j = 1, n
3519 IF( upper )THEN
3520 ibeg = 1
3521 iend = j
3522 ELSE
3523 ibeg = j
3524 iend = n
3525 END IF
3526 DO 30 i = 1, ibeg - 1
3527 IF( aa( i, j ).NE.as( i, j ) )
3528 $ GO TO 70
3529 30 CONTINUE
3530 DO 40 i = iend + 1, lda
3531 IF( aa( i, j ).NE.as( i, j ) )
3532 $ GO TO 70
3533 40 CONTINUE
3534 50 CONTINUE
3535 END IF
3536
3538 GO TO 80
3539 70 CONTINUE
3541 80 RETURN
3542
3543
3544
logical function lzeres(type, uplo, m, n, aa, as, lda)