function [fp,fm] = fluvl(gam,r,u,a)
% Copyright G. Allaire, Octobre 2000
%
% calcul de la decomposition de flux de Van Leer pour les equations d'Euler
%
% r = densite, u = vitesse, a = vitesse du son
% gam = gamma de la loi d'etat du gaz parfait
% fp = flux(+) ,  fm = flux(-)
%
m = u./a  ;
aux = (2/gam)*a.*(((gam-1)/2)*m+1) ;
coef = gam^2/(2*(gam^2-1)) ;
fp(1,:) = 0.25*r.*a.*((1.+m).^2)  ; 
fp(2,:) = fp(1,:).*aux ;
fp(3,:) = fp(2,:).*aux*coef ;
aux = (2/gam)*a.*(((gam-1)/2)*m-1) ;
fm(1,:) = -0.25*r.*a.*(1.-m).^2  ; 
fm(2,:) = fm(1,:).*aux ;
fm(3,:) = fm(2,:).*aux*coef ;


