ScaLAPACK 2.1  2.1 ScaLAPACK: Scalable Linear Algebra PACKage
BI_cvvamx.c
Go to the documentation of this file.
1 #include "Bdef.h"
2 void BI_cvvamx(int N, char *vec1, char *vec2)
3 {
4  SCOMPLEX *v1=(SCOMPLEX*)vec1, *v2=(SCOMPLEX*)vec2;
5  float diff;
6  BI_DistType *dist1, *dist2;
7  int i, k;
8
9  k = N * sizeof(SCOMPLEX);
10  i = k % sizeof(BI_DistType);
11  if (i) k += sizeof(BI_DistType) - i;
12  dist1 = (BI_DistType *) &vec1[k];
13  dist2 = (BI_DistType *) &vec2[k];
14
15  for (k=0; k < N; k++)
16  {
17  diff = Cabs(v1[k]) - Cabs(v2[k]);
18  if (diff < 0)
19  {
20  v1[k].r = v2[k].r;
21  v1[k].i = v2[k].i;
22  dist1[k] = dist2[k];
23  }
24  else if (diff == 0)
25  {
26  if (dist1[k] > dist2[k])
27  {
28  v1[k].r = v2[k].r;
29  v1[k].i = v2[k].i;
30  dist1[k] = dist2[k];
31  }
32  }
33  }
34 }
SCOMPLEX::r
float r
Definition: Bdef.h:38
BI_cvvamx
void BI_cvvamx(int N, char *vec1, char *vec2)
Definition: BI_cvvamx.c:2
SCOMPLEX
Definition: Bdef.h:38
Bdef.h
Cabs
#define Cabs(z)
Definition: Bdef.h:270
BI_DistType
#define BI_DistType
Definition: Bdef.h:72
SCOMPLEX::i
float i
Definition: Bdef.h:38