self.q_mu = Param(np.zeros((self.num_data, self.num_latent)))
q_sqrt = np.array([np.eye(self.num_data)
for _ in range(self.num_latent)]).swapaxes(0, 2)
self.q_sqrt = Param(q_sqrt, transforms.LowerTriangular(self.num_data, self.num_latent))
def _compile(self, optimizer=None):
Before calling the standard compile function, check to see if the size