% Creating a redistribution matrix function redist = redist_ab(dist,alpha,xmax,interval,adv) C = normalize_lbm_abs_advective(xmax,interval,alpha,adv); if alpha > 0 redist = exp(-((dist./alpha).^2))./C; else redist = zeros(length(dist),length(dist)); for z = 1:length(dist) redist(z,z) = 1; end end lambda = max(eig(redist)); redist = redist./lambda;