#include <iostream>

#include <blitz/array.h>
#include <blitz/tinyvec.h>
#include <blitz/tinyvec-et.h>

#include "ternary.h"

using namespace blitz;

float gauss(float val, float mean, float var) {
    return std::exp( -(val - mean)*(val - mean)/(2*var));
}

BZ_DECLARE_FUNCTION3_RET(gauss,double)

int main(int argc, char** argv) {
    Array<float,1> a (10);
    a = tensor::i;

    Array<int,1> m (10);
    m = tensor::i -1;

    Array<short,1> v(10);
    v = tensor::i +1;

    std::cout << a << std::endl;
    std::cout << m << std::endl;
    std::cout << v << std::endl;

    Array<double,1> r(10);

    r = gauss(a+2,m,v);
    std::cout << r << std::endl; 

    return 0;
}