Package evaluation to test RayTraceHeatTransfer on Julia 1.14.0-DEV.2590 (cfde0ec785*) started at 2026-07-03T23:05:22.504 ################################################################################ # Set-up # Installing PkgEval dependencies (TestEnv)... Activating project at `~/.julia/environments/v1.14` Set-up completed after 14.85s ################################################################################ # 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.3 [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.2.0 [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.5+2 [deac9b47] + LibCURL_jll v8.21.0+0 [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.7+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.76s ################################################################################ # Precompilation # Precompiling PkgEval dependencies... Precompiling package dependencies... Precompiling project... 5.4 s ✓ StatsBase 22.7 s ✓ GeometryBasics 8.4 s ✓ RayTraceHeatTransfer 3 dependencies successfully precompiled in 37 seconds. 58 already precompiled. Precompilation completed after 64.65s ################################################################################ # Testing # Testing RayTraceHeatTransfer Status `/tmp/jl_UJtbwO/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_UJtbwO/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.3 [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.2.0 [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.5+2 [deac9b47] LibCURL_jll v8.21.0+0 [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.7+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:05:01 Bin 1 progress: 62%|████████████████████▍ | ETA: 0:00:04 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:06 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.761739987160342e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.09090909090909094) ≥ min gas E (0.03305785123966937); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012262321932357965 Iteration 10: norm(E*F-F'*E) = 2.3916391837725985e-5 Iteration 20: norm(E*F-F'*E) = 3.0698096829047203e-7 Converged after 29 iterations. norm(E*F-F'*E) = 5.304595436661992e-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: 35%|███████████▍ | ETA: 0:00:02 Bin 1 progress: 76%|█████████████████████████▎ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.761739987160345e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.09090909090909094) ≥ min gas E (0.03305785123966936); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012334104976940832 Iteration 10: norm(E*F-F'*E) = 2.7093160663800822e-5 Iteration 20: norm(E*F-F'*E) = 3.026833894519392e-7 Converged after 29 iterations. norm(E*F-F'*E) = 4.132702383203631e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 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.012212001222722483 Iteration 10: norm(E*F-F'*E) = 2.612573748840539e-5 Iteration 20: norm(E*F-F'*E) = 2.5843184626258423e-7 Converged after 29 iterations. norm(E*F-F'*E) = 3.3836489014025146e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 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.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.012407286758703921 Iteration 10: norm(E*F-F'*E) = 3.612259072703983e-5 Iteration 20: norm(E*F-F'*E) = 4.38483344069849e-7 Converged after 29 iterations. norm(E*F-F'*E) = 6.4521715823129145e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▋ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.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.011981195146509815 Iteration 10: norm(E*F-F'*E) = 3.616920098367904e-5 Iteration 20: norm(E*F-F'*E) = 4.159236774512808e-7 Iteration 30: norm(E*F-F'*E) = 5.9055496973091166e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.097881806035392e-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.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.011856893728376195 Iteration 10: norm(E*F-F'*E) = 4.197490229569032e-5 Iteration 20: norm(E*F-F'*E) = 4.796050270186311e-7 Iteration 30: norm(E*F-F'*E) = 6.74575962010266e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.2495003426299833e-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: 77%|█████████████████████████▎ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.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.01176172412143844 Iteration 10: norm(E*F-F'*E) = 3.830269673465132e-5 Iteration 20: norm(E*F-F'*E) = 4.476953562940827e-7 Iteration 30: norm(E*F-F'*E) = 6.439648804840707e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.1990557581172448e-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.011805203463183233 Iteration 10: norm(E*F-F'*E) = 3.3736699390047284e-5 Iteration 20: norm(E*F-F'*E) = 3.782349654154685e-7 Iteration 30: norm(E*F-F'*E) = 5.383781842894097e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.0012287312573679e-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.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.011945024231774534 Iteration 10: norm(E*F-F'*E) = 4.066478027687779e-5 Iteration 20: norm(E*F-F'*E) = 4.761554526978528e-7 Iteration 30: norm(E*F-F'*E) = 6.773362469086272e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.2595048072275385e-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: 75%|████████████████████████▉ | 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.011746666412431411 Iteration 10: norm(E*F-F'*E) = 2.9021286622116312e-5 Iteration 20: norm(E*F-F'*E) = 2.2391289874197802e-7 Iteration 30: norm(E*F-F'*E) = 2.4065512921661473e-9 Converged after 31 iterations. norm(E*F-F'*E) = 9.983669682780785e-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.035796692584650944 Iteration 10: norm(E*F-F'*E) = 0.00012550034891200594 Iteration 20: norm(E*F-F'*E) = 9.857115918507984e-7 Iteration 30: norm(E*F-F'*E) = 1.1882337112550307e-8 Converged after 33 iterations. norm(E*F-F'*E) = 2.1513356892569624e-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.027446802128943428 Iteration 10: norm(E*F-F'*E) = 7.090783821907047e-5 Iteration 20: norm(E*F-F'*E) = 6.868762800311171e-7 Iteration 30: norm(E*F-F'*E) = 8.8795437431508e-9 Converged after 31 iterations. norm(E*F-F'*E) = 3.7839754200302155e-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.025152967470968728 Iteration 10: norm(E*F-F'*E) = 6.087808357841551e-5 Iteration 20: norm(E*F-F'*E) = 5.783448111944265e-7 Converged after 29 iterations. norm(E*F-F'*E) = 7.609690541097877e-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.021364592270443697 Iteration 10: norm(E*F-F'*E) = 5.2270565843891504e-5 Iteration 20: norm(E*F-F'*E) = 6.157613790849603e-7 Converged after 29 iterations. norm(E*F-F'*E) = 9.492690640076806e-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: 75%|████████████████████████▉ | 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.011981195146509815 Iteration 10: norm(E*F-F'*E) = 3.616920098367904e-5 Iteration 20: norm(E*F-F'*E) = 4.159236774512808e-7 Iteration 30: norm(E*F-F'*E) = 5.9055496973091166e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.097881806035392e-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: 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: 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.010314200523744462 Iteration 10: norm(E*F-F'*E) = 4.035201148869848e-5 Iteration 20: norm(E*F-F'*E) = 4.477302195577377e-7 Iteration 30: norm(E*F-F'*E) = 6.1603925786228965e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.774644590291477e-10 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▌ | ETA: 0:00:02 Bin 1 progress: 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.01128495585245405 Iteration 10: norm(E*F-F'*E) = 3.432462752680574e-5 Iteration 20: norm(E*F-F'*E) = 3.723991000767515e-7 Iteration 30: norm(E*F-F'*E) = 4.846340791199859e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.641001889515946e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 84%|███████████████████████████▉ | 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.010314200523744462 Iteration 10: norm(E*F-F'*E) = 4.035201148869848e-5 Iteration 20: norm(E*F-F'*E) = 4.477302195577377e-7 Iteration 30: norm(E*F-F'*E) = 6.1603925786228965e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.774644590291477e-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: 10:53:30 Bin 4 ray tracing: 9%|██▋ | ETA: 0:00:54 Bin 4 ray tracing: 16%|████▉ | ETA: 0:00:31 Bin 4 ray tracing: 24%|███████▎ | ETA: 0:00:23 Bin 4 ray tracing: 32%|█████████▌ | ETA: 0:00:17 Bin 4 ray tracing: 40%|████████████ | ETA: 0:00:14 Bin 4 ray tracing: 48%|██████████████▍ | ETA: 0:00:11 Bin 4 ray tracing: 57%|█████████████████ | ETA: 0:00:09 Bin 4 ray tracing: 66%|███████████████████▋ | ETA: 0:00:06 Bin 4 ray tracing: 74%|██████████████████████▏ | ETA: 0:00:05 Bin 4 ray tracing: 82%|████████████████████████▌ | ETA: 0:00:03 Bin 4 ray tracing: 90%|██████████████████████████▉ | ETA: 0:00:02 Bin 4 ray tracing: 98%|█████████████████████████████▎| 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: 8%|██▍ | ETA: 0:00:12 Bin 5 ray tracing: 16%|████▊ | ETA: 0:00:11 Bin 5 ray tracing: 24%|███████▎ | ETA: 0:00:10 Bin 5 ray tracing: 32%|█████████▋ | ETA: 0:00:09 Bin 5 ray tracing: 40%|████████████ | ETA: 0:00:08 Bin 5 ray tracing: 48%|██████████████▍ | ETA: 0:00:07 Bin 5 ray tracing: 56%|████████████████▊ | ETA: 0:00:06 Bin 5 ray tracing: 64%|███████████████████▏ | ETA: 0:00:05 Bin 5 ray tracing: 72%|█████████████████████▋ | ETA: 0:00:04 Bin 5 ray tracing: 81%|████████████████████████▏ | ETA: 0:00:02 Bin 5 ray tracing: 89%|██████████████████████████▊ | ETA: 0:00:01 Bin 5 ray tracing: 97%|█████████████████████████████▏| ETA: 0:00:00 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: 8%|██▍ | ETA: 0:00:12 Bin 6 ray tracing: 16%|████▊ | ETA: 0:00:11 Bin 6 ray tracing: 24%|███████▏ | ETA: 0:00:10 Bin 6 ray tracing: 32%|█████████▌ | ETA: 0:00:09 Bin 6 ray tracing: 40%|████████████ | ETA: 0:00:07 Bin 6 ray tracing: 49%|██████████████▋ | ETA: 0:00:06 Bin 6 ray tracing: 57%|█████████████████▎ | ETA: 0:00:05 Bin 6 ray tracing: 65%|███████████████████▋ | ETA: 0:00:04 Bin 6 ray tracing: 73%|██████████████████████ | ETA: 0:00:03 Bin 6 ray tracing: 82%|████████████████████████▌ | ETA: 0:00:02 Bin 6 ray tracing: 90%|███████████████████████████▏ | ETA: 0:00:01 Bin 6 ray tracing: 99%|█████████████████████████████▊| ETA: 0:00:00 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: 9%|██▊ | ETA: 0:00:11 Bin 7 ray tracing: 18%|█████▍ | 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: 44%|█████████████▏ | ETA: 0:00:07 Bin 7 ray tracing: 53%|███████████████▊ | ETA: 0:00:06 Bin 7 ray tracing: 61%|██████████████████▍ | ETA: 0:00:05 Bin 7 ray tracing: 69%|████████████████████▊ | ETA: 0:00:04 Bin 7 ray tracing: 77%|███████████████████████▏ | ETA: 0:00:03 Bin 7 ray tracing: 85%|█████████████████████████▌ | ETA: 0:00:02 Bin 7 ray tracing: 93%|███████████████████████████▉ | ETA: 0:00:01 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:11 Bin 8 ray tracing: 17%|█████▏ | ETA: 0:00:10 Bin 8 ray tracing: 25%|███████▌ | ETA: 0:00:09 Bin 8 ray tracing: 33%|██████████ | ETA: 0:00:08 Bin 8 ray tracing: 42%|████████████▋ | ETA: 0:00:07 Bin 8 ray tracing: 51%|███████████████▎ | ETA: 0:00:06 Bin 8 ray tracing: 59%|█████████████████▊ | ETA: 0:00:05 Bin 8 ray tracing: 67%|████████████████████▏ | ETA: 0:00:04 Bin 8 ray tracing: 75%|██████████████████████▋ | ETA: 0:00:03 Bin 8 ray tracing: 83%|█████████████████████████ | ETA: 0:00:02 Bin 8 ray tracing: 91%|███████████████████████████▍ | ETA: 0:00:01 Bin 8 ray tracing: 99%|█████████████████████████████▉| ETA: 0:00:00 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:11 Bin 9 ray tracing: 17%|█████ | ETA: 0:00:11 Bin 9 ray tracing: 25%|███████▍ | 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:01 Bin 9 ray tracing: 96%|████████████████████████████▉ | ETA: 0:00:00 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:11 Bin 10 ray tracing: 17%|████▉ | ETA: 0:00:10 Bin 10 ray tracing: 25%|███████▍ | ETA: 0:00:09 Bin 10 ray tracing: 34%|█████████▉ | ETA: 0:00:08 Bin 10 ray tracing: 43%|████████████▍ | ETA: 0:00:07 Bin 10 ray tracing: 51%|██████████████▉ | ETA: 0:00:06 Bin 10 ray tracing: 60%|█████████████████▍ | ETA: 0:00:05 Bin 10 ray tracing: 69%|███████████████████▉ | ETA: 0:00:04 Bin 10 ray tracing: 77%|██████████████████████▍ | ETA: 0:00:03 Bin 10 ray tracing: 85%|████████████████████████▊ | ETA: 0:00:02 Bin 10 ray tracing: 94%|███████████████████████████▏ | ETA: 0:00:01 Bin 10 ray tracing: 100%|█████████████████████████████| Time: 0:00:11 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: 20%|██████▋ | ETA: 0:00:04 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:03 Bin 1 progress: 62%|████████████████████▌ | ETA: 0:00:02 Bin 1 progress: 84%|███████████████████████████▉ | ETA: 0:00:01 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: 67%|██████████████████████ | ETA: 0:00:02 Bin 2 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 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: 20%|██████▋ | ETA: 0:00:04 Bin 3 progress: 42%|█████████████▉ | ETA: 0:00:03 Bin 3 progress: 64%|█████████████████████▎ | ETA: 0:00:02 Bin 3 progress: 87%|████████████████████████████▋ | 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: 42%|█████████████▉ | ETA: 0:00:03 Bin 5 progress: 64%|█████████████████████▎ | ETA: 0:00:02 Bin 5 progress: 87%|████████████████████████████▋ | 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: 47%|███████████████▍ | ETA: 0:00:02 Bin 7 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 7 progress: 93%|██████████████████████████████▊ | 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: 67%|██████████████████████ | ETA: 0:00:02 Bin 8 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 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: 22%|███████▍ | ETA: 0:00:04 Bin 9 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 9 progress: 67%|██████████████████████ | ETA: 0:00:02 Bin 9 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 Bin 9 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 10/10 Using 1 threads for spectral bin 10 Bin 10 progress: 22%|███████▏ | ETA: 0:00:04 Bin 10 progress: 44%|██████████████▎ | ETA: 0:00:03 Bin 10 progress: 69%|██████████████████████ | ETA: 0:00:01 Bin 10 progress: 91%|█████████████████████████████▏ | ETA: 0:00:00 Bin 10 progress: 100%|████████████████████████████████| Time: 0:00:04 Smoothing F matrix for spectral bin 1/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010314200523744462 Iteration 10: norm(E*F-F'*E) = 4.035201148869848e-5 Iteration 20: norm(E*F-F'*E) = 4.477302195577377e-7 Iteration 30: norm(E*F-F'*E) = 6.1603925786228965e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.774644590291477e-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.011282465710458186 Iteration 10: norm(E*F-F'*E) = 3.347492431637114e-5 Iteration 20: norm(E*F-F'*E) = 3.6198538193098657e-7 Iteration 30: norm(E*F-F'*E) = 4.708296943205918e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.5362759859931507e-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.010113632434334953 Iteration 10: norm(E*F-F'*E) = 3.830603516608284e-5 Iteration 20: norm(E*F-F'*E) = 3.242897419797805e-7 Iteration 30: norm(E*F-F'*E) = 3.636277171060052e-9 Converged after 33 iterations. norm(E*F-F'*E) = 6.241590663247264e-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.010123067464447844 Iteration 10: norm(E*F-F'*E) = 3.090369712048221e-5 Iteration 20: norm(E*F-F'*E) = 2.867290124702505e-7 Iteration 30: norm(E*F-F'*E) = 3.465065186919473e-9 Converged after 33 iterations. norm(E*F-F'*E) = 6.043589700325013e-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.010765185059851723 Iteration 10: norm(E*F-F'*E) = 5.499422015672268e-5 Iteration 20: norm(E*F-F'*E) = 6.117747306641847e-7 Iteration 30: norm(E*F-F'*E) = 7.630436646726819e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.592529085036651e-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.01037057188770565 Iteration 10: norm(E*F-F'*E) = 3.404625861475412e-5 Iteration 20: norm(E*F-F'*E) = 3.402747737987872e-7 Iteration 30: norm(E*F-F'*E) = 4.301246406185086e-9 Converged after 33 iterations. norm(E*F-F'*E) = 7.676987009527532e-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.01093140000795055 Iteration 10: norm(E*F-F'*E) = 2.9115461745118325e-5 Iteration 20: norm(E*F-F'*E) = 3.264579342251619e-7 Iteration 30: norm(E*F-F'*E) = 4.336940152532764e-9 Converged after 33 iterations. norm(E*F-F'*E) = 7.773087455622984e-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.011161942050631337 Iteration 10: norm(E*F-F'*E) = 3.688192985390553e-5 Iteration 20: norm(E*F-F'*E) = 3.576037242058343e-7 Iteration 30: norm(E*F-F'*E) = 4.5111970795258635e-9 Converged after 33 iterations. norm(E*F-F'*E) = 8.044748716259133e-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.010706222380516243 Iteration 10: norm(E*F-F'*E) = 4.614272142804852e-5 Iteration 20: norm(E*F-F'*E) = 5.399210222113785e-7 Iteration 30: norm(E*F-F'*E) = 7.091724619929032e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.330245038344739e-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.010906509441903623 Iteration 10: norm(E*F-F'*E) = 4.51583878593809e-5 Iteration 20: norm(E*F-F'*E) = 5.161060619974198e-7 Iteration 30: norm(E*F-F'*E) = 6.9596733183521815e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.315304212555203e-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: 12%|███▋ | ETA: 0:00:15 Bin 4 ray tracing: 18%|█████▌ | ETA: 0:00:14 Bin 4 ray tracing: 24%|███████▎ | ETA: 0:00:13 Bin 4 ray tracing: 30%|█████████▏ | ETA: 0:00:12 Bin 4 ray tracing: 36%|██████████▉ | ETA: 0:00:11 Bin 4 ray tracing: 42%|████████████▌ | ETA: 0:00:10 Bin 4 ray tracing: 48%|██████████████▎ | ETA: 0:00:09 Bin 4 ray tracing: 53%|████████████████ | ETA: 0:00:08 Bin 4 ray tracing: 60%|█████████████████▉ | ETA: 0:00:07 Bin 4 ray tracing: 66%|███████████████████▊ | ETA: 0:00:06 Bin 4 ray tracing: 72%|█████████████████████▋ | ETA: 0:00:05 Bin 4 ray tracing: 78%|███████████████████████▍ | ETA: 0:00:04 Bin 4 ray tracing: 84%|█████████████████████████▎ | ETA: 0:00:03 Bin 4 ray tracing: 90%|███████████████████████████▏ | ETA: 0:00:02 Bin 4 ray tracing: 97%|█████████████████████████████ | ETA: 0:00:01 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 Updating spectral results for spectral bin 4 Energy per ray: 0.0001853335835185918 Processing spectral bin 5/10 Bin 5 ray tracing: 6%|█▉ | ETA: 0:00:15 Bin 5 ray tracing: 13%|███▊ | ETA: 0:00:14 Bin 5 ray tracing: 19%|█████▋ | ETA: 0:00:13 Bin 5 ray tracing: 25%|███████▌ | 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: 56%|████████████████▉ | ETA: 0:00:07 Bin 5 ray tracing: 62%|██████████████████▋ | ETA: 0:00:06 Bin 5 ray tracing: 68%|████████████████████▌ | ETA: 0:00:05 Bin 5 ray tracing: 74%|██████████████████████▎ | ETA: 0:00:04 Bin 5 ray tracing: 80%|████████████████████████▏ | ETA: 0:00:03 Bin 5 ray tracing: 86%|█████████████████████████▉ | ETA: 0:00:02 Bin 5 ray tracing: 92%|███████████████████████████▋ | ETA: 0:00:01 Bin 5 ray tracing: 98%|█████████████████████████████▌| ETA: 0:00:00 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: 13%|███▊ | ETA: 0:00:14 Bin 6 ray tracing: 19%|█████▋ | ETA: 0:00:13 Bin 6 ray tracing: 25%|███████▍ | ETA: 0:00:13 Bin 6 ray tracing: 31%|█████████▏ | ETA: 0:00:12 Bin 6 ray tracing: 36%|██████████▉ | ETA: 0:00:11 Bin 6 ray tracing: 42%|████████████▋ | ETA: 0:00:10 Bin 6 ray tracing: 48%|██████████████▍ | ETA: 0:00:09 Bin 6 ray tracing: 54%|████████████████▎ | ETA: 0:00:08 Bin 6 ray tracing: 60%|██████████████████▏ | ETA: 0:00:07 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:04 Bin 6 ray tracing: 85%|█████████████████████████▌ | ETA: 0:00:03 Bin 6 ray tracing: 91%|███████████████████████████▍ | ETA: 0:00:01 Bin 6 ray tracing: 98%|█████████████████████████████▎| ETA: 0:00:00 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: 6%|█▉ | ETA: 0:00:15 Bin 7 ray tracing: 13%|███▊ | ETA: 0:00:14 Bin 7 ray tracing: 19%|█████▊ | ETA: 0:00:13 Bin 7 ray tracing: 25%|███████▋ | ETA: 0:00:12 Bin 7 ray tracing: 32%|█████████▌ | ETA: 0:00:11 Bin 7 ray tracing: 38%|███████████▍ | ETA: 0:00:10 Bin 7 ray tracing: 44%|█████████████▏ | ETA: 0:00:09 Bin 7 ray tracing: 50%|███████████████ | ETA: 0:00:08 Bin 7 ray tracing: 56%|████████████████▊ | ETA: 0:00:07 Bin 7 ray tracing: 62%|██████████████████▌ | ETA: 0:00:06 Bin 7 ray tracing: 68%|████████████████████▍ | ETA: 0:00:05 Bin 7 ray tracing: 74%|██████████████████████▏ | ETA: 0:00:04 Bin 7 ray tracing: 80%|███████████████████████▉ | ETA: 0:00:03 Bin 7 ray tracing: 86%|█████████████████████████▊ | ETA: 0:00:02 Bin 7 ray tracing: 92%|███████████████████████████▌ | ETA: 0:00:01 Bin 7 ray tracing: 98%|█████████████████████████████▍| ETA: 0:00:00 Bin 7 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 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: 13%|███▉ | ETA: 0:00:14 Bin 8 ray tracing: 19%|█████▊ | ETA: 0:00:13 Bin 8 ray tracing: 26%|███████▊ | ETA: 0:00:12 Bin 8 ray tracing: 32%|█████████▋ | ETA: 0:00:11 Bin 8 ray tracing: 39%|███████████▋ | ETA: 0:00:10 Bin 8 ray tracing: 45%|█████████████▋ | ETA: 0:00:09 Bin 8 ray tracing: 51%|███████████████▍ | ETA: 0:00:08 Bin 8 ray tracing: 57%|█████████████████▏ | ETA: 0:00:07 Bin 8 ray tracing: 63%|██████████████████▉ | ETA: 0:00:06 Bin 8 ray tracing: 69%|████████████████████▋ | ETA: 0:00:05 Bin 8 ray tracing: 74%|██████████████████████▍ | ETA: 0:00:04 Bin 8 ray tracing: 80%|████████████████████████▏ | ETA: 0:00:03 Bin 8 ray tracing: 86%|█████████████████████████▉ | ETA: 0:00:02 Bin 8 ray tracing: 92%|███████████████████████████▊ | 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:12 Bin 9 ray tracing: 31%|█████████▎ | ETA: 0:00:11 Bin 9 ray tracing: 37%|███████████▏ | ETA: 0:00:10 Bin 9 ray tracing: 44%|█████████████▏ | ETA: 0:00:09 Bin 9 ray tracing: 50%|███████████████ | ETA: 0:00:08 Bin 9 ray tracing: 56%|████████████████▉ | ETA: 0:00:07 Bin 9 ray tracing: 63%|██████████████████▊ | ETA: 0:00:06 Bin 9 ray tracing: 69%|████████████████████▋ | ETA: 0:00:05 Bin 9 ray tracing: 75%|██████████████████████▌ | ETA: 0:00:04 Bin 9 ray tracing: 81%|████████████████████████▏ | ETA: 0:00:03 Bin 9 ray tracing: 86%|█████████████████████████▉ | ETA: 0:00:02 Bin 9 ray tracing: 92%|███████████████████████████▋ | ETA: 0:00:01 Bin 9 ray tracing: 98%|█████████████████████████████▌| ETA: 0:00:00 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:16 Bin 10 ray tracing: 12%|███▌ | ETA: 0:00:15 Bin 10 ray tracing: 18%|█████▎ | ETA: 0:00:14 Bin 10 ray tracing: 24%|██████▉ | ETA: 0:00:13 Bin 10 ray tracing: 31%|████████▉ | ETA: 0:00:11 Bin 10 ray tracing: 37%|██████████▊ | ETA: 0:00:10 Bin 10 ray tracing: 43%|████████████▌ | ETA: 0:00:09 Bin 10 ray tracing: 49%|██████████████▎ | ETA: 0:00:09 Bin 10 ray tracing: 55%|███████████████▉ | ETA: 0:00:08 Bin 10 ray tracing: 61%|█████████████████▋ | ETA: 0:00:07 Bin 10 ray tracing: 66%|███████████████████▎ | ETA: 0:00:06 Bin 10 ray tracing: 72%|█████████████████████ | ETA: 0:00:05 Bin 10 ray tracing: 78%|██████████████████████▊ | ETA: 0:00:04 Bin 10 ray tracing: 84%|████████████████████████▌ | ETA: 0:00:03 Bin 10 ray tracing: 91%|██████████████████████████▎ | ETA: 0:00:02 Bin 10 ray tracing: 97%|████████████████████████████ | ETA: 0:00:01 Bin 10 ray tracing: 100%|█████████████████████████████| Time: 0:00:16 Updating spectral results for spectral bin 10 Energy per ray: 1.5017824710260812e-5 Iter 1: T = 967.2626173438663 K, relative_change = 0.03273738265613369 Iter 2: T = 936.5978662919207 K, relative_change = 0.03170261157838595 Iter 3: T = 907.974541254363 K, relative_change = 0.030560954778682318 Iter 20: T = 676.5602745782883 K, relative_change = 0.006173635188761231 Iter 40: T = 653.5953884211347 K, relative_change = 0.0002246327536905711 Iter 60: T = 652.8461692653391 K, relative_change = 6.911914835889483e-6 Iter 80: T = 652.8232018841973 K, relative_change = 2.1147646266859555e-7 Converged in 96 iterations to T = 652.8225216180018 K Iter 1: T = 970.3104317874231 K, relative_change = 0.029689568212576945 Iter 2: T = 942.7845253174851 K, relative_change = 0.028368144429027936 Iter 3: T = 917.3777239991916 K, relative_change = 0.026948683008705107 Iter 20: T = 729.1066021744392 K, relative_change = 0.004104008789190963 Iter 40: T = 713.0903047245724 K, relative_change = 0.00014044660878469823 Iter 60: T = 712.5797441037951 K, relative_change = 4.312044306623697e-6 Iter 80: T = 712.564105228309 K, relative_change = 1.319221004076264e-7 Converged in 93 iterations to T = 712.5636628722667 K Iter 1: T = 974.5140895921744 K, relative_change = 0.025485910407825618 Iter 2: T = 951.2166735915711 K, relative_change = 0.023906700015341068 Iter 3: T = 930.0321975475116 K, relative_change = 0.02227092589122924 Iter 20: T = 787.0527308783047 K, relative_change = 0.0025577306196192563 Iter 40: T = 776.4847413725167 K, relative_change = 8.379232550525793e-5 Iter 60: T = 776.1532775880188 K, relative_change = 2.568832032581758e-6 Iter 80: T = 776.1431299947772 K, relative_change = 7.858694292839729e-8 Converged in 90 iterations to T = 776.1428657795669 K Iter 1: T = 970.314706494954 K, relative_change = 0.029685293505046 Iter 2: T = 942.7931589235975 K, relative_change = 0.02836352720116129 Iter 3: T = 917.3907750150996 K, relative_change = 0.026943750777212037 Iter 20: T = 729.1728143925685 K, relative_change = 0.004101863982076467 Iter 40: T = 713.1638704018527 K, relative_change = 0.00014036453873343742 Iter 60: T = 712.6535559847384 K, relative_change = 4.30951534817257e-6 Iter 80: T = 712.6379246627332 K, relative_change = 1.3184472115035634e-7 Converged in 93 iterations to T = 712.6374825203566 K Iter 1: T = 969.2997141144024 K, relative_change = 0.030700285885597613 Iter 2: T = 940.73977905622 K, relative_change = 0.02946450374668286 Iter 3: T = 914.2812253370042 K, relative_change = 0.028125263019875456 Iter 20: T = 712.9263946900978 K, relative_change = 0.004659589824528275 Iter 40: T = 695.0186675154208 K, relative_change = 0.000162051120449021 Iter 60: T = 694.4443504974213 K, relative_change = 4.978155598614803e-6 Iter 80: T = 694.4267551362896 K, relative_change = 1.5230365322831394e-7 Converged in 94 iterations to T = 694.4262482268907 K Iter 1: T = 963.5341961006557 K, relative_change = 0.036465803899344305 Iter 2: T = 928.9442685765783 K, relative_change = 0.035899013926085785 Iter 3: T = 896.1965581214129 K, relative_change = 0.03525261047720847 Iter 20: T = 589.4077577109745 K, relative_change = 0.012077762739883816 Iter 40: T = 546.8863248772645 K, relative_change = 0.0005350877693330973 Iter 60: T = 545.3874714232672 K, relative_change = 1.6599029692322428e-5 Iter 80: T = 545.341388497251 K, relative_change = 5.079905180942958e-7 Iter 100: T = 545.3399785810977 K, relative_change = 1.5539858323850487e-8 Converged in 101 iterations to T = 545.3399714625123 K Iter 1: T = 966.90494421338 K, relative_change = 0.03309505578662003 Iter 2: T = 935.867748917612 K, relative_change = 0.032099531067159966 Iter 3: T = 906.8580003123492 K, relative_change = 0.03099770094526102 Iter 20: T = 669.5454354574506 K, relative_change = 0.00651554782958734 Iter 40: T = 645.4510005260042 K, relative_change = 0.00023956951364186436 Iter 60: T = 644.6617789797351 K, relative_change = 7.374391591997487e-6 Iter 80: T = 644.6375819119324 K, relative_change = 2.2562907852496315e-7 Converged in 96 iterations to T = 644.6368652207963 K Iter 1: T = 965.2268500524079 K, relative_change = 0.03477314994759201 Iter 2: T = 932.4306712056308 K, relative_change = 0.033977690161640685 Iter 3: T = 901.5820449779151 K, relative_change = 0.03308409641633567 Iter 20: T = 633.3492778055372 K, relative_change = 0.008601692143456194 Iter 40: T = 602.3954473851472 K, relative_change = 0.0003379210005742153 Iter 60: T = 601.3552629467367 K, relative_change = 1.04286097169622e-5 Iter 80: T = 601.3233418428292 K, relative_change = 3.1910214985596465e-7 Converged in 98 iterations to T = 601.3223781245952 K Iter 1: T = 980.1469649807428 K, relative_change = 0.01985303501925713 Iter 2: T = 962.3374231948101 K, relative_change = 0.018170276930136404 Iter 3: T = 946.450437659655 K, relative_change = 0.01650874750605958 Iter 20: T = 849.0136553200753 K, relative_change = 0.0014478772016091864 Iter 40: T = 842.6476576374607 K, relative_change = 4.6018646825910974e-5 Iter 60: T = 842.4501955050696 K, relative_change = 1.4094145127481e-6 Iter 80: T = 842.4441524459459 K, relative_change = 4.3116187564761604e-8 Converged in 87 iterations to T = 842.44401802269 K Iter 1: T = 976.4158593200098 K, relative_change = 0.023584140679990293 Iter 2: T = 954.9938017331245 K, relative_change = 0.021939481402733337 Iter 3: T = 935.6423655544755 K, relative_change = 0.020263415472990435 Iter 20: T = 809.5708652745603 K, relative_change = 0.0021020894958886233 Iter 40: T = 800.687068331195 K, relative_change = 6.80077662581225e-5 Iter 60: T = 800.4097115679524 K, relative_change = 2.084067004940499e-6 Iter 80: T = 800.4012216914175 K, relative_change = 6.375597597822363e-8 Converged in 89 iterations to T = 800.4010095648391 K Iter 1: T = 980.8221892464974 K, relative_change = 0.019177810753502597 Iter 2: T = 963.6570795372891 K, relative_change = 0.01750073550272675 Iter 3: T = 948.3790764430715 K, relative_change = 0.01585419068529382 Iter 20: T = 855.5976549314953 K, relative_change = 0.0013537228241409602 Iter 40: T = 849.606331624509 K, relative_change = 4.2917386180962435e-5 Iter 60: T = 849.4206628102328 K, relative_change = 1.3143261653966151e-6 Iter 80: T = 849.4149808347732 K, relative_change = 4.0207187488526696e-8 Converged in 86 iterations to T = 849.4148645217659 K Iter 1: T = 967.2814259577909 K, relative_change = 0.03271857404220906 Iter 2: T = 936.6362363572233 K, relative_change = 0.03168177200365767 Iter 3: T = 908.033178801964 K, relative_change = 0.030538064239861766 Iter 20: T = 676.9232820783668 K, relative_change = 0.0061564261513133425 Iter 40: T = 654.0153332096502 K, relative_change = 0.00022388913627075172 Iter 60: T = 653.2681210747728 K, relative_change = 6.88890013522603e-6 Iter 80: T = 653.2452153797032 K, relative_change = 2.1077218055093068e-7 Converged in 96 iterations to T = 653.2445369407211 K Iter 1: T = 973.49617422996 K, relative_change = 0.026503825770039912 Iter 2: T = 949.1854069427104 K, relative_change = 0.024972637726573103 Iter 3: T = 927.0004231874186 K, relative_change = 0.02337265574567649 Iter 20: T = 774.1696131259263 K, relative_change = 0.002850715763816267 Iter 40: T = 762.5416523541847 K, relative_change = 9.415347626160514e-5 Iter 60: T = 762.1758450774786 K, relative_change = 2.88725336205801e-6 Iter 80: T = 762.1646449861232 K, relative_change = 8.83289713437024e-8 Converged in 91 iterations to T = 762.1643434752629 K Iter 1: T = 969.9614194116155 K, relative_change = 0.03003858058838455 Iter 2: T = 942.0792180570697 K, relative_change = 0.028745680803942932 Iter 3: T = 916.3108801234479 K, relative_change = 0.027352623261094847 Iter 20: T = 723.6396451962448 K, relative_change = 0.004284613744492852 Iter 40: T = 707.005656244955 K, relative_change = 0.0001473937871830108 Iter 60: T = 706.4743688760182 K, relative_change = 4.526158248778429e-6 Iter 80: T = 706.4580940601893 K, relative_change = 1.3847344321002284e-7 Converged in 93 iterations to T = 706.4576337152283 K Iter 1: T = 973.5263224431268 K, relative_change = 0.026473677556873187 Iter 2: T = 949.2456646813134 K, relative_change = 0.024940936061060494 Iter 3: T = 927.0905110913704 K, relative_change = 0.023339746931982066 Iter 20: T = 774.5604013492384 K, relative_change = 0.002841450886501656 Iter 40: T = 762.9657181792329 K, relative_change = 9.382324857768369e-5 Iter 60: T = 762.6009926229814 K, relative_change = 2.877102078396407e-6 Iter 80: T = 762.5898256858011 K, relative_change = 8.801839262218687e-8 Converged in 91 iterations to T = 762.5895250674926 K Iter 1: T = 964.3047228034371 K, relative_change = 0.035695277196562826 Iter 2: T = 930.5337840874322 K, relative_change = 0.035021023870779905 Iter 3: T = 898.6561705736384 K, relative_change = 0.03425734138718669 Iter 20: T = 610.5587540298401 K, relative_change = 0.010251683391039441 Iter 40: T = 574.1575442080426 K, relative_change = 0.0004256813040493007 Iter 60: T = 572.9073302332882 K, relative_change = 1.3167229101577754e-5 Iter 80: T = 572.8689318213875 K, relative_change = 4.0292906123155314e-7 Converged in 99 iterations to T = 572.8677641133308 K Iter 1: T = 963.5540307369008 K, relative_change = 0.036445969263099204 Iter 2: T = 928.9852366956765 K, relative_change = 0.035876342102774475 Iter 3: T = 896.2600425395907 K, relative_change = 0.03522681832112489 Iter 20: T = 589.9804936231928 K, relative_change = 0.012024061170834146 Iter 40: T = 547.6409623787464 K, relative_change = 0.0005316640643669937 Iter 60: T = 546.1497053114671 K, relative_change = 1.6491337673701353e-5 Iter 80: T = 546.1038574355229 K, relative_change = 5.046933401480099e-7 Iter 100: T = 546.1024547121878 K, relative_change = 1.5438993506312744e-8 Converged in 101 iterations to T = 546.1024476299186 K Iter 1: T = 969.3785523775979 K, relative_change = 0.030621447622402143 Iter 2: T = 940.8995181848136 K, relative_change = 0.029378651016090357 Iter 3: T = 914.5235288500037 K, relative_change = 0.028032737635677134 Iter 20: T = 714.2279014078534 K, relative_change = 0.004612486177231605 Iter 40: T = 696.4795835012567 K, relative_change = 0.00016019234645382543 Iter 60: T = 695.9106732976134 K, relative_change = 4.920816307281796e-6 Iter 80: T = 695.8932438901159 K, relative_change = 1.505491682188998e-7 Converged in 94 iterations to T = 695.8927417619882 K Iter 1: T = 966.489537960977 K, relative_change = 0.03351046203902296 Iter 2: T = 935.0186945860703 K, relative_change = 0.03256201142259789 Iter 3: T = 905.5577331277352 K, relative_change = 0.03150841970210804 Iter 20: T = 661.1172531479817 K, relative_change = 0.006950989432015177 Iter 40: T = 635.5880583103229 K, relative_change = 0.00025905057455788373 Iter 60: T = 634.747505071246 K, relative_change = 7.978112643764596e-6 Iter 80: T = 634.7217294641758 K, relative_change = 2.4410453382458816e-7 Converged in 97 iterations to T = 634.7209580174195 K Iter 1: T = 966.4554954536471 K, relative_change = 0.03354450454635283 Iter 2: T = 934.9490628055798 K, relative_change = 0.03259998292345416 Iter 3: T = 905.4510088415402 K, relative_change = 0.031550439630927456 Iter 20: T = 660.4124162442524 K, relative_change = 0.006988681194027 Iter 40: T = 634.759181826833 K, relative_change = 0.00026076159355554573 Iter 60: T = 633.9141629454627 K, relative_change = 8.031166731698602e-6 Iter 80: T = 633.8882499805544 K, relative_change = 2.457281564372904e-7 Converged in 97 iterations to T = 633.8874744223845 K Iter 1: T = 976.4100304007469 K, relative_change = 0.023589969599253063 Iter 2: T = 954.9822602269791 K, relative_change = 0.02194546297826669 Iter 3: T = 935.6252766189094 K, relative_change = 0.020269469302465253 Iter 20: T = 809.5046357872181 K, relative_change = 0.0021033333818056715 Iter 40: T = 800.6161732913149 K, relative_change = 6.805032353198313e-5 Iter 60: T = 800.338667399256 K, relative_change = 2.0853734610412604e-6 Iter 80: T = 800.330172954503 K, relative_change = 6.379594540992622e-8 Converged in 89 iterations to T = 800.3299607137815 K Iter 1: T = 965.2233250743828 K, relative_change = 0.03477667492561719 Iter 2: T = 932.4234311427147 K, relative_change = 0.03398166318570936 Iter 3: T = 901.5708967007823 K, relative_change = 0.033088544765677526 Iter 20: T = 633.2666205644388 K, relative_change = 0.008607155748173475 Iter 40: T = 602.294825406989 K, relative_change = 0.0003381961113169906 Iter 60: T = 601.253964598463 K, relative_change = 1.043717508430517e-5 Iter 80: T = 601.2220226549618 K, relative_change = 3.193643101828147e-7 Converged in 98 iterations to T = 601.2210583074901 K Iter 1: T = 964.5859571961693 K, relative_change = 0.035414042803830786 Iter 2: T = 931.1129231067077 K, relative_change = 0.034701971182288426 Iter 3: T = 899.5505458417018 K, relative_change = 0.03389747524897019 Iter 20: T = 617.7722599635895 K, relative_change = 0.009696801644519174 Iter 40: T = 583.2091892187314 K, relative_change = 0.00039508339382235977 Iter 60: T = 582.0309782857526 K, relative_change = 1.2210970767034858e-5 Iter 80: T = 581.9948018435938 K, relative_change = 3.736574169063233e-7 Converged in 99 iterations to T = 581.9937017162039 K Iter 1: T = 964.3170774014563 K, relative_change = 0.03568292259854367 Iter 2: T = 930.5592370185868 K, relative_change = 0.03500699217506005 Iter 3: T = 898.6954979562732 K, relative_change = 0.03424149457094382 Iter 20: T = 610.8808993946352 K, relative_change = 0.010226213807919259 Iter 40: T = 574.5642399819378 K, relative_change = 0.00042425146960290366 Iter 60: T = 573.3173641030423 K, relative_change = 1.3122509163319998e-5 Iter 80: T = 573.279068735548 K, relative_change = 4.015601245319704e-7 Converged in 99 iterations to T = 573.2779041615927 K Iter 1: T = 980.1931982407768 K, relative_change = 0.019806801759223218 Iter 2: T = 962.4278718362493 K, relative_change = 0.01812431104032575 Iter 3: T = 946.582755870941 K, relative_change = 0.01646369190771342 Iter 20: T = 849.4693271185221 K, relative_change = 0.0014412352911711028 Iter 40: T = 843.1296432753894 K, relative_change = 4.579936111135971e-5 Iter 60: T = 842.9330101837132 K, relative_change = 1.402690447812578e-6 Iter 80: T = 842.9269925086252 K, relative_change = 4.291048045682194e-8 Converged in 87 iterations to T = 842.9268586500262 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: 73%|████████████████████████▎ | 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.010314200523744462 Iteration 10: norm(E*F-F'*E) = 4.035201148869848e-5 Iteration 20: norm(E*F-F'*E) = 4.477302195577377e-7 Iteration 30: norm(E*F-F'*E) = 6.1603925786228965e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.774644590291477e-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: 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.010864094255850327 Iteration 10: norm(E*F-F'*E) = 4.8521793359283174e-5 Iteration 20: norm(E*F-F'*E) = 4.915495347875414e-7 Iteration 30: norm(E*F-F'*E) = 5.819282279530765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.1519256203271324e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (5 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00: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.010864094255850327 Iteration 10: norm(E*F-F'*E) = 4.8521793359283174e-5 Iteration 20: norm(E*F-F'*E) = 4.915495347875414e-7 Iteration 30: norm(E*F-F'*E) = 5.819282279530765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.1519256203271324e-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: 44%|██████████████▌ | ETA: 0:00:01 Bin 1 progress: 91%|█████████████████████████████▉ | 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.010864094255850327 Iteration 10: norm(E*F-F'*E) = 4.8521793359283174e-5 Iteration 20: norm(E*F-F'*E) = 4.915495347875414e-7 Iteration 30: norm(E*F-F'*E) = 5.819282279530765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.1519256203271324e-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: 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.010864094255850327 Iteration 10: norm(E*F-F'*E) = 4.8521793359283174e-5 Iteration 20: norm(E*F-F'*E) = 4.915495347875414e-7 Iteration 30: norm(E*F-F'*E) = 5.819282279530765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.1519256203271324e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (50 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00: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.010864094255850327 Iteration 10: norm(E*F-F'*E) = 4.8521793359283174e-5 Iteration 20: norm(E*F-F'*E) = 4.915495347875414e-7 Iteration 30: norm(E*F-F'*E) = 5.819282279530765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.1519256203271324e-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: 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 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.010314200523744462 Iteration 10: norm(E*F-F'*E) = 4.035201148869848e-5 Iteration 20: norm(E*F-F'*E) = 4.477302195577377e-7 Iteration 30: norm(E*F-F'*E) = 6.1603925786228965e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.774644590291477e-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 9m22.3s Testing RayTraceHeatTransfer tests passed Testing completed after 578.36s PkgEval succeeded after 675.27s