Code

Erick/production 1.0#106

Code

Open

erickurquilla1999 wants to merge 74 commits intodevelopmentAMReX-Astro/Emu:developmentfrom erick/production_1.0AMReX-Astro/Emu:erick/production_1.0Copy head branch name to clipboard

+930-215

Commits

Commits on Nov 4, 2024

• Make IMFP calculation energy dependent

  shankar-1729

  committed
  7a53b35

• Ensure that rho, temp and Ye are uniformally spaced.

  shankar-1729

  committed
  de0bac5

Commits on Nov 7, 2024

• Implement function to calculate maximum of IMFP over the grid (UNTESTED)

  

  shankar-1729

  authored and

  erickurquilla1999

  committed
  7daa907

• Move multifab for IMFB inside compute_max_IMFP function, correct temperature units

  

  shankar-1729

  authored and

  erickurquilla1999

  committed
  2c67c2e

• do not input ba, dm and ngrow inside compute_max_IMFP function

  

  shankar-1729

  authored and

  erickurquilla1999

  committed
  795615c

• calculate max_IMFP_abs inside function compute_dt

  Show description for a9d22c4
  

  shankar-1729

  authored and

  erickurquilla1999

  committed
  a9d22c4

• Remove dt calculation from compute_dt function

  

  shankar-1729

  authored and

  erickurquilla1999

  committed
  66eeb3c

• Adding basics comments and documentation to new functions

  erickurquilla1999

  committed
  3480a61

• Change the position of the calculation of the maximum IMFP

  erickurquilla1999

  committed
  8458eae

• Compute maximum of IMFP when IMFP method is 1 too

  erickurquilla1999

  committed
  792c8fd

• Compute the maximum IMFP and pass it as a parameter to the compute_dt function.

  erickurquilla1999

  committed
  4b43fff

• Cleaning up the compute_dt function and adding documentation.

  erickurquilla1999

  committed
  db60725

• Include a generated file to compute the trace of the neutrino number density for mesh quantities

  erickurquilla1999

  committed
  d02f865

• Solving issue in generated file to compute trace of mesh neutrino and antineutrino number densities

  erickurquilla1999

  committed
  0959355

• Compute maximum trace of the neutrino and antineutrino from mesh quantities to compute the correct time step

  erickurquilla1999

  committed
  b4c008b

• Create a function that computes the minimum value in a MultiFab in parallel.

  Show description for 0399808
  erickurquilla1999

  committed
  0399808

• Set some if statements to handle simulations with empty periodic boundary conditions and black holes.

  erickurquilla1999

  committed
  be783b9

• Generate a new file to compute the trace on mesh quantities

  erickurquilla1999

  committed
  d8e9f1e

• Solve syntaxis mistakes and compute the minumin time step rather than the minimum potential

  erickurquilla1999

  committed
  c5c70ed

• Make function compute_max_IMFP return a multifab with the maximum inverse mean free paths in the grid center

  erickurquilla1999

  committed
  9703bac

• Compute the minimum combination of TrN/Opacity*c. This quantity provides a better representation of the time derivative of N than the previous 1/Opacity*c.

  erickurquilla1999

  committed
  aada239

• Delete commented lines and handle periodic empty boundary conditions and black holes in the function that computes the time step.

  erickurquilla1999

  committed
  8a737ea

• Cleaning up and addiing comments to the compute_dt function

  erickurquilla1999

  committed
  3420c23

• Add documentation for header file of Evolve.cpp file

  erickurquilla1999

  committed
  4594186

• Delete unnecessary generation of two files that compute the trace.

  erickurquilla1999

  committed
  acbf75d

• Setting the min(TrN, TrNbar) to one in case it is zero

  erickurquilla1999

  committed
  2296774

• Cleaning space

  erickurquilla1999

  committed
  97e981e

• Cleaning up compute_max_IMFP function

  erickurquilla1999

  committed
  c1be327

• Creating documentation for compute_max_IMFP function

  erickurquilla1999

  committed
  af7a2c8

• Set the new time step to agree with the new function that computes the time step.

  erickurquilla1999

  committed
  08fc11f

• Solving issue when looping over the MultiFabs. Now it does not loop over ghost cells since the MultiFab containing the maximum IMFP for every cell does not have ghost cells.

  erickurquilla1999

  committed
  dac02e9

• Delete unnecessary lines

  erickurquilla1999

  committed
  f76d6fb

• Update empty periodic boundary condition test

  erickurquilla1999

  committed
  1021c4e

• Setting the right cfl factors for particles to equilibrium test

  erickurquilla1999

  committed
  aa6f913

