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