2373
2374
2375
2376
2377
2378
2379
2380
2381
2382
2383
2384
2385
2386
2387 INTEGER LDA, M, N
2388 CHARACTER*1 UPLO
2389 CHARACTER*2 TYPE
2390
2391 REAL AA( LDA, * ), AS( LDA, * )
2392
2393 INTEGER I, IBEG, IEND, J
2394 LOGICAL UPPER
2395
2396 upper = uplo.EQ.'U'
2397 IF( type.EQ.'GE' )THEN
2398 DO 20 j = 1, n
2399 DO 10 i = m + 1, lda
2400 IF( aa( i, j ).NE.as( i, j ) )
2401 $ GO TO 70
2402 10 CONTINUE
2403 20 CONTINUE
2404 ELSE IF( type.EQ.'SY' )THEN
2405 DO 50 j = 1, n
2406 IF( upper )THEN
2407 ibeg = 1
2408 iend = j
2409 ELSE
2410 ibeg = j
2411 iend = n
2412 END IF
2413 DO 30 i = 1, ibeg - 1
2414 IF( aa( i, j ).NE.as( i, j ) )
2415 $ GO TO 70
2416 30 CONTINUE
2417 DO 40 i = iend + 1, lda
2418 IF( aa( i, j ).NE.as( i, j ) )
2419 $ GO TO 70
2420 40 CONTINUE
2421 50 CONTINUE
2422 END IF
2423
2424 60 CONTINUE
2426 GO TO 80
2427 70 CONTINUE
2429 80 RETURN
2430
2431
2432
logical function lseres(type, uplo, m, n, aa, as, lda)