public void Calculate(double sec)
{
var matrix = BurnUp.Matrix;
var density = DensityArray.InitialDensity;
DensityArray.Density = MatExp.Calculate(matrix * sec, density);
DensityArray.Normolize();
}
public Reactor(List <Isotope> isotopes, NeutronSpectra spectra, DensityArray densityArray, Endf[] nuclearData)
{
initialize(nuclearData);
NeutronSpectra = spectra;
BurnUp = new BurnUp(isotopes, NeutronSpectra);
DensityArray = densityArray;
Calculate(1.0E6);
}
public void SetDensityArray(List <NuclideDensity> initDensities)
{
var densities = new List <NuclideDensity>();
foreach (var iso in BurnUp.Isotopes)
{
var dens = initDensities.FirstOrDefault(x => x.NuclideName == iso.Name);
var weight = dens == null ? 0.0 : dens.Density;
densities.Add(new NuclideDensity(iso, weight));
}
DensityArray = new DensityArray(densities);
}
