Skip to content

AsterSeeds: add code to perturb neutron star#145

Open
jaykalinani wants to merge 7 commits into
EinsteinToolkit:devfrom
jaykalinani:perturbNS
Open

AsterSeeds: add code to perturb neutron star#145
jaykalinani wants to merge 7 commits into
EinsteinToolkit:devfrom
jaykalinani:perturbNS

Conversation

@jaykalinani

@jaykalinani jaykalinani commented Jun 16, 2026

Copy link
Copy Markdown
Collaborator

Different ways to perturb initial data of a neutron star are added. Code for external dipolar field for single NS is added as well. Minor changes to function names as well as schedule.

Scripts for CompOSE EOS table are also added.
CompOSE tables also have a bug in the spelling. EOSX now includes the spelling error for consistency.

@jaykalinani jaykalinani requested a review from acwen11 June 16, 2026 03:07

@acwen11 acwen11 left a comment

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Minor suggestions.

Comment thread AsterSeeds/schedule.ccl
} "Set initial entropy"
if (set_entropy_postinitial && !CCTK_Equals(evolution_eos, "Tabulated3d")) {
if (set_temperature_postinitial) {
SCHEDULE SetEntropy IN HydroBaseX_PostInitial AFTER SetTemp

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Could this not be simplified to just say SCHEDULE SetEntropy IN HydroBaseX_PostInitial AFTER SetTemp AFTER SetYe?

Comment thread AsterSeeds/schedule.ccl
WRITES: HydroBaseX::rho(everywhere) HydroBaseX::eps(everywhere) HydroBaseX::press(everywhere) HydroBaseX::temperature(everywhere) HydroBaseX::entropy(everywhere)
} "Apply neutron-star pressure or density perturbation after hydro initial data"
} else {
SCHEDULE AsterSeeds_PerturbNS IN HydroBaseX_PostInitial

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Same comment as above

Comment thread EOSX/scripts/CompOSE/generate_beta_eq_slice.py
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants