#include"pluto.h"

/* ******************************************************************** */
void VelocityLimiter (double *v, double *vp, double *vm)
/*
 *
 *  Check whether the total reconstructed velocity is > 1
 *  If a superluminal value occurs, flatten distribution.
 *
 *
 ********************************************************************** */
{
  #if USE_FOUR_VELOCITY == NO && COMPONENTS > 1
   int    nv;
   double v2m, v2p;

   v2m = EXPAND(vm[VX1]*vm[VX1],  + vm[VX2]*vm[VX2], + vm[VX3]*vm[VX3]);
   v2p = EXPAND(vp[VX1]*vp[VX1],  + vp[VX2]*vp[VX2], + vp[VX3]*vp[VX3]);
   if (v2m >= 1.0 || v2p >= 1.0){
     for (nv = NVAR; nv--;  ) vm[nv] = vp[nv] = v[nv];  
   }
  #endif
}