Mercurial > nand-ecoscentric
diff packages/kernel/current/tests/fpinttestf2.c @ 3293:f35f92c0b860
Add single precision floating point tests.
| author | vae |
|---|---|
| date | Sat, 09 Mar 2013 18:06:46 +0000 |
| parents | |
| children |
line wrap: on
line diff
new file mode 100644 --- /dev/null +++ b/packages/kernel/current/tests/fpinttestf2.c @@ -0,0 +1,400 @@ +//========================================================================== +// +// fpinttestf2.c +// +// Basic FPU integrity test +// +//========================================================================== +// ####ECOSGPLCOPYRIGHTBEGIN#### +// ------------------------------------------- +// This file is part of eCos, the Embedded Configurable Operating System. +// Copyright (C) 2012 Free Software Foundation, Inc. +// +// eCos is free software; you can redistribute it and/or modify it under +// the terms of the GNU General Public License as published by the Free +// Software Foundation; either version 2 or (at your option) any later +// version. +// +// eCos is distributed in the hope that it will be useful, but WITHOUT +// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or +// FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License +// for more details. +// +// You should have received a copy of the GNU General Public License +// along with eCos; if not, write to the Free Software Foundation, Inc., +// 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. +// +// As a special exception, if other files instantiate templates or use +// macros or inline functions from this file, or you compile this file +// and link it with other works to produce a work based on this file, +// this file does not by itself cause the resulting work to be covered by +// the GNU General Public License. However the source code for this file +// must still be made available in accordance with section (3) of the GNU +// General Public License v2. +// +// This exception does not invalidate any other reasons why a work based +// on this file might be covered by the GNU General Public License. +// ------------------------------------------- +// ####ECOSGPLCOPYRIGHTEND#### +//========================================================================== +//#####DESCRIPTIONBEGIN#### +// +// Author(s): ilijak +// Original code: nickg@calivar.com +// Contributors: +// Date: 2012-12-18 +// Description: Simple FPU test with a mix of threads that do and do not use +// floating point. +// This is a modification of the original FPU test, and presents +// both "integer" and "float" threads in order to enforce LAZY +// context switching. Single precision floating point is used. +// This is not very sophisticated as far +// as checking FPU performance or accuracy. It is more +// concerned with checking that several threads doing FP +// operations do not interfere with eachother's use of the +// FPU. +// +//####DESCRIPTIONEND#### +//========================================================================== + +#include <pkgconf/kernel.h> +#include <pkgconf/hal.h> + +#include <cyg/hal/hal_arch.h> + +#include <cyg/kernel/kapi.h> + +#include <cyg/infra/testcase.h> +#include <cyg/infra/diag.h> + +#include <cyg/kernel/test/stackmon.h> +#include CYGHWR_MEMORY_LAYOUT_H + +//========================================================================== + +#if defined(CYGFUN_KERNEL_API_C) && \ + defined(CYGSEM_KERNEL_SCHED_MLQUEUE) && \ + (CYGNUM_KERNEL_SCHED_PRIORITIES > 12) && \ + (CYGMEM_REGION_ram_SIZE >= (49152-4096)) && \ + (!defined(CYGHWR_HAL_CORTEXM_FPU_SWITCH_NONE)) + + +//========================================================================== +// Base priority for all threads. + +#define BASE_PRI 5 + +//========================================================================== +// Runtime +// +// This is the number of ticks that the program will run for. 3000 +// ticks is equal to 30 seconds in the default configuration. For +// simulators we reduce the run time to 3 simulated seconds. + +#define RUN_TICKS 3000 +#define RUN_TICKS_SIM 300 + +//========================================================================== +// Thread parameters + +#define STACK_SIZE CYGNUM_HAL_STACK_SIZE_TYPICAL +#define THREADS_N 2 +static cyg_uint8 stacks[THREADS_N][STACK_SIZE] CYGBLD_ATTRIB_ALIGN_MAX; +static cyg_handle_t thread[THREADS_N]; +static cyg_thread thread_struct[THREADS_N]; +static unsigned long iter_n[THREADS_N]; + +//========================================================================== +// Alarm parameters. + +static cyg_alarm alarm_struct; +static cyg_handle_t alarm; + +volatile static cyg_count32 alarm_ticks = 0; +static cyg_count32 run_ticks = RUN_TICKS; + +//========================================================================== + +static int errors = 0; + +//========================================================================== +// Random number generator. Ripped out of the C library. + +static int CYGBLD_ATTRIB_NO_INLINE rand( unsigned int *seed ) +{ +// This is the code supplied in Knuth Vol 2 section 3.6 p.185 bottom + +#define RAND_MAX 0x7fffffff +#define MM 2147483647 // a Mersenne prime +#define AA 48271 // this does well in the spectral test +#define QQ 44488 // (long)(MM/AA) +#define RR 3399 // MM % AA; it is important that RR<QQ + + *seed = AA*(*seed % QQ) - RR*(unsigned int)(*seed/QQ); + if (*seed < 0) + *seed += MM; + + return (int)( *seed & RAND_MAX ); +} + + +//========================================================================== +// Test calculation. +// Generates an array of random FP values and then thread repeatedly +// alternates following cases: +// - Floatin point calculation occur before it yields; +// - No floating point instructions occur before it yields; +// In addition the re-calculation is called either directly or by means of +// proxy in order to enforce context stacking in different stack locations +// Floating point calculation checks whether the same result is reached +// on every repetition. +// + + +#define V(__i) (values[(__i)%count]) +#define CALC ((V(i-1)*V(i+1))*(V(i-2)*V(i+2))*(V(i-3)*sum)) + +volatile unsigned int fint; + +void CYGBLD_ATTRIB_NO_INLINE +recalc (float *values, float last_sum, cyg_uint32 iter, int thread_i, + int count, const char *name, unsigned int loop_i) +{ + volatile float sum = 1.0; + unsigned int i; + + iter++; + if(thread_i < 3) + iter_n[thread_i]++; + sum = 1.0; + for( i = 0; i < count; i++ ) + sum += CALC; + + if( sum != last_sum ) + { + union float_int_union { + float d; + cyg_uint32 i; + } diu_sum, diu_lastsum; + + diu_sum.d = sum; + diu_lastsum.d = last_sum; + + errors++; + if (sizeof(float) != sizeof(cyg_uint32)) { + diag_printf("Warning: sizeof(float) != sizeof(cyg_uint32), therefore next line may\n" + "have invalid sum/last_sum values\n"); + } + diag_printf("%s: Sum mismatch! %d sum=[%08x] last_sum=[%08x]\n", + name, loop_i, diu_sum.i, diu_lastsum.i ); + + } +} + +void CYGBLD_ATTRIB_NO_INLINE +recalc_proxy (float *values, float last_sum, cyg_uint32 iter, int thread_i, + int count, const char *name, unsigned int loop_i) +{ + int ticks; + + ticks = alarm_ticks; + while(ticks == alarm_ticks) + CYG_EMPTY_STATEMENT; + cyg_thread_yield(); + ticks = alarm_ticks; + while(ticks == alarm_ticks); + CYG_EMPTY_STATEMENT; + cyg_thread_yield(); + recalc(values, last_sum, iter, thread_i, count, name, loop_i); + cyg_thread_yield(); + ticks = alarm_ticks; + while(ticks == alarm_ticks); + CYG_EMPTY_STATEMENT; +} + + +static void do_test1( float *values, + int count, + int loops, + int test, + const char *name) +{ + unsigned int i, j; + // volatiles necessary to force + // values to 32 bits for comparison + volatile float sum = 1.0; + volatile float last_sum; + unsigned int seed; + + cyg_uint32 iter = 0; + int thread_i = name[__builtin_strlen(name)-1]-'1'; + + seed = ((unsigned int)&i)*count; + + // Set up an array of values... + for( i = 0; i < count; i++ ) + values[i] = (float)rand( &seed )/(float)0x7fffffff; + + // Now calculate something from them... + for( i = 0; i < count; i++ ) + sum += CALC; + + last_sum = sum; + + // Now recalculate the sum in a loop and look for errors + for( j = 0; j < loops ; j++ ) + { + if(j % 2) { + recalc_proxy(values, last_sum, iter, thread_i, count, name, j); + } else { + recalc(values, last_sum, iter, thread_i, count, name, j); + } + } +#if 0 + if(thread_i < 3) + diag_printf("INFO:<%s [%d]: %2d calculations done>\n",name,thread_i,j*count); +#endif + if(thread_i < 3) + iter_n[thread_i] = iter; +} + + +//========================================================================== +// Alarm handler +// +// This is called every tick. It lowers the priority of the currently +// running thread and raises the priority of the next. Thus we +// implement a form of timelslicing between the threads at one tick +// granularity. + +static void alarm_fn(cyg_handle_t alarm, cyg_addrword_t data) +{ + alarm_ticks++; + unsigned long iter_sum; + + if( alarm_ticks >= run_ticks ) + { + if( errors ) + CYG_TEST_FAIL("Errors detected"); + else + CYG_TEST_PASS("OK"); + + iter_sum = + iter_n[0] + iter_n[1]; + diag_printf("Iterations = %lu+%lu", + iter_n[0], iter_n[1]); + diag_printf("=%lu\n", iter_sum); + CYG_TEST_FINISH("FP Test done"); + } +} + + +//========================================================================== +// Floating point threads +// +#define FP1_COUNT 10 + +static float fpt1_values[FP1_COUNT]; + +void fptest1( CYG_ADDRWORD id ) +{ + while(1) + do_test1( fpt1_values, FP1_COUNT, 2000000000, id, "fptest1" ); +} + +//========================================================================== +#define FP2_COUNT 11 + +static float fpt2_values[FP2_COUNT]; + +void fptest2( CYG_ADDRWORD id ) +{ + while(1) + do_test1( fpt2_values, FP2_COUNT, 2000000000, id, "fptest2" ); +} + +//====================================================================================== +// Main + + +void fptest_main( void ) +{ + CYG_TEST_INIT(); + + if( cyg_test_is_simulator ) + { + run_ticks = RUN_TICKS_SIM; + } + CYG_TEST_INFO("Run fptest in cyg_start"); + CYG_TEST_INFO( "cyg_start run done"); + + cyg_thread_create( BASE_PRI, + fptest1, + 0, + "thread1", + &stacks[0][0], + STACK_SIZE, + &thread[0], + &thread_struct[0]); + cyg_thread_resume( thread[0] ); + + cyg_thread_create( BASE_PRI, + fptest2, + 1, + "thread2", + &stacks[1][0], + STACK_SIZE, + &thread[1], + &thread_struct[1]); + + cyg_thread_resume( thread[1] ); + + cyg_alarm_create( cyg_real_time_clock(), + alarm_fn, + 0, + &alarm, + &alarm_struct ); + + cyg_alarm_initialize( alarm, cyg_current_time()+1, 1 ); + + cyg_scheduler_start(); + +} + +//========================================================================== + +#ifdef CYGSEM_HAL_STOP_CONSTRUCTORS_ON_FLAG +externC void +cyg_hal_invoke_constructors(); +#endif + +externC void +cyg_start( void ) +{ +#ifdef CYGSEM_HAL_STOP_CONSTRUCTORS_ON_FLAG + cyg_hal_invoke_constructors(); +#endif + fptest_main(); +} + +//========================================================================== + +#else // CYGFUN_KERNEL_API_C... + +externC void +cyg_start( void ) +{ + CYG_TEST_INIT(); + CYG_TEST_INFO("FP test 2 requires:\n" + "CYGFUN_KERNEL_API_C && \n" + "CYGSEM_KERNEL_SCHED_MLQUEUE && \n" + "(CYGNUM_KERNEL_SCHED_PRIORITIES > 12) &&\n" + "(CYGMEM_REGION_ram_SIZE >= (49152-4096))\n" + "(!defined(CYGHWR_HAL_CORTEXM_FPU_SWITCH_NONE))\n" + ); + CYG_TEST_NA("FP test requirements"); +} + +#endif // CYGFUN_KERNEL_API_C, etc. + +//========================================================================== +// EOF fpinttestf2.c
