BI_dvvamx2.c

Go to the documentation of this file.

#include "Bdef.h"
void BI_dvvamx2(Int N, char *vec1, char *vec2)
{
   Int k;
   double *v1=(double*)vec1, *v2=(double*)vec2;
   double diff;
 
   for (k=0; k != N; k++)
   {
      diff = Rabs(v1[k]) - Rabs(v2[k]);
      if (diff < 0) v1[k] = v2[k];
      else if (diff == 0) if (v1[k] < v2[k]) v1[k] = v2[k];
   }
}