// ----------------------------------------------------------------- void setup_bc() { register int i, dir; register site *s; for (dir = 0; dir < NUMLINK; dir++) { FORALLSITES(i, s) { s->bc[dir] = 1.0; s->bc[OPP_LDIR(dir)] = 1.0; if (s->t + offset[dir][TUP] < 0) s->bc[dir] = PBC; else if (s->t + offset[dir][TUP] > nt - 1) s->bc[dir] = PBC; if (s->t - offset[dir][TUP] < 0) s->bc[OPP_LDIR(dir)] = PBC; else if (s->t - offset[dir][TUP] > nt - 1) s->bc[OPP_LDIR(dir)] = PBC; } }
// ----------------------------------------------------------------- void setup_bc() { register int i, dir; register site *s; // Single-offset terms only FORALLDIR(dir) { FORALLSITES(i, s) { s->bc[dir] = 1.0; s->bc[OPP_LDIR(dir)] = 1.0; if (s->t + offset[dir][TUP] < 0) s->bc[dir] = PBC; else if (s->t + offset[dir][TUP] > nt - 1) s->bc[dir] = PBC; if (s->t - offset[dir][TUP] < 0) s->bc[OPP_LDIR(dir)] = PBC; else if (s->t - offset[dir][TUP] > nt - 1) s->bc[OPP_LDIR(dir)] = PBC; } }