Function FOR_MAGNETIC_COMPONENTS

src/monitor.cpp:218–223

  complex<double> tr(0.0, 0.0);
  int nc = 0;
  FOR_MAGNETIC_COMPONENTS(c) {
    if (gv.has_field(c)) {
      tr += get_chi1inv(c, component_direction(c), loc, frequency, false);
      ++nc;
    }
  }
  return complex<double>(nc, 0) / sum_to_all(tr);
}