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