2670
2671
2672
2673
2674
2675
2676
2677
2678
2679
2680
2681
2682
2683
2684 INTEGER LDA, M, N
2685 CHARACTER*1 UPLO
2686 CHARACTER*2 TYPE
2687
2688 COMPLEX AA( LDA, * ), AS( LDA, * )
2689
2690 INTEGER I, IBEG, IEND, J
2691 LOGICAL UPPER
2692
2693 upper = uplo.EQ.'U'
2694 IF( type.EQ.'ge' )THEN
2695 DO 20 j = 1, n
2696 DO 10 i = m + 1, lda
2697 IF( aa( i, j ).NE.as( i, j ) )
2698 $ GO TO 70
2699 10 CONTINUE
2700 20 CONTINUE
2701 ELSE IF( type.EQ.'he'.OR.type.EQ.'sy' )THEN
2702 DO 50 j = 1, n
2703 IF( upper )THEN
2704 ibeg = 1
2705 iend = j
2706 ELSE
2707 ibeg = j
2708 iend = n
2709 END IF
2710 DO 30 i = 1, ibeg - 1
2711 IF( aa( i, j ).NE.as( i, j ) )
2712 $ GO TO 70
2713 30 CONTINUE
2714 DO 40 i = iend + 1, lda
2715 IF( aa( i, j ).NE.as( i, j ) )
2716 $ GO TO 70
2717 40 CONTINUE
2718 50 CONTINUE
2719 END IF
2720
2721 60 CONTINUE
2723 GO TO 80
2724 70 CONTINUE
2726 80 RETURN
2727
2728
2729
logical function lceres(type, uplo, m, n, aa, as, lda)