@property<EOL><INDENT>def precision(self):<DEDENT>
if self._precision is None:<EOL><INDENT>cov = np.atleast_3d(self.covariance)<EOL>self._precision = np.zeros(cov.shape)  <EOL>for p in range(cov.shape[-<NUM_LIT:1>]):<EOL><INDENT>self._precision[:, :, p] = pdinv(cov[:, :, p])[<NUM_LIT:0>]<EOL><DEDENT><DEDENT>return self._precision<EOL>
Inverse of posterior covariance