local angular distributions #36

Yurlungur · 2024-11-17T22:25:59Z

PR Summary

Add local angular distributions code in support of neutrino oscillations. This machinery is intended to find electron lepton number crossings. See, e.g.,
https://arxiv.org/abs/2404.17938

PR Checklist

Adds a test for any bugs fixed. Adds tests for new features.

…hlight into jmm/local-angular-distributions

Yurlungur · 2024-11-21T02:24:55Z

Four species transport now supported. To use do the following:

Clone my fork of NuLib
Modify the file make_example_table.F90 and set mytable_neutrino_scheme to 2 and number_output_species to 4
build with make -j
run the resulting executable make_example_table
NuLib/scripts/table2bhlight.py on the generated hdf5 file
Copy the generated table, which should end with _bhlight.h5 to the location you wish to use for nubhlight
Modify the nublight build.py file or param_template.dat file to point to the file
In the nubhlight build.py file, either add this line bhl.config.set_cparm("RAD_NUM_TYPES", 4) or change RAD_NUM_TYPES to be 4.
build nubhlight and run

Below is the fornax 1 zone equilibrium test run with nue, nuebar, nux, nuxbar:

Note that with four species, at fixed particle number, fewer monte carlo packets are carrying electron number, so to get good Ye values, you may need to add roughly 25% more particles.

Yurlungur · 2024-11-21T04:16:09Z

important paper also for this MR https://arxiv.org/abs/2304.05044

…e to pull out the ability to get the bin indices

Yurlungur · 2024-11-23T04:59:19Z

core/constants.h

@@ -28,6 +28,11 @@
 #define S4THW (S2THW * S2THW)
 #define NUSIGMA0 (1.7611737037e-44) // Fundamental neutrino cross section

+// Frequency scale of neutrino oscillations
+#define ROOT2 (1.4142135623730951)


Square roots are slow in C

Yurlungur · 2024-11-23T04:59:58Z

core/decs.h

+#define LOCAL_NUM_BASES (2)
+#define MOMENTS_A (0)
+#define MOMENTS_B (1)


might want to move these up...

Yurlungur · 2024-11-23T05:00:26Z

core/decs.h

