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