Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Erick/production 1.0 #106

Open
wants to merge 74 commits into
base: development
Choose a base branch
from
Open

Erick/production 1.0 #106

wants to merge 74 commits into from

Conversation

erickurquilla1999
Copy link
Collaborator

No description provided.

shankar-1729 and others added 30 commits November 4, 2024 14:07
@shankar-1729
Make IMFP calculation energy dependent
7a53b35
@shankar-1729
Ensure that rho, temp and Ye are uniformally spaced.
de0bac5
@shankar-1729 @erickurquilla1999
Implement function to calculate maximum of IMFP over the grid (UNTESTED)
7daa907
@shankar-1729 @erickurquilla1999
Move multifab for IMFB inside compute_max_IMFP function, correct temp…
2c67c2e 
…erature units
@shankar-1729 @erickurquilla1999
do not input ba, dm and ngrow inside compute_max_IMFP function
795615c
@shankar-1729 @erickurquilla1999
calculate max_IMFP_abs inside function compute_dt
a9d22c4 
Currently compute_dt runs every iteration. We may need to run compute_dt just once initially during production, since (rho,T,Ye) is static.
@shankar-1729 @erickurquilla1999
Remove dt calculation from compute_dt function
66eeb3c
@erickurquilla1999
Adding basics comments and documentation to new functions
3480a61
@erickurquilla1999
Change the position of the calculation of the maximum IMFP
8458eae
@erickurquilla1999
Compute maximum of IMFP when IMFP method is 1 too
792c8fd
@erickurquilla1999
Compute the maximum IMFP and pass it as a parameter to the compute_dt…
4b43fff 
… function.
@erickurquilla1999
Cleaning up the compute_dt function and adding documentation.
db60725
@erickurquilla1999
Include a generated file to compute the trace of the neutrino number …
d02f865 
…density for mesh quantities
@erickurquilla1999
Solving issue in generated file to compute trace of mesh neutrino and…
0959355 
… antineutrino number densities
@erickurquilla1999
Compute maximum trace of the neutrino and antineutrino from mesh quan…
b4c008b 
…tities to compute the correct time step
@erickurquilla1999
Create a function that computes the minimum value in a MultiFab in pa…
0399808 
…rallel.

This function will be used to compute the minimum trace of mesh quantities, and based on that, it will compute the optimum time step
@erickurquilla1999
Set some if statements to handle simulations with empty periodic boun…
be783b9 
…dary conditions and black holes.
@erickurquilla1999
Generate a new file to compute the trace on mesh quantities
d8e9f1e
@erickurquilla1999
Solve syntaxis mistakes and compute the minumin time step rather than…
c5c70ed 
… the minimum potential
@erickurquilla1999
Make function compute_max_IMFP return a multifab with the maximum inv…
9703bac 
…erse mean free paths in the grid center
@erickurquilla1999
Compute the minimum combination of TrN/Opacity*c. This quantity provi…
aada239 
…des a better representation of the time derivative of N than the previous 1/Opacity*c.
@erickurquilla1999
Delete commented lines and handle periodic empty boundary conditions …
8a737ea 
…and black holes in the function that computes the time step.
@erickurquilla1999
Cleaning up and addiing comments to the compute_dt function
3420c23
@erickurquilla1999
Add documentation for header file of Evolve.cpp file
4594186
@erickurquilla1999
Delete unnecessary generation of two files that compute the trace.
acbf75d
@erickurquilla1999
Setting the min(TrN, TrNbar) to one in case it is zero
2296774
@erickurquilla1999
Cleaning space
97e981e
@erickurquilla1999
Cleaning up compute_max_IMFP function
c1be327
@erickurquilla1999
Creating documentation for compute_max_IMFP function
af7a2c8
@erickurquilla1999
Set the new time step to agree with the new function that computes th…
08fc11f 
…e time step.
@erickurquilla1999
Including the minimum time step barrier in the function that computes…
981c3e2 
… time steps.
@erickurquilla1999
Including a new neutrino container to store the minimum dt for every …
94af9fb 
…particle variable
@erickurquilla1999
Delete documentation in Evolve.H header file and input the time step …
9255937 
…neutrino container in the function that compute the time step
@erickurquilla1999
Add a flag to choose between two different methods for computing the …
bd4be86 
…time step.
@erickurquilla1999
Add a new method to compute the time step based on the time derivativ…
2835ede 
…e of N and Nbar for all particles.
@erickurquilla1999
Delete unnecessary variables and comments.
e49ff1b
@erickurquilla1999
Add time_step_method flag to input parameters.
d437789
@erickurquilla1999
Update input parameters.
4873d89
@erickurquilla1999
Improve script that compute dE and antineutrino chemical potential fo…
f5e9e71 
…r collisional flavor instability test
@erickurquilla1999
Setting asserts at the end of the collisional flavor instability test…
f37885e 
…. In this way the test will always produce a plot that can be convinient for visualization of error
@erickurquilla1999
Fix issue with Method 1 for computing the time step
5ed145b
@erickurquilla1999
Set collisional flavor instability test to use Method 1 for computing…
c7dbd0d 
… the time step
@erickurquilla1999
Fix issue in method 0 to compute time step
e8f4b84
@erickurquilla1999
Adding comments in st9 initial condition to avoid future confusion
d6c37a4
@erickurquilla1999
Adding comments to avoid future confusion in fermi-dirac test
653bc78
@erickurquilla1999
Remove 'FIXME' comments.
40a853b
@erickurquilla1999
Return mf_IMFP at the end of the function instead of in every if stat…
34f27cf 
…ement.
@erickurquilla1999
Delete print line in main
a5f4b30
@erickurquilla1999
Create a script to plot the polarization vector components for single…
c2aa1b4 
… particles
@erickurquilla1999
Add an if statement to compute the IMFP MultiFab only when the time s…
501fe69 
…tep method is zero
@erickurquilla1999
Set the maximum grid size for every direction in the domain.
1652a0e
@erickurquilla1999
Update input file to define max grid size in every direction in the d…
61c025d 
…omain.
@erickurquilla1999
Making convertToHDF5.py script works in 3D simulations
86e8500
@erickurquilla1999
Adding a line to stop execution when a NaN appears in N and N bar
0de65c0
@erickurquilla1999
Improve section the leave the code is N or Nbar is NaN
e02dc74
@erickurquilla1999
Plotting best figures in collisional instability test
9ea0ae6
@erickurquilla1999
Adding units to labels and legends
8b44a36
@erickurquilla1999
Updating plots and legens in CFI test
104bb0f
@erickurquilla1999
Delete lines that import unused packages in the collisional instabili…
8c92437 
…ty test.
@erickurquilla1999
Setting optimal max_grid_size in tests input file
4e61f3d
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Reviewers

@srichers srichers srichers requested changes

@shankar-1729 shankar-1729 Awaiting requested review from shankar-1729

Requested changes must be addressed to merge this pull request.

Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

Successfully merging this pull request may close these issues.
3 participants
@erickurquilla1999 @srichers @shankar-1729