thermoDynamics/upDateTemperature [ Subroutines ]

upDateTemperature/evaluatePSquared [ Subroutines ]

Calculation of pSquared(x,y,z).
This is stored in the field psquared(-gridPoints(1:3):gridPoints(1:3)).
Usage of the smearing weigths to smear particles' momentum over
gridpoints in their neighborhood.

Evaluates temperature as a function of rho(0) and <p**2>. At first call it
initializes a field "temSave" which holds this information. Therefore no
further calculation is necessary after initializing this field. It is
generated in the subroutine "initTempe".

thermoDynamics/integral [ Functions ]

Evaluates for a gas of degenerate neutrons and protons of given "temperature" and
given chemical potential "mu" the following:
Switch=1 : rho=<1>=4*Integral(1/1+exp((E(p)-mu)/T) (dp)**3)/(2pi)**3 over the full p-space.
Switch=2 : <p**2>=4*Integral(p**2/1+exp((E(p)-mu)/T) (dp)**3)/(2pi)**3 over the full p-space.
The first Factor of 4 is due to spin&isospin degeneracy.