@@ -1065,6 +1104,7 @@ void   set_cooling_time(
 void record_lepton_flux(const struct of_photon *ph);
 void check_nu_type(const char *location);
 int  get_lepton_sign(const struct of_photon *ph);
+int  nu_is_heavy(const int radtype);


needed for generalizing to 4 species transport

Yurlungur · 2024-11-23T05:01:39Z

core/io.c

+    double *local_angles_Xharm = safe_malloc(
+        (NDIM - 1) * LOCAL_ANGLES_NX1 * LOCAL_ANGLES_NX2 * sizeof(double));
+#if METRIC == MKS
+    double *local_angles_Xbl = safe_malloc(
+        (NDIM - 1) * LOCAL_ANGLES_NX1 * LOCAL_ANGLES_NX2 * sizeof(double));
+#endif // MKS
+    double *local_angles_Xcart = safe_malloc(
+        (NDIM - 1) * LOCAL_ANGLES_NX1 * LOCAL_ANGLES_NX2 * sizeof(double));


these are diagnostics.

Yurlungur · 2024-11-23T05:01:55Z

core/io.c

+  }
+
+  {
+    double *local_angles_mu = safe_malloc(LOCAL_ANGLES_NMU * sizeof(double));


diagnostics

Yurlungur · 2024-11-23T05:10:20Z

prob/multiscatt/build.py

The changes here are maybe best left as intermediate tests. Should probably undo these changes before merge.

Yurlungur · 2024-11-23T05:10:38Z

prob/oscillations/build.py

this test isn't ready yet.

Yurlungur · 2024-11-23T05:10:48Z

prob/oscillations/problem.c

this test isn't ready yet.

Yurlungur · 2024-11-23T05:11:15Z

prob/torus_cbc/build.py

we should think about the changes here. We may wish to make a new build file rather than continuing to use the old one.

Yurlungur · 2024-11-23T05:12:22Z

script/machine/darwin.py

this is my configuration. May or may not be wise to keep.

Yurlungur · 2024-11-23T23:04:38Z

This MR is ready. @kelslund @payelmuk150 please review.

Yurlungur · 2024-11-23T23:05:26Z

One-zone oscillation test

jonahm-LANL added 3 commits July 10, 2024 14:40

add ability to generate histograms

6c5351f

atomic!

1e4ad19

oscillations in situ machinery seems to run

bafaacc

Yurlungur added the enhancement New feature or request label Nov 17, 2024

Yurlungur self-assigned this Nov 17, 2024

jonahm-LANL added 20 commits November 18, 2024 08:05

loop over photons broken with tracers. Now fixed.

021287c

add local angular distributions to torus_cbc

14d9b3c

darwin

5538e6e

fundamental flaw... coords only specified for Minkowski. now fixed.

e951146

Merge branch 'jmm/local-angular-distributions' of github.com:lanl/nub…

437fbe3

…hlight into jmm/local-angular-distributions

testing

deaee4c

whee

2902ca7

infinity...

558dd53

why is G = 0

7445dfb

pass in struct pointer???

9a8f349

fabs

722b6aa

progress

8f97909

try normalizing kspace

993a836

cleanup...

2b22cde

one more try

426e855

i > mu

31d8c80

should be gcov not gcon

534e4d5

cleanup

8269258

Support 4 species transport

cdfff36

fornax 1zone test modified for 4 species NuLib table

4ec9374

jonahm-LANL added 3 commits November 21, 2024 15:38

output local angles grid positions and mu

1bf5e6e

Define Gnu, moment types but don't set them yet. Refactor binning cod…

4da3b69

…e to pull out the ability to get the bin indices

compute Gnu and moments

e9592fd

jonahm-LANL added 2 commits November 22, 2024 07:37

Better include guards

52b3dfe

Cleanup...

90a3010

Yurlungur requested a review from kelslund November 23, 2024 02:46

jonahm-LANL added 8 commits November 22, 2024 19:50

dt and oscillations

cebaaed

OSCILLATIONS

ada6bdc

oscillations doesn't depend on time. Make sure dt is in code units.

37c4ae9

thread it through step.c

7b6f7e5

better formatting

547ea4b

multiscatt getting crazier...

fa23e77

starting oscillations setup. not done yet.

5f824d9

oops didn't mean to commit that

fa0112e

Yurlungur commented Nov 23, 2024

View reviewed changes

jonahm-LANL added 6 commits November 22, 2024 22:22

NuLib opacities

039644e

debugging, performance, cleanup

4bcc147

oscillations problem test works

420f5f8

simplify/cleanup

903c927

automated test

b55f1f8

revert intermediate fornax and multiscatt tests

be16433

Yurlungur requested a review from payelmuk150 November 23, 2024 23:04

jonahm-LANL added 3 commits November 25, 2024 11:23

can do it badly

f4631f1

fix

c09f8a1

NOW the tests look right

cbb3f38

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

local angular distributions #36

local angular distributions #36

Yurlungur commented Nov 17, 2024

Yurlungur commented Nov 21, 2024 •

edited

Loading

Yurlungur commented Nov 21, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur Nov 23, 2024

Yurlungur commented Nov 23, 2024

Yurlungur commented Nov 23, 2024

local angular distributions #36

Are you sure you want to change the base?

local angular distributions #36

Conversation

Yurlungur commented Nov 17, 2024

PR Summary

PR Checklist

Yurlungur commented Nov 21, 2024 • edited Loading

Yurlungur commented Nov 21, 2024

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Yurlungur commented Nov 23, 2024

Yurlungur commented Nov 23, 2024

Yurlungur commented Nov 21, 2024 •

edited

Loading