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

using namespace blitz;

template<typename T> void sink(T&) { }

int main()
{
   TinyMatrix<complex<double>,4,4> gamma1, gamma3;
   TinyVector<complex<double>,4> Psi;

   complex<double> i(0,1);

   gamma1 = 0, i, 0, 0,
            i, 0, 0, 0,
            0, 0, 0, -i,
            0, 0, -i, 0;

   gamma3 = 0, 1, 0, 0,
            -1, 0, 0, 0,
            0, 0, 0, 1,
            0, 0, -1, 0;

   sink(Psi);

   TinyVector<complex<double>,4> tmp, result;
   tmp = product(gamma3,Psi);
   result = product(gamma1,tmp);

   sink(result);
   return 0;
}