def real_xspec_functions(k, Te, F, alpha=np.pi / 2.):
"""
Calculate analytical expressions for the real component of admittance
and coherence functions.
:type k: np.ndarray
:param k: Wavenumbers (rad/m)
:type Te: float
:param Te: Effective elastic thickness (km)
:type F: float
:param F: Subsurface-to-surface load ratio [0, 1[
:type alpha: float, optional
:param alpha: Phase difference between initial applied loads (rad)
:return:
(tuple): tuple containing:
* admittance (:class:`~numpy.ndarray`): Real admittance function (shape: ``len(k)``)
* coherence (:class:`~numpy.ndarray`): Coherence functions (shape: ``len(k)``)
"""
admittance, coherence = flex.real_xspec_functions(k, Te, F, alpha)
return admittance, coherence
def pred_coh(k, Te, F, alpha):
return flex.real_xspec_functions(k, Te, F, alpha)[1]
def pred_joint(k, Te, F, alpha):
admittance, coherence = flex.real_xspec_functions(k, Te, F, alpha)
return np.array([admittance, coherence]).flatten()
def pred_admit(k, Te, F, alpha):
return flex.real_xspec_functions(k, Te, F, alpha)[0]