RelativisticEuler::Solutions::RotatingStar Class Reference

A solution obtained by reading in rotating neutron star initial data from the RotNS code based on [43] and [44]. More...

#include <RotatingStar.hpp>

Classes
struct	PolytropicConstant
	The polytropic constant of the fluid. More...
struct	RotNsFilename
	The path to the RotNS data file. More...

Public Types
using	options = implementation defined
Public Types inherited from RelativisticEuler::AnalyticSolution< 3 >
using	tags = implementation defined

Public Member Functions
	RotatingStar (const RotatingStar &)
RotatingStar &	operator= (const RotatingStar &)
	RotatingStar (RotatingStar &&)=default
RotatingStar &	operator= (RotatingStar &&)=default
	RotatingStar (std::string rot_ns_filename, std::unique_ptr< EquationsOfState::EquationOfState< true, 1 > > equation_of_state)
	RotatingStar (std::string rot_ns_filename, double polytropic_constant)
auto	get_clone () const -> std::unique_ptr< evolution::initial_data::InitialData > override
template<typename DataType , typename... Tags>
tuples::TaggedTuple< Tags... >	variables (const tnsr::I< DataType, 3 > &x, const double, tmpl::list< Tags... >) const
	Retrieve a collection of variables at (x, t)
void	pup (PUP::er &p) override
const EquationsOfState::EquationOfState< true, 1 > &	equation_of_state () const
double	equatorial_radius () const
virtual auto	get_clone () const -> std::unique_ptr< InitialData >=0
Public Member Functions inherited from hydro::TemperatureInitialization< RotatingStar >
auto	variables (const tnsr::I< DataType, Dim > &x, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >
auto	variables (const tnsr::I< DataType, Dim > &x, const double t, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >
auto	variables (ExtraVars &extra_variables, const tnsr::I< DataType, Dim > &x, Args &... extra_args, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >

Static Public Attributes
static constexpr Options::String	help
Static Public Attributes inherited from RelativisticEuler::AnalyticSolution< 3 >
static constexpr size_t	volume_dim

Protected Types
template<typename DataType >
using	DerivLapse = ::Tags::deriv< gr::Tags::Lapse< DataType >, tmpl::size_t< 3 >, Frame::Inertial >
template<typename DataType >
using	DerivShift = ::Tags::deriv< gr::Tags::Shift< DataType, 3 >, tmpl::size_t< 3 >, Frame::Inertial >
template<typename DataType >
using	DerivSpatialMetric = ::Tags::deriv< gr::Tags::SpatialMetric< DataType, 3 >, tmpl::size_t< 3 >, Frame::Inertial >

Protected Member Functions
template<typename DataType >
void	interpolate_vars_if_necessary (gsl::not_null< IntermediateVariables< DataType > * > vars) const
template<typename DataType >
void	interpolate_deriv_vars_if_necessary (gsl::not_null< IntermediateVariables< DataType > * > vars) const
template<typename DataType >
Scalar< DataType >	lapse (const DataType &gamma, const DataType &rho) const
template<typename DataType >
tnsr::I< DataType, 3, Frame::Inertial >	shift (const DataType &omega, const DataType &phi, const DataType &radius, const DataType &sin_theta) const
template<typename DataType >
tnsr::ii< DataType, 3, Frame::Inertial >	spatial_metric (const DataType &gamma, const DataType &rho, const DataType &alpha, const DataType &phi) const
template<typename DataType >
tnsr::II< DataType, 3, Frame::Inertial >	inverse_spatial_metric (const DataType &gamma, const DataType &rho, const DataType &alpha, const DataType &phi) const
template<typename DataType >
Scalar< DataType >	sqrt_det_spatial_metric (const DataType &gamma, const DataType &rho, const DataType &alpha) const
template<typename DataType >
auto	make_metric_data (size_t num_points) const -> typename IntermediateVariables< DataType >::MetricData
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::RestMassDensity< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::RestMassDensity< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::ElectronFraction< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::ElectronFraction< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::SpecificEnthalpy< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::SpecificEnthalpy< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::Pressure< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Pressure< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::SpecificInternalEnergy< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::SpecificInternalEnergy< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::SpatialVelocity< DataType, 3 > >) const -> tuples::TaggedTuple< hydro::Tags::SpatialVelocity< DataType, 3 > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::LorentzFactor< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::LorentzFactor< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::MagneticField< DataType, 3 > >) const -> tuples::TaggedTuple< hydro::Tags::MagneticField< DataType, 3 > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::DivergenceCleaningField< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::DivergenceCleaningField< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::Lapse< DataType > >) const -> tuples::TaggedTuple< gr::Tags::Lapse< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::Shift< DataType, 3 > >) const -> tuples::TaggedTuple< gr::Tags::Shift< DataType, 3 > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::SpatialMetric< DataType, 3 > >) const -> tuples::TaggedTuple< gr::Tags::SpatialMetric< DataType, 3 > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::SqrtDetSpatialMetric< DataType > >) const -> tuples::TaggedTuple< gr::Tags::SqrtDetSpatialMetric< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::InverseSpatialMetric< DataType, 3 > >) const -> tuples::TaggedTuple< gr::Tags::InverseSpatialMetric< DataType, 3 > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list<::Tags::dt< gr::Tags::Lapse< DataType > > >) const -> tuples::TaggedTuple<::Tags::dt< gr::Tags::Lapse< DataType > > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< DerivLapse< DataType > >) const -> tuples::TaggedTuple< DerivLapse< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list<::Tags::dt< gr::Tags::Shift< DataType, 3 > > >) const -> tuples::TaggedTuple<::Tags::dt< gr::Tags::Shift< DataType, 3 > > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< DerivShift< DataType > >) const -> tuples::TaggedTuple< DerivShift< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list<::Tags::dt< gr::Tags::SpatialMetric< DataType, 3 > > >) const -> tuples::TaggedTuple<::Tags::dt< gr::Tags::SpatialMetric< DataType, 3 > > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< DerivSpatialMetric< DataType > >) const -> tuples::TaggedTuple< DerivSpatialMetric< DataType > >
template<typename DataType >
auto	variables (gsl::not_null< IntermediateVariables< DataType > * > vars, const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::ExtrinsicCurvature< DataType, 3 > >) const -> tuples::TaggedTuple< gr::Tags::ExtrinsicCurvature< DataType, 3 > >

