arm_dotproduct_example_f32.c
Go to the documentation of this file.
00001 /* ---------------------------------------------------------------------- 00002 * Copyright (C) 2010 ARM Limited. All rights reserved. 00003 * 00004 * $Date: 29. November 2010 00005 * $Revision: V1.0.3 00006 * 00007 * Project: CMSIS DSP Library 00008 * Title: arm_dotproduct_example_f32.c 00009 * 00010 * Description: Example code computing dot product of two vectors. 00011 * 00012 * Target Processor: Cortex-M4/Cortex-M3 00013 * 00014 * 00015 * Version 1.0.3 2010/11/29 00016 * Re-organized the CMSIS folders and updated documentation. 00017 * 00018 * Version 1.0.1 2010/10/05 KK 00019 * Production release and review comments incorporated. 00020 * 00021 * Version 1.0.0 2010/09/20 KK 00022 * Production release and review comments incorporated. 00023 * ------------------------------------------------------------------- */ 00024 00070 #include <math.h> 00071 #include "arm_math.h" 00072 00073 /* ---------------------------------------------------------------------- 00074 * Defines each of the tests performed 00075 * ------------------------------------------------------------------- */ 00076 #define MAX_BLOCKSIZE 32 00077 #define DELTA (0.000001f) 00078 00079 /* ---------------------------------------------------------------------- 00080 * Test input data for Floating point Dot Product example for 32-blockSize 00081 * Generated by the MATLAB randn() function 00082 * ------------------------------------------------------------------- */ 00083 /* ---------------------------------------------------------------------- 00084 ** Test input data of srcA for blockSize 32 00085 ** ------------------------------------------------------------------- */ 00086 float32_t srcA_buf_f32[MAX_BLOCKSIZE] = 00087 { 00088 -0.4325648115282207, -1.6655843782380970, 0.1253323064748307, 00089 0.2876764203585489, -1.1464713506814637, 1.1909154656429988, 00090 1.1891642016521031, -0.0376332765933176, 0.3272923614086541, 00091 0.1746391428209245, -0.1867085776814394, 0.7257905482933027, 00092 -0.5883165430141887, 2.1831858181971011, -0.1363958830865957, 00093 0.1139313135208096, 1.0667682113591888, 0.0592814605236053, 00094 -0.0956484054836690, -0.8323494636500225, 0.2944108163926404, 00095 -1.3361818579378040, 0.7143245518189522, 1.6235620644462707, 00096 -0.6917757017022868, 0.8579966728282626, 1.2540014216025324, 00097 -1.5937295764474768, -1.4409644319010200, 0.5711476236581780, 00098 -0.3998855777153632, 0.6899973754643451 00099 }; 00100 00101 /* ---------------------------------------------------------------------- 00102 ** Test input data of srcB for blockSize 32 00103 ** ------------------------------------------------------------------- */ 00104 float32_t srcB_buf_f32[MAX_BLOCKSIZE] = 00105 { 00106 1.7491401329284098, 0.1325982188803279, 0.3252281811989881, 00107 -0.7938091410349637, 0.3149236145048914, -0.5272704888029532, 00108 0.9322666565031119, 1.1646643544607362, -2.0456694357357357, 00109 -0.6443728590041911, 1.7410657940825480, 0.4867684246821860, 00110 1.0488288293660140, 1.4885752747099299, 1.2705014969484090, 00111 -1.8561241921210170, 2.1343209047321410, 1.4358467535865909, 00112 -0.9173023332875400, -1.1060770780029008, 0.8105708062681296, 00113 0.6985430696369063, -0.4015827425012831, 1.2687512030669628, 00114 -0.7836083053674872, 0.2132664971465569, 0.7878984786088954, 00115 0.8966819356782295, -0.1869172943544062, 1.0131816724341454, 00116 0.2484350696132857, 0.0596083377937976 00117 }; 00118 00119 /* Reference dot product output */ 00120 float32_t refDotProdOut = 5.9273644806352142; 00121 00122 /* ---------------------------------------------------------------------- 00123 * Declare Global variables 00124 * ------------------------------------------------------------------- */ 00125 float32_t multOutput[MAX_BLOCKSIZE]; /* Intermediate output */ 00126 float32_t testOutput; /* Final ouput */ 00127 00128 arm_status status; /* Status of the example */ 00129 00130 int32_t main(void) 00131 { 00132 uint32_t i; /* Loop counter */ 00133 float32_t diff; /* Difference between reference and test outputs */ 00134 00135 /* Multiplication of two input buffers */ 00136 arm_mult_f32(srcA_buf_f32, srcB_buf_f32, multOutput, MAX_BLOCKSIZE); 00137 00138 /* Accumulate the multiplication output values to 00139 get the dot product of the two inputs */ 00140 for(i=0; i< MAX_BLOCKSIZE; i++) 00141 { 00142 arm_add_f32(&testOutput, &multOutput[i], &testOutput, 1); 00143 } 00144 00145 /* absolute value of difference between ref and test */ 00146 diff = fabsf(refDotProdOut - testOutput); 00147 00148 /* Comparison of dot product value with reference */ 00149 if(diff > DELTA) 00150 { 00151 status = ARM_MATH_TEST_FAILURE; 00152 } 00153 00154 if( status == ARM_MATH_TEST_FAILURE) 00155 { 00156 while(1); 00157 } 00158 } 00159
