Package evaluation to test RayTraceHeatTransfer on Julia 1.14.0-DEV.2360 (7c0e770819*) started at 2026-06-13T11:20:03.533 ################################################################################ # Set-up # Installing PkgEval dependencies (TestEnv)... Activating project at `~/.julia/environments/v1.14` Set-up completed after 14.73s ################################################################################ # Installation # Installing RayTraceHeatTransfer... Resolving package versions... Updating `~/.julia/environments/v1.14/Project.toml` [7cf1493d] + RayTraceHeatTransfer v0.9.0 Updating `~/.julia/environments/v1.14/Manifest.toml` [66dad0bd] + AliasTables v1.1.3 [49dc2e85] + Calculus v0.5.2 [9a962f9c] + DataAPI v1.16.0 [864edb3b] + DataStructures v0.19.5 [ffbed154] + DocStringExtensions v0.9.5 [411431e0] + Extents v0.1.6 ⌅ [5c1252a2] + GeometryBasics v0.5.10 [92d709cd] + IrrationalConstants v0.2.6 [c8e1da08] + IterTools v1.10.0 [692b3bcd] + JLLWrappers v1.8.0 ⌅ [2ab3a3ac] + LogExpFunctions v0.3.29 [eff96d63] + Measurements v2.14.1 [e1d29d7a] + Missings v1.2.0 ⌅ [bac558e1] + OrderedCollections v1.8.2 [aea7be01] + PrecompileTools v1.3.4 [21216c6a] + Preferences v1.5.2 [92933f4c] + ProgressMeter v1.11.0 [43287f4e] + PtrArrays v1.4.0 [7cf1493d] + RayTraceHeatTransfer v0.9.0 [a2af1166] + SortingAlgorithms v1.2.2 [90137ffa] + StaticArrays v1.9.18 [1e83bf80] + StaticArraysCore v1.4.4 [10745b16] + Statistics v1.11.1 [82ae8749] + StatsAPI v1.8.0 ⌅ [2913bbd2] + StatsBase v0.34.10 [5ae413db] + EarCut_jll v2.2.4+0 [0dad84c5] + ArgTools v1.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.2+0 [deac9b47] + LibCURL_jll v8.20.0+1 [e37daf67] + LibGit2_jll v1.9.4+0 [29816b5a] + LibSSH2_jll v1.11.101+0 [14a3606d] + MozillaCACerts_jll v2026.5.14 [4536629a] + OpenBLAS_jll v0.3.33+0 [458c3c95] + OpenSSL_jll v3.5.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.54s ################################################################################ # Precompilation # Precompiling PkgEval dependencies... Precompiling package dependencies... Precompiling project... 5.3 s ✓ StatsBase 22.6 s ✓ GeometryBasics 12.1 s ✓ RayTraceHeatTransfer 3 dependencies successfully precompiled in 41 seconds. 58 already precompiled. Precompilation completed after 67.98s ################################################################################ # Testing # Testing RayTraceHeatTransfer Status `/tmp/jl_nboJQj/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_nboJQj/Manifest.toml` [66dad0bd] AliasTables v1.1.3 [49dc2e85] Calculus v0.5.2 [9a962f9c] DataAPI v1.16.0 [864edb3b] DataStructures v0.19.5 [ffbed154] DocStringExtensions v0.9.5 [411431e0] Extents v0.1.6 ⌅ [5c1252a2] GeometryBasics v0.5.10 [92d709cd] IrrationalConstants v0.2.6 [c8e1da08] IterTools v1.10.0 [692b3bcd] JLLWrappers v1.8.0 ⌅ [2ab3a3ac] LogExpFunctions v0.3.29 [eff96d63] Measurements v2.14.1 [e1d29d7a] Missings v1.2.0 ⌅ [bac558e1] OrderedCollections v1.8.2 [aea7be01] PrecompileTools v1.3.4 [21216c6a] Preferences v1.5.2 [92933f4c] ProgressMeter v1.11.0 [43287f4e] PtrArrays v1.4.0 [7cf1493d] RayTraceHeatTransfer v0.9.0 [a2af1166] SortingAlgorithms v1.2.2 [90137ffa] StaticArrays v1.9.18 [1e83bf80] StaticArraysCore v1.4.4 [10745b16] Statistics v1.11.1 [82ae8749] StatsAPI v1.8.0 ⌅ [2913bbd2] StatsBase v0.34.10 [5ae413db] EarCut_jll v2.2.4+0 [0dad84c5] ArgTools v1.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.2+0 [deac9b47] LibCURL_jll v8.20.0+1 [e37daf67] LibGit2_jll v1.9.4+0 [29816b5a] LibSSH2_jll v1.11.101+0 [14a3606d] MozillaCACerts_jll v2026.5.14 [4536629a] OpenBLAS_jll v0.3.33+0 [458c3c95] OpenSSL_jll v3.5.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:02 Bin 1 progress: 62%|████████████████████▍ | ETA: 0:00:04 Bin 1 progress: 99%|████████████████████████████████▊| ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:06 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.761739987160342e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.09090909090909094) ≥ min gas E (0.03305785123966937); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.012397999407206336 Iteration 10: norm(E*F-F'*E) = 2.3875264757591455e-5 Iteration 20: norm(E*F-F'*E) = 2.7976751506059103e-7 Converged after 29 iterations. norm(E*F-F'*E) = 4.474329329647801e-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: 32%|██████████▍ | ETA: 0:00:02 Bin 1 progress: 67%|██████████████████████▎ | 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.012493616981982685 Iteration 10: norm(E*F-F'*E) = 2.1068088289139033e-5 Iteration 20: norm(E*F-F'*E) = 2.5898803366381694e-7 Converged after 29 iterations. norm(E*F-F'*E) = 4.262670329743436e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▍ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 165×165 Strategy: Serial Tolerance: 6.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.012391126851555913 Iteration 10: norm(E*F-F'*E) = 2.420465548012659e-5 Iteration 20: norm(E*F-F'*E) = 3.653511772626872e-7 Converged after 29 iterations. norm(E*F-F'*E) = 6.448729502266808e-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: 37%|████████████▎ | 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: 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.012609976696353966 Iteration 10: norm(E*F-F'*E) = 1.6989482239372132e-5 Iteration 20: norm(E*F-F'*E) = 2.013179731771388e-7 Converged after 29 iterations. norm(E*F-F'*E) = 3.385692630232166e-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: 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.01172710358958271 Iteration 10: norm(E*F-F'*E) = 2.6625721071367455e-5 Iteration 20: norm(E*F-F'*E) = 2.887437035768144e-7 Iteration 30: norm(E*F-F'*E) = 4.114897297004715e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.7737736006962738e-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: 34%|███████████▏ | ETA: 0:00:02 Bin 1 progress: 71%|███████████████████████▋ | 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.01162350304284491 Iteration 10: norm(E*F-F'*E) = 2.0338166285111414e-5 Iteration 20: norm(E*F-F'*E) = 1.6794753817382022e-7 Iteration 30: norm(E*F-F'*E) = 2.1576331607412292e-9 Converged after 31 iterations. norm(E*F-F'*E) = 9.229365551394967e-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 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.011729227009971232 Iteration 10: norm(E*F-F'*E) = 2.459824906024309e-5 Iteration 20: norm(E*F-F'*E) = 2.1735050797968644e-7 Iteration 30: norm(E*F-F'*E) = 2.6852833241290763e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.136216390582262e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▋ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374645e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.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.011833437014263027 Iteration 10: norm(E*F-F'*E) = 2.3570821091211986e-5 Iteration 20: norm(E*F-F'*E) = 2.616573840121272e-7 Iteration 30: norm(E*F-F'*E) = 3.802550734472625e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.6410381690805934e-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: 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.01150266169596017 Iteration 10: norm(E*F-F'*E) = 2.222265019710346e-5 Iteration 20: norm(E*F-F'*E) = 1.7029632400647818e-7 Converged after 29 iterations. norm(E*F-F'*E) = 1.7040153595729286e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▋ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374645e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.14285714285714296) ≥ min gas E (0.08163265306122436); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011376146651182664 Iteration 10: norm(E*F-F'*E) = 2.3450161210075634e-5 Iteration 20: norm(E*F-F'*E) = 2.5725860131816783e-7 Iteration 30: norm(E*F-F'*E) = 3.872141492913948e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.683577793990048e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Smoothing single F matrix for grey extinction Matrix size: 45×45 Strategy: Serial Tolerance: 2.793967723846438e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.03578726030307435 Iteration 10: norm(E*F-F'*E) = 0.00014507377259518427 Iteration 20: norm(E*F-F'*E) = 1.582683383191307e-6 Iteration 30: norm(E*F-F'*E) = 1.9840349236065187e-8 Converged after 35 iterations. norm(E*F-F'*E) = 1.4561982779522945e-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.028971988079772217 Iteration 10: norm(E*F-F'*E) = 0.00011336538251825414 Iteration 20: norm(E*F-F'*E) = 1.439374033526798e-6 Iteration 30: norm(E*F-F'*E) = 2.1032849186247916e-8 Converged after 33 iterations. norm(E*F-F'*E) = 3.933254635126081e-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.02352111226421881 Iteration 10: norm(E*F-F'*E) = 6.99775887567126e-5 Iteration 20: norm(E*F-F'*E) = 1.0326986273877616e-6 Iteration 30: norm(E*F-F'*E) = 1.7201959511884828e-8 Converged after 31 iterations. norm(E*F-F'*E) = 7.60032467736039e-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.02146122506774393 Iteration 10: norm(E*F-F'*E) = 5.992028146051813e-5 Iteration 20: norm(E*F-F'*E) = 9.216359475885624e-7 Iteration 30: norm(E*F-F'*E) = 1.6059119316770023e-8 Converged after 31 iterations. norm(E*F-F'*E) = 7.169446898984006e-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: 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.01172710358958271 Iteration 10: norm(E*F-F'*E) = 2.6625721071367455e-5 Iteration 20: norm(E*F-F'*E) = 2.887437035768144e-7 Iteration 30: norm(E*F-F'*E) = 4.114897297004715e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.7737736006962738e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === ✓ 2D Grey Participating Media tests complete ------------------------------------------------------------ Testing 2D Spectral Participating Media ------------------------------------------------------------ Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 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.011181253400046422 Iteration 10: norm(E*F-F'*E) = 1.242557605064265e-5 Iteration 20: norm(E*F-F'*E) = 6.791419359410676e-8 Converged after 29 iterations. norm(E*F-F'*E) = 6.849306478933814e-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: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 82%|███████████████████████████▏ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010483843835571001 Iteration 10: norm(E*F-F'*E) = 2.6521569784630582e-5 Iteration 20: norm(E*F-F'*E) = 2.9129983472977463e-7 Iteration 30: norm(E*F-F'*E) = 3.7932596304931105e-9 Converged after 33 iterations. norm(E*F-F'*E) = 6.750419206957128e-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: 80%|██████████████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011181253400046422 Iteration 10: norm(E*F-F'*E) = 1.242557605064265e-5 Iteration 20: norm(E*F-F'*E) = 6.791419359410676e-8 Converged after 29 iterations. norm(E*F-F'*E) = 6.849306478933814e-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:57:10 Bin 4 ray tracing: 8%|██▎ | ETA: 0:01:02 Bin 4 ray tracing: 15%|████▋ | ETA: 0:00:34 Bin 4 ray tracing: 23%|██████▉ | ETA: 0:00:24 Bin 4 ray tracing: 30%|█████████▏ | ETA: 0:00:19 Bin 4 ray tracing: 38%|███████████▌ | ETA: 0:00:15 Bin 4 ray tracing: 46%|█████████████▉ | ETA: 0:00:12 Bin 4 ray tracing: 54%|████████████████▎ | ETA: 0:00:09 Bin 4 ray tracing: 62%|██████████████████▊ | ETA: 0:00:07 Bin 4 ray tracing: 70%|█████████████████████▏ | ETA: 0:00:06 Bin 4 ray tracing: 78%|███████████████████████▍ | ETA: 0:00:04 Bin 4 ray tracing: 85%|█████████████████████████▋ | ETA: 0:00:03 Bin 4 ray tracing: 94%|████████████████████████████▏ | ETA: 0:00:01 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:17 Updating spectral results for spectral bin 4 Energy per ray: 0.0001853335835185918 Processing spectral bin 5/10 Bin 5 ray tracing: 8%|██▍ | ETA: 0:00:11 Bin 5 ray tracing: 16%|████▉ | ETA: 0:00:10 Bin 5 ray tracing: 25%|███████▍ | ETA: 0:00:09 Bin 5 ray tracing: 33%|█████████▉ | ETA: 0:00:08 Bin 5 ray tracing: 41%|████████████▍ | ETA: 0:00:07 Bin 5 ray tracing: 50%|██████████████▉ | ETA: 0:00:06 Bin 5 ray tracing: 58%|█████████████████▍ | ETA: 0:00:05 Bin 5 ray tracing: 66%|███████████████████▊ | ETA: 0:00:04 Bin 5 ray tracing: 74%|██████████████████████▎ | ETA: 0:00:03 Bin 5 ray tracing: 82%|████████████████████████▋ | ETA: 0:00:02 Bin 5 ray tracing: 90%|███████████████████████████ | ETA: 0:00:01 Bin 5 ray tracing: 98%|█████████████████████████████▎| 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: 15%|████▋ | ETA: 0:00:11 Bin 6 ray tracing: 23%|██████▉ | ETA: 0:00:10 Bin 6 ray tracing: 31%|█████████▎ | ETA: 0:00:09 Bin 6 ray tracing: 38%|███████████▌ | ETA: 0:00:08 Bin 6 ray tracing: 46%|█████████████▉ | ETA: 0:00:07 Bin 6 ray tracing: 54%|████████████████▏ | ETA: 0:00:06 Bin 6 ray tracing: 62%|██████████████████▌ | ETA: 0:00:05 Bin 6 ray tracing: 69%|████████████████████▊ | ETA: 0:00:04 Bin 6 ray tracing: 78%|███████████████████████▎ | ETA: 0:00:03 Bin 6 ray tracing: 86%|█████████████████████████▉ | ETA: 0:00:02 Bin 6 ray tracing: 94%|████████████████████████████▍ | ETA: 0:00:01 Bin 6 ray tracing: 100%|██████████████████████████████| Time: 0:00:12 Updating spectral results for spectral bin 6 Energy per ray: 0.013246116789219256 Processing spectral bin 7/10 Bin 7 ray tracing: 7%|██▎ | ETA: 0:00:13 Bin 7 ray tracing: 15%|████▍ | ETA: 0:00:12 Bin 7 ray tracing: 22%|██████▋ | ETA: 0:00:11 Bin 7 ray tracing: 30%|████████▉ | ETA: 0:00:10 Bin 7 ray tracing: 37%|███████████▏ | ETA: 0:00:08 Bin 7 ray tracing: 45%|█████████████▌ | 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: 8%|██▌ | ETA: 0:00:11 Bin 8 ray tracing: 17%|█████▏ | ETA: 0:00:10 Bin 8 ray tracing: 25%|███████▋ | ETA: 0:00:09 Bin 8 ray tracing: 34%|██████████▏ | ETA: 0:00:08 Bin 8 ray tracing: 42%|████████████▋ | ETA: 0:00:07 Bin 8 ray tracing: 50%|███████████████▏ | ETA: 0:00:06 Bin 8 ray tracing: 58%|█████████████████▌ | ETA: 0:00:05 Bin 8 ray tracing: 66%|███████████████████▉ | ETA: 0:00:04 Bin 8 ray tracing: 74%|██████████████████████▏ | ETA: 0:00:03 Bin 8 ray tracing: 81%|████████████████████████▍ | ETA: 0:00:02 Bin 8 ray tracing: 89%|██████████████████████████▋ | ETA: 0:00:01 Bin 8 ray tracing: 97%|█████████████████████████████ | 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:12 Bin 9 ray tracing: 16%|████▊ | ETA: 0:00:11 Bin 9 ray tracing: 24%|███████▏ | ETA: 0:00:10 Bin 9 ray tracing: 32%|█████████▋ | ETA: 0:00:09 Bin 9 ray tracing: 41%|████████████▏ | 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: 71%|█████████████████████▍ | ETA: 0:00:04 Bin 9 ray tracing: 79%|███████████████████████▋ | ETA: 0:00:03 Bin 9 ray tracing: 86%|█████████████████████████▋ | ETA: 0:00:02 Bin 9 ray tracing: 93%|███████████████████████████▉ | ETA: 0:00:01 Bin 9 ray tracing: 100%|██████████████████████████████| Time: 0:00:13 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: 16%|████▊ | ETA: 0:00:11 Bin 10 ray tracing: 24%|███████ | ETA: 0:00:10 Bin 10 ray tracing: 32%|█████████▎ | ETA: 0:00:10 Bin 10 ray tracing: 39%|███████████▍ | ETA: 0:00:09 Bin 10 ray tracing: 47%|█████████████▊ | ETA: 0:00:07 Bin 10 ray tracing: 56%|████████████████▏ | ETA: 0:00:06 Bin 10 ray tracing: 64%|██████████████████▌ | ETA: 0:00:05 Bin 10 ray tracing: 72%|█████████████████████ | ETA: 0:00:04 Bin 10 ray tracing: 81%|███████████████████████▍ | ETA: 0:00:03 Bin 10 ray tracing: 89%|█████████████████████████▉ | ETA: 0:00:01 Bin 10 ray tracing: 97%|████████████████████████████▎| ETA: 0:00:00 Bin 10 ray tracing: 100%|█████████████████████████████| Time: 0:00:12 Updating spectral results for spectral bin 10 Energy per ray: 1.5017824710260812e-5 Building intermediate mesh... Optimizing mesh... Extinction variation detected across the spectrum, ray tracing each spectral bin separately No spectral variation detected across walls Spectral variation detected across volumes, using spectral solver Building spatial acceleration structures... Computing 10 separate F matrices for variable spectral extinction Computing F matrix for spectral bin 1/10 Using 1 threads for spectral bin 1 Bin 1 progress: 22%|███████▍ | ETA: 0:00:04 Bin 1 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 1 progress: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 91%|██████████████████████████████▏ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 2/10 Using 1 threads for spectral bin 2 Bin 2 progress: 22%|███████▍ | ETA: 0:00:04 Bin 2 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 2 progress: 69%|██████████████████████▊ | ETA: 0:00:02 Bin 2 progress: 91%|██████████████████████████████▏ | ETA: 0:00:00 Bin 2 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 3/10 Using 1 threads for spectral bin 3 Bin 3 progress: 22%|███████▍ | ETA: 0:00:04 Bin 3 progress: 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: 18%|█████▉ | ETA: 0:00:05 Bin 4 progress: 38%|████████████▌ | ETA: 0:00:03 Bin 4 progress: 60%|███████████████████▊ | ETA: 0:00:02 Bin 4 progress: 82%|███████████████████████████▏ | ETA: 0:00:01 Bin 4 progress: 100%|█████████████████████████████████| Time: 0:00:05 Computing F matrix for spectral bin 5/10 Using 1 threads for spectral bin 5 Bin 5 progress: 22%|███████▍ | ETA: 0:00:04 Bin 5 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 5 progress: 67%|██████████████████████ | ETA: 0:00:02 Bin 5 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 Bin 5 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 6/10 Using 1 threads for spectral bin 6 Bin 6 progress: 22%|███████▍ | ETA: 0:00:04 Bin 6 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 6 progress: 67%|██████████████████████ | ETA: 0:00:02 Bin 6 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 Bin 6 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 7/10 Using 1 threads for spectral bin 7 Bin 7 progress: 22%|███████▍ | ETA: 0:00:04 Bin 7 progress: 44%|██████████████▋ | ETA: 0:00:03 Bin 7 progress: 67%|██████████████████████ | ETA: 0:00:02 Bin 7 progress: 89%|█████████████████████████████▍ | ETA: 0:00:01 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: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 9 progress: 91%|██████████████████████████████▏ | ETA: 0:00:00 Bin 9 progress: 100%|█████████████████████████████████| Time: 0:00:04 Computing F matrix for spectral bin 10/10 Using 1 threads for spectral bin 10 Bin 10 progress: 22%|███████▏ | ETA: 0:00:04 Bin 10 progress: 44%|██████████████▎ | ETA: 0:00:03 Bin 10 progress: 69%|██████████████████████ | ETA: 0:00:01 Bin 10 progress: 91%|█████████████████████████████▏ | ETA: 0:00:00 Bin 10 progress: 100%|████████████████████████████████| Time: 0:00:04 Smoothing F matrix for spectral bin 1/10 Matrix size: 45×45 Strategy: Serial Tolerance: 8.381945081369496e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.20000000000000007) ≥ min gas E (0.15999999999999992); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011181253400046422 Iteration 10: norm(E*F-F'*E) = 1.242557605064265e-5 Iteration 20: norm(E*F-F'*E) = 6.791419359410676e-8 Converged after 29 iterations. norm(E*F-F'*E) = 6.849306478933814e-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.010503140838683141 Iteration 10: norm(E*F-F'*E) = 2.638128619415717e-5 Iteration 20: norm(E*F-F'*E) = 2.909166684458091e-7 Iteration 30: norm(E*F-F'*E) = 3.783370876008843e-9 Converged after 33 iterations. norm(E*F-F'*E) = 6.726905662746092e-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.010506129350782285 Iteration 10: norm(E*F-F'*E) = 2.8810074222394544e-5 Iteration 20: norm(E*F-F'*E) = 2.3548982391921868e-7 Iteration 30: norm(E*F-F'*E) = 2.677904331308118e-9 Converged after 33 iterations. norm(E*F-F'*E) = 4.61980167531073e-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.01101980842714011 Iteration 10: norm(E*F-F'*E) = 3.271568588613115e-5 Iteration 20: norm(E*F-F'*E) = 2.6834842350277323e-7 Iteration 30: norm(E*F-F'*E) = 3.1123394138477966e-9 Converged after 33 iterations. norm(E*F-F'*E) = 5.481892006984367e-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.01090226838970252 Iteration 10: norm(E*F-F'*E) = 4.221058321453362e-5 Iteration 20: norm(E*F-F'*E) = 4.701996459772053e-7 Iteration 30: norm(E*F-F'*E) = 6.386671719137411e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.914024279716156e-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.010090899856316874 Iteration 10: norm(E*F-F'*E) = 4.40396011936025e-5 Iteration 20: norm(E*F-F'*E) = 5.115800729385937e-7 Iteration 30: norm(E*F-F'*E) = 6.688369676721898e-9 Converged after 35 iterations. norm(E*F-F'*E) = 5.019821727002163e-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.009977320849951664 Iteration 10: norm(E*F-F'*E) = 3.68460271388428e-5 Iteration 20: norm(E*F-F'*E) = 3.1291509393486424e-7 Iteration 30: norm(E*F-F'*E) = 3.2777766249262726e-9 Converged after 33 iterations. norm(E*F-F'*E) = 5.507271014709659e-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.010735221042844917 Iteration 10: norm(E*F-F'*E) = 2.967417862900675e-5 Iteration 20: norm(E*F-F'*E) = 2.926560533625325e-7 Iteration 30: norm(E*F-F'*E) = 3.4221372087392503e-9 Converged after 33 iterations. norm(E*F-F'*E) = 5.923763771015613e-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.010477184123791937 Iteration 10: norm(E*F-F'*E) = 2.0477542027026246e-5 Iteration 20: norm(E*F-F'*E) = 1.9584024344269946e-7 Iteration 30: norm(E*F-F'*E) = 2.6238127217182816e-9 Converged after 33 iterations. norm(E*F-F'*E) = 4.745750483961158e-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.01046242645808177 Iteration 10: norm(E*F-F'*E) = 3.9871017642155116e-5 Iteration 20: norm(E*F-F'*E) = 4.6922210066958437e-7 Iteration 30: norm(E*F-F'*E) = 6.393878361737152e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.927465097125283e-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: 29%|████████▊ | ETA: 0:00:12 Bin 4 ray tracing: 35%|██████████▌ | ETA: 0:00:11 Bin 4 ray tracing: 41%|████████████▎ | ETA: 0:00:10 Bin 4 ray tracing: 46%|██████████████ | ETA: 0:00:09 Bin 4 ray tracing: 52%|███████████████▋ | ETA: 0:00:08 Bin 4 ray tracing: 58%|█████████████████▍ | ETA: 0:00:07 Bin 4 ray tracing: 64%|███████████████████▏ | ETA: 0:00:06 Bin 4 ray tracing: 70%|█████████████████████ | ETA: 0:00:05 Bin 4 ray tracing: 76%|██████████████████████▊ | ETA: 0:00:04 Bin 4 ray tracing: 81%|████████████████████████▌ | ETA: 0:00:03 Bin 4 ray tracing: 87%|██████████████████████████▏ | ETA: 0:00:02 Bin 4 ray tracing: 93%|███████████████████████████▉ | ETA: 0:00:01 Bin 4 ray tracing: 99%|█████████████████████████████▊| ETA: 0:00:00 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:17 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: 12%|███▊ | ETA: 0:00:14 Bin 5 ray tracing: 18%|█████▌ | ETA: 0:00:13 Bin 5 ray tracing: 25%|███████▍ | ETA: 0:00:12 Bin 5 ray tracing: 31%|█████████▎ | ETA: 0:00:11 Bin 5 ray tracing: 37%|███████████▏ | ETA: 0:00:10 Bin 5 ray tracing: 43%|█████████████ | ETA: 0:00:09 Bin 5 ray tracing: 49%|██████████████▉ | 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: 73%|██████████████████████ | ETA: 0:00:04 Bin 5 ray tracing: 79%|███████████████████████▊ | ETA: 0:00:03 Bin 5 ray tracing: 85%|█████████████████████████▌ | ETA: 0:00:02 Bin 5 ray tracing: 91%|███████████████████████████▎ | ETA: 0:00:02 Bin 5 ray tracing: 96%|████████████████████████████▉ | ETA: 0:00:01 Bin 5 ray tracing: 100%|██████████████████████████████| Time: 0:00:16 Updating spectral results for spectral bin 5 Energy per ray: 0.04303963948070305 Processing spectral bin 6/10 Bin 6 ray tracing: 6%|█▊ | ETA: 0:00:17 Bin 6 ray tracing: 12%|███▌ | ETA: 0:00:16 Bin 6 ray tracing: 17%|█████▏ | ETA: 0:00:15 Bin 6 ray tracing: 23%|██████▉ | ETA: 0:00:14 Bin 6 ray tracing: 29%|████████▋ | ETA: 0:00:13 Bin 6 ray tracing: 35%|██████████▍ | ETA: 0:00:11 Bin 6 ray tracing: 41%|████████████▍ | ETA: 0:00:10 Bin 6 ray tracing: 48%|██████████████▎ | ETA: 0:00:09 Bin 6 ray tracing: 53%|████████████████ | ETA: 0:00:08 Bin 6 ray tracing: 59%|█████████████████▊ | ETA: 0:00:07 Bin 6 ray tracing: 65%|███████████████████▍ | ETA: 0:00:06 Bin 6 ray tracing: 70%|█████████████████████▏ | ETA: 0:00:05 Bin 6 ray tracing: 76%|██████████████████████▉ | ETA: 0:00:04 Bin 6 ray tracing: 82%|████████████████████████▌ | ETA: 0:00:03 Bin 6 ray tracing: 88%|██████████████████████████▍ | ETA: 0:00:02 Bin 6 ray tracing: 94%|████████████████████████████▎ | ETA: 0:00:01 Bin 6 ray tracing: 99%|█████████████████████████████▉| ETA: 0:00:00 Bin 6 ray tracing: 100%|██████████████████████████████| Time: 0:00:17 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:16 Bin 7 ray tracing: 11%|███▍ | ETA: 0:00:16 Bin 7 ray tracing: 18%|█████▎ | ETA: 0:00:14 Bin 7 ray tracing: 24%|███████▏ | ETA: 0:00:13 Bin 7 ray tracing: 30%|████████▉ | ETA: 0:00:12 Bin 7 ray tracing: 36%|██████████▊ | ETA: 0:00:11 Bin 7 ray tracing: 42%|████████████▋ | ETA: 0:00:10 Bin 7 ray tracing: 48%|██████████████▌ | ETA: 0:00:09 Bin 7 ray tracing: 54%|████████████████▍ | ETA: 0:00:08 Bin 7 ray tracing: 61%|██████████████████▏ | ETA: 0:00:07 Bin 7 ray tracing: 67%|████████████████████ | ETA: 0:00:06 Bin 7 ray tracing: 73%|█████████████████████▉ | ETA: 0:00:05 Bin 7 ray tracing: 79%|███████████████████████▋ | ETA: 0:00:04 Bin 7 ray tracing: 84%|█████████████████████████▍ | ETA: 0:00:03 Bin 7 ray tracing: 90%|███████████████████████████ | ETA: 0:00:02 Bin 7 ray tracing: 96%|████████████████████████████▊ | ETA: 0:00:01 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:17 Bin 8 ray tracing: 11%|███▍ | ETA: 0:00:16 Bin 8 ray tracing: 17%|█████▏ | ETA: 0:00:15 Bin 8 ray tracing: 22%|██████▊ | ETA: 0:00:14 Bin 8 ray tracing: 29%|████████▋ | ETA: 0:00:13 Bin 8 ray tracing: 35%|██████████▌ | ETA: 0:00:11 Bin 8 ray tracing: 41%|████████████▎ | ETA: 0:00:10 Bin 8 ray tracing: 47%|██████████████ | ETA: 0:00:09 Bin 8 ray tracing: 52%|███████████████▊ | ETA: 0:00:08 Bin 8 ray tracing: 59%|█████████████████▋ | ETA: 0:00:07 Bin 8 ray tracing: 65%|███████████████████▍ | ETA: 0:00:06 Bin 8 ray tracing: 71%|█████████████████████▎ | ETA: 0:00:05 Bin 8 ray tracing: 77%|███████████████████████ | ETA: 0:00:04 Bin 8 ray tracing: 83%|████████████████████████▉ | ETA: 0:00:03 Bin 8 ray tracing: 89%|██████████████████████████▊ | ETA: 0:00:02 Bin 8 ray tracing: 95%|████████████████████████████▋ | ETA: 0:00:01 Bin 8 ray tracing: 100%|██████████████████████████████| Time: 0:00:17 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:17 Bin 9 ray tracing: 11%|███▍ | ETA: 0:00:16 Bin 9 ray tracing: 17%|█████▏ | ETA: 0:00:15 Bin 9 ray tracing: 23%|██████▉ | ETA: 0:00:14 Bin 9 ray tracing: 29%|████████▋ | ETA: 0:00:12 Bin 9 ray tracing: 35%|██████████▌ | ETA: 0:00:11 Bin 9 ray tracing: 40%|████████████▏ | ETA: 0:00:11 Bin 9 ray tracing: 46%|█████████████▊ | ETA: 0:00:10 Bin 9 ray tracing: 52%|███████████████▌ | ETA: 0:00:09 Bin 9 ray tracing: 57%|█████████████████▏ | ETA: 0:00:08 Bin 9 ray tracing: 63%|██████████████████▉ | ETA: 0:00:07 Bin 9 ray tracing: 69%|████████████████████▋ | ETA: 0:00:06 Bin 9 ray tracing: 75%|██████████████████████▍ | ETA: 0:00:05 Bin 9 ray tracing: 81%|████████████████████████▎ | ETA: 0:00:03 Bin 9 ray tracing: 87%|██████████████████████████ | 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:17 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: 30%|████████▊ | ETA: 0:00:12 Bin 10 ray tracing: 36%|██████████▋ | ETA: 0:00:10 Bin 10 ray tracing: 43%|████████████▍ | ETA: 0:00:09 Bin 10 ray tracing: 49%|██████████████▏ | ETA: 0:00:08 Bin 10 ray tracing: 55%|████████████████ | ETA: 0:00:07 Bin 10 ray tracing: 61%|█████████████████▊ | ETA: 0:00:06 Bin 10 ray tracing: 67%|███████████████████▌ | ETA: 0:00:05 Bin 10 ray tracing: 73%|█████████████████████▎ | ETA: 0:00:04 Bin 10 ray tracing: 80%|███████████████████████▏ | ETA: 0:00:03 Bin 10 ray tracing: 85%|████████████████████████▊ | ETA: 0:00:02 Bin 10 ray tracing: 91%|██████████████████████████▍ | ETA: 0:00:02 Bin 10 ray tracing: 96%|████████████████████████████ | 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.2140796171072 K, relative_change = 0.032785920382892735 Iter 2: T = 936.498837023698 K, relative_change = 0.03175640557834784 Iter 3: T = 907.8231850650819 K, relative_change = 0.030620061472527476 Iter 20: T = 675.6208359086096 K, relative_change = 0.00621838849136509 Iter 40: T = 652.5079214660881 K, relative_change = 0.0002265702187189432 Iter 60: T = 651.7534801755946 K, relative_change = 6.9718829003940875e-6 Iter 80: T = 651.7303522867633 K, relative_change = 2.1331157348413175e-7 Converged in 96 iterations to T = 651.7296672661389 K Iter 1: T = 970.2998143159992 K, relative_change = 0.029700185684000808 Iter 2: T = 942.7630807381582 K, relative_change = 0.028379613364403918 Iter 3: T = 917.3453063886737 K, relative_change = 0.026960935222010412 Iter 20: T = 728.9420680635882 K, relative_change = 0.004109342876995156 Iter 40: T = 712.9074846026475 K, relative_change = 0.00014065075898346357 Iter 60: T = 712.396311686617 K, relative_change = 4.318335166725537e-6 Iter 80: T = 712.3806540259741 K, relative_change = 1.3211458373007097e-7 Converged in 93 iterations to T = 712.3802111385551 K Iter 1: T = 974.4066281068623 K, relative_change = 0.02559337189313767 Iter 2: T = 951.0025488862697 K, relative_change = 0.024018801335602125 Iter 3: T = 929.7130946643356 K, relative_change = 0.022386327194250384 Iter 20: T = 785.721974760731 K, relative_change = 0.0025868387989715735 Iter 40: T = 775.0479402333016 K, relative_change = 8.481421075744954e-5 Iter 60: T = 774.7130508480476 K, relative_change = 2.600229222110077e-6 Iter 80: T = 774.7027982834243 K, relative_change = 7.954752568494474e-8 Converged in 90 iterations to T = 774.7025313349718 K Iter 1: T = 970.3297627176911 K, relative_change = 0.029670237282308942 Iter 2: T = 942.8235669399096 K, relative_change = 0.028347265882829813 Iter 3: T = 917.4367398184258 K, relative_change = 0.026926381575166747 Iter 20: T = 729.4058833029088 K, relative_change = 0.0040943221767982185 Iter 40: T = 713.4227997985546 K, relative_change = 0.00014007603521906686 Iter 60: T = 712.9133511242563 K, relative_change = 4.300625304452493e-6 Iter 80: T = 712.8977463629669 K, relative_change = 1.315727100042858e-7 Converged in 93 iterations to T = 712.8973049719145 K Iter 1: T = 969.3454602335104 K, relative_change = 0.030654539766489672 Iter 2: T = 940.8324731689715 K, relative_change = 0.029414680559467565 Iter 3: T = 914.4218385485917 K, relative_change = 0.028071559362127352 Iter 20: T = 713.682450713961 K, relative_change = 0.004632173328356019 Iter 40: T = 695.8674874761005 K, relative_change = 0.00016096860876905228 Iter 60: T = 695.2963175680335 K, relative_change = 4.944761701218847e-6 Iter 80: T = 695.2788188032795 K, relative_change = 1.512818560236824e-7 Converged in 94 iterations to T = 695.2783146769085 K Iter 1: T = 963.5176058702684 K, relative_change = 0.03648239412973159 Iter 2: T = 928.9099996375905 K, relative_change = 0.0359179801405077 Iter 3: T = 896.1434511009339 K, relative_change = 0.03527419077137753 Iter 20: T = 588.9274605371637 K, relative_change = 0.012122991810052216 Iter 40: T = 546.2527213735309 K, relative_change = 0.0005379816516687375 Iter 60: T = 544.747455518197 K, relative_change = 1.669007135244164e-5 Iter 80: T = 544.701174163692 K, relative_change = 5.107779307746259e-7 Iter 100: T = 544.6997581753724 K, relative_change = 1.5625128726844872e-8 Converged in 101 iterations to T = 544.6997510261286 K Iter 1: T = 966.9865235862322 K, relative_change = 0.033013476413767756 Iter 2: T = 936.0343528333542 K, relative_change = 0.03200889567528477 Iter 3: T = 907.1129100005462 K, relative_change = 0.030897843380708754 Iter 20: T = 671.1644594059592 K, relative_change = 0.006435028074055572 Iter 40: T = 647.3357609887731 K, relative_change = 0.00023602386667000524 Iter 60: T = 646.5559824552105 K, relative_change = 7.264577491099237e-6 Iter 80: T = 646.532075704521 K, relative_change = 2.2226853983983747e-7 Converged in 96 iterations to T = 646.5313676129933 K Iter 1: T = 965.1952495334094 K, relative_change = 0.03480475046659062 Iter 2: T = 932.3657628757287 K, relative_change = 0.034013311476150596 Iter 3: T = 901.4820936340823 K, relative_change = 0.03312398467570368 Iter 20: T = 632.6071465839321 K, relative_change = 0.008650873632837773 Iter 40: T = 601.4915931363231 K, relative_change = 0.0003404009726646833 Iter 60: T = 600.4453160902693 K, relative_change = 1.0505826248766266e-5 Iter 80: T = 600.4132072629524 K, relative_change = 3.214655217417955e-7 Converged in 98 iterations to T = 600.4122378765295 K Iter 1: T = 980.1891698640486 K, relative_change = 0.0198108301359514 Iter 2: T = 962.419991433484 K, relative_change = 0.018128315407758683 Iter 3: T = 946.5712283156332 K, relative_change = 0.016467616278673408 Iter 20: T = 849.4296526047108 K, relative_change = 0.0014418128369832685 Iter 40: T = 843.0876799584363 K, relative_change = 4.581842594534066e-5 Iter 60: T = 842.8909747607117 K, relative_change = 1.4032750390749696e-6 Iter 80: T = 842.8849548778446 K, relative_change = 4.292836464551886e-8 Converged in 87 iterations to T = 842.8848209701347 K Iter 1: T = 976.4487465936817 K, relative_change = 0.023551253406318272 Iter 2: T = 955.0589158627129 K, relative_change = 0.021905738325320945 Iter 3: T = 935.7387706989597 K, relative_change = 0.020229270511862747 Iter 20: T = 809.944236848119 K, relative_change = 0.0020950870561847763 Iter 40: T = 801.0867121105754 K, relative_change = 6.776824453371909e-5 Iter 60: T = 800.8101950299227 K, relative_change = 2.076714035800398e-6 Iter 80: T = 800.8017308750561 K, relative_change = 6.353102113611668e-8 Converged in 89 iterations to T = 800.8015193911693 K Iter 1: T = 980.8380263786943 K, relative_change = 0.01916197362130572 Iter 2: T = 963.6879975705118 K, relative_change = 0.017485077400089467 Iter 3: T = 948.4242134573326 K, relative_change = 0.01583892727901528 Iter 20: T = 855.7502823131472 K, relative_change = 0.0013515857292188762 Iter 40: T = 849.767505473426 K, relative_change = 4.2847176088858316e-5 Iter 60: T = 849.5821053915431 K, relative_change = 1.3121736177305295e-6 Iter 80: T = 849.5764316437682 K, relative_change = 4.014133561986407e-8 Converged in 86 iterations to T = 849.5763154991881 K Iter 1: T = 967.302149125338 K, relative_change = 0.03269785087466196 Iter 2: T = 936.6785093973795 K, relative_change = 0.0316588149376587 Iter 3: T = 908.0977762609944 K, relative_change = 0.0305128524351143 Iter 20: T = 677.3225784198927 K, relative_change = 0.006137550616372462 Iter 40: T = 654.4770908554584 K, relative_change = 0.00022307439903081794 Iter 60: T = 653.7320794343334 K, relative_change = 6.86368532624002e-6 Iter 80: T = 653.7092413773756 K, relative_change = 2.1000057298656856e-7 Converged in 96 iterations to T = 653.7085649419123 K Iter 1: T = 973.4753682899469 K, relative_change = 0.026524631710053054 Iter 2: T = 949.1438183195409 K, relative_change = 0.024994520419297353 Iter 3: T = 926.9382410406424 K, relative_change = 0.023395376812560827 Iter 20: T = 773.8995803284134 K, relative_change = 0.002857132026551781 Iter 40: T = 762.2485823764112 K, relative_change = 9.438227057225075e-5 Iter 60: T = 761.8820261191008 K, relative_change = 2.8942866601243485e-6 Iter 80: T = 761.8708030717397 K, relative_change = 8.854415534439618e-8 Converged in 91 iterations to T = 761.8705009428749 K Iter 1: T = 970.049748845606 K, relative_change = 0.029950251154393998 Iter 2: T = 942.2577963672387 K, relative_change = 0.028650028012935112 Iter 3: T = 916.5811208751273 K, relative_change = 0.027250159766366192 Iter 20: T = 725.0347851657016 K, relative_change = 0.004237856854861892 Iter 40: T = 708.560430240022 K, relative_change = 0.000145588298766831 Iter 60: T = 708.034508091071 K, relative_change = 4.470505148377388e-6 Iter 80: T = 708.0183979005839 K, relative_change = 1.3677059267124082e-7 Converged in 93 iterations to T = 708.0179422124071 K Iter 1: T = 973.5355782347461 K, relative_change = 0.026464421765253883 Iter 2: T = 949.2641632029333 K, relative_change = 0.024931204954853996 Iter 3: T = 927.1181653276636 K, relative_change = 0.023329647040025558 Iter 20: T = 774.6802599468588 K, relative_change = 0.0028386141443445365 Iter 40: T = 763.0957687913559 K, relative_change = 9.372217267816607e-5 Iter 60: T = 762.7313744903025 K, relative_change = 2.8739950141769632e-6 Iter 80: T = 762.7202177058772 K, relative_change = 8.792333200116339e-8 Converged in 91 iterations to T = 762.7199173608934 K Iter 1: T = 964.345779792381 K, relative_change = 0.03565422020761897 Iter 2: T = 930.6183656063188 K, relative_change = 0.034974399113690914 Iter 3: T = 898.7868505894048 K, relative_change = 0.03420469248549057 Iter 20: T = 611.6274050380309 K, relative_change = 0.010167445598875809 Iter 40: T = 575.5057667306658 K, relative_change = 0.00042096174295864233 Iter 60: T = 574.266580699216 K, relative_change = 1.3019631398649746e-5 Iter 80: T = 574.2285226990898 K, relative_change = 3.9841091164972605e-7 Converged in 99 iterations to T = 574.2273653446646 K Iter 1: T = 963.559700982637 K, relative_change = 0.036440299017363016 Iter 2: T = 928.9969479966397 K, relative_change = 0.035869861463436244 Iter 3: T = 896.2781895523846 K, relative_change = 0.035219446645989985 Iter 20: T = 590.1439284131304 K, relative_change = 0.012008783177863618 Iter 40: T = 547.8561236625711 K, relative_change = 0.0005306924514238825 Iter 60: T = 546.3670243550856 K, relative_change = 1.6460779295159343e-5 Iter 80: T = 546.3212432420964 K, relative_change = 5.037577458017121e-7 Iter 100: T = 546.3198425617858 K, relative_change = 1.5410372539199803e-8 Converged in 101 iterations to T = 546.3198354898318 K Iter 1: T = 969.2967172857939 K, relative_change = 0.030703282714206134 Iter 2: T = 940.7337061774077 K, relative_change = 0.02946776833039088 Iter 3: T = 914.2720122265666 K, relative_change = 0.028128782648137417 Iter 20: T = 712.8767827727422 K, relative_change = 0.004661394091892909 Iter 40: T = 694.9629527449954 K, relative_change = 0.00016212242064876151 Iter 60: T = 694.3884286073188 K, relative_change = 4.980355171306166e-6 Iter 80: T = 694.3708268888438 K, relative_change = 1.523709565249044e-7 Converged in 94 iterations to T = 694.3703197962841 K Iter 1: T = 966.4653587276174 K, relative_change = 0.033534641272382586 Iter 2: T = 934.9692383131858 K, relative_change = 0.03258898017399952 Iter 3: T = 905.4819331190309 K, relative_change = 0.03153826242171773 Iter 20: T = 660.6168588844178 K, relative_change = 0.006977727635884665 Iter 40: T = 634.9996693532614 K, relative_change = 0.0002602639443998096 Iter 60: T = 634.1559486125775 K, relative_change = 8.0157354835911e-6 Iter 80: T = 634.1300755774541 K, relative_change = 2.4525591101134855e-7 Converged in 97 iterations to T = 634.1293012144687 K Iter 1: T = 966.3816300261892 K, relative_change = 0.03361836997381082 Iter 2: T = 934.7979487381789 K, relative_change = 0.03268241066125652 Iter 3: T = 905.2193510198554 K, relative_change = 0.03164170156583005 Iter 20: T = 658.8754162126852 K, relative_change = 0.00707157811245731 Iter 40: T = 632.9494461097837 K, relative_change = 0.00026453880290414994 Iter 60: T = 632.0945924841914 K, relative_change = 8.148304603046227e-6 Iter 80: T = 632.0683769963214 K, relative_change = 2.4931296095170535e-7 Converged in 97 iterations to T = 632.067592382976 K Iter 1: T = 976.4484258430059 K, relative_change = 0.023551574156994055 Iter 2: T = 955.0582808357449 K, relative_change = 0.02190606737759248 Iter 3: T = 935.7378305558581 K, relative_change = 0.020229603436326436 Iter 20: T = 809.9405978090467 K, relative_change = 0.0020951552229054086 Iter 40: T = 801.0828172573928 K, relative_change = 6.777057576715533e-5 Iter 60: T = 800.8062920013294 K, relative_change = 2.0767856008459564e-6 Iter 80: T = 800.7978275960294 K, relative_change = 6.353321060668334e-8 Converged in 89 iterations to T = 800.7976161058851 K Iter 1: T = 965.2079622031933 K, relative_change = 0.03479203779680666 Iter 2: T = 932.3918758678507 K, relative_change = 0.0339989801373337 Iter 3: T = 901.5223060503937 K, relative_change = 0.03310793521095858 Iter 20: T = 632.9060071512208 K, relative_change = 0.008631033556542681 Iter 40: T = 601.8556965510418 K, relative_change = 0.000339399584929492 Iter 60: T = 600.811878422619 K, relative_change = 1.0474645785256428e-5 Iter 80: T = 600.7798453598742 K, relative_change = 3.2051117793472376e-7 Converged in 98 iterations to T = 600.7788782611154 K Iter 1: T = 964.5741753372275 K, relative_change = 0.035425824662772484 Iter 2: T = 931.0886719181135 K, relative_change = 0.03471532234149553 Iter 3: T = 899.5131132979062 K, relative_change = 0.03391251507244658 Iter 20: T = 617.4749658953516 K, relative_change = 0.009719044411648614 Iter 40: T = 582.8383580935313 K, relative_change = 0.00039628801293263826 Iter 60: T = 581.6572892147882 K, relative_change = 1.2248589023146486e-5 Iter 80: T = 581.6210246063047 K, relative_change = 3.748089068958943e-7 Converged in 99 iterations to T = 581.6199217973841 K Iter 1: T = 964.3067097350693 K, relative_change = 0.035693290264930726 Iter 2: T = 930.5378776329198 K, relative_change = 0.03501876712174587 Iter 3: T = 898.6624956440853 K, relative_change = 0.03425479258288593 Iter 20: T = 610.6105966142761 K, relative_change = 0.01024758014659245 Iter 40: T = 574.2230094202415 K, relative_change = 0.0004254507845011699 Iter 60: T = 572.9733334462097 K, relative_change = 1.3160019075218858e-5 Iter 80: T = 572.9349516421871 K, relative_change = 4.0270835250968745e-7 Converged in 99 iterations to T = 572.9337844392602 K Iter 1: T = 980.2656084435611 K, relative_change = 0.01973439155643889 Iter 2: T = 962.5695050905395 K, relative_change = 0.018052355607087907 Iter 3: T = 946.7899139855649 K, relative_change = 0.016393196565571963 Iter 20: T = 850.1815403416696 K, relative_change = 0.0014308919041491168 Iter 40: T = 843.8828694603924 K, relative_change = 4.545802518300546e-5 Iter 60: T = 843.6875282817296 K, relative_change = 1.3922240500376576e-6 Iter 80: T = 843.681550162932 K, relative_change = 4.25902854808916e-8 Converged in 87 iterations to T = 843.6814171842457 K Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▌ | ETA: 0:00:02 Bin 1 progress: 76%|████████████████████████▉ | 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.011181253400046422 Iteration 10: norm(E*F-F'*E) = 1.242557605064265e-5 Iteration 20: norm(E*F-F'*E) = 6.791419359410676e-8 Converged after 29 iterations. norm(E*F-F'*E) = 6.849306478933814e-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.010317550096801084 Iteration 10: norm(E*F-F'*E) = 4.993945826378113e-5 Iteration 20: norm(E*F-F'*E) = 4.638904610562472e-7 Iteration 30: norm(E*F-F'*E) = 5.2424806248976576e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.723769409943514e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (5 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010317550096801084 Iteration 10: norm(E*F-F'*E) = 4.993945826378113e-5 Iteration 20: norm(E*F-F'*E) = 4.638904610562472e-7 Iteration 30: norm(E*F-F'*E) = 5.2424806248976576e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.723769409943514e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (10 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010317550096801084 Iteration 10: norm(E*F-F'*E) = 4.993945826378113e-5 Iteration 20: norm(E*F-F'*E) = 4.638904610562472e-7 Iteration 30: norm(E*F-F'*E) = 5.2424806248976576e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.723769409943514e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (20 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 38%|████████████▍ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010317550096801084 Iteration 10: norm(E*F-F'*E) = 4.993945826378113e-5 Iteration 20: norm(E*F-F'*E) = 4.638904610562472e-7 Iteration 30: norm(E*F-F'*E) = 5.2424806248976576e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.723769409943514e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (50 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 41%|█████████████▍ | ETA: 0:00:02 Bin 1 progress: 81%|██████████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for uniform spectral extinction Matrix size: 32×32 Strategy: Serial Tolerance: 4.76837158203125e-10 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.25) ≥ min gas E (0.25); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.010317550096801084 Iteration 10: norm(E*F-F'*E) = 4.993945826378113e-5 Iteration 20: norm(E*F-F'*E) = 4.638904610562472e-7 Iteration 30: norm(E*F-F'*E) = 5.2424806248976576e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.723769409943514e-10 === Using DIRECT spectral solver === Starting spectral steady-state direct solve... Writing spectral results to mesh... ✓ 2D Spectral Participating Media tests complete ------------------------------------------------------------ Testing Spectral Consistency ------------------------------------------------------------ Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-16 === 3D Surface-Only Grey Solver === Found 96 surfaces Populating workspace... Computing emissive powers... Computing B matrix... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures... Writing results to domain... Grey results written: 96 surfaces Computing energy conservation error... === 3D Grey Solution Complete === Computing view factors (geometry only, wavelength-independent)... Matrix size: 96×96 Strategy: Serial Tolerance: 1.0e-15 Iteration 1: norm(E*F-F'*E) = 8.179645123339839e-15 Converged after 3 iterations. norm(E*F-F'*E) = 6.099047402282222e-16 === Using DIRECT solver === === 3D Spectral Surface Radiation Solver (DIRECT) === Spectral mode: spectral_uniform Number of spectral bins: 20 Using optimized direct emission solver Computing GERT matrices for each spectral band... (Using same view factor matrix F for all bands) Setting up boundary conditions... Starting spectral direct solve... Energy conservation errors by band: [0.0, 0.0, 0.0, 0.0, 0.0, 3.1424922343134543e-19, 5.551115123125783e-15, 2.5863755581667647e-12, 1.0103917702508625e-11, 4.177991286269389e-12, 6.323830348264892e-13, 1.0125233984581428e-13, 8.104628079763643e-15, 4.371503159461554e-16, 2.4069288229178198e-17, 1.1993986533120893e-18, 5.815410578492806e-20, 1.4756902127946014e-21, 8.075350730048933e-23, 7.106459654317833e-15] === 3D Spectral Solution Complete (DIRECT) === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for uniform spectral extinction (20 bins) Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 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.011181253400046422 Iteration 10: norm(E*F-F'*E) = 1.242557605064265e-5 Iteration 20: norm(E*F-F'*E) = 6.791419359410676e-8 Converged after 29 iterations. norm(E*F-F'*E) = 6.849306478933814e-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 9m26.5s Testing RayTraceHeatTransfer tests passed Testing completed after 584.65s PkgEval succeeded after 685.21s