3271
3272
3273
3274
3275
3276
3277
3278
3279
3280
3281
3282
3283
3284
3285 INTEGER LDA, M, N
3286 CHARACTER*1 UPLO
3287 CHARACTER*2 TYPE
3288
3289 COMPLEX AA( LDA, * ), AS( LDA, * )
3290
3291 INTEGER I, IBEG, IEND, J
3292 LOGICAL UPPER
3293
3294 upper = uplo.EQ.'U'
3295 IF( type.EQ.'GE' )THEN
3296 DO 20 j = 1, n
3297 DO 10 i = m + 1, lda
3298 IF( aa( i, j ).NE.as( i, j ) )
3299 $ GO TO 70
3300 10 CONTINUE
3301 20 CONTINUE
3302 ELSE IF( type.EQ.'HE'.OR.type.EQ.'SY' )THEN
3303 DO 50 j = 1, n
3304 IF( upper )THEN
3305 ibeg = 1
3306 iend = j
3307 ELSE
3308 ibeg = j
3309 iend = n
3310 END IF
3311 DO 30 i = 1, ibeg - 1
3312 IF( aa( i, j ).NE.as( i, j ) )
3313 $ GO TO 70
3314 30 CONTINUE
3315 DO 40 i = iend + 1, lda
3316 IF( aa( i, j ).NE.as( i, j ) )
3317 $ GO TO 70
3318 40 CONTINUE
3319 50 CONTINUE
3320 END IF
3321
3323 GO TO 80
3324 70 CONTINUE
3326 80 RETURN
3327
3328
3329
logical function lceres(type, uplo, m, n, aa, as, lda)