H_handle->refresh(gauge_state);
H_handle->refreshInternalFields(gauge_state);
int main(int argc, char *argv[]) {  // Initialise QDP  QDP_initialize(&argc, &argv);  // Setup a small lattice  const int nrow_arr[] = {2, 2, 2, 2};  multi1d<int> nrow(Nd);  nrow=nrow_arr;  Layout::setLattSize(nrow);  Layout::create();  multi1d<LatticeColorMatrix> u(Nd);  {    XMLReader file_xml;    XMLReader config_xml;    Cfg_t foo; foo.cfg_type=CFG_TYPE_UNIT;;    // Cfg_t foo; foo.cfg_type=CFG_TYPE_SZIN; foo.cfg_file="./CFGIN";    gaugeStartup(file_xml, config_xml, u, foo);  }  // Dump output  XMLFileWriter xml_out("./XMLDAT");  push(xml_out, "t_gauge_ferm_monomials");  // Read Parameters  multi1d<int> boundary(Nd);           // Ferm BC's  std::string monomial_name;           // String for Factory  XMLReader param_in("DATA");  // Snarf it all  XMLReader paramtop(param_in, "/HamiltonianTest");      Handle< ExactLatColMatHamiltonian > H_handle;  try {     read(paramtop, "./Hamiltonian", H_handle);  }  catch( const std::string& e) {     QDPIO::cerr << "Error Reading Hamiltonian " << e <<  endl;    QDP_abort(1);  }  // Fictitious momenta for now  multi1d<LatticeColorMatrix> p(Nd);  for(int mu=0; mu<Nd; mu++) {     gaussian(p[mu]);  }  // Create a field state  GaugeFieldState gauge_state(p,u);  // Refresh Pseudofermions  H_handle->refresh(gauge_state);  // Compute Force from Monomial  multi1d<LatticeColorMatrix> dsdq(Nd);  H_handle->dsdq(dsdq, gauge_state);  // Compute action from monomial  Double KE, PE;  H_handle->mesE(gauge_state, KE, PE);  QDPIO::cout << "KE = " << KE << "  PE = " << PE << endl;  pop(xml_out);  xml_out.close();    // Finish  QDP_finalize();  exit(0);}