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