shiftleft.m

来自「intlab 工具用于快速计算 各各层的倒数等等」· M 代码 · 共 25 行

function m = shiftleft(m,r)
%SHIFTLEFT    Shift array m left by r bits
%
%   m = shiftleft(m,r)
%
%first r bits of m(:,1) are zero, otherwise lost
%

% written  12/30/98     S.M. Rump
%

global INTLAB_LONG_LOGBETA

  Ones = ones(1,size(m,2));

  % first r bits of m
  factor = 2.^(r-INTLAB_LONG_LOGBETA)*Ones ;
  ms = floor( m .* factor );

  % array m shifted left by r bits, last r bits zero
  m =  ( m - ms ./ factor ) .* ( 2.^r*Ones );

  % add last r bits
  m(:,1:end-1) = m(:,1:end-1) + ms(:,2:end);