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