BI_dvvamx.c File Reference

#include "Bdef.h"

Go to the source code of this file.

Functions
void	BI_dvvamx (Int N, char *vec1, char *vec2)

Function Documentation

BI_dvvamx()

void BI_dvvamx	(	Int	N,
		char *	vec1,
		char *	vec2 )

Definition at line 2 of file BI_dvvamx.c.

{
   double *v1=(double*)vec1, *v2=(double*)vec2;
   double diff;
   BI_DistType *dist1, *dist2;
   Int i, k;
 
   k = N * sizeof(double);
   i = k % sizeof(BI_DistType);
   if (i) k += sizeof(BI_DistType) - i;
   dist1 = (BI_DistType *) &vec1[k];
   dist2 = (BI_DistType *) &vec2[k];
 
   for (k=0; k < N; k++)
   {
      diff = Rabs(v1[k]) - Rabs(v2[k]);
      if (diff < 0)
      {
         v1[k] = v2[k];
         dist1[k] = dist2[k];
      }
      else if (diff == 0)
      {
         if (dist1[k] > dist2[k])
         {
            v1[k] = v2[k];
            dist1[k] = dist2[k];
         }
      }
   }
}