Gaussian(const result_type& mu_=result_type::Zero(), const matrix_type& sigma_=matrix_type::Identity(), const std::random_device::result_type seed=std::random_device()(), const int size_=n) : base_type(seed, size_), mu(mu_), A(sigma_.llt().matrixU()) { }