private void UpdateOutputParameters(
QQState.QQState state
)
{
Energy_MeV = state.Param.Energy_MeV;
GammaDamp_MeV = state.Param.GammaDamp_MeV;
SoftScale_MeV = state.Param.SoftScale_MeV;
AlphaHard = state.AlphaHard;
AlphaSoft = state.AlphaSoft;
AlphaThermal = state.AlphaThermal;
DebyeMass_MeV = state.DebyeMass_MeV;
Sigma_MeV2 = state.Param.Sigma_MeV;
SigmaEff_MeV2 = state.SigmaEff_MeV;
StepSize = state.StepSize_fm;
StepsPerPeriod = state.WaveVector_fm * state.StepSize_fm;
Tchem_MeV = state.Param.Tchem_MeV;
Tcrit_MeV = state.Param.Tcrit_MeV;
Trials = state.Trials;
WaveVector_per_fm = state.WaveVector_fm;
if (state is QQBoundState)
{
QQBoundState boundState = state as QQBoundState;
NumberExtrema = boundState.NumberExtrema;
BoundMass_MeV = boundState.BoundMass_MeV;
AvgInvDisplacement_per_fm = boundState.RadiusExpectationValue(-1);
DisplacementRMS_fm = Math.Sqrt(boundState.RadiusExpectationValue(2));
AlphaUltraSoft = boundState.AlphaUltraSoft;
UltraSoftScale_MeV = boundState.UltraSoftScale_MeV;
}
}
/********************************************************************************************
* Private/protected static members, functions and properties
********************************************************************************************/
private static List <string> CreateWaveFunctionDataList(
QQState.QQState state
)
{
List <string> dataList = new List <string>
{
string.Format("{0,-10}{1,-18}{2,-18}", "#r (fm)", "Re(Psi) (fm^-1/2)", "Im(Psi) (fm^-1/2)")
};
for (int j = 0; j <= state.Param.StepNumber; j++)
{
dataList.Add(string.Format("{0,-10}{1,-18}{2,-18}",
state.Radius_fm[j].ToString("G10"),
state.WaveFunction_fm[j].Re.ToString("G10"),
state.WaveFunction_fm[j].Im.ToString("G10")));
}
return(dataList);
}
