def _setup_args(info): """Helper to test_as_meg_type_evoked.""" coils = _create_coils(info["chs"], FIFF.FWD_COIL_ACCURACY_NORMAL, info["dev_head_t"], "meg") my_origin, int_rad, noise, lut_fun, n_fact = _setup_dots("fast", coils, "meg") args_dict = dict( intrad=int_rad, volume=False, coils1=coils, r0=my_origin, ch_type="meg", lut=lut_fun, n_fact=n_fact ) return args_dict
def _setup_args(info): """Helper to test_as_meg_type_evoked.""" coils = _create_coils(info['chs'], FIFF.FWD_COIL_ACCURACY_NORMAL, info['dev_head_t'], 'meg') my_origin, int_rad, noise, lut_fun, n_fact = _setup_dots( 'fast', coils, 'meg') args_dict = dict(intrad=int_rad, volume=False, coils1=coils, r0=my_origin, ch_type='meg', lut=lut_fun, n_fact=n_fact) return args_dict
def test_magnetic_dipole(): """Test basic magnetic dipole forward calculation """ trans = Transform('mri', 'head', np.eye(4)) info = read_info(fname_raw) picks = pick_types(info, meg=True, eeg=False, exclude=[]) info = pick_info(info, picks[:12]) coils = _create_coils(info['chs'], FIFF.FWD_COIL_ACCURACY_NORMAL, trans) # magnetic dipole at device origin r0 = np.array([0., 13., -6.]) for ch, coil in zip(info['chs'], coils): rr = (ch['coil_trans'][:3, 3] + r0) / 2. far_fwd = _magnetic_dipole_field_vec(r0[np.newaxis, :], [coil]) near_fwd = _magnetic_dipole_field_vec(rr[np.newaxis, :], [coil]) ratio = 8. if ch['ch_name'][-1] == '1' else 16. # grad vs mag assert_allclose(np.median(near_fwd / far_fwd), ratio, atol=1e-1)
def test_magnetic_dipole(): """Test basic magnetic dipole forward calculation """ trans = {'to': FIFF.FIFFV_COORD_HEAD, 'from': FIFF.FIFFV_COORD_MRI, 'trans': np.eye(4)} info = read_info(fname_raw) picks = pick_types(info, meg=True, eeg=False, exclude=[]) info = pick_info(info, picks[:12]) coils = _create_coils(info['chs'], FIFF.FWD_COIL_ACCURACY_NORMAL, trans) # magnetic dipole at device origin r0 = np.array([0., 13., -6.]) for ch, coil in zip(info['chs'], coils): rr = (ch['coil_trans'][:3, 3] + r0) / 2. far_fwd = _magnetic_dipole_field_vec(r0[np.newaxis, :], [coil]) near_fwd = _magnetic_dipole_field_vec(rr[np.newaxis, :], [coil]) ratio = 8. if ch['ch_name'][-1] == '1' else 16. # grad vs mag assert_allclose(np.median(near_fwd / far_fwd), ratio, atol=1e-1)