code: def getModelPosterior(self,min):
          
          Sigma = self.getLaplaceCovar(min)
          n_params = self.vd.getNumberScales()
          ModCompl = 0.5*n_params*SP.log(2*SP.pi)+0.5*SP.log(SP.linalg.det(Sigma))
          RV = min['LML']+ModCompl
          return RV
docs: USES LAPLACE APPROXIMATION TO CALCULATE THE BAYESIAN MODEL POSTERIOR