Protected Attributes
std::string	rot_ns_filename_ {}
detail::CstSolution	cst_solution_ {}
double	polytropic_constant_ = std::numeric_limits<double>::signaling_NaN()
double	polytropic_exponent_ = std::numeric_limits<double>::signaling_NaN()
bool	is_polytrope_ {}
std::unique_ptr< EquationsOfState::EquationOfState< true, 1 > >	equation_of_state_

Static Protected Attributes
static constexpr double	atmosphere_floor_ = 1.e-50

Friends
bool	operator== (const RotatingStar &lhs, const RotatingStar &rhs)

Detailed Description

A solution obtained by reading in rotating neutron star initial data from the RotNS code based on [43] and [44].

The code that generates the initial data is part of a private SXS repository called RotNS.

The metric in spherical coordinates is given by [43]

$$\begin{array}{}\text{(1)}& d{s}^2=-e^{\gamma +\rho }d{t}^2+e^{2\alpha }(d{r}^2+r^{2}d{\theta }^2)+e^{\gamma -\rho }{r}^2{\sin }^2(\theta )(d\varphi -\omega dt{)}^2.\end{array}$$

We use rotation about the $z$-axis. That is,

$$\begin{array}{}\text{(2)}& g_{tt}& =-e^{\gamma +\rho }+e^{\gamma -\rho }{r}^2{\sin }^2(\theta ){\omega }^2\text{(3)}& g_{rr}& =e^{2\alpha }\text{(4)}& g_{\theta \theta }& =e^{2\alpha }{r}^2\text{(5)}& g_{\varphi \varphi }& =e^{\gamma -\rho }{r}^2{\sin }^2(\theta )\text{(6)}& g_{t\varphi }& =-e^{\gamma -\rho }{r}^2{\sin }^2(\theta )\omega .\end{array}$$

We can transform from spherical to Cartesian coordinates using

