Package evaluation to test RayTraceHeatTransfer on Julia 1.14.0-DEV.2302 (b8ba835509*) started at 2026-06-06T12:24:50.947 ################################################################################ # Set-up # Installing PkgEval dependencies (TestEnv)... Activating project at `~/.julia/environments/v1.14` Set-up completed after 16.17s ################################################################################ # Installation # Installing RayTraceHeatTransfer... Resolving package versions... Updating `~/.julia/environments/v1.14/Project.toml` [7cf1493d] + RayTraceHeatTransfer v0.9.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.5 [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.8.0 ⌅ [2ab3a3ac] + LogExpFunctions v0.3.29 [eff96d63] + Measurements v2.14.1 [e1d29d7a] + Missings v1.2.0 ⌅ [bac558e1] + OrderedCollections v1.8.2 [aea7be01] + PrecompileTools v1.3.4 [21216c6a] + Preferences v1.5.2 [92933f4c] + ProgressMeter v1.11.0 [43287f4e] + PtrArrays v1.4.0 [7cf1493d] + RayTraceHeatTransfer v0.9.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.14.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.13.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.5.2+0 [deac9b47] + LibCURL_jll v8.20.0+1 [e37daf67] + LibGit2_jll v1.9.4+0 [29816b5a] + LibSSH2_jll v1.11.101+0 [14a3606d] + MozillaCACerts_jll v2026.5.14 [4536629a] + OpenBLAS_jll v0.3.33+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 Info Packages marked with ⌅ have new versions available but compatibility constraints restrict them from upgrading. To see why use `status --outdated -m` Installation completed after 5.88s ################################################################################ # Precompilation # Precompiling PkgEval dependencies... Precompiling package dependencies... Precompiling project... 5.5 s ✓ StatsBase 12.4 s ✓ RayTraceHeatTransfer 2 dependencies successfully precompiled in 18 seconds. 59 already precompiled. Precompilation completed after 47.38s ################################################################################ # Testing # Testing RayTraceHeatTransfer Status `/tmp/jl_1bIiEz/Project.toml` [5c1252a2] GeometryBasics v0.5.10 [eff96d63] Measurements v2.14.1 [92933f4c] ProgressMeter v1.11.0 [7cf1493d] RayTraceHeatTransfer v0.9.0 [90137ffa] StaticArrays v1.9.18 ⌅ [2913bbd2] StatsBase v0.34.10 [37e2e46d] LinearAlgebra v1.14.0 [9a3f8284] Random v1.11.0 [2f01184e] SparseArrays v1.13.0 [8dfed614] Test v1.11.0 Status `/tmp/jl_1bIiEz/Manifest.toml` [66dad0bd] AliasTables v1.1.3 [49dc2e85] Calculus v0.5.2 [9a962f9c] DataAPI v1.16.0 [864edb3b] DataStructures v0.19.5 [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.8.0 ⌅ [2ab3a3ac] LogExpFunctions v0.3.29 [eff96d63] Measurements v2.14.1 [e1d29d7a] Missings v1.2.0 ⌅ [bac558e1] OrderedCollections v1.8.2 [aea7be01] PrecompileTools v1.3.4 [21216c6a] Preferences v1.5.2 [92933f4c] ProgressMeter v1.11.0 [43287f4e] PtrArrays v1.4.0 [7cf1493d] RayTraceHeatTransfer v0.9.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.14.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.13.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.5.2+0 [deac9b47] LibCURL_jll v8.20.0+1 [e37daf67] LibGit2_jll v1.9.4+0 [29816b5a] LibSSH2_jll v1.11.101+0 [14a3606d] MozillaCACerts_jll v2026.5.14 [4536629a] OpenBLAS_jll v0.3.33+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 Info Packages marked with ⌅ have new versions available but compatibility constraints restrict them from upgrading. Testing Running tests... ================================================================================ STARTING TEST SUITE ================================================================================ ------------------------------------------------------------ Testing 3D View Factors ------------------------------------------------------------ Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 1.8906416838689214e-15 Converged after 3 iterations. norm(E*F-F'*E) = 8.345155887641791e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 1.461852203983726e-15 Converged after 3 iterations. norm(E*F-F'*E) = 7.119724442439612e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 1.9794221072983377e-15 Converged after 3 iterations. norm(E*F-F'*E) = 4.440892098500626e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 2.8318849355912135e-15 Converged after 3 iterations. norm(E*F-F'*E) = 2.633125101432526e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 1.962615573354719e-15 Converged after 3 iterations. norm(E*F-F'*E) = 8.233634315563857e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 2.211058699092487e-15 Converged after 3 iterations. norm(E*F-F'*E) = 4.388541835720876e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 6×6 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 2.0206364052201327e-15 Converged after 3 iterations. norm(E*F-F'*E) = 7.529898907871222e-16 ✓ 3D View Factor tests complete ------------------------------------------------------------ Testing 3D Heat Transfer ------------------------------------------------------------ Computing view factors (geometry only, wavelength-independent)... Matrix size: 150×150 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 1.0431932984236912e-14 Converged after 3 iterations. norm(E*F-F'*E) = 7.560517450615101e-16 === 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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 1.0431932984236912e-14 Converged after 3 iterations. norm(E*F-F'*E) = 7.560517450615101e-16 === 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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 1.0431932984236912e-14 Converged after 3 iterations. norm(E*F-F'*E) = 7.560517450615101e-16 === 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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.604788474939644e-15 Converged after 3 iterations. norm(E*F-F'*E) = 9.862222410840723e-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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.836202745879726e-15 Converged after 3 iterations. norm(E*F-F'*E) = 8.677919832028415e-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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.57995603044491e-15 Converged after 3 iterations. norm(E*F-F'*E) = 8.051026671140432e-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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-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:06:35 Bin 1 progress: 36%|███████████▊ | ETA: 0:00:12 Bin 1 progress: 76%|█████████████████████████ | ETA: 0:00:02 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:08 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012724132541177872 Iteration 10: norm(E*F-F'*E) = 2.2562390371091414e-5 Iteration 20: norm(E*F-F'*E) = 2.81747934726781e-7 Converged after 29 iterations. norm(E*F-F'*E) = 4.51230822153629e-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: 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012615028774925617 Iteration 10: norm(E*F-F'*E) = 2.929368274539419e-5 Iteration 20: norm(E*F-F'*E) = 4.31266418626085e-7 Iteration 30: norm(E*F-F'*E) = 7.361062995897725e-9 Converged after 31 iterations. norm(E*F-F'*E) = 3.266672926912583e-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: 42%|██████████████ | ETA: 0:00:01 Bin 1 progress: 85%|████████████████████████████▎ | ETA: 0:00:00 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012528556424676958 Iteration 10: norm(E*F-F'*E) = 4.019144026273758e-5 Iteration 20: norm(E*F-F'*E) = 6.59631236691233e-7 Iteration 30: norm(E*F-F'*E) = 1.1539917271474175e-8 Converged after 31 iterations. norm(E*F-F'*E) = 5.138695639811133e-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: 41%|█████████████▋ | ETA: 0:00:01 Bin 1 progress: 83%|███████████████████████████▍ | ETA: 0:00:00 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.01275138904213884 Iteration 10: norm(E*F-F'*E) = 4.377660723551281e-5 Iteration 20: norm(E*F-F'*E) = 6.784475087231262e-7 Iteration 30: norm(E*F-F'*E) = 1.1471219797925177e-8 Converged after 31 iterations. norm(E*F-F'*E) = 5.083601712094415e-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: 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011200057001712635 Iteration 10: norm(E*F-F'*E) = 3.1072252836633964e-5 Iteration 20: norm(E*F-F'*E) = 3.015890404132241e-7 Iteration 30: norm(E*F-F'*E) = 3.795079769888138e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.5991191034806917e-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: 42%|█████████████▊ | ETA: 0:00:01 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.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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010941659664567473 Iteration 10: norm(E*F-F'*E) = 3.3723607708257e-5 Iteration 20: norm(E*F-F'*E) = 3.5014924661516674e-7 Iteration 30: norm(E*F-F'*E) = 4.386648472399213e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.8447360422845952e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011106413130397605 Iteration 10: norm(E*F-F'*E) = 2.2801204185807325e-5 Iteration 20: norm(E*F-F'*E) = 2.3018303265855515e-7 Iteration 30: norm(E*F-F'*E) = 3.2307290316027986e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.3870990665504233e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011276241084731062 Iteration 10: norm(E*F-F'*E) = 2.668275423291254e-5 Iteration 20: norm(E*F-F'*E) = 2.7998816089086246e-7 Iteration 30: norm(E*F-F'*E) = 3.83088859799256e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.6391301094927685e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011509825487293577 Iteration 10: norm(E*F-F'*E) = 3.32781048006822e-5 Iteration 20: norm(E*F-F'*E) = 3.991620602778427e-7 Iteration 30: norm(E*F-F'*E) = 5.7834117129448705e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.0742383324176481e-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: 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.14285714285714296) ≥ min gas E (0.08163265306122436); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011482657368106119 Iteration 10: norm(E*F-F'*E) = 2.5302006279955488e-5 Iteration 20: norm(E*F-F'*E) = 2.0229725209180965e-7 Iteration 30: norm(E*F-F'*E) = 2.204377269365422e-9 Converged after 31 iterations. norm(E*F-F'*E) = 9.077873481705735e-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 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.03522342402436201 Iteration 10: norm(E*F-F'*E) = 0.00011925112490717111 Iteration 20: norm(E*F-F'*E) = 1.2558424326117652e-6 Iteration 30: norm(E*F-F'*E) = 1.655667585206054e-8 Converged after 35 iterations. norm(E*F-F'*E) = 1.261109740580851e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.028327750958833154 Iteration 10: norm(E*F-F'*E) = 6.60750929422855e-5 Iteration 20: norm(E*F-F'*E) = 8.394578057822155e-7 Iteration 30: norm(E*F-F'*E) = 1.2285095251431934e-8 Converged after 33 iterations. norm(E*F-F'*E) = 2.284426945864392e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.024442576933276162 Iteration 10: norm(E*F-F'*E) = 8.482844031890458e-5 Iteration 20: norm(E*F-F'*E) = 1.1653785657557274e-6 Iteration 30: norm(E*F-F'*E) = 1.893396438659578e-8 Converged after 33 iterations. norm(E*F-F'*E) = 3.685388544672953e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.021783316432619257 Iteration 10: norm(E*F-F'*E) = 2.81990825998982e-5 Iteration 20: norm(E*F-F'*E) = 2.209346593518276e-7 Converged after 27 iterations. norm(E*F-F'*E) = 7.371114547988138e-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: 43%|██████████████▏ | ETA: 0:00:01 Bin 1 progress: 88%|█████████████████████████████▏ | 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011200057001712635 Iteration 10: norm(E*F-F'*E) = 3.1072252836633964e-5 Iteration 20: norm(E*F-F'*E) = 3.015890404132241e-7 Iteration 30: norm(E*F-F'*E) = 3.795079769888138e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.5991191034806917e-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: 42%|█████████████▉ | ETA: 0:00:01 Bin 1 progress: 89%|█████████████████████████████▍ | ETA: 0:00:00 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011165403270466182 Iteration 10: norm(E*F-F'*E) = 2.2924394455790618e-5 Iteration 20: norm(E*F-F'*E) = 1.563534137504874e-7 Iteration 30: norm(E*F-F'*E) = 1.8311879498287595e-9 Converged after 31 iterations. norm(E*F-F'*E) = 7.739436362847707e-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: 44%|██████████████▋ | ETA: 0:00:01 Bin 1 progress: 89%|█████████████████████████████▍ | 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.009973478231560627 Iteration 10: norm(E*F-F'*E) = 2.6154762054020128e-5 Iteration 20: norm(E*F-F'*E) = 2.6243570446441545e-7 Iteration 30: norm(E*F-F'*E) = 3.45748879999977e-9 Converged after 33 iterations. norm(E*F-F'*E) = 6.2302945385852e-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: 38%|████████████▌ | 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011165403270466182 Iteration 10: norm(E*F-F'*E) = 2.2924394455790618e-5 Iteration 20: norm(E*F-F'*E) = 1.563534137504874e-7 Iteration 30: norm(E*F-F'*E) = 1.8311879498287595e-9 Converged after 31 iterations. norm(E*F-F'*E) = 7.739436362847707e-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/VDrcq/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/VDrcq/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/VDrcq/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 4/10 Bin 4 ray tracing: 0%| | ETA: 13:22:27 Bin 4 ray tracing: 8%|██▌ | ETA: 0:01:05 Bin 4 ray tracing: 17%|█████▏ | ETA: 0:00:33 Bin 4 ray tracing: 26%|███████▊ | ETA: 0:00:23 Bin 4 ray tracing: 34%|██████████▎ | ETA: 0:00:17 Bin 4 ray tracing: 42%|████████████▋ | ETA: 0:00:14 Bin 4 ray tracing: 50%|███████████████ | ETA: 0:00:11 Bin 4 ray tracing: 58%|█████████████████▍ | ETA: 0:00:09 Bin 4 ray tracing: 65%|███████████████████▋ | ETA: 0:00:07 Bin 4 ray tracing: 73%|██████████████████████ | ETA: 0:00:05 Bin 4 ray tracing: 81%|████████████████████████▍ | ETA: 0:00:03 Bin 4 ray tracing: 89%|██████████████████████████▊ | ETA: 0:00:02 Bin 4 ray tracing: 97%|█████████████████████████████▏| ETA: 0:00:01 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:17 Updating spectral results for spectral bin 4 Energy per ray: 0.0001853335835185918 Processing spectral bin 5/10 Bin 5 ray tracing: 8%|██▍ | ETA: 0:00:12 Bin 5 ray tracing: 16%|████▊ | ETA: 0:00:11 Bin 5 ray tracing: 25%|███████▍ | ETA: 0:00:09 Bin 5 ray tracing: 33%|██████████ | ETA: 0:00:08 Bin 5 ray tracing: 42%|████████████▌ | ETA: 0:00:07 Bin 5 ray tracing: 50%|██████████████▉ | ETA: 0:00:06 Bin 5 ray tracing: 58%|█████████████████▎ | ETA: 0:00:05 Bin 5 ray tracing: 66%|███████████████████▉ | ETA: 0:00:04 Bin 5 ray tracing: 75%|██████████████████████▌ | ETA: 0:00:03 Bin 5 ray tracing: 84%|█████████████████████████▏ | ETA: 0:00:02 Bin 5 ray tracing: 93%|███████████████████████████▊ | ETA: 0:00:01 Bin 5 ray tracing: 100%|██████████████████████████████| Time: 0:00:12 Updating spectral results for spectral bin 5 Energy per ray: 0.04303963948070305 Processing spectral bin 6/10 Bin 6 ray tracing: 9%|██▋ | ETA: 0:00:10 Bin 6 ray tracing: 18%|█████▎ | ETA: 0:00:09 Bin 6 ray tracing: 26%|███████▉ | ETA: 0:00:08 Bin 6 ray tracing: 35%|██████████▌ | ETA: 0:00:07 Bin 6 ray tracing: 44%|█████████████ | ETA: 0:00:07 Bin 6 ray tracing: 52%|███████████████▌ | ETA: 0:00:06 Bin 6 ray tracing: 60%|██████████████████ | ETA: 0:00:05 Bin 6 ray tracing: 68%|████████████████████▌ | ETA: 0:00:04 Bin 6 ray tracing: 76%|██████████████████████▉ | ETA: 0:00:03 Bin 6 ray tracing: 85%|█████████████████████████▍ | ETA: 0:00:02 Bin 6 ray tracing: 93%|███████████████████████████▉ | 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:11 Bin 7 ray tracing: 17%|█████ | ETA: 0:00:10 Bin 7 ray tracing: 26%|███████▊ | ETA: 0:00:09 Bin 7 ray tracing: 35%|██████████▌ | ETA: 0:00:08 Bin 7 ray tracing: 43%|████████████▉ | ETA: 0:00:07 Bin 7 ray tracing: 51%|███████████████▎ | ETA: 0:00:06 Bin 7 ray tracing: 58%|█████████████████▌ | ETA: 0:00:05 Bin 7 ray tracing: 66%|███████████████████▉ | ETA: 0:00:04 Bin 7 ray tracing: 74%|██████████████████████▎ | ETA: 0:00:03 Bin 7 ray tracing: 83%|████████████████████████▊ | ETA: 0:00:02 Bin 7 ray tracing: 91%|███████████████████████████▍ | ETA: 0:00:01 Bin 7 ray tracing: 99%|██████████████████████████████| ETA: 0:00:00 Bin 7 ray tracing: 100%|██████████████████████████████| Time: 0:00:12 Updating spectral results for spectral bin 7 Energy per ray: 0.00021661482457376898 Processing spectral bin 8/10 Bin 8 ray tracing: 9%|██▋ | ETA: 0:00:10 Bin 8 ray tracing: 18%|█████▍ | ETA: 0:00:09 Bin 8 ray tracing: 27%|████████ | ETA: 0:00:08 Bin 8 ray tracing: 35%|██████████▋ | ETA: 0:00:08 Bin 8 ray tracing: 44%|█████████████ | ETA: 0:00:07 Bin 8 ray tracing: 52%|███████████████▌ | ETA: 0:00:06 Bin 8 ray tracing: 60%|██████████████████ | ETA: 0:00:05 Bin 8 ray tracing: 68%|████████████████████▌ | ETA: 0:00:04 Bin 8 ray tracing: 77%|███████████████████████ | ETA: 0:00:03 Bin 8 ray tracing: 85%|█████████████████████████▌ | ETA: 0:00:02 Bin 8 ray tracing: 93%|███████████████████████████▉ | ETA: 0:00:01 Bin 8 ray tracing: 100%|██████████████████████████████| Time: 0:00:12 Updating spectral results for spectral bin 8 Energy per ray: 1.0195075180910974e-6 Processing spectral bin 9/10 Bin 9 ray tracing: 8%|██▍ | ETA: 0:00:12 Bin 9 ray tracing: 16%|████▉ | ETA: 0:00:11 Bin 9 ray tracing: 24%|███████▏ | ETA: 0:00:10 Bin 9 ray tracing: 32%|█████████▌ | ETA: 0:00:09 Bin 9 ray tracing: 40%|███████████▉ | ETA: 0:00:08 Bin 9 ray tracing: 48%|██████████████▎ | ETA: 0:00:07 Bin 9 ray tracing: 56%|████████████████▋ | ETA: 0:00:06 Bin 9 ray tracing: 64%|███████████████████▏ | ETA: 0:00:05 Bin 9 ray tracing: 72%|█████████████████████▌ | ETA: 0:00:04 Bin 9 ray tracing: 80%|████████████████████████ | ETA: 0:00:03 Bin 9 ray tracing: 88%|██████████████████████████▍ | ETA: 0:00:02 Bin 9 ray tracing: 96%|████████████████████████████▊ | ETA: 0:00:01 Bin 9 ray tracing: 100%|██████████████████████████████| Time: 0:00:12 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:12 Bin 10 ray tracing: 16%|████▋ | ETA: 0:00:11 Bin 10 ray tracing: 24%|██████▉ | ETA: 0:00:10 Bin 10 ray tracing: 32%|█████████▏ | ETA: 0:00:09 Bin 10 ray tracing: 40%|███████████▌ | ETA: 0:00:08 Bin 10 ray tracing: 48%|██████████████ | ETA: 0:00:06 Bin 10 ray tracing: 57%|████████████████▌ | ETA: 0:00:05 Bin 10 ray tracing: 65%|██████████████████▉ | ETA: 0:00:04 Bin 10 ray tracing: 73%|█████████████████████▎ | ETA: 0:00:03 Bin 10 ray tracing: 82%|███████████████████████▉ | ETA: 0:00:02 Bin 10 ray tracing: 91%|██████████████████████████▍ | ETA: 0:00:01 Bin 10 ray tracing: 99%|████████████████████████████▉| ETA: 0:00:00 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: 24%|████████▏ | ETA: 0:00:03 Bin 1 progress: 47%|███████████████▍ | ETA: 0:00:02 Bin 1 progress: 71%|███████████████████████▌ | 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: 22%|███████▍ | ETA: 0:00:04 Bin 2 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 2 progress: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 2 progress: 91%|██████████████████████████████▏ | 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: 67%|██████████████████████ | ETA: 0:00:02 Bin 3 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 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: 22%|███████▍ | ETA: 0:00:04 Bin 4 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 4 progress: 67%|██████████████████████ | ETA: 0:00:02 Bin 4 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 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: 67%|██████████████████████ | ETA: 0:00:02 Bin 5 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 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: 22%|███████▍ | ETA: 0:00:04 Bin 6 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 6 progress: 67%|██████████████████████ | ETA: 0:00:02 Bin 6 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 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: 49%|████████████████▏ | ETA: 0:00:02 Bin 7 progress: 73%|████████████████████████▎ | 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: 22%|███████▍ | ETA: 0:00:04 Bin 8 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 8 progress: 69%|██████████████████████▊ | 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: 98%|████████████████████████████████▎| 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: 24%|███████▉ | ETA: 0:00:03 Bin 10 progress: 49%|███████████████▋ | ETA: 0:00:02 Bin 10 progress: 73%|███████████████████████▌ | ETA: 0:00:01 Bin 10 progress: 98%|███████████████████████████████▎| 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011165403270466182 Iteration 10: norm(E*F-F'*E) = 2.2924394455790618e-5 Iteration 20: norm(E*F-F'*E) = 1.563534137504874e-7 Iteration 30: norm(E*F-F'*E) = 1.8311879498287595e-9 Converged after 31 iterations. norm(E*F-F'*E) = 7.739436362847707e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.009991465127794058 Iteration 10: norm(E*F-F'*E) = 2.6355781212038057e-5 Iteration 20: norm(E*F-F'*E) = 2.67975051901387e-7 Iteration 30: norm(E*F-F'*E) = 3.5435082314192637e-9 Converged after 33 iterations. norm(E*F-F'*E) = 6.386069733597781e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010721910152575703 Iteration 10: norm(E*F-F'*E) = 3.6940844228855025e-5 Iteration 20: norm(E*F-F'*E) = 3.838874943073066e-7 Iteration 30: norm(E*F-F'*E) = 4.87260052301876e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.6345596357894995e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010405506392554593 Iteration 10: norm(E*F-F'*E) = 2.5567548963130732e-5 Iteration 20: norm(E*F-F'*E) = 2.2800740030811457e-7 Iteration 30: norm(E*F-F'*E) = 2.8960766124378582e-9 Converged after 33 iterations. norm(E*F-F'*E) = 5.173308438692455e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011083879712045256 Iteration 10: norm(E*F-F'*E) = 5.2729891561736435e-5 Iteration 20: norm(E*F-F'*E) = 6.134766235863196e-7 Iteration 30: norm(E*F-F'*E) = 7.98202150371361e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.984161257437805e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010826922109238838 Iteration 10: norm(E*F-F'*E) = 2.3417898886011587e-5 Iteration 20: norm(E*F-F'*E) = 1.8596397500744126e-7 Iteration 30: norm(E*F-F'*E) = 2.270200412202148e-9 Converged after 33 iterations. norm(E*F-F'*E) = 4.019271607317107e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010436725003758788 Iteration 10: norm(E*F-F'*E) = 3.474615102453393e-5 Iteration 20: norm(E*F-F'*E) = 3.084311695458338e-7 Iteration 30: norm(E*F-F'*E) = 3.6297133653161557e-9 Converged after 33 iterations. norm(E*F-F'*E) = 6.351178400396532e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011246339829203767 Iteration 10: norm(E*F-F'*E) = 3.771739455462875e-5 Iteration 20: norm(E*F-F'*E) = 4.6932802009431876e-7 Iteration 30: norm(E*F-F'*E) = 6.4517353899695215e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.949036050956758e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010945033755620977 Iteration 10: norm(E*F-F'*E) = 4.492756334815193e-5 Iteration 20: norm(E*F-F'*E) = 5.761615955419631e-7 Iteration 30: norm(E*F-F'*E) = 7.966875368232075e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.146771540616354e-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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.00985595211116036 Iteration 10: norm(E*F-F'*E) = 4.756274758097569e-5 Iteration 20: norm(E*F-F'*E) = 5.753161224317444e-7 Iteration 30: norm(E*F-F'*E) = 7.840831799907705e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.014953234516988e-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/VDrcq/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/VDrcq/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/VDrcq/src/RayTracing/RayTracing2D/DirectTracing2D/directRayTracing.jl:25 Processing spectral bin 4/10 Bin 4 ray tracing: 6%|█▊ | ETA: 0:00:16 Bin 4 ray tracing: 13%|███▊ | ETA: 0:00:14 Bin 4 ray tracing: 19%|█████▉ | ETA: 0:00:13 Bin 4 ray tracing: 26%|███████▋ | ETA: 0:00:12 Bin 4 ray tracing: 31%|█████████▌ | ETA: 0:00:11 Bin 4 ray tracing: 37%|███████████▎ | ETA: 0:00:10 Bin 4 ray tracing: 43%|█████████████ | ETA: 0:00:09 Bin 4 ray tracing: 49%|██████████████▊ | ETA: 0:00:08 Bin 4 ray tracing: 55%|████████████████▌ | ETA: 0:00:08 Bin 4 ray tracing: 61%|██████████████████▎ | ETA: 0:00:07 Bin 4 ray tracing: 67%|████████████████████▎ | ETA: 0:00:05 Bin 4 ray tracing: 74%|██████████████████████▏ | ETA: 0:00:04 Bin 4 ray tracing: 80%|████████████████████████ | ETA: 0:00:03 Bin 4 ray tracing: 86%|█████████████████████████▊ | ETA: 0:00:02 Bin 4 ray tracing: 93%|███████████████████████████▊ | ETA: 0:00:01 Bin 4 ray tracing: 99%|█████████████████████████████▊| ETA: 0:00:00 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:15 Bin 5 ray tracing: 13%|███▉ | ETA: 0:00:13 Bin 5 ray tracing: 20%|█████▉ | ETA: 0:00:12 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: 44%|█████████████▎ | ETA: 0:00:09 Bin 5 ray tracing: 50%|███████████████▏ | ETA: 0:00:08 Bin 5 ray tracing: 57%|█████████████████ | ETA: 0:00:07 Bin 5 ray tracing: 63%|██████████████████▉ | ETA: 0:00:06 Bin 5 ray tracing: 69%|████████████████████▉ | ETA: 0:00:05 Bin 5 ray tracing: 76%|██████████████████████▋ | ETA: 0:00:04 Bin 5 ray tracing: 82%|████████████████████████▌ | ETA: 0:00:03 Bin 5 ray tracing: 88%|██████████████████████████▍ | ETA: 0:00:02 Bin 5 ray tracing: 94%|████████████████████████████▎ | ETA: 0:00:01 Bin 5 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 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: 12%|███▊ | ETA: 0:00:15 Bin 6 ray tracing: 18%|█████▌ | ETA: 0:00:14 Bin 6 ray tracing: 24%|███████▍ | ETA: 0:00:13 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: 50%|██████████████▉ | ETA: 0:00:08 Bin 6 ray tracing: 55%|████████████████▋ | ETA: 0:00:07 Bin 6 ray tracing: 61%|██████████████████▍ | ETA: 0:00:06 Bin 6 ray tracing: 67%|████████████████████▏ | ETA: 0:00:06 Bin 6 ray tracing: 73%|██████████████████████ | ETA: 0:00:05 Bin 6 ray tracing: 79%|███████████████████████▊ | ETA: 0:00:03 Bin 6 ray tracing: 85%|█████████████████████████▋ | ETA: 0:00:02 Bin 6 ray tracing: 91%|███████████████████████████▍ | ETA: 0:00:01 Bin 6 ray tracing: 97%|█████████████████████████████▏| ETA: 0:00:01 Bin 6 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 Updating spectral results for spectral bin 6 Energy per ray: 0.013246116789219256 Processing spectral bin 7/10 Bin 7 ray tracing: 7%|██ | ETA: 0:00:15 Bin 7 ray tracing: 13%|████ | ETA: 0:00:14 Bin 7 ray tracing: 20%|█████▉ | ETA: 0:00:13 Bin 7 ray tracing: 26%|███████▊ | 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: 45%|█████████████▍ | ETA: 0:00:09 Bin 7 ray tracing: 51%|███████████████▍ | ETA: 0:00:08 Bin 7 ray tracing: 58%|█████████████████▎ | ETA: 0:00:07 Bin 7 ray tracing: 64%|███████████████████▎ | ETA: 0:00:06 Bin 7 ray tracing: 71%|█████████████████████▎ | ETA: 0:00:05 Bin 7 ray tracing: 78%|███████████████████████▎ | ETA: 0:00:04 Bin 7 ray tracing: 84%|█████████████████████████▎ | ETA: 0:00:02 Bin 7 ray tracing: 91%|███████████████████████████▎ | ETA: 0:00:01 Bin 7 ray tracing: 97%|█████████████████████████████▏| ETA: 0:00:00 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: 6%|█▉ | ETA: 0:00:15 Bin 8 ray tracing: 12%|███▊ | ETA: 0:00:14 Bin 8 ray tracing: 19%|█████▋ | ETA: 0:00:13 Bin 8 ray tracing: 25%|███████▍ | ETA: 0:00:12 Bin 8 ray tracing: 31%|█████████▎ | ETA: 0:00:11 Bin 8 ray tracing: 37%|███████████▎ | ETA: 0:00:10 Bin 8 ray tracing: 44%|█████████████▏ | ETA: 0:00:09 Bin 8 ray tracing: 50%|██████████████▉ | ETA: 0:00:08 Bin 8 ray tracing: 56%|████████████████▊ | 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: 88%|██████████████████████████▎ | ETA: 0:00:02 Bin 8 ray tracing: 94%|████████████████████████████▏ | ETA: 0:00:01 Bin 8 ray tracing: 99%|█████████████████████████████▉| ETA: 0:00:00 Bin 8 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 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:16 Bin 9 ray tracing: 12%|███▌ | ETA: 0:00:15 Bin 9 ray tracing: 18%|█████▎ | ETA: 0:00:14 Bin 9 ray tracing: 24%|███████▏ | ETA: 0:00:13 Bin 9 ray tracing: 30%|█████████ | ETA: 0:00:12 Bin 9 ray tracing: 36%|██████████▉ | ETA: 0:00:11 Bin 9 ray tracing: 43%|████████████▉ | ETA: 0:00:09 Bin 9 ray tracing: 49%|██████████████▊ | ETA: 0:00:08 Bin 9 ray tracing: 55%|████████████████▋ | ETA: 0:00:07 Bin 9 ray tracing: 61%|██████████████████▍ | ETA: 0:00:06 Bin 9 ray tracing: 67%|████████████████████▏ | ETA: 0:00:05 Bin 9 ray tracing: 73%|█████████████████████▉ | ETA: 0:00:04 Bin 9 ray tracing: 79%|███████████████████████▊ | ETA: 0:00:04 Bin 9 ray tracing: 85%|█████████████████████████▍ | ETA: 0:00:03 Bin 9 ray tracing: 91%|███████████████████████████▎ | ETA: 0:00:02 Bin 9 ray tracing: 97%|█████████████████████████████ | ETA: 0:00:01 Bin 9 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 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: 12%|███▋ | 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: 31%|█████████ | ETA: 0:00:11 Bin 10 ray tracing: 37%|██████████▊ | ETA: 0:00:10 Bin 10 ray tracing: 44%|████████████▋ | ETA: 0:00:09 Bin 10 ray tracing: 50%|██████████████▋ | ETA: 0:00:08 Bin 10 ray tracing: 57%|████████████████▌ | ETA: 0:00:07 Bin 10 ray tracing: 63%|██████████████████▍ | ETA: 0:00:06 Bin 10 ray tracing: 70%|████████████████████▎ | ETA: 0:00:05 Bin 10 ray tracing: 76%|██████████████████████▏ | ETA: 0:00:04 Bin 10 ray tracing: 83%|████████████████████████ | ETA: 0:00:03 Bin 10 ray tracing: 89%|█████████████████████████▉ | 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.2859184941946 K, relative_change = 0.03271408150580549 Iter 2: T = 936.6454008959283 K, relative_change = 0.031676794846724705 Iter 3: T = 908.0471835537227 K, relative_change = 0.03053259783782707 Iter 20: T = 677.0099036280781 K, relative_change = 0.006152326593061095 Iter 40: T = 654.1155197464349 K, relative_change = 0.00022371210512771091 Iter 60: T = 653.3687856412091 K, relative_change = 6.883421214514681e-6 Iter 80: T = 653.3458946383026 K, relative_change = 2.1060451801942432e-7 Converged in 96 iterations to T = 653.3452166345195 K Iter 1: T = 970.298397280664 K, relative_change = 0.02970160271933601 Iter 2: T = 942.7602186323768 K, relative_change = 0.02838114411552677 Iter 3: T = 917.3409796730942 K, relative_change = 0.02696257060587174 Iter 20: T = 728.9201006094092 K, relative_change = 0.004110055516895025 Iter 40: T = 712.8830743203098 K, relative_change = 0.00014067803840286578 Iter 60: T = 712.3718196018052 K, relative_change = 4.319175783286015e-6 Iter 80: T = 712.3561594314646 K, relative_change = 1.3214030430548678e-7 Converged in 93 iterations to T = 712.3557164730535 K Iter 1: T = 974.4306972309998 K, relative_change = 0.02556930276900016 Iter 2: T = 951.0505148190365 K, relative_change = 0.02399368418749718 Iter 3: T = 929.7845866832826 K, relative_change = 0.02236046119989789 Iter 20: T = 786.0206172294432 K, relative_change = 0.0025802839982264417 Iter 40: T = 775.3704478719526 K, relative_change = 8.458395159306299e-5 Iter 60: T = 775.036329607292 K, relative_change = 2.593154416534455e-6 Iter 80: T = 775.0261006724717 K, relative_change = 7.93310750512283e-8 Converged in 90 iterations to T = 775.0258343392918 K Iter 1: T = 970.3776445499944 K, relative_change = 0.0296223554500056 Iter 2: T = 942.9202606006995 K, relative_change = 0.028295565240508094 Iter 3: T = 917.5828859217571 K, relative_change = 0.02687117430566277 Iter 20: T = 730.1456301590321 K, relative_change = 0.00407046700651691 Iter 40: T = 714.2443804376179 K, relative_change = 0.0001391643029619674 Iter 60: T = 713.7376703142886 K, relative_change = 4.272531760084931e-6 Iter 80: T = 713.7221495697991 K, relative_change = 1.30713124809035e-7 Converged in 93 iterations to T = 713.7217105553328 K Iter 1: T = 969.2839399503046 K, relative_change = 0.03071606004969541 Iter 2: T = 940.7078130639812 K, relative_change = 0.02948168819116967 Iter 3: T = 914.2327289109844 K, relative_change = 0.028143791074472676 Iter 20: T = 712.6651420023269 K, relative_change = 0.00466909819615821 Iter 40: T = 694.7252558129303 K, relative_change = 0.00016242695213545253 Iter 60: T = 694.1498472735486 K, relative_change = 4.98974989541231e-6 Iter 80: T = 694.1322184090874 K, relative_change = 1.5265841967934237e-7 Converged in 94 iterations to T = 694.1317105344273 K Iter 1: T = 963.5306861461851 K, relative_change = 0.03646931385381493 Iter 2: T = 928.9370185395686 K, relative_change = 0.03590302634260677 Iter 3: T = 896.1853229249781 K, relative_change = 0.035257175632941416 Iter 20: T = 589.306237261358 K, relative_change = 0.012087307966267717 Iter 40: T = 546.7524582785368 K, relative_change = 0.0005356977110969616 Iter 60: T = 545.2522526460749 K, relative_change = 1.661821733397025e-5 Iter 80: T = 545.2061278786306 K, relative_change = 5.085779829491705e-7 Iter 100: T = 545.2047166820769 K, relative_change = 1.5557829608328555e-8 Converged in 101 iterations to T = 545.2047095570266 K Iter 1: T = 966.928095424127 K, relative_change = 0.03307190457587294 Iter 2: T = 935.9150336102338 K, relative_change = 0.03207380358545678 Iter 3: T = 906.9303552482596 K, relative_change = 0.030969347986822608 Iter 20: T = 670.0060643055942 K, relative_change = 0.0064925391683114475 Iter 40: T = 645.9875484913892 K, relative_change = 0.00023855455743247185 Iter 60: T = 645.2010267790273 K, relative_change = 7.342954771772872e-6 Iter 80: T = 645.1769127175107 K, relative_change = 2.2466704480053109e-7 Converged in 96 iterations to T = 645.1761984851418 K Iter 1: T = 965.1907385927628 K, relative_change = 0.034809261407237176 Iter 2: T = 932.3564967234271 K, relative_change = 0.034018397148327106 Iter 3: T = 901.4678238649153 K, relative_change = 0.033129680510688334 Iter 20: T = 632.5010000262616 K, relative_change = 0.008657931441971278 Iter 40: T = 601.3622364436253 K, relative_change = 0.0003407575128744026 Iter 60: T = 600.3150843057265 K, relative_change = 1.0516928325988e-5 Iter 80: T = 600.2829485145389 K, relative_change = 3.2180532456863786e-7 Converged in 98 iterations to T = 600.2819783139583 K Iter 1: T = 980.116463160623 K, relative_change = 0.01988353683937708 Iter 2: T = 962.2777435310346 K, relative_change = 0.01820061217221375 Iter 3: T = 946.3631211700479 K, relative_change = 0.01653849158205485 Iter 20: T = 848.7126330449323 K, relative_change = 0.0014522753410618398 Iter 40: T = 842.3292203993134 K, relative_change = 4.6163896478266084e-5 Iter 60: T = 842.1312095278472 K, relative_change = 1.4138684159722342e-6 Iter 80: T = 842.1251496670578 K, relative_change = 4.325244440391165e-8 Converged in 87 iterations to T = 842.1250148700561 K Iter 1: T = 976.3663958244659 K, relative_change = 0.0236336041755341 Iter 2: T = 954.8958550406616 K, relative_change = 0.02199024963950561 Iter 3: T = 935.4973303714698 K, relative_change = 0.020314806653303187 Iter 20: T = 809.0083383619148 K, relative_change = 0.0021126716518719527 Iter 40: T = 800.0848628499982 K, relative_change = 6.836990758572455e-5 Iter 60: T = 799.8062376768 K, relative_change = 2.0951843918433434e-6 Iter 80: T = 799.7977089453135 K, relative_change = 6.4096098644661e-8 Converged in 89 iterations to T = 799.79749584789 K Iter 1: T = 980.847165434604 K, relative_change = 0.019152834565395974 Iter 2: T = 963.7058385793945 K, relative_change = 0.01747604260814108 Iter 3: T = 948.4502584043486 K, relative_change = 0.01583012114727275 Iter 20: T = 855.8383214323095 K, relative_change = 0.0013503539252558456 Iter 40: T = 849.8604715764335 K, relative_change = 4.2806711221779404e-5 Iter 60: T = 849.675226411464 K, relative_change = 1.3109330225028275e-6 Iter 80: T = 849.6695574067201 K, relative_change = 4.010338268040067e-8 Converged in 86 iterations to T = 849.6694413592336 K Iter 1: T = 967.3255768158936 K, relative_change = 0.03267442318410643 Iter 2: T = 936.726295884264 K, relative_change = 0.03163286660149315 Iter 3: T = 908.170792944027 K, relative_change = 0.030484361403862154 Iter 20: T = 677.7731531563445 K, relative_change = 0.006116318550904184 Iter 40: T = 654.997937971381 K, relative_change = 0.000222159057384075 Iter 60: T = 654.2554011907714 K, relative_change = 6.835358251276064e-6 Iter 80: T = 654.2326391895432 K, relative_change = 2.0913372704145574e-7 Converged in 96 iterations to T = 654.2319650069367 K Iter 1: T = 973.5954896673522 K, relative_change = 0.026404510332647714 Iter 2: T = 949.3838880709552 K, relative_change = 0.024868235168868102 Iter 3: T = 927.2971262144985 K, relative_change = 0.02326431081670732 Iter 20: T = 775.4547617476566 K, relative_change = 0.002820338819594777 Iter 40: T = 763.935965141352 K, relative_change = 9.307138676361909e-5 Iter 60: T = 763.573705490621 K, relative_change = 2.8539903041457368e-6 Iter 80: T = 763.5626141313087 K, relative_change = 8.731128810896569e-8 Converged in 91 iterations to T = 763.5623155476503 K Iter 1: T = 970.0755772722218 K, relative_change = 0.029924422727778132 Iter 2: T = 942.3100047026486 K, relative_change = 0.028622071537609397 Iter 3: T = 916.6601112990506 K, relative_change = 0.02722022824292519 Iter 20: T = 725.4412432982162 K, relative_change = 0.004224321658010246 Iter 40: T = 709.0131351228247 K, relative_change = 0.0001450665522495178 Iter 60: T = 708.4887661710085 K, relative_change = 4.454423599004185e-6 Iter 80: T = 708.4727036373262 K, relative_change = 1.362785370274584e-7 Converged in 93 iterations to T = 708.4722492972243 K Iter 1: T = 973.5584842423605 K, relative_change = 0.02644151575763955 Iter 2: T = 949.3099405020713 K, relative_change = 0.0249071259023127 Iter 3: T = 927.1865960533501 K, relative_change = 0.023304659000010716 Iter 20: T = 774.9766470186395 K, relative_change = 0.002831609261014197 Iter 40: T = 763.4173294638633 K, relative_change = 9.347265006791077e-5 Iter 60: T = 763.0537532487629 K, relative_change = 2.866324781248207e-6 Iter 80: T = 763.042621538043 K, relative_change = 8.768866123131812e-8 Converged in 91 iterations to T = 763.0423218680746 K Iter 1: T = 964.1897405037683 K, relative_change = 0.035810259496231615 Iter 2: T = 930.2968473154758 K, relative_change = 0.035151684118298304 Iter 3: T = 898.2899909009913 K, relative_change = 0.034404992886781935 Iter 20: T = 607.5365403587965 K, relative_change = 0.01049386344774024 Iter 40: T = 570.3304282804482 K, relative_change = 0.00043940218800799305 Iter 60: T = 569.0483082577508 K, relative_change = 1.3596538167978905e-5 Iter 80: T = 569.0089247687113 K, relative_change = 4.1607094346811834e-7 Converged in 100 iterations to T = 569.0077198570356 K Iter 1: T = 963.5784558451621 K, relative_change = 0.03642154415483787 Iter 2: T = 929.0356826262386 K, relative_change = 0.03584842833438598 Iter 3: T = 896.3382072310643 K, relative_change = 0.035195069475419526 Iter 20: T = 590.6835690023773 K, relative_change = 0.011958482275150948 Iter 40: T = 548.565992331675 K, relative_change = 0.0005275010932172247 Iter 60: T = 547.0839866014978 K, relative_change = 1.636041827792307e-5 Iter 80: T = 547.0384249599363 K, relative_change = 5.006850406485224e-7 Iter 100: T = 547.0370309956991 K, relative_change = 1.531637470896521e-8 Converged in 101 iterations to T = 547.0370239576544 K Iter 1: T = 969.2614200652861 K, relative_change = 0.03073857993471394 Iter 2: T = 940.6621741441237 K, relative_change = 0.02950622538884926 Iter 3: T = 914.1634842135094 K, relative_change = 0.02817025140266171 Iter 20: T = 712.2916756269317 K, relative_change = 0.004682721679133156 Iter 40: T = 694.3057236276834 K, relative_change = 0.00016296580211019887 Iter 60: T = 693.7287511069829 K, relative_change = 5.006373654577154e-6 Iter 80: T = 693.7110742369615 K, relative_change = 1.5316707977030368e-7 Converged in 94 iterations to T = 693.7105649792172 K Iter 1: T = 966.4933249228687 K, relative_change = 0.03350667507713127 Iter 2: T = 935.0264400895309 K, relative_change = 0.032557788059063016 Iter 3: T = 905.5696037959059 K, relative_change = 0.03150374687886307 Iter 20: T = 661.1955244080307 K, relative_change = 0.006946816257240968 Iter 40: T = 635.680064318198 K, relative_change = 0.0002588613786442809 Iter 60: T = 634.8400052789109 K, relative_change = 7.972246479487973e-6 Iter 80: T = 634.8142448725254 K, relative_change = 2.439250109242795e-7 Converged in 97 iterations to T = 634.813473880758 K Iter 1: T = 966.5220873398653 K, relative_change = 0.03347791266013471 Iter 2: T = 935.0852649163236 K, relative_change = 0.03252571548578317 Iter 3: T = 905.6597526381612 K, relative_change = 0.03146826645888339 Iter 20: T = 661.7891171172017 K, relative_change = 0.006915248454289256 Iter 40: T = 636.3775614091257 K, relative_change = 0.0002574317962585567 Iter 60: T = 635.541239252935 K, relative_change = 7.92792305598637e-6 Iter 80: T = 635.5155937842509 K, relative_change = 2.4256857807831673e-7 Converged in 97 iterations to T = 635.5148262328105 K Iter 1: T = 976.3781353304792 K, relative_change = 0.023621864669520803 Iter 2: T = 954.9191028043325 K, relative_change = 0.02197819855816749 Iter 3: T = 935.5317567733352 K, relative_change = 0.02030260571189952 Iter 20: T = 809.1419517595394 K, relative_change = 0.0021101546303871005 Iter 40: T = 800.227911366316 K, relative_change = 6.82837514048723e-5 Iter 60: T = 799.9495878338805 K, relative_change = 2.0925394618145147e-6 Iter 80: T = 799.9410683425227 K, relative_change = 6.401518028666913e-8 Converged in 89 iterations to T = 799.9408554759767 K Iter 1: T = 965.1884062089327 K, relative_change = 0.034811593791067304 Iter 2: T = 932.351705600764 K, relative_change = 0.03402102677252896 Iter 3: T = 901.4604454961155 K, relative_change = 0.03313262572383415 Iter 20: T = 632.4460964449721 K, relative_change = 0.008661584350748602 Iter 40: T = 601.2953197948284 K, relative_change = 0.00034094211158906183 Iter 60: T = 600.2477146596705 K, relative_change = 1.0522676508533114e-5 Iter 80: T = 600.2155649103696 K, relative_change = 3.219812600963819e-7 Converged in 98 iterations to T = 600.2145942883323 K Iter 1: T = 964.5741525900285 K, relative_change = 0.035425847409971525 Iter 2: T = 931.0886250954961 K, relative_change = 0.0347153481198087 Iter 3: T = 899.5130410239792 K, relative_change = 0.03391254411284248 Iter 20: T = 617.4743915061176 K, relative_change = 0.009719087436953447 Iter 40: T = 582.8376414469781 K, relative_change = 0.0003962903448311163 Iter 60: T = 581.6565670377352 K, relative_change = 1.22486618467899e-5 Iter 80: T = 581.6203022586359 K, relative_change = 3.7481113601644456e-7 Converged in 99 iterations to T = 581.6191994445263 K Iter 1: T = 964.2981705710226 K, relative_change = 0.03570182942897736 Iter 2: T = 930.5202847585182 K, relative_change = 0.03502846613563773 Iter 3: T = 898.6353119844159 K, relative_change = 0.0342657471270245 Iter 20: T = 610.3877035539883 K, relative_change = 0.010265233832674915 Iter 40: T = 573.9415033336871 K, relative_change = 0.00042644302420211474 Iter 60: T = 572.6895120757591 K, relative_change = 1.3191054262278713e-5 Iter 80: T = 572.651058799208 K, relative_change = 4.036583823295927e-7 Converged in 99 iterations to T = 572.6498894224387 K Iter 1: T = 980.1937122761686 K, relative_change = 0.019806287723831355 Iter 2: T = 962.4288773967744 K, relative_change = 0.018123800078396136 Iter 3: T = 946.5842268074242 K, relative_change = 0.0164631911629746 Iter 20: T = 849.4743893332686 K, relative_change = 0.0014411616103356803 Iter 40: T = 843.1349974953512 K, relative_change = 4.579692894355515e-5 Iter 60: T = 842.938373602891 K, relative_change = 1.4026158694934848e-6 Iter 80: T = 842.9323562094689 K, relative_change = 4.29081988924435e-8 Converged in 87 iterations to T = 842.9322223571356 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: 36%|███████████▊ | 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: 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011165403270466182 Iteration 10: norm(E*F-F'*E) = 2.2924394455790618e-5 Iteration 20: norm(E*F-F'*E) = 1.563534137504874e-7 Iteration 30: norm(E*F-F'*E) = 1.8311879498287595e-9 Converged after 31 iterations. norm(E*F-F'*E) = 7.739436362847707e-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: 44%|██████████████▌ | ETA: 0:00:01 Bin 1 progress: 88%|████████████████████████████▉ | 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010631671928723167 Iteration 10: norm(E*F-F'*E) = 5.028049455359402e-5 Iteration 20: norm(E*F-F'*E) = 5.616760671256716e-7 Iteration 30: norm(E*F-F'*E) = 6.9909631569948065e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.1378970037722763e-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: 44%|██████████████▌ | ETA: 0:00:01 Bin 1 progress: 88%|████████████████████████████▉ | 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010631671928723167 Iteration 10: norm(E*F-F'*E) = 5.028049455359402e-5 Iteration 20: norm(E*F-F'*E) = 5.616760671256716e-7 Iteration 30: norm(E*F-F'*E) = 6.9909631569948065e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.1378970037722763e-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: 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010631671928723167 Iteration 10: norm(E*F-F'*E) = 5.028049455359402e-5 Iteration 20: norm(E*F-F'*E) = 5.616760671256716e-7 Iteration 30: norm(E*F-F'*E) = 6.9909631569948065e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.1378970037722763e-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: 41%|█████████████▍ | ETA: 0:00:01 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010631671928723167 Iteration 10: norm(E*F-F'*E) = 5.028049455359402e-5 Iteration 20: norm(E*F-F'*E) = 5.616760671256716e-7 Iteration 30: norm(E*F-F'*E) = 6.9909631569948065e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.1378970037722763e-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: 44%|██████████████▌ | ETA: 0:00:01 Bin 1 progress: 88%|████████████████████████████▉ | 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010631671928723167 Iteration 10: norm(E*F-F'*E) = 5.028049455359402e-5 Iteration 20: norm(E*F-F'*E) = 5.616760671256716e-7 Iteration 30: norm(E*F-F'*E) = 6.9909631569948065e-9 Converged after 37 iterations. norm(E*F-F'*E) = 2.1378970037722763e-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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-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, 3.1424922343134543e-19, 5.551115123125783e-15, 2.5863755581667647e-12, 1.0103917702508625e-11, 4.177991286269389e-12, 6.323830348264892e-13, 1.0125233984581428e-13, 8.104628079763643e-15, 4.371503159461554e-16, 2.4069288229178198e-17, 1.1993986533120893e-18, 5.815410578492806e-20, 1.4756902127946014e-21, 8.075350730048933e-23, 7.106459654317833e-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: 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 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/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011165403270466182 Iteration 10: norm(E*F-F'*E) = 2.2924394455790618e-5 Iteration 20: norm(E*F-F'*E) = 1.563534137504874e-7 Iteration 30: norm(E*F-F'*E) = 1.8311879498287595e-9 Converged after 31 iterations. norm(E*F-F'*E) = 7.739436362847707e-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: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-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.8973538018496328e-19, 2.6645352591003757e-14, 3.311129148642067e-12, 9.627854069549358e-12, 6.281197784119286e-12, 8.242295734817162e-13, 1.0746958878371515e-13, 8.770761894538737e-15, 4.3368086899420177e-16, 2.970713952610282e-17, 1.1926223897340549e-18, 5.651297944962285e-20, 1.616310916928616e-21, 7.935764001680609e-23, 6.600592449666959e-15] === 3D Spectral Solution Complete (DIRECT) === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-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.1845193777663408 Iteration 3: convergence error = 0.4842939367584489 Iteration 4: convergence error = 0.2209144278951865 Iteration 5: convergence error = 0.0915991099786516 Iteration 6: convergence error = 0.03576093742294521 Iteration 7: convergence error = 0.01359208191288714 Iteration 8: convergence error = 0.005109475888597641 Iteration 9: convergence error = 0.0019123763701812358 Iteration 10: convergence error = 0.0007145628006044933 Converged after 10 iterations Energy conservation errors by band: [-4.579669976578771e-16, -5.2562121322097255e-16, -4.0245584642661925e-16, -4.597017211338539e-16, -4.926614671774132e-16, -5.134781488891349e-16, 5.928590951498336e-14, 5.5138116294983774e-12, 2.1167068098293385e-11, 7.858602657506708e-12, 1.9895196601282805e-13, 3.197442310920451e-14, 3.164135620181696e-15, 5.551115123125783e-17, 1.1796119636642288e-16, -1.5265566588595902e-16, -6.418476861114186e-17, -5.377642775528102e-17, -1.231653667943533e-16, 1.0269562977782698e-15] Writing spectral results to mesh... === 3D Spectral Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 54×54 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 6.233122709088313e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.800123372539454e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 54×54 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 6.233122709088313e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.800123372539454e-16 Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-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.877406409464993e-19, 8.060219158778636e-14, 1.4352963262354024e-12, 1.3169909607313457e-11, 4.661160346586257e-12, 4.121147867408581e-13, 8.08242361927114e-14, 6.938893903907228e-15, 3.7816971776294395e-16, 3.133344278483108e-17, 1.1807639284724947e-18, 5.590417451878382e-20, 1.5948042210022374e-21, 7.153819828876575e-23, 4.494952646004891e-15] === 3D Spectral Solution Complete (DIRECT) === ✓ Spectral Consistency tests complete ================================================================================ TEST SUITE COMPLETE ================================================================================ Test Summary: | Pass Total Time RayTraceHeatTransfer.jl | 1394 1394 9m39.6s Testing RayTraceHeatTransfer tests passed Testing completed after 601.34s PkgEval succeeded after 682.84s