vfp9.c

三星s3c2460开发板完整功能测试代码

#include <stdio.h>
#include <math.h>


/*
void FFT(int dir, long m, complex <double> x[])
{
   long i, i1, i2,j, k, l, l1, l2, n;
   complex <double> tx, t1, u, c;

   //Calculate the number of points
   n = 1;
   for(i = 0; i < m; i++) 
      n <<= 1;   

   // Do the bit reversal
   i2 = n >> 1;
   j = 0;

   for (i = 0; i < n-1 ; i++)
   {
      if (i < j)
         swap(x[i], x[j]);

      k = i2;

      while (k <= j) 
	  {
         j -= k;
         k >>= 1;
      }

      j += k;
   }

   // Compute the FFT 
   c.real(-1.0);
   c.imag(0.0);
   l2 = 1;
   for (l = 0; l < m; l++) 
   {
      l1 = l2;
      l2 <<= 1;
      u.real(1.0);
      u.imag(0.0);

      for (j = 0; j < l1; j++) 
	  {
         for (i = j; i < n; i += l2) 
		 {
            i1 = i + l1;
            t1 = u * x[i1];
            x[i1] = x[i] - t1; 
            x[i] += t1;
         }

         u = u * c;
      }

      c.imag(sqrt((1.0 - c.real()) / 2.0));
      if (dir == 1)
         c.imag(-c.imag());
      c.real(sqrt((1.0 + c.real()) / 2.0));
   }

   // Scaling for forward transform
   if (dir == 1) 
   {
      for (i = 0; i < n; i++)
         x[i] /= n;      
   }   
   return;
}
*/

void VFP9( void)
{
	double	a,b,c,d;
	int 	A,B,C,D;
	
	a=1.1111;
	b=2.2222;
	c=3.3333;

	printf("%f %f %f\n",a,b,c);
	
	d=a*b;
	d=b*c;
	
	d=c-(a+b);
	
	printf("d: %f\n",d);
	
	d=b/(c-(a+b));
	
	printf("d: %f\n",d);


	A=(int)(a);
	B=(int)(b);
	C=(int)(c);
	
	D=C-(A+B);
	
	D=B/D;
	printf("D: %d\n",D);
}