FormFac(u, quark_prop_1, quark_prop_2, t_source, t_sink, j_decay);
FormFac(u, quark_prop_1, quark_prop_2, t_source, t_sink);
int main(int argc, char **argv){  // Setup the geometry  const int foo[Nd] = {LX0,LX1};  multi1d<int> nrow(Nd);  nrow = foo;  geom.Init(nrow);  //! Test out propagators  multi1d<LatticeGauge> u(Nd);  for(int m=0; m < u.size(); ++m)    gaussian(u[m]);  WRITE_NAMELIST(cerr,u);  LatticePropagator quark_prop_1, quark_prop_2;  gaussian(quark_prop_1);  gaussian(quark_prop_2);  int j_decay = Nd-1;  int length = geom.LattSize()[j_decay];  multi1d<int> t_source(Nd);  t_source = 0;  int t_sink = 3;  FormFac(u, quark_prop_1, quark_prop_2, t_source, t_sink, j_decay);}