Example #1
0
        private void AdvanceFields()
        {
            if (Param.ExpansionMode == ExpansionMode.Transverse)
            {
                Solver.SolveUntil(CurrentTime);
                VelocityX.SetValues(Solver.VX);
                VelocityY.SetValues(Solver.VY);
                Temperature.Advance(Solver);
            }
            else
            {
                Temperature.Advance(CurrentTime);
            }

            if (Param.UseElectricField)
            {
                ElectricField.Advance(CurrentTime, QGPConductivity);
            }
            if (Param.UseMagneticField)
            {
                MagneticField.Advance(CurrentTime, QGPConductivity);
            }

            DecayWidth.Advance(CurrentTime);
            DampingFactor.SetValues((x, y, pT, state) => DampingFactor[x, y, pT, state]
                                    * Math.Exp(-TimeStep * DecayWidth[x, y, pT, state] / Constants.HbarC_MeV_fm));
        }