#include <iostream>
#include <complex>
#include <blitz/array.h>

struct Mean { 
  double operator() ( double i1, double i2 ) { return (i1+i2)/2.0; }
  BZ_DECLARE_FUNCTOR2( Mean );
};

int main(int argc, char** argv) 
{
    const int N = 100;
    blitz::Array<double,1> A(N);
    Mean mean;

    A = blitz::tensor::i;
    A = mean( A, A );
    std::cout << A << std::endl;
    return 0;
}