Compute multipolar moments for point-like mags (in fine cal).
(exp, trans, cal, ignore_ref=False, mag_scale=100.0)
| 2330 | |
| 2331 | |
| 2332 | def _sss_basis_point(exp, trans, cal, ignore_ref=False, mag_scale=100.0): |
| 2333 | """Compute multipolar moments for point-like mags (in fine cal).""" |
| 2334 | # Loop over all coordinate directions desired and create point mags |
| 2335 | S_tot = 0.0 |
| 2336 | # These are magnetometers, so use a uniform coil_scale of 100. |
| 2337 | this_cs = np.array([mag_scale], float) |
| 2338 | for imb, coils in zip(cal["grad_imbalances"], cal["grad_coilsets"]): |
| 2339 | S_add = _trans_sss_basis(exp, coils, trans, this_cs) |
| 2340 | # Scale spaces by gradiometer imbalance |
| 2341 | S_add *= imb[:, np.newaxis] |
| 2342 | S_tot += S_add |
| 2343 | |
| 2344 | # Return point-like mag bases |
| 2345 | return S_tot |
| 2346 | |
| 2347 | |
| 2348 | def _regularize_out(int_order, ext_order, mag_or_fine, extended_remove): |
no test coverage detected