$$\begin{array}{}\text{(7)}& \frac{\partial (r,\theta ,\varphi )}{\partial (x,y,z)}=\left(\begin{array}{ccc}\cos \left(\varphi \right)\sin \left(\theta \right)& \sin \left(\theta \right)\sin \left(\varphi \right)& \cos \left(\theta \right)\\ \frac{\cos \left(\varphi \right)\cos \left(\theta \right)}{r}& \frac{\sin \left(\varphi \right)\cos \left(\theta \right)}{r}& -\frac{\sin \left(\theta \right)}{r}\\ -\frac{\sin \left(\varphi \right)}{r\sin \left(\theta \right)}& \frac{\cos \left(\varphi \right)}{r\sin \left(\theta \right)}& 0\end{array}\right)\end{array}$$

and

$$\begin{array}{}\text{(8)}& \frac{\partial (x,y,z)}{\partial (r,\theta ,\varphi )}=\left(\begin{array}{ccc}\cos \left(\varphi \right)\sin \left(\theta \right)& r\cos \left(\varphi \right)\cos \left(\theta \right)& -r\sin \left(\theta \right)\sin \left(\varphi \right)\\ \sin \left(\theta \right)\sin \left(\varphi \right)& r\sin \left(\varphi \right)\cos \left(\theta \right)& r\cos \left(\varphi \right)\sin \left(\theta \right)\\ \cos \left(\theta \right)& -r\sin \left(\theta \right)& 0\end{array}\right)\end{array}$$

We denote the lapse as $N$ since $\alpha$ is already being used,

$$\begin{array}{}\text{(9)}& N=e^{(\gamma +\rho )/2}.\end{array}$$

The shift is

$$\begin{array}{}\text{(10)}& {\beta }^{\varphi }=& -\omega \end{array}$$

so

$$\begin{array}{}\text{(11)}& {\beta }^x& ={\partial }_{\varphi }x{\beta }^{\varphi }=\sin \left(\varphi \right)\omega r\sin \left(\theta \right),\text{(12)}& {\beta }^y& ={\partial }_{\varphi }y{\beta }^{\varphi }=-\cos \left(\varphi \right)\omega r\sin \left(\theta \right),\text{(13)}& {\beta }^z& =0.\end{array}$$

The spatial metric is

$$\begin{array}{rl}{\gamma }_{xx}& ={\sin }^2(\theta ){\cos }^2(\varphi )e^{(2\alpha )}+{\cos }^2(\theta ){\cos }^2(\varphi )e^{(2\alpha )}+{\sin }^2(\varphi )e^{(\gamma -\rho )}\\ \text{(14)}& & ={\cos }^2(\varphi )e^{2\alpha }+{\sin }^2(\varphi )e^{(\gamma -\rho )}{\gamma }_{xy}& ={\sin }^2(\theta )\cos \left(\varphi \right)\sin \left(\varphi \right)e^{(2\alpha )}+{\cos }^2(\theta )\cos \left(\varphi \right)\sin \left(\varphi \right)e^{(2\alpha )}-\sin \left(\varphi \right)\cos \left(\varphi \right)e^{(\gamma -\rho )}\\ \text{(15)}& & =\cos \left(\varphi \right)\sin \left(\varphi \right)e^{2\alpha }-\sin \left(\varphi \right)\cos \left(\varphi \right)e^{(\gamma -\rho )}{\gamma }_{yy}& ={\sin }^2(\theta ){\sin }^2(\varphi )e^{(2\alpha )}+{\cos }^2(\theta ){\sin }^2(\varphi )e^{(2\alpha )}+{\cos }^2(\varphi )e^{(\gamma -\rho )}\\ \text{(16)}& & ={\sin }^2(\varphi )e^{2\alpha }+{\cos }^2(\varphi )e^{(\gamma -\rho )}\text{(17)}& {\gamma }_{xz}& =\sin \left(\theta \right)\cos \left(\varphi \right)\cos \left(\theta \right)e^{2\alpha }-\cos \left(\theta \right)\cos \left(\varphi \right)\sin \left(\theta \right)e^{2\alpha }=0\text{(18)}& {\gamma }_{yz}& =\sin \left(\theta \right)\sin \left(\varphi \right)\cos \left(\theta \right)e^{2\alpha }-\cos \left(\theta \right)\sin \left(\varphi \right)\sin \left(\theta \right)e^{2\alpha }=0\text{(19)}& {\gamma }_{zz}& ={\cos }^2(\theta )e^{2\alpha }+{\sin }^2(\theta )e^{2\alpha }=e^{2\alpha }\end{array}$$

