Function connect_phase_from_phase

src/boundaries.cpp:362–366

  return fabs(thephase.imag()) < 1e-13 && fabs(thephase.real() - realphase) < 1e-13;
}
static connect_phase connect_phase_from_phase(std::complex<double> thephase) {
  return phase_isclose(thephase, 1.0)    ? CONNECT_COPY
         : phase_isclose(thephase, -1.0) ? CONNECT_NEGATE
                                         : CONNECT_PHASE;
}

void fields::connect_the_chunks() {
  /* For some of the chunks, H==B, and we definitely don't need to