sumsq.c 1021 Bytes
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50
/* Compute the sum of the squares of a vector of signed shorts

 * Copyright 2004 Phil Karn, KA9Q
 * May be used under the terms of the GNU Lesser General Public License (LGPL)
 */

#include <stdlib.h>
#include "fec.h"

unsigned long long sumsq_port(signed short *,int);

#ifdef __i386__
unsigned long long sumsq_mmx(signed short *,int);
unsigned long long sumsq_sse(signed short *,int);
unsigned long long sumsq_sse2(signed short *,int);
#endif

#ifdef __x86_64__
unsigned long long sumsq_sse2(signed short *,int);
#endif

#ifdef __VEC__
unsigned long long sumsq_av(signed short *,int);
#endif

unsigned long long sumsq(signed short *in,int cnt){
  switch(Cpu_mode){
  case PORT:
  default:
    return sumsq_port(in,cnt);
#ifdef __i386__
  case SSE:
  case MMX:
    return sumsq_mmx(in,cnt);
  case SSE2:
    return sumsq_sse2(in,cnt);
#endif

#ifdef __x86_64__
  case SSE2:
    return sumsq_port(in,cnt);
    //return sumsq_sse2(in,cnt);
#endif

#ifdef __VEC__
  case ALTIVEC:
    return sumsq_av(in,cnt);
#endif
  }
}