@property<EOL><INDENT>def branchScale(self):<DEDENT>
bs = -(self.Phi_x * scipy.diagonal(self.Pxy[<NUM_LIT:0>])).sum() * self.mu<EOL>assert bs > <NUM_LIT:0><EOL>return bs<EOL>
See docs for `Model` abstract base class.