fix xray_source turnover bug, improve box summary printing

src/py21cmfast/drivers/single_field.py

@@ -715,8 +715,11 @@ def compute_xray_source_field(
         R_outer = R_range[i].to("Mpc").value
 
         if zpp_avg[i] >= z_max:
+            box.filtered_sfr[i] = 0
+            box.filtered_sfr_mini[i] = 0
+            box.filtered_xray[i] = 0
+            box.mean_log10_Mcrit_LW[i] = inputs.astro_params.M_TURN  # minimum
             logger.debug(f"ignoring Radius {i} which is above Z_HEAT_MAX")
-            box.filtered_sfr[i, ...] = 0
             continue
 
         hbox_interp = interp_halo_boxes(
@@ -730,6 +733,7 @@ def compute_xray_source_field(
             box.filtered_sfr[i] = 0
             box.filtered_sfr_mini[i] = 0
             box.filtered_xray[i] = 0
+            box.mean_log10_Mcrit_LW[i] = hbox_interp.log10_Mcrit_MCG_ave
             logger.debug(f"ignoring Radius {i} due to no stars")
             continue
 

src/py21cmfast/wrapper/structs.py

@@ -1082,26 +1082,33 @@ def ensure_input_computed(self, input_box, load=False) -> bool:
     def summarize(self, indent=0) -> str:
         """Generate a string summary of the struct."""
         indent = indent * "    "
-        out = f"\n{indent}{self.__class__.__name__}\n"
 
-        out += "".join(
-            f"{indent}    {fieldname:>15}: {getattr(self, fieldname, 'non-existent')}\n"
-            for fieldname in self.struct.primitive_fields
+        # print array type and column headings
+        out = (
+            f"\n{indent}{self.__class__.__name__:>25}    "
+            + "   1st:         End:         Min:         Max:         Mean:         \n"
         )
 
+        # print array extrema and means
         for fieldname, state in self._array_state.items():
             if not state.initialized:
-                out += f"{indent}    {fieldname:>15}: uninitialized\n"
+                out += f"{indent}    {fieldname:>25}:  uninitialized\n"
             elif not state.computed:
-                out += f"{indent}    {fieldname:>15}: initialized\n"
+                out += f"{indent}    {fieldname:>25}:  initialized\n"
             elif not state.computed_in_mem:
-                out += f"{indent}    {fieldname:>15}: computed on disk\n"
+                out += f"{indent}    {fieldname:>25}:  computed on disk\n"
             else:
                 x = getattr(self, fieldname).flatten()
                 if len(x) > 0:
-                    out += f"{indent}    {fieldname:>15}: {x[0]:1.4e}, {x[-1]:1.4e}, {x.min():1.4e}, {x.max():1.4e}, {np.mean(x):1.4e}\n"
+                    out += f"{indent}    {fieldname:>25}: {x[0]:11.4e}, {x[-1]:11.4e}, {x.min():11.4e}, {x.max():11.4e}, {np.mean(x):11.4e}\n"
                 else:
-                    out += f"{indent}    {fieldname:>15}: size zero\n"
+                    out += f"{indent}    {fieldname:>25}: size zero\n"
+
+        # print primitive fields
+        out += "".join(
+            f"{indent}    {fieldname:>25}: {getattr(self, fieldname, 'non-existent')}\n"
+            for fieldname in self.struct.primitive_fields
+        )
 
         return out
 

tests/test_segfaults.py

@@ -1,10 +1,14 @@
 """
-This file contains tests which run a lightcone under various flag options to test the C backend for segfaults.
-These will not test the outputs of the run, just that the run completes without error
+This file contains tests which run a lightcone under various flag options
+to test the C backend for segfaults.
+These will not test the outputs of the run past the fact that they are finite,
+just that the run completes without error
 """
 
 import pytest
 
+import numpy as np
+
 import py21cmfast as p21c
 
 from . import produce_integration_test_data as prd
@@ -131,3 +135,4 @@ def test_lc_runs(name, max_redshift):
         )
 
     assert isinstance(lightcone, p21c.LightCone)
+    assert np.all(np.isfinite(lightcone.brightness_temp))