arm_dotproduct_example_f32.c

/* ---------------------------------------------------------------------- 
* Copyright (C) 2010 ARM Limited. All rights reserved.   
*  
* $Date:        29. November 2010  
* $Revision:    V1.0.3 
*  
* Project:          CMSIS DSP Library  
* Title:            arm_dotproduct_example_f32.c                  
*  
* Description:  Example code computing dot product of two vectors. 
* 
* Target Processor: Cortex-M4/Cortex-M3  
*
*
* Version 1.0.3 2010/11/29 
*    Re-organized the CMSIS folders and updated documentation. 
* 
* Version 1.0.1 2010/10/05 KK 
*    Production release and review comments incorporated.  
*
* Version 1.0.0 2010/09/20 KK
*    Production release and review comments incorporated.
* ------------------------------------------------------------------- */ 
 
#include <math.h>     
#include "arm_math.h" 
 
/* ---------------------------------------------------------------------- 
* Defines each of the tests performed 
* ------------------------------------------------------------------- */ 
#define MAX_BLOCKSIZE   32 
#define DELTA           (0.000001f) 
 
/* ---------------------------------------------------------------------- 
* Test input data for Floating point Dot Product example for 32-blockSize 
* Generated by the MATLAB randn() function 
* ------------------------------------------------------------------- */  
/* ----------------------------------------------------------------------  
** Test input data of srcA for blockSize 32   
** ------------------------------------------------------------------- */  
float32_t srcA_buf_f32[MAX_BLOCKSIZE] =   
{   
-0.4325648115282207,    -1.6655843782380970,    0.1253323064748307,      
 0.2876764203585489,    -1.1464713506814637,    1.1909154656429988,      
 1.1891642016521031,    -0.0376332765933176,    0.3272923614086541,      
 0.1746391428209245,    -0.1867085776814394,    0.7257905482933027,      
-0.5883165430141887,     2.1831858181971011,   -0.1363958830865957,      
 0.1139313135208096,     1.0667682113591888,    0.0592814605236053,      
-0.0956484054836690,    -0.8323494636500225,    0.2944108163926404,      
-1.3361818579378040,     0.7143245518189522,    1.6235620644462707,      
-0.6917757017022868,     0.8579966728282626,    1.2540014216025324,      
-1.5937295764474768,    -1.4409644319010200,    0.5711476236581780,      
-0.3998855777153632,     0.6899973754643451 
};   
  
/* ----------------------------------------------------------------------  
** Test input data of srcB for blockSize 32   
** ------------------------------------------------------------------- */   
float32_t srcB_buf_f32[MAX_BLOCKSIZE] =   
{   
 1.7491401329284098,    0.1325982188803279,      0.3252281811989881,     
-0.7938091410349637,    0.3149236145048914,     -0.5272704888029532,     
 0.9322666565031119,    1.1646643544607362,     -2.0456694357357357,     
-0.6443728590041911,    1.7410657940825480,      0.4867684246821860,     
 1.0488288293660140,    1.4885752747099299,      1.2705014969484090,     
-1.8561241921210170,    2.1343209047321410,  1.4358467535865909,         
-0.9173023332875400,   -1.1060770780029008,      0.8105708062681296,     
 0.6985430696369063,   -0.4015827425012831,      1.2687512030669628,     
-0.7836083053674872,    0.2132664971465569,      0.7878984786088954,     
 0.8966819356782295,   -0.1869172943544062,      1.0131816724341454,     
 0.2484350696132857,    0.0596083377937976 
};   
 
/* Reference dot product output */ 
float32_t  refDotProdOut = 5.9273644806352142;   
 
/* ---------------------------------------------------------------------- 
* Declare Global variables  
* ------------------------------------------------------------------- */ 
float32_t multOutput[MAX_BLOCKSIZE];  /* Intermediate output */ 
float32_t testOutput;  /* Final ouput */ 
 
arm_status status;       /* Status of the example */ 

int32_t main(void) 
{ 
        uint32_t i;                      /* Loop counter */ 
        float32_t diff;          /* Difference between reference and test outputs */ 
 
        /* Multiplication of two input buffers */ 
        arm_mult_f32(srcA_buf_f32, srcB_buf_f32, multOutput, MAX_BLOCKSIZE); 
         
        /* Accumulate the multiplication output values to  
           get the dot product of the two inputs */ 
        for(i=0; i< MAX_BLOCKSIZE; i++) 
    {          
                arm_add_f32(&testOutput, &multOutput[i], &testOutput, 1);        
    } 
 
        /* absolute value of difference between ref and test */ 
        diff = fabsf(refDotProdOut - testOutput); 
         
        /* Comparison of dot product value with reference */ 
        if(diff > DELTA) 
        { 
                status = ARM_MATH_TEST_FAILURE; 
        } 
                 
        if( status == ARM_MATH_TEST_FAILURE) 
        { 
          while(1); 
        } 

    while(1);                             /* main function does not return */
}