Package evaluation to test RayTraceHeatTransfer on Julia 1.14.0-DEV.2064 (1d5dcac2d2*) started at 2026-04-21T15:56:43.922 ################################################################################ # Set-up # Installing PkgEval dependencies (TestEnv)... Activating project at `~/.julia/environments/v1.14` Set-up completed after 14.41s ################################################################################ # Installation # Installing RayTraceHeatTransfer... Resolving package versions... Updating `~/.julia/environments/v1.14/Project.toml` [7cf1493d] + RayTraceHeatTransfer v0.8.0 Updating `~/.julia/environments/v1.14/Manifest.toml` [66dad0bd] + AliasTables v1.1.3 [49dc2e85] + Calculus v0.5.2 [9a962f9c] + DataAPI v1.16.0 [864edb3b] + DataStructures v0.19.4 [ffbed154] + DocStringExtensions v0.9.5 [411431e0] + Extents v0.1.6 [5c1252a2] + GeometryBasics v0.5.10 [92d709cd] + IrrationalConstants v0.2.6 [c8e1da08] + IterTools v1.10.0 [692b3bcd] + JLLWrappers v1.7.1 [2ab3a3ac] + LogExpFunctions v0.3.29 [eff96d63] + Measurements v2.14.1 [e1d29d7a] + Missings v1.2.0 [bac558e1] + OrderedCollections v1.8.1 [aea7be01] + PrecompileTools v1.3.3 [21216c6a] + Preferences v1.5.2 [92933f4c] + ProgressMeter v1.11.0 [43287f4e] + PtrArrays v1.4.0 [7cf1493d] + RayTraceHeatTransfer v0.8.0 [a2af1166] + SortingAlgorithms v1.2.2 [90137ffa] + StaticArrays v1.9.18 [1e83bf80] + StaticArraysCore v1.4.4 [10745b16] + Statistics v1.11.1 [82ae8749] + StatsAPI v1.8.0 [2913bbd2] + StatsBase v0.34.10 [5ae413db] + EarCut_jll v2.2.4+0 [0dad84c5] + ArgTools v1.1.2 [56f22d72] + Artifacts v1.11.0 [2a0f44e3] + Base64 v1.11.0 [ade2ca70] + Dates v1.11.0 [8ba89e20] + Distributed v1.11.0 [f43a241f] + Downloads v1.7.0 [7b1f6079] + FileWatching v1.11.0 [ac6e5ff7] + JuliaSyntaxHighlighting v1.13.0 [b27032c2] + LibCURL v1.0.0 [76f85450] + LibGit2 v1.11.0 [8f399da3] + Libdl v1.11.0 [37e2e46d] + LinearAlgebra v1.13.0 [56ddb016] + Logging v1.11.0 [d6f4376e] + Markdown v1.11.0 [ca575930] + NetworkOptions v1.3.0 [44cfe95a] + Pkg v1.14.0 [de0858da] + Printf v1.11.0 [9a3f8284] + Random v1.11.0 [ea8e919c] + SHA v1.0.0 [9e88b42a] + Serialization v1.11.0 [6462fe0b] + Sockets v1.11.0 [2f01184e] + SparseArrays v1.13.0 [f489334b] + StyledStrings v1.13.0 [fa267f1f] + TOML v1.0.3 [a4e569a6] + Tar v1.10.0 [cf7118a7] + UUIDs v1.11.0 [4ec0a83e] + Unicode v1.11.0 [e66e0078] + CompilerSupportLibraries_jll v1.3.0+1 [deac9b47] + LibCURL_jll v8.19.0+0 [e37daf67] + LibGit2_jll v1.9.2+0 [29816b5a] + LibSSH2_jll v1.11.3+1 [14a3606d] + MozillaCACerts_jll v2026.3.19 [4536629a] + OpenBLAS_jll v0.3.30+0 [458c3c95] + OpenSSL_jll v3.5.6+0 [efcefdf7] + PCRE2_jll v10.47.0+0 [bea87d4a] + SuiteSparse_jll v7.10.1+0 [83775a58] + Zlib_jll v1.3.2+0 [3161d3a3] + Zstd_jll v1.5.7+1 [8e850b90] + libblastrampoline_jll v5.15.0+0 [8e850ede] + nghttp2_jll v1.69.0+0 [3f19e933] + p7zip_jll v17.8.0+0 Installation completed after 5.93s ################################################################################ # Precompilation # Precompiling PkgEval dependencies... Precompiling package dependencies... Precompiling project... 7.4 s ✓ RayTraceHeatTransfer 1 dependency successfully precompiled in 8 seconds. 60 already precompiled. Precompilation completed after 34.23s ################################################################################ # Testing # Testing RayTraceHeatTransfer Status `/tmp/jl_oGHYlP/Project.toml` [5c1252a2] GeometryBasics v0.5.10 [eff96d63] Measurements v2.14.1 [92933f4c] ProgressMeter v1.11.0 [7cf1493d] RayTraceHeatTransfer v0.8.0 [90137ffa] StaticArrays v1.9.18 [2913bbd2] StatsBase v0.34.10 [37e2e46d] LinearAlgebra v1.13.0 [9a3f8284] Random v1.11.0 [2f01184e] SparseArrays v1.13.0 [8dfed614] Test v1.11.0 Status `/tmp/jl_oGHYlP/Manifest.toml` [66dad0bd] AliasTables v1.1.3 [49dc2e85] Calculus v0.5.2 [9a962f9c] DataAPI v1.16.0 [864edb3b] DataStructures v0.19.4 [ffbed154] DocStringExtensions v0.9.5 [411431e0] Extents v0.1.6 [5c1252a2] GeometryBasics v0.5.10 [92d709cd] IrrationalConstants v0.2.6 [c8e1da08] IterTools v1.10.0 [692b3bcd] JLLWrappers v1.7.1 [2ab3a3ac] LogExpFunctions v0.3.29 [eff96d63] Measurements v2.14.1 [e1d29d7a] Missings v1.2.0 [bac558e1] OrderedCollections v1.8.1 [aea7be01] PrecompileTools v1.3.3 [21216c6a] Preferences v1.5.2 [92933f4c] ProgressMeter v1.11.0 [43287f4e] PtrArrays v1.4.0 [7cf1493d] RayTraceHeatTransfer v0.8.0 [a2af1166] SortingAlgorithms v1.2.2 [90137ffa] StaticArrays v1.9.18 [1e83bf80] StaticArraysCore v1.4.4 [10745b16] Statistics v1.11.1 [82ae8749] StatsAPI v1.8.0 [2913bbd2] StatsBase v0.34.10 [5ae413db] EarCut_jll v2.2.4+0 [0dad84c5] ArgTools v1.1.2 [56f22d72] Artifacts v1.11.0 [2a0f44e3] Base64 v1.11.0 [ade2ca70] Dates v1.11.0 [8ba89e20] Distributed v1.11.0 [f43a241f] Downloads v1.7.0 [7b1f6079] FileWatching v1.11.0 [b77e0a4c] InteractiveUtils v1.11.0 [ac6e5ff7] JuliaSyntaxHighlighting v1.13.0 [b27032c2] LibCURL v1.0.0 [76f85450] LibGit2 v1.11.0 [8f399da3] Libdl v1.11.0 [37e2e46d] LinearAlgebra v1.13.0 [56ddb016] Logging v1.11.0 [d6f4376e] Markdown v1.11.0 [ca575930] NetworkOptions v1.3.0 [44cfe95a] Pkg v1.14.0 [de0858da] Printf v1.11.0 [9a3f8284] Random v1.11.0 [ea8e919c] SHA v1.0.0 [9e88b42a] Serialization v1.11.0 [6462fe0b] Sockets v1.11.0 [2f01184e] SparseArrays v1.13.0 [f489334b] StyledStrings v1.13.0 [fa267f1f] TOML v1.0.3 [a4e569a6] Tar v1.10.0 [8dfed614] Test v1.11.0 [cf7118a7] UUIDs v1.11.0 [4ec0a83e] Unicode v1.11.0 [e66e0078] CompilerSupportLibraries_jll v1.3.0+1 [deac9b47] LibCURL_jll v8.19.0+0 [e37daf67] LibGit2_jll v1.9.2+0 [29816b5a] LibSSH2_jll v1.11.3+1 [14a3606d] MozillaCACerts_jll v2026.3.19 [4536629a] OpenBLAS_jll v0.3.30+0 [458c3c95] OpenSSL_jll v3.5.6+0 [efcefdf7] PCRE2_jll v10.47.0+0 [bea87d4a] SuiteSparse_jll v7.10.1+0 [83775a58] Zlib_jll v1.3.2+0 [3161d3a3] Zstd_jll v1.5.7+1 [8e850b90] libblastrampoline_jll v5.15.0+0 [8e850ede] nghttp2_jll v1.69.0+0 [3f19e933] p7zip_jll v17.8.0+0 Testing Running tests... ================================================================================ STARTING TEST SUITE ================================================================================ ------------------------------------------------------------ Testing 3D View Factors ------------------------------------------------------------ Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 1.8906416838689214e-15 Converged after 7 iterations. norm(E*F-F'*E) = 1.7554167342883506e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 1.461852203983726e-15 Converged after 7 iterations. norm(E*F-F'*E) = 8.777083671441753e-17 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 1.9794221072983377e-15 ┌ Warning: No convergence progress detected at iteration 6; norm(E*F-F'*E) = 2.5133742693021536e-16. Stopping. └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:331 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 2.8318849355912135e-15 ┌ Warning: No convergence progress detected at iteration 5; norm(E*F-F'*E) = 2.8576114088871287e-16. Stopping. └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:331 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 1.962615573354719e-15 Converged after 9 iterations. norm(E*F-F'*E) = 1.7554167342883506e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 2.211058699092487e-15 ┌ Warning: No convergence progress detected at iteration 5; norm(E*F-F'*E) = 4.791355229691893e-16. Stopping. └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:331 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 2.0206364052201327e-15 Converged after 7 iterations. norm(E*F-F'*E) = 1.3597399555105182e-16 ✓ 3D View Factor tests complete ------------------------------------------------------------ Testing 3D Heat Transfer ------------------------------------------------------------ Computing view factors (geometry only, wavelength-independent)... Matrix size: 150×150 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 1.0431932984236912e-14 Converged after 7 iterations. norm(E*F-F'*E) = 9.838939700371186e-17 === 3D Surface-Only Grey Solver === Found 150 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 150 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 150×150 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 1.0431932984236912e-14 Converged after 7 iterations. norm(E*F-F'*E) = 9.838939700371186e-17 === 3D Surface-Only Grey Solver === Found 150 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 150 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 150×150 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 1.0431932984236912e-14 Converged after 7 iterations. norm(E*F-F'*E) = 9.838939700371186e-17 === 3D Surface-Only Grey Solver === Found 150 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 150 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.8549284161345355e-16 === 3D Surface-Only Grey Solver === Found 96 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 96 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.604788474939644e-15 Converged after 7 iterations. norm(E*F-F'*E) = 1.4226597660905571e-16 === 3D Surface-Only Grey Solver === Found 96 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 96 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.836202745879726e-15 Converged after 7 iterations. norm(E*F-F'*E) = 1.1934330404686084e-16 === 3D Surface-Only Grey Solver === Found 96 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 96 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.57995603044491e-15 Converged after 7 iterations. norm(E*F-F'*E) = 1.0528380661555636e-16 === 3D Surface-Only Grey Solver === Found 96 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 96 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.8549284161345355e-16 === 3D Surface-Only Grey Solver === Found 96 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 96 surfaces Computing energy conservation error... === 3D Grey Solution Complete === ✓ 3D Heat Transfer tests complete ------------------------------------------------------------ Testing 2D Grey Participating Media ------------------------------------------------------------ Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 1%|▍ | ETA: 0:05:00 Bin 1 progress: 62%|████████████████████▍ | ETA: 0:00:04 Bin 1 progress: 99%|████████████████████████████████▊| ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:06 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.761739987160342e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.09090909090909094) ≥ min gas E (0.03305785123966937); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012448274195210366 Iteration 10: norm(E*F-F'*E) = 2.5748496081298536e-5 Iteration 20: norm(E*F-F'*E) = 3.4113198883395137e-7 Converged after 29 iterations. norm(E*F-F'*E) = 5.729751385290669e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 36%|███████████▊ | ETA: 0:00:03 Bin 1 progress: 78%|█████████████████████████▋ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:03 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.761739987160345e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.09090909090909094) ≥ min gas E (0.03305785123966936); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012484856713269935 Iteration 10: norm(E*F-F'*E) = 2.8279972488489214e-5 Iteration 20: norm(E*F-F'*E) = 3.4371870581013447e-7 Converged after 29 iterations. norm(E*F-F'*E) = 5.252758881613633e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 80%|██████████████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.761739987160345e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.09090909090909094) ≥ min gas E (0.03305785123966936); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012204661858603912 Iteration 10: norm(E*F-F'*E) = 2.1533921861620564e-5 Iteration 20: norm(E*F-F'*E) = 2.607523701623098e-7 Converged after 29 iterations. norm(E*F-F'*E) = 4.40546968849867e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▋ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▋ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.761739987160352e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.090909090909091) ≥ min gas E (0.03305785123966934); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012479417358099465 Iteration 10: norm(E*F-F'*E) = 2.857428562864479e-5 Iteration 20: norm(E*F-F'*E) = 3.353069656592337e-7 Converged after 29 iterations. norm(E*F-F'*E) = 5.095512817392086e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▍ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374624e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.1428571428571429) ≥ min gas E (0.08163265306122441); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011350080399587185 Iteration 10: norm(E*F-F'*E) = 3.347183657760137e-5 Iteration 20: norm(E*F-F'*E) = 4.399218463208738e-7 Iteration 30: norm(E*F-F'*E) = 6.5939331244600214e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.2363529294659277e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 79%|██████████████████████████▏ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374636e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.1428571428571429) ≥ min gas E (0.08163265306122437); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011120417905298755 Iteration 10: norm(E*F-F'*E) = 4.27322036375148e-5 Iteration 20: norm(E*F-F'*E) = 5.010784588665617e-7 Iteration 30: norm(E*F-F'*E) = 7.103072414405315e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.3181369074786615e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 79%|██████████████████████████▏ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374653e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.14285714285714293) ≥ min gas E (0.08163265306122433); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011148833619830106 Iteration 10: norm(E*F-F'*E) = 3.2876993895034756e-5 Iteration 20: norm(E*F-F'*E) = 4.4955678867970466e-7 Iteration 30: norm(E*F-F'*E) = 6.8930501289030285e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.2975959002729782e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▍ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374645e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.142857142857143) ≥ min gas E (0.08163265306122439); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011229823965646089 Iteration 10: norm(E*F-F'*E) = 3.479377660455181e-5 Iteration 20: norm(E*F-F'*E) = 4.852668012079107e-7 Iteration 30: norm(E*F-F'*E) = 7.414910349849079e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.3965777369535577e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 82%|███████████████████████████ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374624e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.1428571428571429) ≥ min gas E (0.08163265306122441); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.01105049810782192 Iteration 10: norm(E*F-F'*E) = 3.299360432511707e-5 Iteration 20: norm(E*F-F'*E) = 4.85952895193734e-7 Iteration 30: norm(E*F-F'*E) = 7.526508224095674e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.4197932443921591e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▋ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374645e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.14285714285714296) ≥ min gas E (0.08163265306122436); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011060439811770607 Iteration 10: norm(E*F-F'*E) = 2.802580057911332e-5 Iteration 20: norm(E*F-F'*E) = 3.922263497661051e-7 Iteration 30: norm(E*F-F'*E) = 6.139241518203554e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.161791264634951e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Smoothing single F matrix for grey extinction Matrix size: 45×45 Strategy: Serial Tolerance: 2.793967723846438e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.033579898749102895 Iteration 10: norm(E*F-F'*E) = 8.378900192675343e-5 Iteration 20: norm(E*F-F'*E) = 6.751572886914089e-7 Iteration 30: norm(E*F-F'*E) = 8.87634449311576e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.6308674935409745e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 5.116683173643062e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.1428571428571429) ≥ min gas E (0.08163265306122441); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.029136396700027854 Iteration 10: norm(E*F-F'*E) = 8.311805107405826e-5 Iteration 20: norm(E*F-F'*E) = 7.099063789905807e-7 Iteration 30: norm(E*F-F'*E) = 7.999256673181218e-9 Converged after 31 iterations. norm(E*F-F'*E) = 3.3414072804983947e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Smoothing single F matrix for grey extinction Matrix size: 117×117 Strategy: Serial Tolerance: 8.109247654827728e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.11111111111111116) ≥ min gas E (0.04938271604938266); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.024809618029634713 Iteration 10: norm(E*F-F'*E) = 7.360108378871834e-5 Iteration 20: norm(E*F-F'*E) = 1.1931361484424063e-6 Iteration 30: norm(E*F-F'*E) = 2.0406393606885333e-8 Converged after 33 iterations. norm(E*F-F'*E) = 4.009299317936783e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 1.1770089922830171e-8 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.09090909090909094) ≥ min gas E (0.03305785123966937); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.021700664920722777 Iteration 10: norm(E*F-F'*E) = 5.968496388779386e-5 Iteration 20: norm(E*F-F'*E) = 8.058144648047221e-7 Iteration 30: norm(E*F-F'*E) = 1.3269152672741507e-8 Converged after 31 iterations. norm(E*F-F'*E) = 5.853630510784972e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 36%|████████████ | ETA: 0:00:02 Bin 1 progress: 77%|█████████████████████████▎ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374624e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.1428571428571429) ≥ min gas E (0.08163265306122441); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011350080399587185 Iteration 10: norm(E*F-F'*E) = 3.347183657760137e-5 Iteration 20: norm(E*F-F'*E) = 4.399218463208738e-7 Iteration 30: norm(E*F-F'*E) = 6.5939331244600214e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.2363529294659277e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === ✓ 2D Grey Participating Media tests complete ------------------------------------------------------------ Testing 2D Spectral Participating Media ------------------------------------------------------------ Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 80%|██████████████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010875429008884887 Iteration 10: norm(E*F-F'*E) = 4.7073324895837455e-5 Iteration 20: norm(E*F-F'*E) = 5.270698631735426e-7 Iteration 30: norm(E*F-F'*E) = 6.846201464537942e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.139332014503023e-10 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▌ | ETA: 0:00:02 Bin 1 progress: 80%|██████████████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.01007952574314013 Iteration 10: norm(E*F-F'*E) = 2.8144717144846884e-5 Iteration 20: norm(E*F-F'*E) = 2.896117333355649e-7 Iteration 30: norm(E*F-F'*E) = 3.923584763730424e-9 Converged after 33 iterations. norm(E*F-F'*E) = 7.126999961087163e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 82%|███████████████████████████▏ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010875429008884887 Iteration 10: norm(E*F-F'*E) = 4.7073324895837455e-5 Iteration 20: norm(E*F-F'*E) = 5.270698631735426e-7 Iteration 30: norm(E*F-F'*E) = 6.846201464537942e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.139332014503023e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Running direct ray tracing for 10 spectral bins Processing spectral bin 1/10 ┌ Warning: No emitters found for spectral bin 1, skipping ray tracing └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 2/10 ┌ Warning: No emitters found for spectral bin 2, skipping ray tracing └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 3/10 ┌ Warning: No emitters found for spectral bin 3, skipping ray tracing └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 4/10 Bin 4 ray tracing: 0%| | ETA: 11:06:01 Bin 4 ray tracing: 9%|██▋ | ETA: 0:00:53 Bin 4 ray tracing: 17%|█████▏ | ETA: 0:00:29 Bin 4 ray tracing: 25%|███████▋ | ETA: 0:00:21 Bin 4 ray tracing: 34%|██████████▏ | ETA: 0:00:16 Bin 4 ray tracing: 42%|████████████▋ | ETA: 0:00:12 Bin 4 ray tracing: 51%|███████████████▏ | ETA: 0:00:10 Bin 4 ray tracing: 59%|█████████████████▋ | ETA: 0:00:08 Bin 4 ray tracing: 67%|████████████████████▏ | ETA: 0:00:06 Bin 4 ray tracing: 76%|██████████████████████▊ | ETA: 0:00:04 Bin 4 ray tracing: 84%|█████████████████████████▎ | ETA: 0:00:03 Bin 4 ray tracing: 93%|███████████████████████████▉ | ETA: 0:00:01 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 4 Energy per ray: 0.0001853335835185918 Processing spectral bin 5/10 Bin 5 ray tracing: 9%|██▋ | ETA: 0:00:11 Bin 5 ray tracing: 17%|█████▎ | ETA: 0:00:10 Bin 5 ray tracing: 26%|███████▊ | ETA: 0:00:09 Bin 5 ray tracing: 35%|██████████▍ | ETA: 0:00:08 Bin 5 ray tracing: 43%|█████████████ | ETA: 0:00:07 Bin 5 ray tracing: 52%|███████████████▌ | ETA: 0:00:06 Bin 5 ray tracing: 60%|██████████████████▏ | ETA: 0:00:05 Bin 5 ray tracing: 69%|████████████████████▊ | ETA: 0:00:04 Bin 5 ray tracing: 78%|███████████████████████▎ | ETA: 0:00:03 Bin 5 ray tracing: 86%|█████████████████████████▉ | ETA: 0:00:02 Bin 5 ray tracing: 95%|████████████████████████████▌ | ETA: 0:00:01 Bin 5 ray tracing: 100%|██████████████████████████████| Time: 0:00:11 Updating spectral results for spectral bin 5 Energy per ray: 0.04303963948070305 Processing spectral bin 6/10 Bin 6 ray tracing: 8%|██▌ | ETA: 0:00:11 Bin 6 ray tracing: 17%|█████▏ | ETA: 0:00:10 Bin 6 ray tracing: 25%|███████▋ | ETA: 0:00:09 Bin 6 ray tracing: 34%|██████████▏ | ETA: 0:00:08 Bin 6 ray tracing: 42%|████████████▋ | ETA: 0:00:07 Bin 6 ray tracing: 50%|███████████████▏ | ETA: 0:00:06 Bin 6 ray tracing: 59%|█████████████████▋ | ETA: 0:00:05 Bin 6 ray tracing: 67%|████████████████████▏ | ETA: 0:00:04 Bin 6 ray tracing: 75%|██████████████████████▋ | ETA: 0:00:03 Bin 6 ray tracing: 84%|█████████████████████████▎ | ETA: 0:00:02 Bin 6 ray tracing: 92%|███████████████████████████▋ | ETA: 0:00:01 Bin 6 ray tracing: 100%|██████████████████████████████| Time: 0:00:12 Updating spectral results for spectral bin 6 Energy per ray: 0.013246116789219256 Processing spectral bin 7/10 Bin 7 ray tracing: 8%|██▍ | ETA: 0:00:12 Bin 7 ray tracing: 16%|████▉ | ETA: 0:00:10 Bin 7 ray tracing: 25%|███████▍ | ETA: 0:00:09 Bin 7 ray tracing: 33%|██████████ | ETA: 0:00:08 Bin 7 ray tracing: 42%|████████████▋ | ETA: 0:00:07 Bin 7 ray tracing: 50%|███████████████▏ | ETA: 0:00:06 Bin 7 ray tracing: 59%|█████████████████▋ | ETA: 0:00:05 Bin 7 ray tracing: 67%|████████████████████▏ | ETA: 0:00:04 Bin 7 ray tracing: 76%|██████████████████████▋ | ETA: 0:00:03 Bin 7 ray tracing: 84%|█████████████████████████▏ | ETA: 0:00:02 Bin 7 ray tracing: 92%|███████████████████████████▊ | ETA: 0:00:01 Bin 7 ray tracing: 100%|██████████████████████████████| Time: 0:00:11 Updating spectral results for spectral bin 7 Energy per ray: 0.00021661482457376898 Processing spectral bin 8/10 Bin 8 ray tracing: 8%|██▌ | ETA: 0:00:11 Bin 8 ray tracing: 17%|█████ | ETA: 0:00:10 Bin 8 ray tracing: 25%|███████▌ | ETA: 0:00:09 Bin 8 ray tracing: 33%|██████████ | ETA: 0:00:08 Bin 8 ray tracing: 42%|████████████▌ | ETA: 0:00:07 Bin 8 ray tracing: 50%|███████████████▏ | ETA: 0:00:06 Bin 8 ray tracing: 59%|█████████████████▋ | ETA: 0:00:05 Bin 8 ray tracing: 67%|████████████████████▎ | ETA: 0:00:04 Bin 8 ray tracing: 76%|██████████████████████▊ | ETA: 0:00:03 Bin 8 ray tracing: 84%|█████████████████████████▎ | ETA: 0:00:02 Bin 8 ray tracing: 93%|███████████████████████████▊ | ETA: 0:00:01 Bin 8 ray tracing: 100%|██████████████████████████████| Time: 0:00:11 Updating spectral results for spectral bin 8 Energy per ray: 1.0195075180910974e-6 Processing spectral bin 9/10 Bin 9 ray tracing: 9%|██▌ | ETA: 0:00:11 Bin 9 ray tracing: 17%|█████▏ | ETA: 0:00:10 Bin 9 ray tracing: 26%|███████▊ | ETA: 0:00:09 Bin 9 ray tracing: 34%|██████████▎ | ETA: 0:00:08 Bin 9 ray tracing: 43%|████████████▊ | ETA: 0:00:07 Bin 9 ray tracing: 51%|███████████████▎ | ETA: 0:00:06 Bin 9 ray tracing: 59%|█████████████████▊ | ETA: 0:00:05 Bin 9 ray tracing: 68%|████████████████████▍ | ETA: 0:00:04 Bin 9 ray tracing: 76%|██████████████████████▉ | ETA: 0:00:03 Bin 9 ray tracing: 85%|█████████████████████████▍ | ETA: 0:00:02 Bin 9 ray tracing: 93%|████████████████████████████ | ETA: 0:00:01 Bin 9 ray tracing: 100%|██████████████████████████████| Time: 0:00:11 Updating spectral results for spectral bin 9 Energy per ray: 2.17242363082386e-9 Processing spectral bin 10/10 Bin 10 ray tracing: 8%|██▍ | ETA: 0:00:11 Bin 10 ray tracing: 17%|████▉ | ETA: 0:00:10 Bin 10 ray tracing: 25%|███████▎ | ETA: 0:00:09 Bin 10 ray tracing: 33%|█████████▋ | ETA: 0:00:08 Bin 10 ray tracing: 42%|████████████ | ETA: 0:00:07 Bin 10 ray tracing: 50%|██████████████▍ | ETA: 0:00:06 Bin 10 ray tracing: 58%|████████████████▉ | ETA: 0:00:05 Bin 10 ray tracing: 67%|███████████████████▍ | ETA: 0:00:04 Bin 10 ray tracing: 75%|█████████████████████▉ | ETA: 0:00:03 Bin 10 ray tracing: 84%|████████████████████████▎ | ETA: 0:00:02 Bin 10 ray tracing: 92%|██████████████████████████▋ | ETA: 0:00:01 Bin 10 ray tracing: 100%|█████████████████████████████| Time: 0:00:12 Updating spectral results for spectral bin 10 Energy per ray: 1.5017824710260812e-5 Building intermediate mesh... Optimizing mesh... Extinction variation detected across the spectrum, ray tracing each spectral bin separately No spectral variation detected across walls Spectral variation detected across volumes, using spectral solver Building spatial acceleration structures... Computing 10 separate F matrices for variable spectral extinction Computing F matrix for spectral bin 1/10 Using 1 threads for spectral bin 1 Bin 1 progress: 22%|███████▍ | ETA: 0:00:04 Bin 1 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 1 progress: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 93%|██████████████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 2/10 Using 1 threads for spectral bin 2 Bin 2 progress: 24%|████████▏ | ETA: 0:00:03 Bin 2 progress: 49%|████████████████▏ | ETA: 0:00:02 Bin 2 progress: 73%|████████████████████████▎ | ETA: 0:00:01 Bin 2 progress: 98%|████████████████████████████████▎| ETA: 0:00:00 Bin 2 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 3/10 Using 1 threads for spectral bin 3 Bin 3 progress: 22%|███████▍ | ETA: 0:00:04 Bin 3 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 3 progress: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 3 progress: 93%|██████████████████████████████▊ | ETA: 0:00:00 Bin 3 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 4/10 Using 1 threads for spectral bin 4 Bin 4 progress: 24%|████████▏ | ETA: 0:00:03 Bin 4 progress: 47%|███████████████▍ | ETA: 0:00:02 Bin 4 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 4 progress: 96%|███████████████████████████████▌ | ETA: 0:00:00 Bin 4 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 5/10 Using 1 threads for spectral bin 5 Bin 5 progress: 22%|███████▍ | ETA: 0:00:04 Bin 5 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 5 progress: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 5 progress: 91%|██████████████████████████████▏ | ETA: 0:00:00 Bin 5 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 6/10 Using 1 threads for spectral bin 6 Bin 6 progress: 24%|████████▏ | ETA: 0:00:03 Bin 6 progress: 47%|███████████████▍ | ETA: 0:00:02 Bin 6 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 6 progress: 96%|███████████████████████████████▌ | ETA: 0:00:00 Bin 6 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 7/10 Using 1 threads for spectral bin 7 Bin 7 progress: 24%|████████▏ | ETA: 0:00:03 Bin 7 progress: 47%|███████████████▍ | ETA: 0:00:02 Bin 7 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 7 progress: 96%|███████████████████████████████▌ | ETA: 0:00:00 Bin 7 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 8/10 Using 1 threads for spectral bin 8 Bin 8 progress: 24%|████████▏ | ETA: 0:00:03 Bin 8 progress: 47%|███████████████▍ | ETA: 0:00:02 Bin 8 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 8 progress: 93%|██████████████████████████████▊ | ETA: 0:00:00 Bin 8 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 9/10 Using 1 threads for spectral bin 9 Bin 9 progress: 24%|████████▏ | ETA: 0:00:03 Bin 9 progress: 49%|████████████████▏ | ETA: 0:00:02 Bin 9 progress: 73%|████████████████████████▎ | ETA: 0:00:01 Bin 9 progress: 96%|███████████████████████████████▌ | ETA: 0:00:00 Bin 9 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 10/10 Using 1 threads for spectral bin 10 Bin 10 progress: 22%|███████▏ | ETA: 0:00:04 Bin 10 progress: 44%|██████████████▎ | ETA: 0:00:03 Bin 10 progress: 69%|██████████████████████ | ETA: 0:00:01 Bin 10 progress: 91%|█████████████████████████████▏ | ETA: 0:00:00 Bin 10 progress: 100%|████████████████████████████████| Time: 0:00:04 Smoothing F matrix for spectral bin 1/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010875429008884887 Iteration 10: norm(E*F-F'*E) = 4.7073324895837455e-5 Iteration 20: norm(E*F-F'*E) = 5.270698631735426e-7 Iteration 30: norm(E*F-F'*E) = 6.846201464537942e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.139332014503023e-10 Smoothing F matrix for spectral bin 2/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.372642145652105e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.1601777777777777); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010048318677072478 Iteration 10: norm(E*F-F'*E) = 2.856639769314018e-5 Iteration 20: norm(E*F-F'*E) = 2.9256177520589275e-7 Iteration 30: norm(E*F-F'*E) = 3.943990224498512e-9 Converged after 33 iterations. norm(E*F-F'*E) = 7.154635651711205e-10 Smoothing F matrix for spectral bin 3/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.363359837286635e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.16035555555555547); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010104455960999155 Iteration 10: norm(E*F-F'*E) = 4.0629822915258255e-5 Iteration 20: norm(E*F-F'*E) = 4.635732225057233e-7 Iteration 30: norm(E*F-F'*E) = 5.852654160691107e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.3129943423159126e-10 Smoothing F matrix for spectral bin 4/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.354098087743682e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.16053333333333328); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010701939423257789 Iteration 10: norm(E*F-F'*E) = 3.6014204894019805e-5 Iteration 20: norm(E*F-F'*E) = 3.9911675216502977e-7 Iteration 30: norm(E*F-F'*E) = 5.286467426174762e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.018817906372915e-10 Smoothing F matrix for spectral bin 5/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.344856828797063e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.16071111111111103); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010202461027404105 Iteration 10: norm(E*F-F'*E) = 3.395667994730951e-5 Iteration 20: norm(E*F-F'*E) = 3.726833294829022e-7 Iteration 30: norm(E*F-F'*E) = 4.84588326450314e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.656851892078086e-10 Smoothing F matrix for spectral bin 6/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.335635992522152e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.1608888888888888); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010404806114800018 Iteration 10: norm(E*F-F'*E) = 2.5386597719335338e-5 Iteration 20: norm(E*F-F'*E) = 2.3597460846722936e-7 Iteration 30: norm(E*F-F'*E) = 2.996608618605295e-9 Converged after 33 iterations. norm(E*F-F'*E) = 5.286916304827337e-10 Smoothing F matrix for spectral bin 7/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.326435511294201e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.16106666666666658); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.0104863342307487 Iteration 10: norm(E*F-F'*E) = 3.276783853721903e-5 Iteration 20: norm(E*F-F'*E) = 3.417959021253451e-7 Iteration 30: norm(E*F-F'*E) = 4.560817917658331e-9 Converged after 33 iterations. norm(E*F-F'*E) = 8.220904035298203e-10 Smoothing F matrix for spectral bin 8/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.317255317786709e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.16124444444444438); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011316191122273297 Iteration 10: norm(E*F-F'*E) = 5.301242802888478e-5 Iteration 20: norm(E*F-F'*E) = 6.511976977880824e-7 Iteration 30: norm(E*F-F'*E) = 8.817323854186376e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.735376742077596e-10 Smoothing F matrix for spectral bin 9/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.308095344969763e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.16142222222222216); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010531255481722852 Iteration 10: norm(E*F-F'*E) = 3.832289818813306e-5 Iteration 20: norm(E*F-F'*E) = 4.5872483610221405e-7 Iteration 30: norm(E*F-F'*E) = 6.40298998875529e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.9880748886734e-10 Smoothing F matrix for spectral bin 10/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.298955526108412e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.1615999999999999); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011461542330417978 Iteration 10: norm(E*F-F'*E) = 5.031922469788149e-5 Iteration 20: norm(E*F-F'*E) = 5.93509926102589e-7 Iteration 30: norm(E*F-F'*E) = 7.495065557905907e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.511317621749189e-10 === Using FULL spectral solver === ==== Building and Factorizing Block matrix ==== Starting spectral steady-state iteration... Converged after 8 iterations Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... Extinction variation detected across the spectrum, ray tracing each spectral bin separately No spectral variation detected across walls Spectral variation detected across volumes, using spectral solver Building spatial acceleration structures... Running direct ray tracing for 10 spectral bins Processing spectral bin 1/10 ┌ Warning: No emitters found for spectral bin 1, skipping ray tracing └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 2/10 ┌ Warning: No emitters found for spectral bin 2, skipping ray tracing └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 3/10 ┌ Warning: No emitters found for spectral bin 3, skipping ray tracing └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 4/10 Bin 4 ray tracing: 7%|██ | ETA: 0:00:15 Bin 4 ray tracing: 13%|███▉ | ETA: 0:00:14 Bin 4 ray tracing: 19%|█████▊ | ETA: 0:00:13 Bin 4 ray tracing: 25%|███████▋ | ETA: 0:00:12 Bin 4 ray tracing: 32%|█████████▌ | ETA: 0:00:11 Bin 4 ray tracing: 38%|███████████▌ | ETA: 0:00:10 Bin 4 ray tracing: 44%|█████████████▍ | ETA: 0:00:09 Bin 4 ray tracing: 51%|███████████████▎ | ETA: 0:00:08 Bin 4 ray tracing: 57%|█████████████████▏ | ETA: 0:00:07 Bin 4 ray tracing: 63%|███████████████████ | ETA: 0:00:06 Bin 4 ray tracing: 69%|████████████████████▉ | ETA: 0:00:05 Bin 4 ray tracing: 76%|██████████████████████▊ | ETA: 0:00:04 Bin 4 ray tracing: 82%|████████████████████████▋ | ETA: 0:00:03 Bin 4 ray tracing: 88%|██████████████████████████▌ | ETA: 0:00:02 Bin 4 ray tracing: 94%|████████████████████████████▍ | ETA: 0:00:01 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 Updating spectral results for spectral bin 4 Energy per ray: 0.0001853335835185918 Processing spectral bin 5/10 Bin 5 ray tracing: 6%|██ | ETA: 0:00:14 Bin 5 ray tracing: 13%|███▉ | ETA: 0:00:13 Bin 5 ray tracing: 19%|█████▉ | ETA: 0:00:13 Bin 5 ray tracing: 26%|███████▊ | ETA: 0:00:12 Bin 5 ray tracing: 32%|█████████▋ | ETA: 0:00:11 Bin 5 ray tracing: 38%|███████████▌ | ETA: 0:00:10 Bin 5 ray tracing: 45%|█████████████▍ | ETA: 0:00:09 Bin 5 ray tracing: 51%|███████████████▎ | ETA: 0:00:08 Bin 5 ray tracing: 57%|█████████████████▎ | ETA: 0:00:07 Bin 5 ray tracing: 64%|███████████████████▏ | ETA: 0:00:06 Bin 5 ray tracing: 70%|█████████████████████ | ETA: 0:00:05 Bin 5 ray tracing: 76%|██████████████████████▉ | ETA: 0:00:04 Bin 5 ray tracing: 83%|████████████████████████▊ | ETA: 0:00:03 Bin 5 ray tracing: 89%|██████████████████████████▊ | ETA: 0:00:02 Bin 5 ray tracing: 95%|████████████████████████████▋ | ETA: 0:00:01 Bin 5 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 5 Energy per ray: 0.04303963948070305 Processing spectral bin 6/10 Bin 6 ray tracing: 6%|█▉ | ETA: 0:00:15 Bin 6 ray tracing: 13%|███▊ | ETA: 0:00:14 Bin 6 ray tracing: 19%|█████▋ | ETA: 0:00:13 Bin 6 ray tracing: 25%|███████▌ | ETA: 0:00:12 Bin 6 ray tracing: 31%|█████████▍ | ETA: 0:00:11 Bin 6 ray tracing: 38%|███████████▍ | ETA: 0:00:10 Bin 6 ray tracing: 44%|█████████████▍ | ETA: 0:00:09 Bin 6 ray tracing: 51%|███████████████▎ | ETA: 0:00:08 Bin 6 ray tracing: 57%|█████████████████▏ | ETA: 0:00:07 Bin 6 ray tracing: 64%|███████████████████▏ | ETA: 0:00:06 Bin 6 ray tracing: 70%|█████████████████████ | ETA: 0:00:05 Bin 6 ray tracing: 76%|██████████████████████▉ | ETA: 0:00:04 Bin 6 ray tracing: 83%|████████████████████████▉ | ETA: 0:00:03 Bin 6 ray tracing: 89%|██████████████████████████▋ | ETA: 0:00:02 Bin 6 ray tracing: 95%|████████████████████████████▌ | ETA: 0:00:01 Bin 6 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 6 Energy per ray: 0.013246116789219256 Processing spectral bin 7/10 Bin 7 ray tracing: 6%|█▉ | ETA: 0:00:15 Bin 7 ray tracing: 13%|███▊ | ETA: 0:00:14 Bin 7 ray tracing: 19%|█████▋ | ETA: 0:00:13 Bin 7 ray tracing: 25%|███████▋ | ETA: 0:00:12 Bin 7 ray tracing: 32%|█████████▌ | ETA: 0:00:11 Bin 7 ray tracing: 38%|███████████▍ | ETA: 0:00:10 Bin 7 ray tracing: 44%|█████████████▏ | ETA: 0:00:09 Bin 7 ray tracing: 50%|███████████████▏ | ETA: 0:00:08 Bin 7 ray tracing: 57%|█████████████████ | ETA: 0:00:07 Bin 7 ray tracing: 63%|██████████████████▉ | ETA: 0:00:06 Bin 7 ray tracing: 69%|████████████████████▉ | ETA: 0:00:05 Bin 7 ray tracing: 76%|██████████████████████▊ | ETA: 0:00:04 Bin 7 ray tracing: 82%|████████████████████████▋ | ETA: 0:00:03 Bin 7 ray tracing: 89%|██████████████████████████▋ | ETA: 0:00:02 Bin 7 ray tracing: 95%|████████████████████████████▌ | ETA: 0:00:01 Bin 7 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 7 Energy per ray: 0.00021661482457376898 Processing spectral bin 8/10 Bin 8 ray tracing: 7%|██ | ETA: 0:00:15 Bin 8 ray tracing: 13%|███▉ | ETA: 0:00:14 Bin 8 ray tracing: 19%|█████▊ | ETA: 0:00:13 Bin 8 ray tracing: 26%|███████▋ | ETA: 0:00:12 Bin 8 ray tracing: 32%|█████████▋ | ETA: 0:00:11 Bin 8 ray tracing: 38%|███████████▌ | ETA: 0:00:10 Bin 8 ray tracing: 44%|█████████████▍ | ETA: 0:00:09 Bin 8 ray tracing: 51%|███████████████▎ | ETA: 0:00:08 Bin 8 ray tracing: 57%|█████████████████ | ETA: 0:00:07 Bin 8 ray tracing: 63%|██████████████████▉ | ETA: 0:00:06 Bin 8 ray tracing: 69%|████████████████████▉ | ETA: 0:00:05 Bin 8 ray tracing: 76%|██████████████████████▊ | ETA: 0:00:04 Bin 8 ray tracing: 82%|████████████████████████▋ | ETA: 0:00:03 Bin 8 ray tracing: 89%|██████████████████████████▋ | ETA: 0:00:02 Bin 8 ray tracing: 95%|████████████████████████████▌ | ETA: 0:00:01 Bin 8 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 8 Energy per ray: 1.0195075180910974e-6 Processing spectral bin 9/10 Bin 9 ray tracing: 6%|█▉ | ETA: 0:00:15 Bin 9 ray tracing: 13%|███▉ | ETA: 0:00:14 Bin 9 ray tracing: 19%|█████▊ | ETA: 0:00:13 Bin 9 ray tracing: 26%|███████▊ | ETA: 0:00:12 Bin 9 ray tracing: 32%|█████████▊ | ETA: 0:00:11 Bin 9 ray tracing: 39%|███████████▋ | ETA: 0:00:10 Bin 9 ray tracing: 45%|█████████████▋ | ETA: 0:00:09 Bin 9 ray tracing: 52%|███████████████▌ | ETA: 0:00:08 Bin 9 ray tracing: 58%|█████████████████▌ | ETA: 0:00:07 Bin 9 ray tracing: 65%|███████████████████▍ | ETA: 0:00:06 Bin 9 ray tracing: 71%|█████████████████████▎ | ETA: 0:00:05 Bin 9 ray tracing: 77%|███████████████████████▎ | ETA: 0:00:04 Bin 9 ray tracing: 84%|█████████████████████████▏ | ETA: 0:00:03 Bin 9 ray tracing: 90%|███████████████████████████▏ | ETA: 0:00:02 Bin 9 ray tracing: 96%|████████████████████████████▉ | ETA: 0:00:01 Bin 9 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 9 Energy per ray: 2.17242363082386e-9 Processing spectral bin 10/10 Bin 10 ray tracing: 6%|█▉ | ETA: 0:00:15 Bin 10 ray tracing: 13%|███▊ | ETA: 0:00:14 Bin 10 ray tracing: 19%|█████▌ | ETA: 0:00:13 Bin 10 ray tracing: 25%|███████▍ | ETA: 0:00:12 Bin 10 ray tracing: 32%|█████████▎ | ETA: 0:00:11 Bin 10 ray tracing: 38%|███████████▏ | ETA: 0:00:10 Bin 10 ray tracing: 45%|████████████▉ | ETA: 0:00:09 Bin 10 ray tracing: 51%|██████████████▊ | ETA: 0:00:08 Bin 10 ray tracing: 57%|████████████████▋ | ETA: 0:00:07 Bin 10 ray tracing: 64%|██████████████████▌ | ETA: 0:00:06 Bin 10 ray tracing: 70%|████████████████████▍ | ETA: 0:00:05 Bin 10 ray tracing: 77%|██████████████████████▎ | ETA: 0:00:04 Bin 10 ray tracing: 83%|████████████████████████▎ | ETA: 0:00:03 Bin 10 ray tracing: 90%|██████████████████████████ | ETA: 0:00:02 Bin 10 ray tracing: 96%|████████████████████████████ | ETA: 0:00:01 Bin 10 ray tracing: 100%|█████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 10 Energy per ray: 1.5017824710260812e-5 Iter 1: T = 967.286944167656 K, relative_change = 0.032713055832344 Iter 2: T = 936.6474931969096 K, relative_change = 0.031675658557669656 Iter 3: T = 908.0503808626047 K, relative_change = 0.030531349885642708 Iter 20: T = 677.0296752837695 K, relative_change = 0.006151391227235582 Iter 40: T = 654.1383864960015 K, relative_change = 0.00022367171940294693 Iter 60: T = 653.3917614546724 K, relative_change = 6.882171327491049e-6 Iter 80: T = 653.3688738038192 K, relative_change = 2.1056626975524984e-7 Converged in 96 iterations to T = 653.3681958993276 K Iter 1: T = 970.2879224485687 K, relative_change = 0.029712077551431328 Iter 2: T = 942.7390613179441 K, relative_change = 0.028392460107205763 Iter 3: T = 917.3089949720064 K, relative_change = 0.026974660740572913 Iter 20: T = 728.7576548078575 K, relative_change = 0.004115328810681001 Iter 40: T = 712.7025539002323 K, relative_change = 0.00014087993153765138 Iter 60: T = 712.1906938771575 K, relative_change = 4.3253971668921165e-6 Iter 80: T = 712.175015136042 K, relative_change = 1.3233066191499306e-7 Converged in 93 iterations to T = 712.1745716523184 K Iter 1: T = 974.5084625778018 K, relative_change = 0.025491537422198212 Iter 2: T = 951.2054632118934 K, relative_change = 0.023912567474546644 Iter 3: T = 930.0154939454912 K, relative_change = 0.022276963375347902 Iter 20: T = 786.9832137750134 K, relative_change = 0.002559244819925182 Iter 40: T = 776.4097036197868 K, relative_change = 8.384544324835509e-5 Iter 60: T = 776.0780615561412 K, relative_change = 2.570464021890355e-6 Iter 80: T = 776.0679084999113 K, relative_change = 7.863687285009992e-8 Converged in 90 iterations to T = 776.067644142456 K Iter 1: T = 970.4459632083435 K, relative_change = 0.029554036791656496 Iter 2: T = 943.0581985726596 K, relative_change = 0.02822183374861857 Iter 3: T = 917.7913274994804 K, relative_change = 0.026792483339226836 Iter 20: T = 731.197294095624 K, relative_change = 0.004036766895857317 Iter 40: T = 715.4117437039116 K, relative_change = 0.00013787842499342142 Iter 60: T = 714.9089027578826 K, relative_change = 4.232911792685695e-6 Iter 80: T = 714.8935007142167 K, relative_change = 1.295008649214207e-7 Converged in 93 iterations to T = 714.8930650574428 K Iter 1: T = 969.3540847470745 K, relative_change = 0.030645915252925513 Iter 2: T = 940.8499472229784 K, relative_change = 0.029405289535178975 Iter 3: T = 914.4483434149098 K, relative_change = 0.02806143943143729 Iter 20: T = 713.8247253026699 K, relative_change = 0.004627030722893739 Iter 40: T = 696.0271681813962 K, relative_change = 0.00016076575060992857 Iter 60: T = 695.4565885576128 K, relative_change = 4.93850403416667e-6 Iter 80: T = 695.4391079106045 K, relative_change = 1.5109038210908124e-7 Converged in 94 iterations to T = 695.4386043062218 K Iter 1: T = 963.560077441969 K, relative_change = 0.03643992255803092 Iter 2: T = 928.99772552625 K, relative_change = 0.03586943121125801 Iter 3: T = 896.2793943439965 K, relative_change = 0.03521895725171929 Iter 20: T = 590.1547745305548 K, relative_change = 0.012007769999465526 Iter 40: T = 547.8703997078031 K, relative_change = 0.0005306280556283289 Iter 60: T = 546.3814434416892 K, relative_change = 1.6458754026211648e-5 Iter 80: T = 546.3356667544736 K, relative_change = 5.03695738936158e-7 Iter 100: T = 546.3342662095969 K, relative_change = 1.540847567357114e-8 Converged in 101 iterations to T = 546.3342591383267 K Iter 1: T = 966.8734899744658 K, relative_change = 0.033126510025534225 Iter 2: T = 935.8035000734723 K, relative_change = 0.03213449352284341 Iter 3: T = 906.7596770337203 K, relative_change = 0.031036241088510146 Iter 20: T = 668.9181081627411 K, relative_change = 0.006547011863425353 Iter 40: T = 644.7198743847567 K, relative_change = 0.00024095976751950776 Iter 60: T = 643.9269589126915 K, relative_change = 7.417455425592291e-6 Iter 80: T = 643.9026482723381 K, relative_change = 2.2694692609432572e-7 Converged in 96 iterations to T = 643.9019282170151 K Iter 1: T = 965.193308076543 K, relative_change = 0.03480669192345707 Iter 2: T = 932.3617748476245 K, relative_change = 0.03401550026734631 Iter 3: T = 901.4759521456529 K, relative_change = 0.03312643604143817 Iter 20: T = 632.5614687199453 K, relative_change = 0.008653910087742045 Iter 40: T = 601.4359297256046 K, relative_change = 0.00034055434557507015 Iter 60: T = 600.3892762158067 K, relative_change = 1.051060200418661e-5 Iter 80: T = 600.3571557886818 K, relative_change = 3.216116939020975e-7 Converged in 98 iterations to T = 600.3561860520099 K Iter 1: T = 980.1115211345509 K, relative_change = 0.019888478865449134 Iter 2: T = 962.2680734488563 K, relative_change = 0.01820552794343187 Iter 3: T = 946.3489722190885 K, relative_change = 0.016543312273379567 Iter 20: T = 848.663830334824 K, relative_change = 0.0014529891637382927 Iter 40: T = 842.2775919367575 K, relative_change = 4.618747388756374e-5 Iter 60: T = 842.0794920219631 K, relative_change = 1.414591391619414e-6 Iter 80: T = 842.0734294347799 K, relative_change = 4.327456216845732e-8 Converged in 87 iterations to T = 842.0732945771307 K Iter 1: T = 976.409675824731 K, relative_change = 0.02359032417526895 Iter 2: T = 954.9815581445108 K, relative_change = 0.02194582684990373 Iter 3: T = 935.6242370696938 K, relative_change = 0.020269837579301054 Iter 20: T = 809.5006064960725 K, relative_change = 0.0021034090750447687 Iter 40: T = 800.6118601038842 K, relative_change = 6.805291333047397e-5 Iter 60: T = 800.3343451372863 K, relative_change = 2.085452964785743e-6 Iter 80: T = 800.3258504145555 K, relative_change = 6.379837771938795e-8 Converged in 89 iterations to T = 800.3256381668882 K Iter 1: T = 980.8765766520859 K, relative_change = 0.019123423347914073 Iter 2: T = 963.7632508160249 K, relative_change = 0.01744697166127867 Iter 3: T = 948.5340658215957 K, relative_change = 0.01580179051394047 Iter 20: T = 856.1214654389715 K, relative_change = 0.0013463968757483691 Iter 40: T = 850.1594473493172 K, relative_change = 4.267673987515924e-5 Iter 60: T = 849.9746999476826 K, relative_change = 1.3069483036889187e-6 Iter 80: T = 849.9690461827382 K, relative_change = 3.9981480099062065e-8 Converged in 86 iterations to T = 849.9689304472226 K Iter 1: T = 967.3244364033394 K, relative_change = 0.032675563596660545 Iter 2: T = 936.7239698205686 K, relative_change = 0.031634129596216996 Iter 3: T = 908.1672389151189 K, relative_change = 0.030485748017017037 Iter 20: T = 677.7512404940857 K, relative_change = 0.006117349471517494 Iter 40: T = 654.9726129194307 K, relative_change = 0.0002222034745316533 Iter 60: T = 654.2299560029285 K, relative_change = 6.8367327973769716e-6 Iter 80: T = 654.2071903094721 K, relative_change = 2.0917578992912073e-7 Converged in 96 iterations to T = 654.2065160174977 K Iter 1: T = 973.4990409920373 K, relative_change = 0.02650095900796266 Iter 2: T = 949.191137042029 K, relative_change = 0.0249696229030051 Iter 3: T = 927.0089903278689 K, relative_change = 0.023369525745137624 Iter 20: T = 774.2067979261365 K, relative_change = 0.0028498331287511976 Iter 40: T = 762.5820067593266 K, relative_change = 9.412200920299757e-5 Iter 60: T = 762.216302523598 K, relative_change = 2.8862860490876026e-6 Iter 80: T = 762.2051055904057 K, relative_change = 8.829937638702002e-8 Converged in 91 iterations to T = 762.2048041645669 K Iter 1: T = 970.030713779482 K, relative_change = 0.029969286220518067 Iter 2: T = 942.2193169701299 K, relative_change = 0.0286706352843117 Iter 3: T = 916.5228975176335 K, relative_change = 0.027272227378151916 Iter 20: T = 724.7348016353941 K, relative_change = 0.004247871481688214 Iter 40: T = 708.2262393847146 K, relative_change = 0.00014597459827979487 Iter 60: T = 707.6991680555624 K, relative_change = 4.482412159592231e-6 Iter 80: T = 707.6830226043941 K, relative_change = 1.3713491806649636e-7 Converged in 93 iterations to T = 707.6825659187931 K Iter 1: T = 973.5623843776241 K, relative_change = 0.02643761562237588 Iter 2: T = 949.3177345222334 K, relative_change = 0.024903026497772623 Iter 3: T = 927.1982465054932 K, relative_change = 0.023300405346237757 Iter 20: T = 775.0270784765041 K, relative_change = 0.002830418743761514 Iter 40: T = 763.4720401661591 K, relative_change = 9.343025202526483e-5 Iter 60: T = 763.1086030038892 K, relative_change = 2.8650214909254846e-6 Iter 80: T = 763.0974755550186 K, relative_change = 8.764878708098425e-8 Converged in 91 iterations to T = 763.0971759997846 K Iter 1: T = 964.3082623802768 K, relative_change = 0.035691737619723246 Iter 2: T = 930.5410764275539 K, relative_change = 0.03501700365956912 Iter 3: T = 898.6674381728003 K, relative_change = 0.034252800937192185 Iter 20: T = 610.6510989564626 K, relative_change = 0.010244375648701284 Iter 40: T = 574.2741502504275 K, relative_change = 0.00042527080125208454 Iter 60: T = 573.0248943780207 K, relative_change = 1.315438974781628e-5 Iter 80: T = 572.9865255423178 K, relative_change = 4.0253603122748547e-7 Converged in 99 iterations to T = 572.9853587338231 K Iter 1: T = 963.5786496758909 K, relative_change = 0.03642135032410912 Iter 2: T = 929.0360829343492 K, relative_change = 0.03584820684140355 Iter 3: T = 896.3388274692396 K, relative_change = 0.035194817581073644 Iter 20: T = 590.6891386882064 K, relative_change = 0.011957964272856659 Iter 40: T = 548.5733144425892 K, relative_change = 0.0005274682885517264 Iter 60: T = 547.0913816783886 K, relative_change = 1.63593867326161e-5 Iter 80: T = 547.0458222942751 K, relative_change = 5.006534584393187e-7 Iter 100: T = 547.0444283991183 K, relative_change = 1.5315408552706643e-8 Converged in 101 iterations to T = 547.0444213614223 K Iter 1: T = 969.2828291676557 K, relative_change = 0.030717170832344327 Iter 2: T = 940.7055620251699 K, relative_change = 0.02948289836829744 Iter 3: T = 914.229313699623 K, relative_change = 0.028145095973014404 Iter 20: T = 712.6467344923897 K, relative_change = 0.004669768816978213 Iter 40: T = 694.7045803875204 K, relative_change = 0.00016245346721614263 Iter 60: T = 694.1290948623358 K, relative_change = 4.99056788637947e-6 Iter 80: T = 694.1114636348501 K, relative_change = 1.526834488764649e-7 Converged in 94 iterations to T = 694.110955692109 K Iter 1: T = 966.4940039421023 K, relative_change = 0.033505996057897704 Iter 2: T = 935.0278288827184 K, relative_change = 0.032557030805199796 Iter 3: T = 905.5717322266336 K, relative_change = 0.031502909053821894 Iter 20: T = 661.2095559179752 K, relative_change = 0.006946068404160954 Iter 40: T = 635.6965571885887 K, relative_change = 0.00025882747898769526 Iter 60: T = 634.8565867074809 K, relative_change = 7.97119540079656e-6 Iter 80: T = 634.8308290249769 K, relative_change = 2.4389284464927896e-7 Converged in 97 iterations to T = 634.8300581147412 K Iter 1: T = 966.4849870831456 K, relative_change = 0.03351501291685436 Iter 2: T = 935.0093865105378 K, relative_change = 0.03256708691109761 Iter 3: T = 905.5434674609077 K, relative_change = 0.03151403555380034 Iter 20: T = 661.0231567728022 K, relative_change = 0.006956009639769221 Iter 40: T = 635.4774397742057 K, relative_change = 0.00025927823652787405 Iter 60: T = 634.636291966781 K, relative_change = 7.98517155444268e-6 Iter 80: T = 634.6104980717896 K, relative_change = 2.443205585073906e-7 Converged in 97 iterations to T = 634.609726077637 K Iter 1: T = 976.4427467073914 K, relative_change = 0.023557253292608653 Iter 2: T = 955.047037089778 K, relative_change = 0.021911893646361504 Iter 3: T = 935.721184279946 K, relative_change = 0.02023549841976551 Iter 20: T = 809.876157805581 K, relative_change = 0.0020963625846471467 Iter 40: T = 801.0138465030531 K, relative_change = 6.781186776408001e-5 Iter 60: T = 800.737176449365 K, relative_change = 2.0780531987903566e-6 Iter 80: T = 800.7287076085449 K, relative_change = 6.357199113537287e-8 Converged in 89 iterations to T = 800.7284960075732 K Iter 1: T = 965.2030329842974 K, relative_change = 0.034796967015702694 Iter 2: T = 932.3817509303726 K, relative_change = 0.03400453679931479 Iter 3: T = 901.5067144906077 K, relative_change = 0.033114157810313634 Iter 20: T = 632.7901758018505 K, relative_change = 0.008638717590082985 Iter 40: T = 601.7145972134712 K, relative_change = 0.0003397872677787193 Iter 60: T = 600.6698269309826 K, relative_change = 1.0486716970476723e-5 Iter 80: T = 600.6377645305063 K, relative_change = 3.208806417794945e-7 Converged in 98 iterations to T = 600.6367965459133 K Iter 1: T = 964.5855598247807 K, relative_change = 0.03541444017521926 Iter 2: T = 931.1121051928616 K, relative_change = 0.034702421460673386 Iter 3: T = 899.5492833906078 K, relative_change = 0.03389798245155041 Iter 20: T = 617.7622398170273 K, relative_change = 0.009697550468463812 Iter 40: T = 583.1966935591897 K, relative_change = 0.00039512391923411064 Iter 60: T = 582.0183864481587 K, relative_change = 1.2212236269755607e-5 Iter 80: T = 581.9822070389886 K, relative_change = 3.736961537337229e-7 Converged in 99 iterations to T = 581.9811068213587 K Iter 1: T = 964.3150842265311 K, relative_change = 0.03568491577346893 Iter 2: T = 930.5551307524878 K, relative_change = 0.03500925581924493 Iter 3: T = 898.6891534792311 K, relative_change = 0.03424405091129684 Iter 20: T = 610.8289610094349 K, relative_change = 0.010230315731479816 Iter 40: T = 574.4986858365222 K, relative_change = 0.00042448157932198187 Iter 60: T = 573.2512725719386 K, relative_change = 1.3129705918732122e-5 Iter 80: T = 573.2129606159854 K, relative_change = 4.01780426591335e-7 Converged in 99 iterations to T = 573.2117955374913 K Iter 1: T = 980.1862371100171 K, relative_change = 0.019813762889982798 Iter 2: T = 962.4142542493108 K, relative_change = 0.018131230768047943 Iter 3: T = 946.5628357965695 K, relative_change = 0.01647047348140691 Iter 20: T = 849.4007651503271 K, relative_change = 0.0014422334448133418 Iter 40: T = 843.0571257258894 K, relative_change = 4.58323106160024e-5 Iter 60: T = 842.8603680177165 K, relative_change = 1.4037007896480831e-6 Iter 80: T = 842.8543465270575 K, relative_change = 4.2941389475393796e-8 Converged in 87 iterations to T = 842.8542125835829 K Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▌ | ETA: 0:00:02 Bin 1 progress: 80%|██████████████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010875429008884887 Iteration 10: norm(E*F-F'*E) = 4.7073324895837455e-5 Iteration 20: norm(E*F-F'*E) = 5.270698631735426e-7 Iteration 30: norm(E*F-F'*E) = 6.846201464537942e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.139332014503023e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only Spectral variation detected across walls, using spectral solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▍ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010809915448327982 Iteration 10: norm(E*F-F'*E) = 5.808705469000852e-5 Iteration 20: norm(E*F-F'*E) = 6.56399608634201e-7 Iteration 30: norm(E*F-F'*E) = 8.56169167453908e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.6761649802195253e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (5 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00:02 Bin 1 progress: 84%|███████████████████████████▉ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010809915448327982 Iteration 10: norm(E*F-F'*E) = 5.808705469000852e-5 Iteration 20: norm(E*F-F'*E) = 6.56399608634201e-7 Iteration 30: norm(E*F-F'*E) = 8.56169167453908e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.6761649802195253e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010809915448327982 Iteration 10: norm(E*F-F'*E) = 5.808705469000852e-5 Iteration 20: norm(E*F-F'*E) = 6.56399608634201e-7 Iteration 30: norm(E*F-F'*E) = 8.56169167453908e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.6761649802195253e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (20 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▍ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010809915448327982 Iteration 10: norm(E*F-F'*E) = 5.808705469000852e-5 Iteration 20: norm(E*F-F'*E) = 6.56399608634201e-7 Iteration 30: norm(E*F-F'*E) = 8.56169167453908e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.6761649802195253e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (50 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010809915448327982 Iteration 10: norm(E*F-F'*E) = 5.808705469000852e-5 Iteration 20: norm(E*F-F'*E) = 6.56399608634201e-7 Iteration 30: norm(E*F-F'*E) = 8.56169167453908e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.6761649802195253e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... ✓ 2D Spectral Participating Media tests complete ------------------------------------------------------------ Testing Spectral Consistency ------------------------------------------------------------ Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.8549284161345355e-16 === 3D Surface-Only Grey Solver === Found 96 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 96 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.8549284161345355e-16 === Using DIRECT solver === === 3D Spectral Surface Radiation Solver (DIRECT) === Spectral mode: spectral_uniform Number of spectral bins: 20 Using optimized direct emission solver Computing GERT matrices for each spectral band... (Using same view factor matrix F for all bands) Setting up boundary conditions... Starting spectral direct solve... Energy conservation errors by band: [0.0, 0.0, 0.0, 0.0, 0.0, -1.5585406229479126e-19, 4.440892098500626e-15, -1.7053025658242404e-13, 5.3859139370615594e-12, 3.325340003357269e-12, 9.876544027065393e-13, 7.638334409421077e-14, 5.662137425588298e-15, 4.926614671774132e-16, 2.0599841277224584e-17, 1.179069862577986e-18, 5.741295195608054e-20, 1.3532674821367533e-21, 9.005928919171089e-23, 3.1669026403744026e-15] === 3D Spectral Solution Complete (DIRECT) === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (20 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▌ | ETA: 0:00:02 Bin 1 progress: 80%|██████████████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/BoIOs/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010875429008884887 Iteration 10: norm(E*F-F'*E) = 4.7073324895837455e-5 Iteration 20: norm(E*F-F'*E) = 5.270698631735426e-7 Iteration 30: norm(E*F-F'*E) = 6.846201464537942e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.139332014503023e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.8549284161345355e-16 === Using DIRECT solver === === 3D Spectral Surface Radiation Solver (DIRECT) === Spectral mode: spectral_uniform Number of spectral bins: 20 Using optimized direct emission solver Computing GERT matrices for each spectral band... (Using same view factor matrix F for all bands) Setting up boundary conditions... Starting spectral direct solve... Energy conservation errors by band: [0.0, 0.0, 0.0, 0.0, 0.0, -4.633270221481023e-19, -2.4424906541753444e-15, -2.2595258997171186e-12, 1.667643800828955e-11, 4.746425474877469e-12, 9.805489753489383e-13, 1.127986593019159e-13, 5.329070518200751e-15, 4.85722573273506e-16, 2.6237692574149207e-17, 1.2841019480375193e-18, 4.438982039204568e-20, 1.4558378780933287e-21, 8.499280793982359e-23, 4.3504392608903414e-15] === 3D Spectral Solution Complete (DIRECT) === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.8549284161345355e-16 === Using FULL solver === === 3D Spectral Surface Radiation Solver === Spectral mode: spectral_uniform Number of spectral bins: 20 Computing GERT matrices for each spectral band... (Using same view factor matrix F for all bands) Building matrices for spectral bin 1... Building matrices for spectral bin 2... Building matrices for spectral bin 3... Building matrices for spectral bin 4... Building matrices for spectral bin 5... Building matrices for spectral bin 6... Building matrices for spectral bin 7... Building matrices for spectral bin 8... Building matrices for spectral bin 9... Building matrices for spectral bin 10... Building matrices for spectral bin 11... Building matrices for spectral bin 12... Building matrices for spectral bin 13... Building matrices for spectral bin 14... Building matrices for spectral bin 15... Building matrices for spectral bin 16... Building matrices for spectral bin 17... Building matrices for spectral bin 18... Building matrices for spectral bin 19... Building matrices for spectral bin 20... Assembling block matrix structure... Setting up boundary conditions... Starting spectral iteration... Iteration 1: convergence error = 1.0 Iteration 2: convergence error = 1.1845193777663423 Iteration 3: convergence error = 0.48429393675844934 Iteration 4: convergence error = 0.2209144278951869 Iteration 5: convergence error = 0.09159910997865119 Iteration 6: convergence error = 0.03576093742294504 Iteration 7: convergence error = 0.013592081912886961 Iteration 8: convergence error = 0.005109475888597634 Iteration 9: convergence error = 0.0019123763701817337 Iteration 10: convergence error = 0.0007145628006044927 Converged after 10 iterations Energy conservation errors by band: [-4.77048955893622e-16, -3.95516952522712e-16, -5.412337245047638e-16, -4.440892098500626e-16, -4.2500725161431774e-16, -3.7816971776294395e-16, 4.618527782440651e-14, 3.751665644813329e-12, 1.6619594589428743e-11, 8.44124770082999e-12, 2.1316282072803006e-13, 3.241851231905457e-14, 3.011479954295737e-15, 1.0495077029659683e-16, -1.734723475976807e-16, 1.3183898417423734e-16, -5.0306980803327406e-17, -3.122502256758253e-17, -1.8735013540549517e-16, 7.355227538141662e-16] Writing spectral results to mesh... === 3D Spectral Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 54×54 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 6.233122709088313e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.993453929734661e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 54×54 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 6.233122709088313e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.993453929734661e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 2.220446049250313e-16 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 5 iterations. norm(E*F-F'*E) = 1.8549284161345355e-16 === Using DIRECT solver === === 3D Spectral Surface Radiation Solver (DIRECT) === Spectral mode: spectral_uniform Number of spectral bins: 20 Using optimized direct emission solver Computing GERT matrices for each spectral band... (Using same view factor matrix F for all bands) Setting up boundary conditions... Starting spectral direct solve... Energy conservation errors by band: [0.0, 0.0, 0.0, 0.0, 0.0, 7.352245982167327e-19, 9.126033262418787e-14, 4.689582056016661e-13, 4.728661906483467e-12, 4.092726157978177e-12, 5.258016244624741e-13, 8.79296635503124e-14, 7.049916206369744e-15, 4.683753385137379e-16, 3.133344278483108e-17, 1.3061248046661311e-18, 4.828087799349512e-20, 1.828069153742191e-21, 8.88055935758102e-23, 4.8870001515755484e-15] === 3D Spectral Solution Complete (DIRECT) === ✓ Spectral Consistency tests complete ================================================================================ TEST SUITE COMPLETE ================================================================================ Test Summary: | Pass Total Time RayTraceHeatTransfer.jl | 1394 1394 8m55.2s Testing RayTraceHeatTransfer tests passed Testing completed after 549.64s PkgEval succeeded after 623.13s