def _forward_prop_random_thru_post_mm(self, mx, vx):
"""Summary
Args:
mx (TYPE): Description
vx (TYPE): Description
Returns:
TYPE: Description
"""
Kuuinv = self.Kuuinv
A = np.einsum('ab,db->da', Kuuinv, self.mu)
Smm = self.Su + np.einsum('da,db->dab', self.mu, self.mu)
B = np.einsum(
'ab,dbc->dac',
Kuuinv, np.einsum('dab,bc->dac', Smm, Kuuinv)) - Kuuinv
psi0 = np.exp(2.0 * self.sf)
psi1, psi2 = compute_psi_weave(
2 * self.ls, 2 * self.sf, mx, vx, self.zu)
mout = np.einsum('nm,dm->nd', psi1, A)
Bpsi2 = np.einsum('dab,nab->nd', B, psi2)
vout = psi0 + Bpsi2 - mout**2
return mout, vout
评论列表
文章目录