function [fp,fm] = fluvl(gam,r,u,a,fp,fm)
// 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 ;


endfunction