2355
2356
2357
2358
2359
2360
2361
2362
2363
2364
2365
2366
2367
2368
2369 INTEGER LDA, M, N
2370 CHARACTER*1 UPLO
2371 CHARACTER*2 TYPE
2372
2373 REAL AA( LDA, * ), AS( LDA, * )
2374
2375 INTEGER I, IBEG, IEND, J
2376 LOGICAL UPPER
2377
2378 upper = uplo.EQ.'U'
2379 IF( type.EQ.'GE' )THEN
2380 DO 20 j = 1, n
2381 DO 10 i = m + 1, lda
2382 IF( aa( i, j ).NE.as( i, j ) )
2383 $ GO TO 70
2384 10 CONTINUE
2385 20 CONTINUE
2386 ELSE IF( type.EQ.'SY' )THEN
2387 DO 50 j = 1, n
2388 IF( upper )THEN
2389 ibeg = 1
2390 iend = j
2391 ELSE
2392 ibeg = j
2393 iend = n
2394 END IF
2395 DO 30 i = 1, ibeg - 1
2396 IF( aa( i, j ).NE.as( i, j ) )
2397 $ GO TO 70
2398 30 CONTINUE
2399 DO 40 i = iend + 1, lda
2400 IF( aa( i, j ).NE.as( i, j ) )
2401 $ GO TO 70
2402 40 CONTINUE
2403 50 CONTINUE
2404 END IF
2405
2406 60 CONTINUE
2408 GO TO 80
2409 70 CONTINUE
2411 80 RETURN
2412
2413
2414
logical function lseres(type, uplo, m, n, aa, as, lda)