function pbeDef = paracetamolGrowth_4u_volume_scaledReduced()
% automatically generated by writePbeDefinitionForProjection.m

dataApprox = load('/data/havena/elk-dev/code-matlab/cs-dynamic-approximation/database/paracetamolGrowth_4u_volume_scaledReduced.matx', '-mat');
dataApprox = dataApprox.dataApprox;
pbeDef = paracetamolGrowth;
pbeDef.dec = dataApprox.dec;

%% initial condition
heMatrix = mapToExtendedEmbedding(pbeDef.seeding.positionMatrix, ...
    dataApprox.dec.properData);
pbeDef.seeding.positionMatrix = mapFromProperEmbedding(projectVectorDec(...
    heMatrix, pbeDef.dec), pbeDef.dec);

%% growth rate mapping
pbeDef.growthRate = @(H, P, y, t)(dataApprox.mappingMatrixOtoC * ...
    pbeDef.growthRate(ones(4, size(H,2)), P,y,t));
  
pbeDef.pivotProperty = @(H, up)([up; ...
    computeMeasureDec(pbeDef.dec, 'hC', H, 'volume')]);

pbeDef.optionStruct.odeSolver = @ode45;