• Set up the Fermi-Dirac test for multi-energy particles with the following time step CFL factors.

  erickurquilla1999

  committed
  ec29a9b

• Delete unnecessary functions

  erickurquilla1999

  committed
  8f7094b

• Solving possible issue with zero inverse mean free path

  erickurquilla1999

  committed
  a114cc8

• Improve documentation

  erickurquilla1999

  committed
  7358916

• Solve syntax issue

  erickurquilla1999

  committed
  2a3525d

• Solving issue with fermi dirac test

  Show description for 94a9bf2
  erickurquilla1999

  committed
  94a9bf2

Commits on Nov 25, 2024

• Fix time step function: now calculates the time step based on a basic equation. A more advanced time step calculation will be implemented later

  erickurquilla1999

  committed
  3f60f26

• Adding comments on the definitions of V_stupid and V_adaptive

  erickurquilla1999

  committed
  f1c9cd7

• Delete unnecessary line that computes the trace of matrices.

  erickurquilla1999

  committed
  216fc2a

• Add input parameter for the minimum allowed time step in the simulation

  erickurquilla1999

  committed
  f09220b

• Including the minimum time step barrier in the function that computes time steps.

  erickurquilla1999

  committed
  981c3e2

• Including a new neutrino container to store the minimum dt for every particle variable

  erickurquilla1999

  committed
  94af9fb

• Delete documentation in Evolve.H header file and input the time step neutrino container in the function that compute the time step

  erickurquilla1999

  committed
  9255937

• Add a flag to choose between two different methods for computing the time step.

  erickurquilla1999

  committed
  bd4be86

• Add a new method to compute the time step based on the time derivative of N and Nbar for all particles.

  erickurquilla1999

  committed
  2835ede

• Delete unnecessary variables and comments.

  erickurquilla1999

  committed
  e49ff1b

• Add time_step_method flag to input parameters.

  erickurquilla1999

  committed
  d437789

• Update input parameters.

  erickurquilla1999

  committed
  4873d89

• Improve script that compute dE and antineutrino chemical potential for collisional flavor instability test

  erickurquilla1999

  committed
  f5e9e71

• Setting asserts at the end of the collisional flavor instability test. In this way the test will always produce a plot that can be convinient for visualization of error

  erickurquilla1999

  committed
  f37885e

• Fix issue with Method 1 for computing the time step

  erickurquilla1999

  committed
  5ed145b

• Set collisional flavor instability test to use Method 1 for computing the time step

  erickurquilla1999

  committed
  c7dbd0d

• Fix issue in method 0 to compute time step

  erickurquilla1999

  committed
  e8f4b84

• Adding comments in st9 initial condition to avoid future confusion

  erickurquilla1999

  committed
  d6c37a4

• Adding comments to avoid future confusion in fermi-dirac test

  erickurquilla1999

  committed
  653bc78

• Remove 'FIXME' comments.

  erickurquilla1999

  committed
  40a853b

• Return mf_IMFP at the end of the function instead of in every if statement.

  erickurquilla1999

  committed
  34f27cf

• Delete print line in main

  erickurquilla1999

  committed
  a5f4b30

• Create a script to plot the polarization vector components for single particles

  erickurquilla1999

  committed
  c2aa1b4

• Add an if statement to compute the IMFP MultiFab only when the time step method is zero

  erickurquilla1999

  committed
  501fe69

• Set the maximum grid size for every direction in the domain.

  erickurquilla1999

  committed
  1652a0e

• Update input file to define max grid size in every direction in the domain.

  erickurquilla1999

  committed
  61c025d

Commits on Nov 27, 2024

• Making convertToHDF5.py script works in 3D simulations

  erickurquilla1999

  committed
  86e8500

Commits on Nov 29, 2024

• Adding a line to stop execution when a NaN appears in N and N bar

  erickurquilla1999

  committed
  0de65c0

• Improve section the leave the code is N or Nbar is NaN

  erickurquilla1999

  committed
  e02dc74

Commits on Dec 2, 2024

• Plotting best figures in collisional instability test

  erickurquilla1999

  committed
  9ea0ae6

• Adding units to labels and legends

  erickurquilla1999

  committed
  8b44a36

• Updating plots and legens in CFI test

  erickurquilla1999

  committed
  104bb0f

Commits on Dec 4, 2024

• Delete lines that import unused packages in the collisional instability test.

  erickurquilla1999

  committed
  8c92437

• Setting optimal max_grid_size in tests input file

  erickurquilla1999

  committed
  4e61f3d