2620
2621
2622
2623
2624
2625
2626
2627
2628
2629
2630
2631
2632 INTEGER LDA, M, N
2633 CHARACTER*1 UPLO
2634 CHARACTER*2 TYPE
2635
2636 COMPLEX*16 AA( LDA, * ), AS( LDA, * )
2637
2638 INTEGER I, IBEG, IEND, J
2639 LOGICAL UPPER
2640
2641 upper = uplo.EQ.'U'
2642 IF( type.EQ.'ge' )THEN
2643 DO 20 j = 1, n
2644 DO 10 i = m + 1, lda
2645 IF( aa( i, j ).NE.as( i, j ) )
2646 $ GO TO 70
2647 10 CONTINUE
2648 20 CONTINUE
2649 ELSE IF( type.EQ.'he' )THEN
2650 DO 50 j = 1, n
2651 IF( upper )THEN
2652 ibeg = 1
2653 iend = j
2654 ELSE
2655 ibeg = j
2656 iend = n
2657 END IF
2658 DO 30 i = 1, ibeg - 1
2659 IF( aa( i, j ).NE.as( i, j ) )
2660 $ GO TO 70
2661 30 CONTINUE
2662 DO 40 i = iend + 1, lda
2663 IF( aa( i, j ).NE.as( i, j ) )
2664 $ GO TO 70
2665 40 CONTINUE
2666 50 CONTINUE
2667 END IF
2668
2669 60 CONTINUE
2671 GO TO 80
2672 70 CONTINUE
2674 80 RETURN
2675
2676
2677
logical function lzeres(type, uplo, m, n, aa, as, lda)