diff --git a/jb5.log b/jb5.log
new file mode 100644
index 00000000000..5e28ff8bdaa
--- /dev/null
+++ b/jb5.log
@@ -0,0 +1,1101 @@
+{'text/html': ' \n '}
+
+============================= test session starts ==============================
+platform linux -- Python 3.12.4, pytest-8.2.2, pluggy-1.5.0 -- /home/aryaatharva18/miniforge3/envs/tardis_numpy/bin/python3.12
+cachedir: .pytest_cache
+metadata: {'Python': '3.12.4', 'Platform': 'Linux-5.15.0-67-generic-x86_64-with-glibc2.31', 'Packages': {'pytest': '8.2.2', 'pluggy': '1.5.0'}, 'Plugins': {'doctestplus': '1.2.1', 'anyio': '4.4.0', 'cov': '5.0.0', 'metadata': '3.1.1', 'html': '4.1.1'}}
+rootdir: /home/aryaatharva18/tardis-main/tardis
+configfile: pyproject.toml
+plugins: doctestplus-1.2.1, anyio-4.4.0, cov-5.0.0, metadata-3.1.1, html-4.1.1
+collecting ... collected 1 item
+
+tardis/tests/test_tardis_full.py::TestTransportSimple::test_j_blue_estimators [[1mtardis.io.model.parse_density_configuration[0m][[1;33mWARNING[0m]
+ Number of density points larger than number of shells. Assuming inner point irrelevant ([1mparse_density_configuration.py[0m:114)
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [9967.48844219 9952.56047019 9937.67714557 9922.83826832 9908.04363964
+ 9893.2930619 9878.58633865 9863.9232746 9849.30367564 9834.72734879] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.40039164 0.33245223 0.28966798 0.2577141 0.23224568 0.21120466
+ 0.19341188 0.17811402 0.16479446 0.1530809 ]
+ Shape: (20,)
+ dtype: float64
+
+Geometry boundaries:
+ v_inner_boundary_index: 0
+ v_outer_boundary_index: 20
+==================================================
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [9967.48844219 9952.56047019 9937.67714557 9922.83826832 9908.04363964
+ 9893.2930619 9878.58633865 9863.9232746 9849.30367564 9834.72734879] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.40039164 0.33245223 0.28966798 0.2577141 0.23224568 0.21120466
+ 0.19341188 0.17811402 0.16479446 0.1530809 ]
+ Shape: (20,)
+ dtype: float64
+==================================================
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[6.95817130e-92 4.17753016e-92 2.63191016e-92 1.69312480e-92
+ 1.10326201e-92 7.25462111e-93 4.80368606e-93 3.19867259e-93
+ 2.13990515e-93 1.43731896e-93 9.68770866e-94 6.54970950e-94
+ 4.44031918e-94 3.01773364e-94 2.05553393e-94 1.40301612e-94
+ 9.59453315e-95 6.57274047e-95 4.50999952e-95 3.09931792e-95]
+ [1.24522312e-88 7.56171680e-89 4.81859727e-89 3.13535917e-89
+ 2.06645401e-89 1.37439185e-89 9.20490286e-90 6.19959240e-90
+ 4.19504504e-90 2.84999446e-90 1.94294580e-90 1.32864944e-90
+ 9.11069100e-91 6.26277469e-91 4.31478628e-91 2.97883232e-91
+ 2.06042115e-91 1.42766854e-91 9.90845994e-92 6.88722935e-92]
+ [1.36934391e-88 8.31665272e-89 5.30043590e-89 3.44938009e-89
+ 2.27374783e-89 1.51248098e-89 1.01312109e-89 6.82445671e-90
+ 4.61853641e-90 3.13815620e-90 2.13970561e-90 1.46341178e-90
+ 1.00362218e-90 6.89999095e-91 4.75448870e-91 3.28286774e-91
+ 2.27104702e-91 1.57383904e-91 1.09245078e-91 7.59456783e-92]
+ [2.56405915e-88 1.55875782e-88 9.94388834e-89 6.47739346e-89
+ 4.27381872e-89 2.84562913e-89 1.90793821e-89 1.28642836e-89
+ 8.71437949e-90 5.92681304e-90 4.04496968e-90 2.76912376e-90
+ 1.90090602e-90 1.30813778e-90 9.02242667e-91 6.23573286e-91
+ 4.31792222e-91 2.99518449e-91 2.08103660e-91 1.44808974e-91]
+ [3.75410809e-86 2.29960754e-86 1.47818117e-86 9.70215145e-87
+ 6.45030446e-87 4.32751677e-87 2.92362128e-87 1.98627218e-87
+ 1.35576984e-87 9.29109870e-88 6.38936152e-88 4.40738413e-88
+ 3.04856549e-88 2.11390136e-88 1.46909857e-88 1.02308431e-88
+ 7.13830469e-89 4.98930658e-89 3.49295247e-89 2.44909021e-89]
+ [4.85853862e-86 2.97730361e-86 1.91455381e-86 1.25712503e-86
+ 8.36105637e-87 5.61164527e-87 3.79265272e-87 2.57769326e-87
+ 1.76014707e-87 1.20670350e-87 8.30159372e-88 5.72869143e-88
+ 3.96406318e-88 2.74979469e-88 1.91177598e-88 1.33188890e-88
+ 9.29655774e-89 6.50036575e-89 4.55261360e-89 3.19332712e-89]
+ [5.26925596e-86 3.22938941e-86 2.07691406e-86 1.36390173e-86
+ 9.07234206e-87 6.08978741e-87 4.11631553e-87 2.79801780e-87
+ 1.91082912e-87 1.31016842e-87 9.01450086e-88 6.22141678e-88
+ 4.30554459e-88 2.98704283e-88 2.07697772e-88 1.44715982e-88
+ 1.01023945e-88 7.06469804e-89 4.94846223e-89 3.47141499e-89]
+ [1.00335597e-85 6.15533956e-86 3.96256119e-86 2.60475207e-86
+ 1.73431740e-86 1.16529841e-86 7.88441604e-87 5.36459824e-87
+ 3.66719865e-87 2.51689772e-87 1.73342869e-87 1.19751082e-87
+ 8.29552922e-88 5.76080681e-88 4.00958620e-88 2.79646937e-88
+ 1.95408635e-88 1.36785129e-88 9.59050323e-89 6.73447161e-89]
+ [1.08451393e-85 6.65401144e-86 4.28409350e-86 2.81644181e-86
+ 1.87548850e-86 1.26030137e-86 8.52821656e-87 5.80333074e-87
+ 3.96758270e-87 2.72338192e-87 1.87585988e-87 1.29606056e-87
+ 8.97927716e-88 6.23637214e-88 4.34109908e-88 3.02804024e-88
+ 2.11615151e-88 1.48147158e-88 1.03883666e-88 7.29559689e-89]
+ [4.12554590e-84 2.54528363e-84 1.64785098e-84 1.08934671e-84
+ 7.29434167e-85 4.92892781e-85 3.35384217e-85 2.29492278e-85
+ 1.57769516e-85 1.08896020e-85 7.54241121e-86 5.24012271e-86
+ 3.65059699e-86 2.54953448e-86 1.78457546e-86 1.25170853e-86
+ 8.79619162e-87 6.19223744e-87 4.36624293e-87 3.08338539e-87]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Initialized j_blue_estimator with shape: (29224, 20)
+[J_BLUE_DEBUG] Initial values: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [10231.64445737 10291.83257817 10330.25756259 10347.40451513
+ 10359.55191309 10387.73833116 10354.85751326 10331.70472781
+ 10302.30271457 10274.61475832] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.4974496 0.38393696 0.31166896 0.26024359 0.22194936 0.19074918
+ 0.1697988 0.1521222 0.13776531 0.1258416 ]
+ Shape: (20,)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Input j_blue_estimator: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] Normalization factor: [2.85766110e+13 2.64164249e+13 2.44922250e+13 2.27708505e+13
+ 2.12247665e+13 1.98309571e+13 1.85700632e+13 1.74257072e+13
+ 1.63839583e+13 1.54329069e+13 1.45623223e+13 1.37633766e+13
+ 1.30284203e+13 1.23507981e+13 1.17246981e+13 1.11450264e+13
+ 1.06073037e+13 1.01075778e+13 9.64235139e+12 9.20852047e+12]
+[J_BLUE_DEBUG] After normalization: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] line_list_nu: [4.48979300e+16 4.33201055e+16 4.33000835e+16 ... 3.01741106e+12
+ 3.00451999e+12 2.99822435e+12]
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[2.29428933e-89 6.06790293e-89 1.07321699e-88 1.26616501e-88
+ 1.37851072e-88 2.08319104e-88 9.59706759e-89 5.39355744e-89
+ 2.69355837e-89 1.40033387e-89 6.15486308e-90 1.57060715e-90
+ 4.15146521e-91 1.37918322e-91 3.97154649e-92 1.05278818e-92
+ 3.35015072e-93 6.28090047e-94 1.03468105e-94 1.53463550e-95]
+ [3.37457263e-86 8.54697091e-86 1.47087381e-85 1.71436277e-85
+ 1.85052880e-85 2.74158063e-85 1.29259782e-85 7.38443501e-86
+ 3.76575926e-86 1.99691649e-86 9.00727877e-87 2.40599872e-87
+ 6.64921822e-88 2.29128084e-88 6.87952351e-89 1.90742246e-89
+ 6.30817772e-90 1.25250713e-90 2.19548372e-91 3.47810725e-92]
+ [3.70171706e-86 9.37040018e-86 1.61202029e-85 1.87858509e-85
+ 2.02757394e-85 3.00311920e-85 1.41632374e-85 8.09294711e-86
+ 4.12816699e-86 2.18964528e-86 9.87984623e-87 2.64060776e-87
+ 7.30170906e-88 2.51729355e-88 7.56212966e-89 2.09787782e-89
+ 6.94143983e-90 1.37924704e-90 2.41954807e-91 3.83627829e-92]
+ [6.81840483e-86 1.71974087e-85 2.95175068e-85 3.43635557e-85
+ 3.70622452e-85 5.48032089e-85 2.58963242e-85 1.48176221e-85
+ 7.57164219e-86 4.02278870e-86 1.81905501e-86 4.88048143e-87
+ 1.35457665e-87 4.68429078e-88 1.41212762e-88 3.93227791e-89
+ 1.30531526e-89 2.60612616e-90 4.59564651e-91 7.32688809e-92]
+ [8.76037092e-84 2.14675323e-83 3.61810988e-83 4.17817104e-83
+ 4.48061360e-83 6.53844017e-83 3.13763567e-83 1.81502902e-83
+ 9.40472471e-84 5.06306008e-84 2.32928809e-84 6.44263879e-85
+ 1.84199270e-85 6.52695393e-86 2.02313654e-86 5.80482465e-87
+ 1.97699946e-87 4.10095336e-88 7.53698365e-89 1.25553893e-89]
+ [1.12612447e-83 2.75548698e-83 4.63968461e-83 5.35563725e-83
+ 5.74161443e-83 8.37286122e-83 4.02113481e-83 2.32742083e-83
+ 1.20684205e-83 6.50150378e-84 2.99372197e-84 8.29346615e-85
+ 2.37479769e-85 8.42551902e-86 2.61538975e-86 7.51575235e-87
+ 2.56310547e-87 5.32725330e-88 9.81172128e-89 1.63818796e-89]
+ [1.21872703e-83 2.98067485e-83 5.01736651e-83 5.79083674e-83
+ 6.20760087e-83 9.05045175e-83 4.34764423e-83 2.51685067e-83
+ 1.30536349e-83 7.03376926e-84 3.23972172e-84 8.97940529e-85
+ 2.57245523e-85 9.13040684e-86 2.83548004e-86 8.15218725e-87
+ 2.78131129e-87 5.78437783e-88 1.06608267e-88 1.78122860e-89]
+ [2.28182362e-83 5.55997313e-83 9.33708076e-83 1.07652133e-82
+ 1.15314630e-82 1.67837747e-82 8.07864393e-83 4.68332573e-83
+ 2.43338093e-83 1.31343019e-83 6.06309560e-84 1.68710712e-84
+ 4.85184916e-85 1.72749474e-85 5.38411591e-86 1.55396275e-86
+ 5.31931905e-87 1.11175073e-87 2.05997576e-88 3.46141128e-89]
+ [2.46136989e-83 5.99476481e-83 1.00643804e-82 1.16022907e-82
+ 1.24270138e-82 1.80834999e-82 8.70634343e-83 5.04807339e-83
+ 2.62346822e-83 1.41632233e-83 6.53982928e-84 1.82062701e-84
+ 5.23825487e-85 1.86578322e-85 5.81764749e-86 1.67987258e-86
+ 5.75262003e-87 1.20302887e-87 2.23054498e-88 3.75058890e-89]
+ [8.51144594e-82 2.02982454e-81 3.36274482e-81 3.85376651e-81
+ 4.11051563e-81 5.92411566e-81 2.88446397e-81 1.68583746e-81
+ 8.85081116e-82 4.82449876e-82 2.25592864e-82 6.42145346e-83
+ 1.88800786e-83 6.84547451e-84 2.17826928e-84 6.42874136e-85
+ 2.24311800e-85 4.82368910e-86 9.21777902e-87 1.60039964e-87]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Planck j_blues: [[2.29428933e-89 6.06790293e-89 1.07321699e-88 ... 6.28090047e-94
+ 1.03468105e-94 1.53463550e-95]
+ [3.37457263e-86 8.54697091e-86 1.47087381e-85 ... 1.25250713e-90
+ 2.19548372e-91 3.47810725e-92]
+ [3.70171706e-86 9.37040018e-86 1.61202029e-85 ... 1.37924704e-90
+ 2.41954807e-91 3.83627829e-92]
+ ...
+ [1.41370335e-08 1.09757554e-08 8.94330146e-09 ... 2.14575993e-09
+ 2.05197282e-09 1.96744670e-09]
+ [1.40169230e-08 1.08825016e-08 8.86731515e-09 ... 2.12753187e-09
+ 2.03454199e-09 1.95073442e-09]
+ [1.39584493e-08 1.08371027e-08 8.83032254e-09 ... 2.11865784e-09
+ 2.02605607e-09 1.94259832e-09]]
+[J_BLUE_DEBUG] Final j_blues after zero handling: [[2.29428933e-099 6.06790293e-099 1.07321699e-098 ... 6.28090047e-104
+ 1.03468105e-104 1.53463550e-105]
+ [3.37457263e-096 8.54697091e-096 1.47087381e-095 ... 1.25250713e-100
+ 2.19548372e-101 3.47810725e-102]
+ [3.70171706e-096 9.37040018e-096 1.61202029e-095 ... 1.37924704e-100
+ 2.41954807e-101 3.83627829e-102]
+ ...
+ [1.41370335e-018 1.09757554e-018 8.94330146e-019 ... 2.14575993e-019
+ 2.05197282e-019 1.96744670e-019]
+ [1.40169230e-018 1.08825016e-018 8.86731515e-019 ... 2.12753187e-019
+ 2.03454199e-019 1.95073442e-019]
+ [1.39584493e-018 1.08371027e-018 8.83032254e-019 ... 2.11865784e-019
+ 2.02605607e-019 1.94259832e-019]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [10099.56644978 10122.19652418 10133.96735408 10135.12139173
+ 10133.79777637 10140.51569653 10116.72192596 10097.81400121
+ 10075.80319511 10054.67105355] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.44892062 0.3581946 0.30066847 0.25897885 0.22709752 0.20097692
+ 0.18160534 0.16511811 0.15127989 0.13946125]
+ Shape: (20,)
+ dtype: float64
+==================================================
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[1.31817623e-90 1.69463539e-90 1.82149631e-90 1.60738294e-90
+ 1.37090558e-90 1.39674791e-90 7.65696987e-91 4.67216486e-91
+ 2.68569814e-91 1.57953690e-91 8.52609218e-92 3.47247271e-92
+ 1.44559780e-92 6.79201860e-93 2.97209224e-93 1.24988238e-93
+ 5.78599261e-94 2.05334101e-94 6.86039616e-95 2.18717679e-95]
+ [2.13588477e-87 2.70022984e-87 2.87725960e-87 2.53688391e-87
+ 2.16577132e-87 2.19570092e-87 1.22500974e-87 7.58032120e-88
+ 4.42936178e-88 2.64650997e-88 1.45589443e-88 6.10544937e-89
+ 2.61532322e-89 1.25911646e-89 5.66051357e-90 2.44939373e-90
+ 1.16284140e-90 4.27246285e-91 1.48118401e-91 4.90836186e-92]
+ [2.34582522e-87 2.96501016e-87 3.15905084e-87 2.78530957e-87
+ 2.37788504e-87 2.41059448e-87 1.34520138e-87 8.32554384e-88
+ 4.86582446e-88 2.90787589e-88 1.60006170e-88 6.71252049e-89
+ 2.87640941e-89 1.38524185e-89 6.22965998e-90 2.69664892e-90
+ 1.28063476e-90 4.70732718e-91 1.63270860e-91 5.41314090e-92]
+ [4.35608680e-87 5.49815778e-87 5.85371387e-87 5.16080432e-87
+ 4.40626162e-87 4.46502072e-87 2.49531629e-87 1.54618276e-87
+ 9.04900144e-88 5.41495380e-88 2.98432134e-88 1.25504215e-88
+ 5.39091027e-89 2.60150158e-89 1.17258735e-89 5.08796033e-90
+ 2.42137568e-90 8.92635443e-91 3.10564897e-91 1.03299879e-91]
+ [5.96945582e-85 7.45085854e-85 7.88691131e-85 6.94939070e-85
+ 5.93720327e-85 5.99657234e-85 3.39067224e-85 2.12068156e-85
+ 1.25474142e-85 7.58783047e-86 4.23501039e-86 1.81604212e-86
+ 7.95039479e-87 3.89945641e-87 1.78951839e-87 7.91391787e-88
+ 3.83006240e-88 1.44497034e-88 5.15256032e-89 1.75855594e-89]
+ [7.69922383e-85 9.60434358e-85 1.01633832e-84 8.95499486e-85
+ 7.65094593e-85 7.72613363e-85 4.37127127e-85 2.73531365e-85
+ 1.61931389e-85 9.79785002e-86 5.47206729e-86 2.34887775e-86
+ 1.02932000e-86 5.05278296e-87 2.32095535e-87 1.02742296e-87
+ 4.97669444e-88 1.87980787e-88 6.71166661e-89 2.29372891e-89]
+ [8.34108886e-85 1.04031427e-84 1.10076409e-84 9.69878400e-85
+ 8.28651030e-85 8.36749471e-85 4.73504042e-85 2.96339172e-85
+ 1.75464843e-85 1.06185260e-85 5.93163170e-86 2.54695254e-86
+ 1.11646545e-86 5.48201644e-87 2.51885763e-87 1.11537397e-87
+ 5.40419457e-88 2.04205254e-88 7.29386577e-89 2.49374221e-89]
+ [1.57476409e-84 1.96124052e-84 2.07365171e-84 1.82695164e-84
+ 1.56105394e-84 1.57563883e-84 8.92979079e-85 5.59539110e-85
+ 3.31775048e-85 2.01051885e-85 1.12493346e-85 4.84246276e-86
+ 2.12793180e-86 1.04704170e-86 4.82210170e-87 2.14052328e-87
+ 1.03937670e-87 3.93918273e-88 1.41148899e-88 4.84190609e-89]
+ [1.70038462e-84 2.11732177e-84 2.23847682e-84 1.97215024e-84
+ 1.68513718e-84 1.70079387e-84 9.64085525e-85 6.04182283e-85
+ 3.58306923e-85 2.17165545e-85 1.21533275e-85 5.23319162e-86
+ 2.30031291e-86 1.13214809e-86 5.21551863e-87 2.31584731e-87
+ 1.12480381e-87 4.26448505e-88 1.52863815e-88 5.24587680e-89]
+ [6.16329149e-83 7.61224104e-83 8.01388835e-83 7.05750196e-83
+ 6.03326431e-83 6.07467960e-83 3.47293083e-83 2.19133266e-83
+ 1.30995079e-83 8.00066914e-84 4.51892807e-84 1.97370286e-84
+ 8.79693047e-85 4.38123447e-85 2.04499770e-85 9.20731076e-86
+ 4.52715823e-86 1.74560339e-86 6.37066044e-87 2.22773458e-87]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Initialized j_blue_estimator with shape: (29224, 20)
+[J_BLUE_DEBUG] Initial values: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [10925.93307493 10943.60134063 10982.73157655 11009.41071627
+ 11024.35730256 10994.26070361 10987.49523571 10980.48709583
+ 10949.05349676 10880.43462949] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.51203941 0.40161425 0.32636519 0.27164137 0.23084568 0.20309009
+ 0.17918147 0.15961881 0.14476607 0.13401689]
+ Shape: (20,)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Input j_blue_estimator: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] Normalization factor: [3.83584691e+13 3.54588448e+13 3.28759858e+13 3.05653796e+13
+ 2.84900666e+13 2.66191521e+13 2.49266505e+13 2.33905780e+13
+ 2.19922355e+13 2.07156364e+13 1.95470480e+13 1.84746209e+13
+ 1.74880869e+13 1.65785126e+13 1.57380967e+13 1.49600018e+13
+ 1.42382149e+13 1.35674314e+13 1.29429566e+13 1.23606241e+13]
+[J_BLUE_DEBUG] After normalization: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] line_list_nu: [4.48979300e+16 4.33201055e+16 4.33000835e+16 ... 3.01741106e+12
+ 3.00451999e+12 2.99822435e+12]
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[1.53151192e-83 1.65160368e-83 2.70687020e-83 3.62443933e-83
+ 4.01613181e-83 2.06920800e-83 1.61807344e-83 1.27183542e-83
+ 6.56665811e-84 1.75718743e-84 8.25923391e-85 2.69763909e-85
+ 8.56807207e-86 2.79690014e-86 8.82079056e-87 1.98210905e-87
+ 5.51803254e-88 1.11475993e-88 3.87350970e-89 8.25538918e-90]
+ [1.40749151e-80 1.50096910e-80 2.40008483e-80 3.16041221e-80
+ 3.46945285e-80 1.82147522e-80 1.43040550e-80 1.12928132e-80
+ 5.94721881e-81 1.66237908e-81 8.00115067e-82 2.71136699e-82
+ 8.94515846e-83 3.03070368e-83 9.93487017e-84 2.34876755e-84
+ 6.82774953e-85 1.45707075e-85 5.24579719e-86 1.17880334e-86]
+ [1.53475278e-80 1.63644996e-80 2.61590343e-80 3.44387013e-80
+ 3.78018134e-80 1.98508238e-80 1.55897020e-80 1.23084971e-80
+ 6.48374473e-81 1.81335331e-81 8.73042779e-82 2.95988141e-82
+ 9.76974950e-83 3.31164634e-83 1.08611474e-83 2.56941026e-84
+ 7.47324669e-85 1.59593231e-85 5.74831870e-86 1.29259515e-86]
+ [2.71773856e-80 2.89510875e-80 4.61835080e-80 6.07161318e-80
+ 6.65933239e-80 3.50252057e-80 2.75165682e-80 2.17330917e-80
+ 1.14673403e-80 3.21888971e-81 1.55282358e-81 5.28081057e-82
+ 1.74860617e-82 5.94574721e-83 1.95634119e-83 4.64782966e-84
+ 1.35674866e-84 2.91070946e-85 1.05151621e-85 2.37500315e-86]
+ [2.55270428e-78 2.69911150e-78 4.23556295e-78 5.50674478e-78
+ 6.00238639e-78 3.19678439e-78 2.51856842e-78 1.99504838e-78
+ 1.06664991e-78 3.08235598e-79 1.51063684e-79 5.26490573e-80
+ 1.78802222e-80 6.23235973e-81 2.10413307e-81 5.17096288e-82
+ 1.55356867e-82 3.45690118e-83 1.27869507e-83 2.99180571e-84]
+ [3.22869487e-78 3.41255677e-78 5.35058529e-78 6.95240528e-78
+ 7.57573090e-78 4.03733981e-78 3.18126006e-78 2.52037215e-78
+ 1.34843308e-78 3.90250051e-79 1.91414646e-79 6.67969604e-80
+ 2.27147283e-80 7.92764157e-81 2.68005152e-81 6.59783390e-82
+ 1.98521486e-82 4.42572351e-83 1.63906125e-83 3.84196987e-84]
+ [3.47642210e-78 3.67394538e-78 5.75887978e-78 7.48157656e-78
+ 8.15152139e-78 4.34508202e-78 3.42390589e-78 2.71273911e-78
+ 1.45166359e-78 4.20324727e-79 2.06219069e-79 7.19919189e-80
+ 2.44913940e-80 8.55116324e-81 2.89205375e-81 7.12366946e-82
+ 2.14443819e-82 4.78352901e-83 1.77225567e-83 4.15656420e-84]
+ [6.25075138e-78 6.59954879e-78 1.03228167e-77 1.33914901e-77
+ 1.45789400e-77 7.78373372e-78 6.13578742e-78 4.86318623e-78
+ 2.60686634e-78 7.57648511e-79 3.72475801e-79 1.30445509e-79
+ 4.45224238e-80 1.55948515e-80 5.29184457e-81 1.30918982e-81
+ 3.95574728e-82 8.86567877e-83 3.29469957e-83 7.76252290e-84]
+ [6.70971904e-78 7.08330348e-78 1.10766529e-77 1.43669249e-77
+ 1.56393528e-77 8.35152274e-78 6.58365620e-78 5.21840274e-78
+ 2.79785252e-78 8.13524610e-79 4.00044256e-79 1.40153939e-79
+ 4.78549056e-80 1.67686027e-80 5.69242331e-81 1.40903562e-81
+ 4.25934755e-82 9.55155203e-83 3.55089559e-83 8.37074252e-84]
+ [1.84593225e-76 1.93813634e-76 2.99468275e-76 3.85281102e-76
+ 4.17506653e-76 2.24999257e-76 1.77737127e-76 1.41181401e-76
+ 7.64267976e-77 2.26987734e-77 1.12914182e-77 4.02737273e-78
+ 1.40076861e-78 4.99799076e-79 1.72885702e-79 4.38641337e-80
+ 1.35413798e-80 3.11867078e-81 1.17957477e-81 2.85321118e-82]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Planck j_blues: [[1.53151192e-83 1.65160368e-83 2.70687020e-83 ... 1.11475993e-88
+ 3.87350970e-89 8.25538918e-90]
+ [1.40749151e-80 1.50096910e-80 2.40008483e-80 ... 1.45707075e-85
+ 5.24579719e-86 1.17880334e-86]
+ [1.53475278e-80 1.63644996e-80 2.61590343e-80 ... 1.59593231e-85
+ 5.74831870e-86 1.29259515e-86]
+ ...
+ [1.55460990e-08 1.22133146e-08 9.96067175e-09 ... 2.42643575e-09
+ 2.30249155e-09 2.19821365e-09]
+ [1.54139871e-08 1.21095244e-08 9.87602377e-09 ... 2.40581919e-09
+ 2.28292844e-09 2.17953701e-09]
+ [1.53496709e-08 1.20589960e-08 9.83481445e-09 ... 2.39578238e-09
+ 2.27340448e-09 2.17044461e-09]]
+[J_BLUE_DEBUG] Final j_blues after zero handling: [[1.53151192e-093 1.65160368e-093 2.70687020e-093 ... 1.11475993e-098
+ 3.87350970e-099 8.25538918e-100]
+ [1.40749151e-090 1.50096910e-090 2.40008483e-090 ... 1.45707075e-095
+ 5.24579719e-096 1.17880334e-096]
+ [1.53475278e-090 1.63644996e-090 2.61590343e-090 ... 1.59593231e-095
+ 5.74831870e-096 1.29259515e-096]
+ ...
+ [1.55460990e-018 1.22133146e-018 9.96067175e-019 ... 2.42643575e-019
+ 2.30249155e-019 2.19821365e-019]
+ [1.54139871e-018 1.21095244e-018 9.87602377e-019 ... 2.40581919e-019
+ 2.28292844e-019 2.17953701e-019]
+ [1.53496709e-018 1.20589960e-018 9.83481445e-019 ... 2.39578238e-019
+ 2.27340448e-019 2.17044461e-019]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [10512.74976236 10532.8989324 10558.34946532 10572.266054
+ 10579.07753947 10567.38820007 10552.10858083 10539.15054852
+ 10512.42834593 10467.55284152] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.48048001 0.37990443 0.31351683 0.26531011 0.2289716 0.2020335
+ 0.1803934 0.16236846 0.14802298 0.13673907]
+ Shape: (20,)
+ dtype: float64
+==================================================
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[6.18306755e-87 7.23582388e-87 9.77762772e-87 1.08241679e-86
+ 1.06515730e-86 7.50250307e-87 4.98624509e-87 3.49153748e-87
+ 1.89293695e-87 7.26219417e-88 3.63532083e-88 1.31311241e-88
+ 4.72011350e-89 1.82153296e-89 6.67812912e-90 2.00548660e-90
+ 7.07030885e-91 1.85660553e-91 6.32795339e-92 1.62728651e-92]
+ [7.46155088e-84 8.61249267e-84 1.14379527e-83 1.25432279e-83
+ 1.22864304e-83 8.72281829e-84 5.85774658e-84 4.13814333e-84
+ 2.28484587e-84 9.04065927e-85 4.62422152e-85 1.72699809e-85
+ 6.42049263e-86 2.55654018e-86 9.69010193e-87 3.03023141e-87
+ 1.10620897e-87 3.03986991e-88 1.07428400e-88 2.89348439e-89]
+ [8.16437470e-84 9.42207922e-84 1.25103852e-83 1.37176485e-83
+ 1.34360207e-83 9.53993469e-84 6.40731888e-84 4.52688970e-84
+ 2.50006872e-84 9.89612850e-85 5.06317194e-85 1.89173294e-85
+ 7.03593664e-86 2.80271388e-86 1.06276674e-86 3.32512916e-87
+ 1.21440083e-87 3.33910694e-88 1.18057593e-88 3.18164000e-89]
+ [1.47911706e-83 1.70500261e-83 2.26057278e-83 2.47676049e-83
+ 2.42497476e-83 1.72293950e-83 1.15818644e-83 8.18885266e-84
+ 4.52938853e-84 1.79753138e-84 9.21336454e-85 3.45199075e-85
+ 1.28752889e-85 5.14224920e-86 1.95534547e-86 6.13856715e-87
+ 2.24848040e-87 6.20598583e-88 2.20085450e-88 5.95432315e-89]
+ [1.66569057e-81 1.90252668e-81 2.49350997e-81 2.71484870e-81
+ 2.64993227e-81 1.89272737e-81 1.28114458e-81 9.11160984e-82
+ 5.10146682e-82 2.06664227e-82 1.07459339e-82 4.11670710e-83
+ 1.57028911e-83 6.40375114e-84 2.48963226e-84 8.02953405e-85
+ 3.01022506e-85 8.56381130e-86 3.11113592e-86 8.68059384e-87]
+ [2.12665440e-81 2.42787966e-81 3.18015432e-81 3.46131758e-81
+ 3.37801506e-81 2.41342247e-81 1.63417553e-81 1.16259300e-81
+ 6.51329808e-82 2.64139467e-82 1.37446848e-82 5.27157017e-83
+ 2.01313820e-83 8.21858407e-84 3.19886504e-84 1.03313453e-84
+ 3.87781648e-85 1.10493164e-85 4.01910144e-86 1.12318688e-86]
+ [2.29659859e-81 2.62150354e-81 3.43312712e-81 3.73627369e-81
+ 3.64617156e-81 2.60523050e-81 1.76425097e-81 1.25525197e-81
+ 7.03380245e-82 2.85343461e-82 1.48515211e-82 5.69814240e-83
+ 2.17683484e-83 8.88988921e-84 3.46140188e-84 1.11841665e-84
+ 4.19950582e-85 1.19718245e-85 4.35636589e-86 1.21805066e-86]
+ [4.22731741e-81 4.81964721e-81 6.30241569e-81 6.85335179e-81
+ 6.68542613e-81 4.78006780e-81 3.23993315e-81 2.30694021e-81
+ 1.29472687e-81 5.26635362e-82 2.74611550e-82 1.05664333e-82
+ 4.04836075e-83 1.65775508e-83 6.47321740e-84 2.09886650e-84
+ 7.90464663e-85 2.26226254e-85 8.25771339e-86 2.31809155e-86]
+ [4.55057678e-81 5.18745861e-81 6.78216298e-81 7.37431341e-81
+ 7.19327801e-81 5.14360423e-81 3.48671437e-81 2.48288432e-81
+ 1.39373672e-81 5.67089994e-82 2.95772697e-82 1.13846146e-82
+ 4.36336084e-83 1.78732542e-83 6.98157997e-84 2.26464991e-84
+ 8.53210311e-85 2.44299040e-85 8.92076178e-86 2.50542630e-86]
+ [1.42922819e-79 1.61837781e-79 2.09813731e-79 2.27087710e-79
+ 2.21016658e-79 1.58649158e-79 1.08088082e-79 7.73003011e-80
+ 4.37786515e-80 1.80823557e-80 9.53046769e-81 3.72839798e-81
+ 1.45256973e-81 6.04132765e-82 2.39835082e-82 7.93431843e-83
+ 3.04038004e-83 8.90001062e-84 3.30761422e-84 9.50099019e-85]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Initialized j_blue_estimator with shape: (29224, 20)
+[J_BLUE_DEBUG] Initial values: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [10961.30468434 11045.76181757 11119.89428668 11130.89377416
+ 11162.91276394 11138.01140074 11125.05312176 11097.91367833
+ 11023.38144177 10959.09272047] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.50051149 0.38390544 0.30827136 0.26014908 0.22024417 0.19383445
+ 0.17137491 0.15398179 0.14205648 0.13087607]
+ Shape: (20,)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Input j_blue_estimator: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] Normalization factor: [3.81033480e+13 3.52230091e+13 3.26573285e+13 3.03620902e+13
+ 2.83005800e+13 2.64421089e+13 2.47608641e+13 2.32350080e+13
+ 2.18459658e+13 2.05778573e+13 1.94170412e+13 1.83517468e+13
+ 1.73717741e+13 1.64682494e+13 1.56334231e+13 1.48605033e+13
+ 1.41435170e+13 1.34771948e+13 1.28568734e+13 1.22784140e+13]
+[J_BLUE_DEBUG] After normalization: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] line_list_nu: [4.48979300e+16 4.33201055e+16 4.33000835e+16 ... 3.01741106e+12
+ 3.00451999e+12 2.99822435e+12]
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[2.82890217e-83 9.75440078e-83 2.87548166e-82 2.93874402e-82
+ 4.33502666e-82 2.47796973e-82 1.74883446e-82 9.78515502e-83
+ 2.42902893e-83 7.10939399e-84 2.85817800e-84 5.88575429e-85
+ 1.95725242e-85 7.39328096e-86 1.89039016e-86 3.52894281e-87
+ 8.92390984e-88 3.05484914e-88 8.81359008e-89 1.20562856e-89]
+ [2.54232116e-80 8.31520647e-80 2.34171767e-79 2.37718603e-79
+ 3.43889290e-79 1.99576480e-79 1.41971520e-79 8.07698280e-80
+ 2.09966399e-80 6.39809412e-81 2.64875435e-81 5.75323625e-82
+ 1.98428535e-82 7.73892445e-83 2.07210255e-83 4.09683384e-84
+ 1.08543342e-84 3.85151069e-85 1.15909041e-85 1.69845016e-86]
+ [2.77140377e-80 9.05839654e-80 2.54953486e-79 2.58792988e-79
+ 3.74283297e-79 2.17257479e-79 1.54564677e-79 8.79528483e-80
+ 2.28773015e-80 6.97473483e-81 2.88855308e-81 6.27833869e-82
+ 2.16639554e-82 8.45260247e-83 2.26451438e-83 4.48052306e-84
+ 1.18779215e-84 4.21663724e-85 1.26965276e-85 1.86208522e-86]
+ [4.89841263e-80 1.59398829e-79 4.46922249e-79 4.53397056e-79
+ 6.54661344e-79 3.80489572e-79 2.70873416e-79 1.54351777e-79
+ 4.03037054e-80 1.23291743e-80 5.11862658e-81 1.11752650e-81
+ 3.86792833e-82 1.51319325e-82 4.06963892e-83 8.09094315e-84
+ 2.15331128e-84 7.66723742e-85 2.31682175e-85 3.41751107e-86]
+ [4.53292425e-78 1.42405891e-77 3.87305615e-77 3.91159458e-77
+ 5.57502712e-77 3.27310755e-77 2.34247367e-77 1.34969338e-77
+ 3.63424444e-78 1.14198513e-78 4.83461855e-79 1.09368392e-79
+ 3.87853908e-80 1.55000690e-80 4.29860910e-81 8.87945430e-82
+ 2.43769058e-82 8.88990915e-83 2.76284164e-83 4.26651021e-84]
+ [5.72888749e-78 1.79650889e-77 4.87833158e-77 4.92573066e-77
+ 7.01571300e-77 4.12108951e-77 2.95015605e-77 1.70080477e-77
+ 4.58694808e-78 1.44335495e-78 6.11665072e-79 1.38624888e-79
+ 4.92224899e-80 1.96928010e-80 5.47005250e-81 1.13216383e-81
+ 3.11314328e-82 1.13672427e-82 3.53789662e-83 5.47635118e-84]
+ [6.16695074e-78 1.93277290e-77 5.24574969e-77 5.29633224e-77
+ 7.54196457e-77 4.43094272e-77 3.17224240e-77 1.82917065e-77
+ 4.93560984e-78 1.55374555e-78 6.58655780e-79 1.49360969e-79
+ 5.30556107e-80 2.12337057e-80 5.90101645e-81 1.22212341e-81
+ 3.36220661e-82 1.22814457e-82 3.82417862e-83 5.92390714e-84]
+ [1.10671138e-77 3.45279399e-77 9.33446387e-77 9.41901622e-77
+ 1.33901534e-76 7.87705700e-77 5.64326141e-77 3.25866785e-77
+ 8.82776618e-78 2.78866217e-78 1.18514188e-78 2.69986311e-79
+ 9.62054424e-80 3.86090448e-80 1.07724022e-80 2.24206035e-81
+ 6.19296624e-82 2.26916282e-82 7.09138549e-83 1.10502470e-83]
+ [1.18769668e-77 3.70342384e-77 1.00072753e-76 1.00972158e-76
+ 1.43513769e-76 8.44384901e-77 6.04981901e-77 3.49403652e-77
+ 9.46992182e-78 2.99277015e-78 1.27227253e-78 2.89996146e-79
+ 1.03374833e-79 4.15000437e-80 1.15845744e-80 2.41253745e-81
+ 6.66708320e-82 2.44379562e-82 7.64048183e-83 1.19143985e-83]
+ [3.23216867e-76 9.82284309e-76 2.59596133e-75 2.61073579e-75
+ 3.67564941e-75 2.17863097e-75 1.56695998e-75 9.12342807e-76
+ 2.52882880e-76 8.14997534e-77 3.51443260e-77 8.22105975e-78
+ 2.98302001e-78 1.21630383e-78 3.47221542e-79 7.43583471e-80
+ 2.10201303e-80 7.84055356e-81 2.50214391e-81 4.03450639e-82]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Planck j_blues: [[2.82890217e-83 9.75440078e-83 2.87548166e-82 ... 3.05484914e-88
+ 8.81359008e-89 1.20562856e-89]
+ [2.54232116e-80 8.31520647e-80 2.34171767e-79 ... 3.85151069e-85
+ 1.15909041e-85 1.69845016e-86]
+ [2.77140377e-80 9.05839654e-80 2.54953486e-79 ... 4.21663724e-85
+ 1.26965276e-85 1.86208522e-86]
+ ...
+ [1.52456208e-08 1.17844886e-08 9.52672553e-09 ... 2.39547303e-09
+ 2.28250174e-09 2.19012662e-09]
+ [1.51160609e-08 1.16843396e-08 9.44576202e-09 ... 2.37511920e-09
+ 2.26310820e-09 2.17151857e-09]
+ [1.50529871e-08 1.16355838e-08 9.40634644e-09 ... 2.36521030e-09
+ 2.25366679e-09 2.16245956e-09]]
+[J_BLUE_DEBUG] Final j_blues after zero handling: [[2.82890217e-93 9.75440078e-93 2.87548166e-92 ... 3.05484914e-98
+ 8.81359008e-99 1.20562856e-99]
+ [2.54232116e-90 8.31520647e-90 2.34171767e-89 ... 3.85151069e-95
+ 1.15909041e-95 1.69845016e-96]
+ [2.77140377e-90 9.05839654e-90 2.54953486e-89 ... 4.21663724e-95
+ 1.26965276e-95 1.86208522e-96]
+ ...
+ [1.52456208e-18 1.17844886e-18 9.52672553e-19 ... 2.39547303e-19
+ 2.28250174e-19 2.19012662e-19]
+ [1.51160609e-18 1.16843396e-18 9.44576202e-19 ... 2.37511920e-19
+ 2.26310820e-19 2.17151857e-19]
+ [1.50529871e-18 1.16355838e-18 9.40634644e-19 ... 2.36521030e-19
+ 2.25366679e-19 2.16245956e-19]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [10737.02722335 10789.33037499 10839.121876 10851.57991408
+ 10870.9951517 10852.69980041 10838.58085129 10818.53211342
+ 10767.90489385 10713.322781 ] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.49049575 0.38190493 0.3108941 0.26272959 0.22460789 0.19793398
+ 0.17588416 0.15817513 0.14503973 0.13380757]
+ Shape: (20,)
+ dtype: float64
+==================================================
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[4.56610647e-85 9.40528356e-85 1.91627247e-84 2.03455898e-84
+ 2.47971392e-84 1.56451722e-84 1.07340085e-84 6.67836010e-85
+ 2.40062578e-85 7.99002290e-86 3.56512841e-86 9.59901909e-87
+ 3.31180415e-87 1.26352247e-87 3.83885570e-88 9.00282323e-89
+ 2.67131751e-89 8.04461131e-90 2.51423264e-90 4.66405884e-91]
+ [4.74054042e-82 9.43638848e-82 1.86161310e-81 1.96073658e-81
+ 2.36014164e-81 1.50666462e-81 1.04314708e-81 6.57470648e-82
+ 2.44243628e-82 8.42571646e-83 3.85713571e-83 1.08502028e-83
+ 3.87742948e-84 1.52692146e-84 4.82825063e-85 1.18952915e-85
+ 3.67714470e-86 1.15315677e-86 3.74853169e-87 7.36842603e-88]
+ [5.17717083e-82 1.03010629e-81 2.03136505e-81 2.13930940e-81
+ 2.57468282e-81 1.64386812e-81 1.13827191e-81 7.17543433e-82
+ 2.66671394e-82 9.20359458e-83 4.21460431e-83 1.18623604e-83
+ 4.24102581e-84 1.67077607e-84 5.28581029e-85 1.30307232e-85
+ 4.03023049e-86 1.26453492e-86 4.11263632e-87 8.09008446e-88]
+ [9.26185724e-82 1.83756947e-81 3.61390498e-81 3.80338732e-81
+ 4.57264270e-81 2.92238701e-81 2.02510460e-81 1.27797008e-81
+ 4.76262330e-82 1.64865885e-82 7.56593492e-83 2.13732410e-83
+ 7.66388996e-84 3.02725391e-84 9.60939405e-85 2.37873686e-85
+ 7.38240721e-86 2.32419837e-86 7.58390504e-87 1.49910811e-87]
+ [9.43550100e-80 1.82987120e-79 3.52230513e-79 3.68723363e-79
+ 4.39634877e-79 2.83178011e-79 1.97423382e-79 1.25665834e-79
+ 4.78699270e-80 1.69711711e-80 7.92242027e-81 2.30428207e-81
+ 8.45833036e-82 3.41229073e-82 1.11237905e-82 2.84553029e-83
+ 9.07536043e-84 2.93535805e-84 9.83224428e-85 2.01998599e-85]
+ [1.19843949e-79 2.32145395e-79 4.46358785e-79 4.67129991e-79
+ 5.56727589e-79 3.58744992e-79 2.50184797e-79 1.59321097e-79
+ 6.07590934e-80 2.15673339e-80 1.00768792e-80 2.93534536e-81
+ 1.07878348e-81 4.35681971e-82 1.42224468e-82 3.64436920e-83
+ 1.16395439e-83 3.76998411e-84 1.26450147e-84 2.60304554e-85]
+ [1.29209901e-79 2.50195082e-79 4.80895780e-79 5.03230394e-79
+ 5.99671171e-79 3.86466246e-79 2.69543856e-79 1.71673308e-79
+ 6.54931289e-80 2.32567823e-80 1.08692582e-80 3.16766538e-81
+ 1.16460842e-81 4.70505130e-82 1.53658736e-82 3.93946666e-83
+ 1.25876281e-83 4.07885501e-84 1.36868402e-84 2.81928103e-85]
+ [2.34775392e-79 4.53270904e-79 8.68810672e-79 9.08534252e-79
+ 1.08148837e-78 6.97684540e-79 4.86986196e-79 3.10508762e-79
+ 1.18794603e-79 4.23145364e-80 1.98197291e-80 5.79794301e-81
+ 2.13810215e-81 8.66156762e-82 2.83846002e-82 7.30811220e-83
+ 2.34337862e-83 7.61993371e-84 2.56557819e-84 5.31110748e-85]
+ [2.52333584e-79 4.86996691e-79 9.33142128e-79 9.75725721e-79
+ 1.16132039e-78 7.49276821e-79 5.23047192e-79 3.33546603e-79
+ 1.27652085e-79 4.54864979e-80 2.13111218e-80 6.23706436e-81
+ 2.30087707e-81 9.32404614e-82 3.05682838e-82 7.87437288e-83
+ 2.52602771e-83 8.21731124e-84 2.76784159e-84 5.73327138e-85]
+ [7.36608535e-78 1.39819534e-77 2.63743174e-77 2.74705570e-77
+ 3.24982460e-77 2.10877356e-77 1.47859063e-77 9.48850441e-78
+ 3.68993582e-78 1.33794892e-78 6.34710673e-79 1.89757468e-79
+ 7.12092351e-80 2.93045788e-80 9.79581827e-81 2.58461306e-81
+ 8.45799773e-82 2.80618517e-82 9.63450399e-83 2.05268514e-83]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Initialized j_blue_estimator with shape: (29224, 20)
+[J_BLUE_DEBUG] Initial values: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [11010.21256748 11106.83751853 11175.59248275 11207.1147306
+ 11217.45325877 11192.39937563 11149.09994681 11139.59575014
+ 11080.81018853 11024.55353733] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.48615902 0.37128549 0.29979087 0.25038224 0.21468554 0.18848197
+ 0.16851345 0.15070096 0.1380795 0.12734655]
+ Shape: (20,)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Input j_blue_estimator: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] Normalization factor: [3.77672353e+13 3.49123041e+13 3.23692556e+13 3.00942637e+13
+ 2.80509383e+13 2.62088609e+13 2.45424465e+13 2.30300501e+13
+ 2.16532607e+13 2.03963383e+13 1.92457619e+13 1.81898645e+13
+ 1.72185363e+13 1.63229817e+13 1.54955194e+13 1.47294175e+13
+ 1.40187559e+13 1.33583114e+13 1.27434619e+13 1.21701051e+13]
+[J_BLUE_DEBUG] After normalization: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+[J_BLUE_DEBUG] line_list_nu: [4.48979300e+16 4.33201055e+16 4.33000835e+16 ... 3.01741106e+12
+ 3.00451999e+12 2.99822435e+12]
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[6.57981450e-83 2.75775324e-82 7.34547465e-82 1.05518726e-81
+ 1.08015696e-81 6.16896527e-82 2.61133103e-82 1.98029890e-82
+ 6.50240803e-83 2.22322967e-83 4.04133868e-84 1.85860796e-84
+ 5.96300914e-85 1.47765031e-85 3.09493294e-86 9.33887301e-87
+ 1.87485379e-87 5.35974703e-88 1.07290837e-88 1.46583760e-89]
+ [5.73454547e-80 2.26389343e-79 5.78234762e-79 8.14961243e-79
+ 8.29067303e-79 4.80704077e-79 2.08900127e-79 1.59339763e-79
+ 5.42412359e-80 1.92036713e-80 3.69829698e-81 1.74324933e-81
+ 5.80739366e-82 1.50834760e-82 3.33221553e-83 1.04682106e-83
+ 2.22025340e-84 6.62175609e-85 1.40073608e-85 2.05019022e-86]
+ [6.24883815e-80 2.46505446e-79 6.29279490e-79 8.86688969e-79
+ 9.01965276e-79 5.23071566e-79 2.27387621e-79 1.73453958e-79
+ 5.90729104e-80 2.09235435e-80 4.03246833e-81 1.90136024e-81
+ 6.33714388e-82 1.64692106e-82 3.64081174e-83 1.14435193e-83
+ 2.42880384e-84 7.24763787e-85 1.53420617e-85 2.24750762e-86]
+ [1.10163666e-79 4.32403091e-79 1.09996739e-78 1.54744120e-78
+ 1.57328029e-78 9.13539152e-79 3.98005279e-79 3.03750539e-79
+ 1.03760789e-79 3.68594495e-80 7.13814437e-81 3.37267693e-81
+ 1.12764682e-81 2.94213410e-82 6.53322224e-83 2.06041417e-83
+ 4.39325039e-84 1.31562133e-84 2.79785858e-85 4.12239172e-86]
+ [9.98813326e-78 3.76726687e-77 9.31930752e-77 1.29450808e-76
+ 1.31067574e-76 7.68754404e-77 3.40839799e-77 2.61128848e-77
+ 9.13699554e-78 3.32205435e-78 6.68569616e-79 3.21114610e-79
+ 1.10089309e-79 2.96368590e-80 6.81896786e-81 2.20902101e-81
+ 4.88565087e-82 1.50494068e-82 3.32029815e-83 5.12184679e-84]
+ [1.26100545e-77 4.74639305e-77 1.17244723e-76 1.62753089e-76
+ 1.64750427e-76 9.66818657e-77 4.29043069e-77 3.28769980e-77
+ 1.15180871e-77 4.19281301e-78 8.45491855e-79 4.06435251e-79
+ 1.39521004e-79 3.76209709e-80 8.67189856e-81 2.81318461e-81
+ 6.23365597e-82 1.92297461e-82 4.25066596e-83 6.57260004e-84]
+ [1.35697751e-77 5.10431705e-77 1.26028781e-76 1.74910519e-76
+ 1.77045101e-76 1.03913876e-76 4.61267800e-77 3.53485571e-77
+ 1.23888133e-77 4.51148016e-78 9.10321670e-79 4.37716365e-79
+ 1.50320504e-79 4.05536519e-80 9.35330695e-81 3.03556079e-81
+ 6.73042060e-82 2.07717153e-82 4.59426062e-83 7.10919143e-84]
+ [2.42878855e-77 9.08905296e-77 2.23605851e-76 3.09825460e-76
+ 3.13438514e-76 1.84207249e-76 8.19536477e-77 6.28352583e-77
+ 2.20906848e-77 8.06866637e-78 1.63619667e-78 7.88412967e-79
+ 2.71634953e-79 7.35790340e-80 1.70483206e-80 5.55214688e-81
+ 1.23685047e-81 3.83115912e-82 8.51404024e-83 1.32530098e-83]
+ [2.60568771e-77 9.74498631e-77 2.39638418e-76 3.31974222e-76
+ 3.35823830e-76 1.97394056e-76 8.78444321e-77 6.73558795e-77
+ 2.36888544e-77 8.65553734e-78 1.75625833e-78 8.46482156e-79
+ 2.91755844e-79 7.90678908e-80 1.83302460e-80 5.97213301e-81
+ 1.33117045e-81 4.12513176e-82 9.17259879e-83 1.42883513e-83]
+ [6.98603036e-76 2.53778098e-75 6.11465033e-75 8.39256786e-75
+ 8.46421988e-75 5.01187891e-75 2.25906916e-75 1.73705883e-75
+ 6.21725606e-76 2.31053653e-76 4.82170683e-77 2.35196462e-77
+ 8.25619445e-78 2.28921656e-78 5.44644107e-79 1.80959569e-79
+ 4.14272770e-80 1.31049262e-80 2.99324406e-81 4.82145316e-82]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Planck j_blues: [[6.57981450e-83 2.75775324e-82 7.34547465e-82 ... 5.35974703e-88
+ 1.07290837e-88 1.46583760e-89]
+ [5.73454547e-80 2.26389343e-79 5.78234762e-79 ... 6.62175609e-85
+ 1.40073608e-85 2.05019022e-86]
+ [6.24883815e-80 2.46505446e-79 6.29279490e-79 ... 7.24763787e-85
+ 1.53420617e-85 2.24750762e-86]
+ ...
+ [1.48749542e-08 1.14605348e-08 9.31135538e-09 ... 2.36596570e-09
+ 2.25960551e-09 2.17108297e-09]
+ [1.47485425e-08 1.13631371e-08 9.23222091e-09 ... 2.34586240e-09
+ 2.24040644e-09 2.15263665e-09]
+ [1.46870013e-08 1.13157208e-08 9.19369578e-09 ... 2.33607546e-09
+ 2.23105971e-09 2.14365638e-09]]
+[J_BLUE_DEBUG] Final j_blues after zero handling: [[6.57981450e-93 2.75775324e-92 7.34547465e-92 ... 5.35974703e-98
+ 1.07290837e-98 1.46583760e-99]
+ [5.73454547e-90 2.26389343e-89 5.78234762e-89 ... 6.62175609e-95
+ 1.40073608e-95 2.05019022e-96]
+ [6.24883815e-90 2.46505446e-89 6.29279490e-89 ... 7.24763787e-95
+ 1.53420617e-95 2.24750762e-96]
+ ...
+ [1.48749542e-18 1.14605348e-18 9.31135538e-19 ... 2.36596570e-19
+ 2.25960551e-19 2.17108297e-19]
+ [1.47485425e-18 1.13631371e-18 9.23222091e-19 ... 2.34586240e-19
+ 2.24040644e-19 2.15263665e-19]
+ [1.46870013e-18 1.13157208e-18 9.19369578e-19 ... 2.33607546e-19
+ 2.23105971e-19 2.14365638e-19]]
+
+=== DilutePlanckianRadiationField.__init__ ===
+Temperature array:
+ First 10 elements: [10873.61989542 10948.08394676 11007.35717937 11029.34732234
+ 11044.22420524 11022.54958802 10993.84039905 10979.06393178
+ 10924.35754119 10868.93815916] K
+ Shape: (20,)
+ dtype: float64
+ unit: K
+
+Dilution factor array:
+ First 10 elements: [0.48832738 0.37659521 0.30534248 0.25655592 0.21964671 0.19320798
+ 0.1721988 0.15443804 0.14155961 0.13057706]
+ Shape: (20,)
+ dtype: float64
+==================================================
+
+=== DilutePlanckianRadiationField.calculate_mean_intensity ===
+Input frequency (nu):
+ First 10 elements: [4.48979300e+16 4.33201055e+16 4.33000835e+16 4.31679038e+16
+ 4.21169214e+16 4.20625564e+16 4.20454487e+16 4.19096722e+16
+ 4.18932740e+16 4.11260505e+16]
+ Shape: (29224,)
+ dtype: float64
+
+Output intensity:
+ First 10 elements: [[5.65544479e-84 1.67873990e-83 3.92802737e-83 4.87610187e-83
+ 5.43132758e-83 3.25529242e-83 1.74137154e-83 1.19964831e-83
+ 4.11532190e-84 1.38851965e-84 3.90585838e-85 1.38383144e-85
+ 4.59909504e-86 1.40615093e-86 3.53104016e-87 9.41887540e-88
+ 2.28841091e-88 6.71004354e-89 1.67081583e-89 2.65233875e-90]
+ [5.37369514e-81 1.52131419e-80 3.42953057e-80 4.19929294e-80
+ 4.63439392e-80 2.81534917e-80 1.53329217e-80 1.06613794e-80
+ 3.78584921e-81 1.32330831e-81 3.88253574e-82 1.42310806e-82
+ 4.90512500e-83 1.56019632e-83 4.10555193e-84 1.14519404e-84
+ 2.91948022e-85 8.92283475e-86 2.32966146e-86 3.94067958e-87]
+ [5.86204875e-81 1.65857183e-80 3.73718662e-80 4.57520639e-80
+ 5.04866433e-80 3.06753918e-80 1.67101996e-80 1.16204035e-80
+ 4.12820700e-81 1.44362360e-81 4.23780030e-82 1.55399690e-82
+ 5.35874557e-83 1.70533687e-83 4.49014491e-84 1.25318216e-84
+ 3.19672863e-85 9.77533325e-86 2.55377584e-86 4.32325683e-87]
+ [1.04095459e-80 2.93354909e-80 6.58944443e-80 8.05778517e-80
+ 8.88474606e-80 5.40442084e-80 2.94844734e-80 2.05196627e-80
+ 7.31083664e-81 2.56416204e-81 7.55377760e-82 2.77785768e-82
+ 9.60834365e-83 3.06784865e-83 8.10935873e-84 2.27178118e-84
+ 5.81846217e-85 1.78542840e-85 4.68293641e-86 7.96997240e-87]
+ [9.99701291e-79 2.72979313e-78 5.98151028e-78 7.24786269e-78
+ 7.94262300e-78 4.87492893e-78 2.69154844e-78 1.88478029e-78
+ 6.87146063e-79 2.46746893e-79 7.47572450e-80 2.81208087e-80
+ 9.96565811e-81 3.26684099e-81 8.90874856e-82 2.57114377e-82
+ 6.79960717e-83 2.14492326e-83 5.80631121e-84 1.03088066e-84]
+ [1.26588921e-78 3.45101146e-78 7.55214268e-78 9.14669180e-78
+ 1.00202752e-77 6.15298329e-78 3.39928908e-78 2.38114229e-78
+ 8.69141558e-79 3.12479855e-79 9.48099271e-80 3.57056339e-80
+ 1.26695188e-80 4.15885445e-81 1.13595891e-81 3.28353088e-82
+ 8.69797984e-83 2.74768372e-83 7.45014321e-84 1.32563181e-84]
+ [1.36350962e-78 3.71523141e-78 8.12707510e-78 9.84155096e-78
+ 1.07804180e-77 6.62071927e-78 3.65840695e-78 2.56290740e-78
+ 9.35838100e-79 3.36588024e-79 1.02171262e-79 3.84921070e-80
+ 1.36636465e-80 4.48710655e-81 1.22624040e-81 3.54621099e-82
+ 9.39871286e-83 2.97037963e-83 8.05810740e-84 1.43479669e-84]
+ [2.45869156e-78 6.67208452e-78 1.45484929e-77 1.75968398e-77
+ 1.92602140e-77 1.18422588e-77 6.55377981e-78 4.59493312e-78
+ 1.68282053e-78 6.07094932e-79 1.84952673e-79 6.98832849e-80
+ 2.48845540e-80 8.19987294e-81 2.24990973e-81 6.53167293e-82
+ 1.73830515e-82 5.51339439e-83 1.50179862e-83 2.68869894e-84]
+ [2.64013818e-78 7.16094532e-78 1.56084119e-77 1.88761527e-77
+ 2.06584706e-77 1.27037653e-77 7.03186836e-78 4.93060171e-78
+ 1.80640207e-78 6.51917674e-79 1.98694957e-79 7.51022416e-80
+ 2.67530867e-80 8.81920831e-81 2.42102220e-81 7.03169211e-82
+ 1.87231344e-82 5.94098881e-83 1.61906891e-83 2.90056379e-84]
+ [7.38209200e-77 1.95668332e-76 4.18835593e-76 5.03155059e-76
+ 5.48193032e-76 3.39324309e-76 1.89470567e-76 1.33453078e-76
+ 4.97206647e-77 1.82548815e-77 5.67893209e-78 2.18226634e-78
+ 7.91266452e-79 2.65905458e-79 7.46754543e-80 2.21655080e-80
+ 6.04164182e-81 1.95609812e-81 5.45516144e-82 1.00793972e-82]]
+ Shape: (29224, 20)
+ dtype: float64
+==================================================
+[J_BLUE_DEBUG] Initialized j_blue_estimator with shape: (29224, 20)
+[J_BLUE_DEBUG] Initial values: [[0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ ...
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]
+ [0. 0. 0. ... 0. 0. 0.]]
+FAILED
+
+=================================== FAILURES ===================================
+__________________ TestTransportSimple.test_j_blue_estimators __________________
+
+self =
+simulation =
+
+ def test_j_blue_estimators(self, simulation):
+ key = "simulation/transport/transport_state/j_blue_estimator"
+ expected = self.get_expected_data(key)
+
+ # np.save("j_blue_estimator.npy", simulation.transport.transport_state.radfield_mc_estimators.j_blue_estimator)
+
+
+> npt.assert_allclose(
+ simulation.transport.transport_state.radfield_mc_estimators.j_blue_estimator,
+ expected.values,
+ )
+
+tardis/tests/test_tardis_full.py:77:
+_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
+
+args = (.compare at 0x7f2fb77d5760>, array([[0., 0., 0., ..., 0., 0., 0.],
+ [0., 0., .....,
+ [0., 0., 0., ..., 0., 0., 0.],
+ [0., 0., 0., ..., 0., 0., 0.],
+ [0., 0., 0., ..., 0., 0., 0.]]))
+kwds = {'equal_nan': True, 'err_msg': '', 'header': 'Not equal to tolerance rtol=1e-07, atol=0', 'strict': False, ...}
+
+ @wraps(func)
+ def inner(*args, **kwds):
+ with self._recreate_cm():
+> return func(*args, **kwds)
+E AssertionError:
+E Not equal to tolerance rtol=1e-07, atol=0
+E
+E Mismatched elements: 490652 / 584480 (83.9%)
+E Max absolute difference among violations: 5.32418652e-19
+E Max relative difference among violations: 8.19556014
+E ACTUAL: array([[0., 0., 0., ..., 0., 0., 0.],
+E [0., 0., 0., ..., 0., 0., 0.],
+E [0., 0., 0., ..., 0., 0., 0.],...
+E DESIRED: array([[0., 0., 0., ..., 0., 0., 0.],
+E [0., 0., 0., ..., 0., 0., 0.],
+E [0., 0., 0., ..., 0., 0., 0.],...
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/contextlib.py:81: AssertionError
+------------------------------ Captured log setup ------------------------------
+WARNING tardis.io.model.parse_density_configuration:parse_density_configuration.py:114 Number of density points larger than number of shells. Assuming inner point irrelevant
+=============================== warnings summary ===============================
+tardis/io/configuration/config_validator.py:6
+tardis/io/configuration/config_validator.py:6
+ /home/aryaatharva18/tardis-main/tardis/tardis/io/configuration/config_validator.py:6: DeprecationWarning:
+
+ jsonschema.RefResolver is deprecated as of v4.18.0, in favor of the https://github.com/python-jsonschema/referencing library, which provides more compliant referencing behavior as well as more flexible APIs for customization. A future release will remove RefResolver. Please file a feature request (on referencing) if you are missing an API for the kind of customization you need.
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/jupyter_client/connect.py:22
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/jupyter_client/connect.py:22: DeprecationWarning:
+
+ Jupyter is migrating its paths to use standard platformdirs
+ given by the platformdirs library. To remove this warning and
+ see the appropriate new directories, set the environment variable
+ `JUPYTER_PLATFORM_DIRS=1` and then run `jupyter --paths`.
+ The use of platformdirs will be the default in `jupyter_core` v6
+
+tardis/plasma/properties/__init__.py:17
+ /home/aryaatharva18/tardis-main/tardis/tardis/plasma/properties/__init__.py:17: DeprecationWarning:
+
+ The numpy.linalg.linalg has been made private and renamed to numpy.linalg._linalg. All public functions exported by it are available from numpy.linalg. Please use numpy.linalg.LinAlgError instead.
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2599
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2599: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int8.
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2598
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2598: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int8.
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2599
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2599: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int16.
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2598
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/indexes/multi.py:2598: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int16.
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/sorting.py:349
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/sorting.py:349: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int8.
+
+../../miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/sorting.py:349
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/pandas/core/sorting.py:349: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int16.
+
+tardis/tests/test_tardis_full.py::TestTransportSimple::test_j_blue_estimators
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/numexpr/expressions.py:285: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int32.
+
+tardis/tests/test_tardis_full.py::TestTransportSimple::test_j_blue_estimators
+ /home/aryaatharva18/miniforge3/envs/tardis_numpy/lib/python3.12/site-packages/numexpr/expressions.py:295: UserWarning:
+
+ result dtype changed due to the removal of value-based promotion from NumPy. Changed from int64 to int32.
+
+tardis/tests/test_tardis_full.py::TestTransportSimple::test_j_blue_estimators
+ /home/aryaatharva18/tardis-main/tardis/tardis/transport/montecarlo/montecarlo_main_loop.py:123: NumbaTypeSafetyWarning:
+
+ [1m[1m[1munsafe cast from uint64 to int64. Precision may be lost.[0m[0m[0m
+
+-- Docs: https://docs.pytest.org/en/stable/how-to/capture-warnings.html
+=========================== short test summary info ============================
+FAILED tardis/tests/test_tardis_full.py::TestTransportSimple::test_j_blue_estimators - AssertionError:
+Not equal to tolerance rtol=1e-07, atol=0
+
+Mismatched elements: 490652 / 584480 (83.9%)
+Max absolute difference among violations: 5.32418652e-19
+Max relative difference among violations: 8.19556014
+ ACTUAL: array([[0., 0., 0., ..., 0., 0., 0.],
+ [0., 0., 0., ..., 0., 0., 0.],
+ [0., 0., 0., ..., 0., 0., 0.],...
+ DESIRED: array([[0., 0., 0., ..., 0., 0., 0.],
+ [0., 0., 0., ..., 0., 0., 0.],
+ [0., 0., 0., ..., 0., 0., 0.],...
+======================= 1 failed, 13 warnings in 55.32s ========================
diff --git a/tardis/plasma/radiation_field/planck_rad_field.py b/tardis/plasma/radiation_field/planck_rad_field.py
index ab3a4b741db..5d16a592405 100644
--- a/tardis/plasma/radiation_field/planck_rad_field.py
+++ b/tardis/plasma/radiation_field/planck_rad_field.py
@@ -27,6 +27,24 @@ def __init__(
dilution_factor: np.ndarray,
geometry=None,
):
+ print("\n=== DilutePlanckianRadiationField.__init__ ===")
+ print("Temperature array:")
+ print(f" First 10 elements: {temperature[:10]}")
+ print(f" Shape: {temperature.shape}")
+ print(f" dtype: {temperature.dtype}")
+ print(f" unit: {temperature.unit}")
+
+ print("\nDilution factor array:")
+ print(f" First 10 elements: {dilution_factor[:10]}")
+ print(f" Shape: {dilution_factor.shape}")
+ print(f" dtype: {dilution_factor.dtype}")
+
+ if geometry is not None:
+ print("\nGeometry boundaries:")
+ print(f" v_inner_boundary_index: {geometry.v_inner_boundary_index}")
+ print(f" v_outer_boundary_index: {geometry.v_outer_boundary_index}")
+ print("="*50)
+
# ensuring that the radiation_field has both
# dilution_factor and t_radiative equal length
assert len(temperature) == len(dilution_factor)
@@ -53,7 +71,13 @@ def __init__(
@property
def temperature_kelvin(self):
- return self.temperature.to(u.K).value
+ print("\n=== DilutePlanckianRadiationField.temperature_kelvin ===")
+ result = self.temperature.to(u.K).value
+ print(f"First 10 elements: {result[:10]}")
+ print(f"Shape: {result.shape}")
+ print(f"dtype: {result.dtype}")
+ print("="*50)
+ return result
def calculate_mean_intensity(self, nu: Union[u.Quantity, np.ndarray]):
"""
@@ -69,11 +93,34 @@ def calculate_mean_intensity(self, nu: Union[u.Quantity, np.ndarray]):
intensity : u.Quantity
Intensity of the radiation field at the given frequency
"""
- return self.dilution_factor * intensity_black_body(
+ print("\n=== DilutePlanckianRadiationField.calculate_mean_intensity ===")
+ print("Input frequency (nu):")
+ print(f" First 10 elements: {nu[:10]}")
+ print(f" Shape: {nu.shape}")
+ print(f" dtype: {nu.dtype}")
+ if isinstance(nu, u.Quantity):
+ print(f" unit: {nu.unit}")
+
+ result = self.dilution_factor * intensity_black_body(
nu[np.newaxis].T, self.temperature
)
+ print("\nOutput intensity:")
+ print(f" First 10 elements: {result[:10]}")
+ print(f" Shape: {result.shape}")
+ print(f" dtype: {result.dtype}")
+ if isinstance(result, u.Quantity):
+ print(f" unit: {result.unit}")
+ print("="*50)
+ return result
def to_planckian_radiation_field(self):
+ print("\n=== DilutePlanckianRadiationField.to_planckian_radiation_field ===")
+ print("Converting with temperature:")
+ print(f" First 10 elements: {self.temperature[:10]}")
+ print(f" Shape: {self.temperature.shape}")
+ print(f" dtype: {self.temperature.dtype}")
+ print(f" unit: {self.temperature.unit}")
+ print("="*50)
return PlanckianRadiationField(self.temperature)