and its determinant is

$$\begin{array}{}\text{(20)}& \gamma =e^{4\alpha +(\gamma -\rho )}=e^{4\alpha }{e}^{(\gamma -\rho )}.\end{array}$$

At $r=0$ we have $2\alpha =\gamma -\rho$ and so the ${\gamma }_{xx}={\gamma }_{yy}={\gamma }_{zz}=e^{2\alpha }$ and all other components are zero. The inverse spatial metric is given by

$$\begin{array}{rl}{\gamma }^{xx}& =\frac{{\gamma }_{yy}}{e^{2\alpha }{e}^{(\gamma -\rho )}}=[{\sin }^2(\varphi )e^{2\alpha }+{\cos }^2(\varphi )e^{(\gamma -\rho )}]e^{-2\alpha }{e}^{-(\gamma -\rho )}\\ \text{(21)}& & ={\sin }^2(\varphi )e^{-(\gamma -\rho )}+{\cos }^2(\varphi )e^{-2\alpha }{\gamma }^{yy}& =\frac{{\gamma }_{xx}}{e^{2\alpha }{e}^{(\gamma -\rho )}}=[{\cos }^2(\varphi )e^{2\alpha }+{\sin }^2(\varphi )e^{(\gamma -\rho )}]e^{-2\alpha }{e}^{-(\gamma -\rho )}\\ \text{(22)}& & ={\cos }^2(\varphi )e^{-(\gamma -\rho )}+{\sin }^2(\varphi )e^{-2\alpha }{\gamma }^{xy}& =\frac{-{\gamma }_{xy}}{e^{2\alpha }{e}^{(\gamma -\rho )}}=-[\cos \left(\varphi \right)\sin \left(\varphi \right)e^{2\alpha }-\sin \left(\varphi \right)\cos \left(\varphi \right)e^{(\gamma -\rho )}]e^{-2\alpha }{e}^{-(\gamma -\rho )}\\ & =-\cos \left(\varphi \right)\sin \left(\varphi \right)e^{-(\gamma -\rho )}-\sin \left(\varphi \right)\cos \left(\varphi \right)e^{-2\alpha }\\ \text{(23)}& & =\cos \left(\varphi \right)\sin \left(\varphi \right)[e^{-2\alpha }-e^{-(\gamma -\rho )}]\text{(24)}& {\gamma }^{xz}& =0\text{(25)}& {\gamma }^{yz}& =0\text{(26)}& {\gamma }^{zz}& =e^{-2\alpha }.\end{array}$$

The 4-velocity in spherical coordinates is given by

$$\begin{array}{}\text{(27)}& u^{\overline{a}}=\frac{e^{-(\rho +\gamma )/2}}{\sqrt{1-v^2}}[1,0,0,\mathrm{\Omega }],\end{array}$$

where

$$\begin{array}{}\text{(28)}& v=(\mathrm{\Omega }-\omega )r\sin \left(\theta \right)e^{-\rho }.\end{array}$$

Transforming to Cartesian coordinates we have

$$\begin{array}{}\text{(29)}& u^t& =\frac{e^{-(\rho +\gamma )/2}}{\sqrt{1-v^2}}\text{(30)}& u^x& ={\partial }_{\varphi }x{u}^{\varphi }=-r\sin \left(\theta \right)\sin \left(\varphi \right)u^t\mathrm{\Omega }\text{(31)}& u^y& ={\partial }_{\varphi }y{u}^{\varphi }=r\sin \left(\theta \right)\cos \left(\varphi \right)u^t\mathrm{\Omega }\text{(32)}& u^z& =0.\end{array}$$

The Lorentz factor is given by

$$\begin{array}{rl}W& =N{u}^t=e^{(\gamma +\rho )/2}\frac{e^{-(\rho +\gamma )/2}}{\sqrt{1-v^2}}\\ \text{(33)}& & =\frac{1}{\sqrt{1-v^2}}.\end{array}$$

Using

$$\begin{array}{}\text{(34)}& v^i=\frac{1}{N}(\frac{u^i}{u^t}+{\beta }^i)\end{array}$$

we get

