[BACK]Return to dlasyf_aa.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dlasyf_aa.f between versions 1.4 and 1.5

version 1.4, 2018/05/29 07:18:01 version 1.5, 2020/05/21 21:46:00
Line 284 Line 284
 *  *
 *              Swap A(I1, I2+1:M) with A(I2, I2+1:M)  *              Swap A(I1, I2+1:M) with A(I2, I2+1:M)
 *  *
                CALL DSWAP( M-I2, A( J1+I1-1, I2+1 ), LDA,                 IF( I2.LT.M )
      $                           A( J1+I2-1, I2+1 ), LDA )       $            CALL DSWAP( M-I2, A( J1+I1-1, I2+1 ), LDA,
        $                              A( J1+I2-1, I2+1 ), LDA )
 *  *
 *              Swap A(I1, I1) with A(I2,I2)  *              Swap A(I1, I1) with A(I2,I2)
 *  *
Line 325 Line 326
 *           Compute L(J+2, J+1) = WORK( 3:M ) / T(J, J+1),  *           Compute L(J+2, J+1) = WORK( 3:M ) / T(J, J+1),
 *            where A(J, J+1) = T(J, J+1) and A(J+2:M, J) = L(J+2:M, J+1)  *            where A(J, J+1) = T(J, J+1) and A(J+2:M, J) = L(J+2:M, J+1)
 *  *
             IF( A( K, J+1 ).NE.ZERO ) THEN              IF( J.LT.(M-1) ) THEN
                ALPHA = ONE / A( K, J+1 )                 IF( A( K, J+1 ).NE.ZERO ) THEN
                CALL DCOPY( M-J-1, WORK( 3 ), 1, A( K, J+2 ), LDA )                    ALPHA = ONE / A( K, J+1 )
                CALL DSCAL( M-J-1, ALPHA, A( K, J+2 ), LDA )                    CALL DCOPY( M-J-1, WORK( 3 ), 1, A( K, J+2 ), LDA )
             ELSE                    CALL DSCAL( M-J-1, ALPHA, A( K, J+2 ), LDA )
                CALL DLASET( 'Full', 1, M-J-1, ZERO, ZERO,                 ELSE
      $                      A( K, J+2 ), LDA)                    CALL DLASET( 'Full', 1, M-J-1, ZERO, ZERO,
        $                         A( K, J+2 ), LDA)
                  END IF
             END IF              END IF
          END IF           END IF
          J = J + 1           J = J + 1
Line 432 Line 435
 *  *
 *              Swap A(I2+1:M, I1) with A(I2+1:M, I2)  *              Swap A(I2+1:M, I1) with A(I2+1:M, I2)
 *  *
                CALL DSWAP( M-I2, A( I2+1, J1+I1-1 ), 1,                 IF( I2.LT.M )
      $                           A( I2+1, J1+I2-1 ), 1 )       $            CALL DSWAP( M-I2, A( I2+1, J1+I1-1 ), 1,
        $                              A( I2+1, J1+I2-1 ), 1 )
 *  *
 *              Swap A(I1, I1) with A(I2, I2)  *              Swap A(I1, I1) with A(I2, I2)
 *  *
Line 473 Line 477
 *           Compute L(J+2, J+1) = WORK( 3:M ) / T(J, J+1),  *           Compute L(J+2, J+1) = WORK( 3:M ) / T(J, J+1),
 *            where A(J, J+1) = T(J, J+1) and A(J+2:M, J) = L(J+2:M, J+1)  *            where A(J, J+1) = T(J, J+1) and A(J+2:M, J) = L(J+2:M, J+1)
 *  *
             IF( A( J+1, K ).NE.ZERO ) THEN              IF( J.LT.(M-1) ) THEN
                ALPHA = ONE / A( J+1, K )                 IF( A( J+1, K ).NE.ZERO ) THEN
                CALL DCOPY( M-J-1, WORK( 3 ), 1, A( J+2, K ), 1 )                    ALPHA = ONE / A( J+1, K )
                CALL DSCAL( M-J-1, ALPHA, A( J+2, K ), 1 )                    CALL DCOPY( M-J-1, WORK( 3 ), 1, A( J+2, K ), 1 )
             ELSE                    CALL DSCAL( M-J-1, ALPHA, A( J+2, K ), 1 )
                CALL DLASET( 'Full', M-J-1, 1, ZERO, ZERO,                 ELSE
      $                      A( J+2, K ), LDA )                    CALL DLASET( 'Full', M-J-1, 1, ZERO, ZERO,
        $                         A( J+2, K ), LDA )
                  END IF
             END IF              END IF
          END IF           END IF
          J = J + 1           J = J + 1
Removed from v.1.4  
changed lines
  Added in v.1.5

CVSweb interface <joel.bertrand@systella.fr>