Skip to content

Feature/adaptive moving puncture#171

Open
kazewong wants to merge 4 commits intomainfrom
feature/adaptive_moving_puncture
Open

Feature/adaptive moving puncture#171
kazewong wants to merge 4 commits intomainfrom
feature/adaptive_moving_puncture

Conversation

@kazewong
Copy link
Contributor

@kazewong kazewong commented Apr 2, 2021

Adding adaptive moving puncture. The initial request contains 3 commits:

  1. Passing cell coordinate to CCZ4RHS and MovingPunctureGauge

    This commit passes the cell coordinate object to the CCZ4RHS and subsequently the gauge. It is expected to be compatible with the original main branch (commit #a48c11813e861578deeb27c289c1583987442c69) since no actual change is made in the computation.

  2. Change MovingPunctureGauge class to use adaptive gaussian around each puncture

    Add gaussian smoothing to advection coefficient according to the puncture masses. This commit is not expected to be compatible with the default behaviour, which assumes eta to be constant throughout the whole simulation domain. I think there should be a way to switch on and off this adaptive behaviour. A better way to define this will be making a new gauge class, however, I failed to compile the code when I create a new gauge class and try to pass it to CCZ4RHS<>{} in BinaryBHLevel.cpp. It complains about the need to use the original MovingPunctureGauge params, which is not defined in the new gauge class.

  3. Update MovingPunctureGauge adaptive parameters in BinaryBHLevel.cpp

    Make corresponding changes in BinaryBHLevel.cpp to update the puncture location in the gauge during the simulation.

@kazewong kazewong requested a review from mirenradia April 2, 2021 15:34
@kazewong kazewong self-assigned this Apr 25, 2021
@kazewong kazewong force-pushed the feature/adaptive_moving_puncture branch from a803fe5 to 90ed830 Compare April 26, 2021 21:07
mirenradia and others added 4 commits May 5, 2021 10:20
* Prepend hdf5_path to restart_file when checking

Fixes #172

* Add more CH_USE_HDF5 guards

The guarded code is unnecessary if building without HDF5.

* Move `pout_prefix` under CH_MPI guards

This parameter is irrelevant when compiling without MPI.
@kazewong kazewong force-pushed the feature/adaptive_moving_puncture branch from c87547d to 6b786ed Compare May 5, 2021 14:22
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