$$\begin{array}{}\text{(35)}& v^x& =-e^{-(\gamma +\rho )/2}r\sin \left(\theta \right)\sin \left(\varphi \right)(\mathrm{\Omega }-\omega )=-e^{-(\gamma -\rho )/2}\sin \left(\varphi \right)v\text{(36)}& v^y& =e^{-(\gamma +\rho )/2}r\sin \left(\theta \right)\cos \left(\varphi \right)(\mathrm{\Omega }-\omega )=e^{-(\gamma -\rho )/2}\cos \left(\varphi \right)v\text{(37)}& v^z& =0.\end{array}$$

Lowering with the spatial metric we get

$$\begin{array}{rl}{v}_x& ={\gamma }_{xx}{v}^x+{\gamma }_{xy}{v}^y\\ & =-[{\cos }^2(\varphi )e^{2\alpha }+{\sin }^2(\varphi )e^{\gamma -\rho }]e^{-(\gamma -\rho )/2}\sin \left(\varphi \right)v\\ & +[\cos \left(\varphi \right)\sin \left(\varphi \right)e^{2\alpha }-\sin \left(\varphi \right)\cos \left(\varphi \right)e^{\gamma -\rho }]e^{-(\gamma -\rho )/2}\cos \left(\varphi \right)v\\ & =-e^{-(\gamma -\rho )/2}v\sin \left(\varphi \right)e^{\gamma -\rho }[{\sin }^2(\varphi )+{\cos }^2(\varphi )]\\ \text{(38)}& & =-e^{(\gamma -\rho )/2}v\sin \left(\varphi \right)v_y& ={\gamma }_{yx}{v}^x+{\gamma }_{yy}{v}^y\\ & =-[\cos \left(\varphi \right)\sin \left(\varphi \right)e^{2\alpha }-\sin \left(\varphi \right)\cos \left(\varphi \right)e^{(\gamma -\rho )}]e^{-(\gamma -\rho )/2}\sin \left(\varphi \right)v\\ & +[{\sin }^2(\varphi )e^{2\alpha }+{\cos }^2(\varphi )e^{(\gamma -\rho )}]e^{-(\gamma -\rho )/2}\cos \left(\varphi \right)v\\ \text{(39)}& & =e^{(\gamma -\rho )/2}v\cos \left(\varphi \right)\text{(40)}& v_z& =0.\end{array}$$

This is consistent with the Lorentz factor read off from $u^t$ since $v^i{v}_i=v^2$. For completeness, $u_i=W{v}_i$ so

$$\begin{array}{}\text{(41)}& u_x& =-\frac{e^{(\gamma -\rho )/2}v\sin \left(\varphi \right)}{\sqrt{1-v^2}}\text{(42)}& u_y& =\frac{e^{(\gamma -\rho )/2}v\cos \left(\varphi \right)}{\sqrt{1-v^2}}\text{(43)}& u_z& =0.\end{array}$$

Warning

Near (within 1e-2) $r=0$ the numerical errors from interpolation and computing the metric derivatives by finite difference no longer cancel out and so the tilde_s time derivative only vanishes to roughly 1e-8 rather than machine precision. Computing the Cartesian derivatives from analytic differentiation of the radial and angular polynomial fits might improve the situation but is a decent about of work to implement.

Member Function Documentation

get_clone()

auto RelativisticEuler::Solutions::RotatingStar::get_clone ( ) const -> std::unique_ptr< evolution::initial_data::InitialData >

overridevirtual

Implements evolution::initial_data::InitialData.

Member Data Documentation

help

constexpr Options::String RelativisticEuler::Solutions::RotatingStar::help

staticconstexpr

Initial value:

= {
      "Rotating neutron star initial data solved by the RotNS solver. The data "
      "is read in from disk. If RotNS used a polytropic equation of state, use "
      "the PolytropicConstant option, so that the polytropic index may be read "
      "directly from the RotNS output. Otherwise, set the equation of state to "
      "whichever one you would like to use to load the initial data. Note that "
      "if a polytrope is put into the EquationOfState option rather than using "
      "the PolytropicConstant option, the generic unit conversion will be used "
      "instead of the specific polytrope one."}

---

The documentation for this class was generated from the following file:

• src/PointwiseFunctions/AnalyticSolutions/RelativisticEuler/RotatingStar.hpp