def opd(p, xy, ray_pkg, fld, wvl, foc):
if ray_pkg[mc.ray] is not None:
fod = opt_model.optical_spec.parax_data.fod
opd = wave_abr_full_calc(fod, fld, wvl, foc, ray_pkg,
fld.chief_ray, fld.ref_sphere)
return convert_to_waves*opd
else:
return None
def wave(p, wi, ray_pkg, fld, wvl, foc):
x = p[0]
y = p[1]
if ray_pkg is not None:
fod = opt_model.optical_spec.parax_data.fod
opd = wave_abr_full_calc(fod, fld, wvl, foc, ray_pkg,
fld.chief_ray, fld.ref_sphere)
opd = convert_to_waves*opd
else:
opd = 0.0
return np.array([x, y, opd])
def wave(p, ray_pkg, fld, wvl, foc):
x = p[0]
y = p[1]
if ray_pkg is not None:
fod = self.opt_model.optical_spec.parax_data.fod
opd = analyses.wave_abr_full_calc(fod, fld, wvl, foc, ray_pkg,
fld.chief_ray,
fld.ref_sphere)
opd = opd/self.opt_model.nm_to_sys_units(wvl)
else:
opd = 0.0
return np.array([x, y, opd])