/Users/mdipierro/fermiqcd/development/Libraries/fermiqcd_check_differences.h

Go to the documentation of this file.

// compares two fields and returns the max distrance between 
// equivalent components.
float check_differences(mdp_field<mdp_complex> &chi,
                        mdp_field<mdp_complex> &psi) {
  begin_function("check_differences");

  mdp_int i;
  mdp_int i_min=psi.physical_local_start();
  mdp_int i_max=psi.physical_local_stop();
  float max=0, tmp;
  if(&chi.lattice()!=&psi.lattice()) 
    error("check_differences()\nFields defined on different lattices");
  for(i=i_min; i<i_max; i++) {
    tmp=abs(chi[i]-psi[i]);
    if(tmp>max) max=tmp;
  }
  mdp << "Fields agree/disagree within precision=" << max << '\n';

  end_function("check_differences");
  return max;
}

---