Compare commits

...

2 commits

Author SHA1 Message Date
8c92adcbcf merge fixes 2025-06-02 18:19:32 +02:00
7d142b35b5 merge fixes 2025-06-02 18:18:53 +02:00
2 changed files with 0 additions and 16 deletions

View file

@ -133,17 +133,10 @@ def main(parsed_args):
case "allsCOLA":
print_message(f"Running ICs, pre_sCOLA, sCOLA and post_sCOLA.", 1, "control", verbose=parsed_args.verbose)
<<<<<<< HEAD:main.py
from parameters_card import parse_arguments_card, main_parameter_card
from timestepping import main_timestepping
from ICs import main_ICs
from scola import main_scola, main_pre_scola, main_post_scola, have_all_tiles
=======
from .parameters_card import parse_arguments_card, main_parameter_card
from .timestepping import main_timestepping
from .ICs import main_ICs
from .scola import main_scola, main_pre_scola, main_post_scola, have_all_tiles
>>>>>>> main:sbmy_control/main.py
card_dict = main_parameter_card(parsed_args)

View file

@ -47,7 +47,6 @@ def register_arguments_card(parser:ArgumentParser):
parser.add_argument("--OutputFCsDensity", type=str, default=None, help="Output FCs density file.")
parser.add_argument("--OutputFCsSnapshot", type=str, default=None, help="Output FCs snapshot file.")
parser.add_argument("--OutputRngStateLPT", type=str, default=None, help="Output RNG state file.")
<<<<<<< HEAD:parameters_card.py
## Tests with phiBCs and density
parser.add_argument("--WriteGravPot", type=bool, default=False, help="Write gravitational potential.")
parser.add_argument("--OutputGravitationalPotentialBase", type=str, default=None, help="Output gravitational potential base.")
@ -57,12 +56,10 @@ def register_arguments_card(parser:ArgumentParser):
parser.add_argument("--MeshDensity", type=int, default=None, help="Mesh for density.")
parser.add_argument("--LoadPhiBCs", type=bool, default=False, help="Load phiBCs.")
parser.add_argument("--InputPhiBCsBase", type=str, default=None, help="Input phiBCs file base.")
=======
parser.add_argument("--WriteReferenceFrame", type=bool, default=False, help="Write reference frame (COCA).")
parser.add_argument("--OutputMomentaBase", type=str, default=None, help="Output momenta base (COCA).")
parser.add_argument("--ReadReferenceFrame", type=bool, default=False, help="Read reference frame (COCA).")
parser.add_argument("--InputMomentaBase", type=str, default=None, help="Read momenta base (COCA).")
>>>>>>> main:sbmy_control/parameters_card.py
def register_arguments_card_for_ICs(parser:ArgumentParser):
@ -233,7 +230,6 @@ def parse_arguments_card(parsed_args):
card_dict["OutputFCsSnapshot"] = main_dict["resultdir"]+"final_particles_"+main_dict["simname"]+".gadget3"
if card_dict["OutputRngStateLPT"] is None:
card_dict["OutputRngStateLPT"] = main_dict["workdir"]+"rng_state.h5"
<<<<<<< HEAD:parameters_card.py
## Tests with phiBCs and density
if card_dict["OutputGravitationalPotentialBase"] is None:
card_dict["OutputGravitationalPotentialBase"] = main_dict["workdir"]+"gravpot_"+main_dict["simname"]
@ -245,12 +241,10 @@ def parse_arguments_card(parsed_args):
card_dict["MeshDensity"] = card_dict["N_PM_mesh"]
if card_dict["InputPhiBCsBase"] is None:
card_dict["InputPhiBCsBase"] = main_dict["workdir"]+"gravpot_tCOLA"
=======
if card_dict["OutputMomentaBase"] is None:
card_dict["OutputMomentaBase"] = main_dict["workdir"]+"momenta_"+main_dict["simname"]+"_"
if card_dict["InputMomentaBase"] is None:
card_dict["InputMomentaBase"] = main_dict["workdir"]+"momenta_"+main_dict["simname"]+"_"
>>>>>>> main:sbmy_control/parameters_card.py
return card_dict
@ -369,7 +363,6 @@ def create_parameter_card_dict(
OutputFCsDensity:str = 'fcs_density.h5',
OutputFCsSnapshot:str = 'fcs_particles.gadget3',
<<<<<<< HEAD:parameters_card.py
## Tests with phiBCs and density
WriteGravPot:bool = True,
OutputGravitationalPotentialBase:str = 'gravitational_potential.h5',
@ -379,13 +372,11 @@ def create_parameter_card_dict(
MeshDensity:int = 128,
LoadPhiBCs:bool = False,
InputPhiBCsBase:str = 'gravitational_potential.h5',
=======
## COCA parameters
WriteReferenceFrame:bool = False,
OutputMomentaBase:str = 'momenta_',
ReadReferenceFrame:bool = False,
InputMomentaBase:str = 'momenta_',
>>>>>>> main:sbmy_control/parameters_card.py
## Cosmological parameters
h:float = 0.6732,