source: doc/theses/mike_brooks_MMath/programs/ar-bchk/control.c

#include <stddef.h>
#include <stdio.h>



#ifdef UNSOUND_BOUND
        #define BND( correct_limit ) 100
#else
        #define BND( correct_limit ) correct_limit
#endif










double sum(
        size_t n, float x[] ) {
        double sum = 0;
        for ( size_t i = 0;
            i < BND( n ); i++ )
        sum += @x[i]@;
        return sum;
}






#ifdef RUNIT

#define EXPSZ 7





int main() {
        float x[EXPSZ];
        for ( size_t i = 0; i < EXPSZ; i++ ) x[i] = 0.1 * (i + 1);
        for ( size_t i = 0; i < EXPSZ; i++ ) printf("elm %zd %g\n", i, x[i]);
        double sum_ret = sum( EXPSZ, x );
        printf( "sum   %2g\n", sum_ret );
}



#endif