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

#include <cstddef>
#include <iostream>
#include <vector>
using namespace std;

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










double sum(
        vector<float> & x ) {
        double sum = 0;
        for ( size_t i = 0;
            i < BND( x.size() ); i++ )
        sum += @x.at(i)@;
        return sum;
}






#ifdef RUNIT

#define EXPSZ 7





int main() {
        vector<float> x( EXPSZ );
        for ( size_t i = 0; i < EXPSZ; i++ ) x.at(i) = 0.1 * (i + 1);
        for ( size_t i = 0; i < EXPSZ; i++ ) cout << "elm " << i << " " <<  x[i] << endl;
        double sum_ret = sum( x );
        cout << "sum   " << sum_ret << endl;
}



#endif