function pbeDefinition = constGrowth_normal(seed)

%% call common definition
pbeDefinition = constGrowth_common;

%% Initial distribution
rng(seed);
% discrete sampling
% pbeDefinition.seeding.nSample = 1e2;
pbeDefinition.seeding.disD = obtainDistributionDiscrete('normal', ...
    'nSample', pbeDefinition.seeding.nSample, ...
    'mean', pbeDefinition.seeding.mean, ...
    'cov', pbeDefinition.seeding.cov);

% printDistributionStatistics(pbeDefinition)

%% testing vector vs. nonvector calls
% thisHandle = @(h)(computeMeasureDec(pbeDefinition.dec, 'hC', h, 'volume', ...
%     [], 0, -1));
% % % thisHandle = @(h)(computePartitionIndexDec(pbeDefinition.dec, h, [], 0));
% % 
% vectorSum = sum(thisHandle(pbeDefinition.seeding.disD.positionMatrix))
% forSum = 0;
% for iSample = 1:pbeDefinition.seeding.nSample
%     forSum = forSum + thisHandle(pbeDefinition.seeding.disD.positionMatrix(:, iSample));
% end
% forSum
