2746
2747
2748
2749
2750
2751
2752
2753
2754
2755
2756
2757
2758
2759
2760 INTEGER LDA, M, N
2761 CHARACTER*1 UPLO
2762 CHARACTER*2 TYPE
2763
2764 DOUBLE PRECISION AA( LDA, * ), AS( LDA, * )
2765
2766 INTEGER I, IBEG, IEND, J
2767 LOGICAL UPPER
2768
2769 upper = uplo.EQ.'U'
2770 IF( type.EQ.'GE' )THEN
2771 DO 20 j = 1, n
2772 DO 10 i = m + 1, lda
2773 IF( aa( i, j ).NE.as( i, j ) )
2774 $ GO TO 70
2775 10 CONTINUE
2776 20 CONTINUE
2777 ELSE IF( type.EQ.'SY' )THEN
2778 DO 50 j = 1, n
2779 IF( upper )THEN
2780 ibeg = 1
2781 iend = j
2782 ELSE
2783 ibeg = j
2784 iend = n
2785 END IF
2786 DO 30 i = 1, ibeg - 1
2787 IF( aa( i, j ).NE.as( i, j ) )
2788 $ GO TO 70
2789 30 CONTINUE
2790 DO 40 i = iend + 1, lda
2791 IF( aa( i, j ).NE.as( i, j ) )
2792 $ GO TO 70
2793 40 CONTINUE
2794 50 CONTINUE
2795 END IF
2796
2798 GO TO 80
2799 70 CONTINUE
2801 80 RETURN
2802
2803
2804
logical function lderes(type, uplo, m, n, aa, as, lda)