Package evaluation to test RayTraceHeatTransfer on Julia 1.14.0-DEV.2114 (cccbcd9611*) started at 2026-05-05T11:22:09.017 ################################################################################ # Set-up # Installing PkgEval dependencies (TestEnv)... Activating project at `~/.julia/environments/v1.14` Set-up completed after 14.75s ################################################################################ # 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.4 [ffbed154] + DocStringExtensions v0.9.5 [411431e0] + Extents v0.1.6 [5c1252a2] + GeometryBasics v0.5.10 [92d709cd] + IrrationalConstants v0.2.6 [c8e1da08] + IterTools v1.10.0 [692b3bcd] + JLLWrappers v1.7.1 [2ab3a3ac] + LogExpFunctions v0.3.29 [eff96d63] + Measurements v2.14.1 [e1d29d7a] + Missings v1.2.0 [bac558e1] + OrderedCollections v1.8.1 [aea7be01] + PrecompileTools v1.3.3 [21216c6a] + Preferences v1.5.2 [92933f4c] + ProgressMeter v1.11.0 [43287f4e] + PtrArrays v1.4.0 [7cf1493d] + RayTraceHeatTransfer v0.9.0 [a2af1166] + SortingAlgorithms v1.2.2 [90137ffa] + StaticArrays v1.9.18 [1e83bf80] + StaticArraysCore v1.4.4 [10745b16] + Statistics v1.11.1 [82ae8749] + StatsAPI v1.8.0 [2913bbd2] + StatsBase v0.34.10 [5ae413db] + EarCut_jll v2.2.4+0 [0dad84c5] + ArgTools v1.1.2 [56f22d72] + Artifacts v1.11.0 [2a0f44e3] + Base64 v1.11.0 [ade2ca70] + Dates v1.11.0 [8ba89e20] + Distributed v1.11.0 [f43a241f] + Downloads v1.7.0 [7b1f6079] + FileWatching v1.11.0 [ac6e5ff7] + JuliaSyntaxHighlighting v1.13.0 [b27032c2] + LibCURL v1.0.0 [76f85450] + LibGit2 v1.11.0 [8f399da3] + Libdl v1.11.0 [37e2e46d] + LinearAlgebra v1.13.0 [56ddb016] + Logging v1.11.0 [d6f4376e] + Markdown v1.11.0 [ca575930] + NetworkOptions v1.3.0 [44cfe95a] + Pkg v1.14.0 [de0858da] + Printf v1.11.0 [9a3f8284] + Random v1.11.0 [ea8e919c] + SHA v1.0.0 [9e88b42a] + Serialization v1.11.0 [6462fe0b] + Sockets v1.11.0 [2f01184e] + SparseArrays v1.13.0 [f489334b] + StyledStrings v1.13.0 [fa267f1f] + TOML v1.0.3 [a4e569a6] + Tar v1.10.0 [cf7118a7] + UUIDs v1.11.0 [4ec0a83e] + Unicode v1.11.0 [e66e0078] + CompilerSupportLibraries_jll v1.5.1+0 [deac9b47] + LibCURL_jll v8.19.0+0 [e37daf67] + LibGit2_jll v1.9.2+0 [29816b5a] + LibSSH2_jll v1.11.3+1 [14a3606d] + MozillaCACerts_jll v2026.3.19 [4536629a] + OpenBLAS_jll v0.3.33+0 [458c3c95] + OpenSSL_jll v3.5.6+0 [efcefdf7] + PCRE2_jll v10.47.0+0 [bea87d4a] + SuiteSparse_jll v7.10.1+0 [83775a58] + Zlib_jll v1.3.2+0 [3161d3a3] + Zstd_jll v1.5.7+1 [8e850b90] + libblastrampoline_jll v5.15.0+0 [8e850ede] + nghttp2_jll v1.69.0+0 [3f19e933] + p7zip_jll v17.8.0+0 Installation completed after 5.97s ################################################################################ # Precompilation # Precompiling PkgEval dependencies... Precompiling package dependencies... Precompiling project... 7.3 s ✓ RayTraceHeatTransfer 1 dependency successfully precompiled in 8 seconds. 60 already precompiled. Precompilation completed after 34.16s ################################################################################ # Testing # Testing RayTraceHeatTransfer Status `/tmp/jl_rDqtd5/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.13.0 [9a3f8284] Random v1.11.0 [2f01184e] SparseArrays v1.13.0 [8dfed614] Test v1.11.0 Status `/tmp/jl_rDqtd5/Manifest.toml` [66dad0bd] AliasTables v1.1.3 [49dc2e85] Calculus v0.5.2 [9a962f9c] DataAPI v1.16.0 [864edb3b] DataStructures v0.19.4 [ffbed154] DocStringExtensions v0.9.5 [411431e0] Extents v0.1.6 [5c1252a2] GeometryBasics v0.5.10 [92d709cd] IrrationalConstants v0.2.6 [c8e1da08] IterTools v1.10.0 [692b3bcd] JLLWrappers v1.7.1 [2ab3a3ac] LogExpFunctions v0.3.29 [eff96d63] Measurements v2.14.1 [e1d29d7a] Missings v1.2.0 [bac558e1] OrderedCollections v1.8.1 [aea7be01] PrecompileTools v1.3.3 [21216c6a] Preferences v1.5.2 [92933f4c] ProgressMeter v1.11.0 [43287f4e] PtrArrays v1.4.0 [7cf1493d] RayTraceHeatTransfer v0.9.0 [a2af1166] SortingAlgorithms v1.2.2 [90137ffa] StaticArrays v1.9.18 [1e83bf80] StaticArraysCore v1.4.4 [10745b16] Statistics v1.11.1 [82ae8749] StatsAPI v1.8.0 [2913bbd2] StatsBase v0.34.10 [5ae413db] EarCut_jll v2.2.4+0 [0dad84c5] ArgTools v1.1.2 [56f22d72] Artifacts v1.11.0 [2a0f44e3] Base64 v1.11.0 [ade2ca70] Dates v1.11.0 [8ba89e20] Distributed v1.11.0 [f43a241f] Downloads v1.7.0 [7b1f6079] FileWatching v1.11.0 [b77e0a4c] InteractiveUtils v1.11.0 [ac6e5ff7] JuliaSyntaxHighlighting v1.13.0 [b27032c2] LibCURL v1.0.0 [76f85450] LibGit2 v1.11.0 [8f399da3] Libdl v1.11.0 [37e2e46d] LinearAlgebra v1.13.0 [56ddb016] Logging v1.11.0 [d6f4376e] Markdown v1.11.0 [ca575930] NetworkOptions v1.3.0 [44cfe95a] Pkg v1.14.0 [de0858da] Printf v1.11.0 [9a3f8284] Random v1.11.0 [ea8e919c] SHA v1.0.0 [9e88b42a] Serialization v1.11.0 [6462fe0b] Sockets v1.11.0 [2f01184e] SparseArrays v1.13.0 [f489334b] StyledStrings v1.13.0 [fa267f1f] TOML v1.0.3 [a4e569a6] Tar v1.10.0 [8dfed614] Test v1.11.0 [cf7118a7] UUIDs v1.11.0 [4ec0a83e] Unicode v1.11.0 [e66e0078] CompilerSupportLibraries_jll v1.5.1+0 [deac9b47] LibCURL_jll v8.19.0+0 [e37daf67] LibGit2_jll v1.9.2+0 [29816b5a] LibSSH2_jll v1.11.3+1 [14a3606d] MozillaCACerts_jll v2026.3.19 [4536629a] OpenBLAS_jll v0.3.33+0 [458c3c95] OpenSSL_jll v3.5.6+0 [efcefdf7] PCRE2_jll v10.47.0+0 [bea87d4a] SuiteSparse_jll v7.10.1+0 [83775a58] Zlib_jll v1.3.2+0 [3161d3a3] Zstd_jll v1.5.7+1 [8e850b90] libblastrampoline_jll v5.15.0+0 [8e850ede] nghttp2_jll v1.69.0+0 [3f19e933] p7zip_jll v17.8.0+0 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:03:27 Bin 1 progress: 42%|█████████████▊ | ETA: 0:00:05 Bin 1 progress: 99%|████████████████████████████████▊| ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:04 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.012754438542014053 Iteration 10: norm(E*F-F'*E) = 2.7403353729333193e-5 Iteration 20: norm(E*F-F'*E) = 3.490820197472742e-7 Converged after 29 iterations. norm(E*F-F'*E) = 5.68971188124946e-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: 70%|███████████████████████ | 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.012673093725116442 Iteration 10: norm(E*F-F'*E) = 1.9835586381639356e-5 Iteration 20: norm(E*F-F'*E) = 1.7112013635505075e-7 Converged after 27 iterations. norm(E*F-F'*E) = 5.214505854118364e-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.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.012754701164499541 Iteration 10: norm(E*F-F'*E) = 2.7060680539228474e-5 Iteration 20: norm(E*F-F'*E) = 4.0054501139416443e-7 Iteration 30: norm(E*F-F'*E) = 6.914075584441387e-9 Converged after 31 iterations. norm(E*F-F'*E) = 3.076534664387112e-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: 74%|████████████████████████▍ | 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.012831672566988986 Iteration 10: norm(E*F-F'*E) = 3.834888899061628e-5 Iteration 20: norm(E*F-F'*E) = 5.802493771502401e-7 Iteration 30: norm(E*F-F'*E) = 1.0012445651975504e-8 Converged after 31 iterations. norm(E*F-F'*E) = 4.466014339503966e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 79%|██████████████████████████▏ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374624e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.1428571428571429) ≥ min gas E (0.08163265306122441); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011406926654678283 Iteration 10: norm(E*F-F'*E) = 2.7767121856409525e-5 Iteration 20: norm(E*F-F'*E) = 3.017390473420065e-7 Iteration 30: norm(E*F-F'*E) = 4.413920115365303e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.910118126780592e-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: 62%|████████████████████▋ | 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.011203202153524649 Iteration 10: norm(E*F-F'*E) = 2.9153672568834708e-5 Iteration 20: norm(E*F-F'*E) = 2.836010964731208e-7 Iteration 30: norm(E*F-F'*E) = 3.389846431858306e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.4145911982007704e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 42%|█████████████▊ | ETA: 0:00:01 Bin 1 progress: 82%|███████████████████████████ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.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.011254082851074055 Iteration 10: norm(E*F-F'*E) = 2.8662986257905484e-5 Iteration 20: norm(E*F-F'*E) = 3.071919100568651e-7 Iteration 30: norm(E*F-F'*E) = 4.1727075180388965e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.7843815067384435e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 39%|████████████▉ | ETA: 0:00:02 Bin 1 progress: 78%|█████████████████████████▊ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.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.011273495263163522 Iteration 10: norm(E*F-F'*E) = 3.645501546381811e-5 Iteration 20: norm(E*F-F'*E) = 3.568255648356876e-7 Iteration 30: norm(E*F-F'*E) = 4.4732540738913e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.8894638352420103e-9 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 40%|█████████████▎ | ETA: 0:00:02 Bin 1 progress: 82%|███████████████████████████ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374624e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.1428571428571429) ≥ min gas E (0.08163265306122441); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011199570170478675 Iteration 10: norm(E*F-F'*E) = 4.917755228110711e-5 Iteration 20: norm(E*F-F'*E) = 5.791184052101956e-7 Iteration 30: norm(E*F-F'*E) = 7.915719425678834e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.4499346598776505e-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: 74%|████████████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:02 Smoothing single F matrix for grey extinction Matrix size: 77×77 Strategy: Serial Tolerance: 2.0079324748374645e-9 ┌ Warning: Algorithm 1 convergence check failed: max surface E (0.14285714285714296) ≥ min gas E (0.08163265306122436); convergence not guaranteed, consider refining the mesh └ @ RayTraceHeatTransfer ~/.julia/packages/RayTraceHeatTransfer/VDrcq/src/HeatTransfer/exchangeFactorSmoothing/smoothExchangeFactors.jl:285 Iteration 1: norm(E*F-F'*E) = 0.011139615020015994 Iteration 10: norm(E*F-F'*E) = 4.033597926749573e-5 Iteration 20: norm(E*F-F'*E) = 4.924491704941697e-7 Iteration 30: norm(E*F-F'*E) = 6.99689826344582e-9 Converged after 33 iterations. norm(E*F-F'*E) = 1.2967754836527434e-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.03524653741858906 Iteration 10: norm(E*F-F'*E) = 0.00010074766501363204 Iteration 20: norm(E*F-F'*E) = 8.5551486825656e-7 Iteration 30: norm(E*F-F'*E) = 1.1142364630863217e-8 Converged after 33 iterations. norm(E*F-F'*E) = 2.0214727000413364e-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.02901465972516473 Iteration 10: norm(E*F-F'*E) = 8.702086010043505e-5 Iteration 20: norm(E*F-F'*E) = 1.142234985593809e-6 Iteration 30: norm(E*F-F'*E) = 1.767428614664339e-8 Converged after 33 iterations. norm(E*F-F'*E) = 3.3589776429763255e-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.02476334875473252 Iteration 10: norm(E*F-F'*E) = 5.2736125736974514e-5 Iteration 20: norm(E*F-F'*E) = 6.316381852390321e-7 Iteration 30: norm(E*F-F'*E) = 1.0345996324596481e-8 Converged after 31 iterations. norm(E*F-F'*E) = 4.5796287588279604e-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.02172289642794673 Iteration 10: norm(E*F-F'*E) = 5.3500681660744296e-5 Iteration 20: norm(E*F-F'*E) = 6.938374828323008e-7 Converged after 29 iterations. norm(E*F-F'*E) = 1.0839948862538186e-8 === Variable Extinction Memory-Optimized Steady State Solver === Allocating workspace... Populating workspace from mesh... Computing emissive powers with variable extinction... Computing B matrix with variable extinction... Computing K matrix... Solving for S_infty... Assembling linear system... Solving linear system... Computing absorbed and reflected energies... Computing temperatures with variable extinction... Writing results to mesh... Computing energy conservation error... === Variable Extinction Steady State Solution Complete === Building intermediate mesh... Optimizing mesh... No extinction variation detected across the spectrum, ray tracing grey domain only No spectral variation detected across walls No spectral variation detected across volumes No spectral variation detected across mesh, using efficient grey solver Building spatial acceleration structures... Computing single F matrix for grey extinction Using 1 threads for spectral bin 1 Bin 1 progress: 43%|██████████████▏ | ETA: 0:00:01 Bin 1 progress: 91%|██████████████████████████████ | 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.011406926654678283 Iteration 10: norm(E*F-F'*E) = 2.7767121856409525e-5 Iteration 20: norm(E*F-F'*E) = 3.017390473420065e-7 Iteration 30: norm(E*F-F'*E) = 4.413920115365303e-9 Converged after 31 iterations. norm(E*F-F'*E) = 1.910118126780592e-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: 58%|███████████████████▏ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010518564707839467 Iteration 10: norm(E*F-F'*E) = 5.683709284498492e-5 Iteration 20: norm(E*F-F'*E) = 6.759036373708808e-7 Iteration 30: norm(E*F-F'*E) = 8.675777522537765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.420214670317196e-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: 56%|██████████████████▍ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.011705586946480452 Iteration 10: norm(E*F-F'*E) = 3.885765472011106e-5 Iteration 20: norm(E*F-F'*E) = 3.460780482047115e-7 Iteration 30: norm(E*F-F'*E) = 4.165512178051003e-9 Converged after 33 iterations. norm(E*F-F'*E) = 7.268318279908246e-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: 62%|████████████████████▌ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010518564707839467 Iteration 10: norm(E*F-F'*E) = 5.683709284498492e-5 Iteration 20: norm(E*F-F'*E) = 6.759036373708808e-7 Iteration 30: norm(E*F-F'*E) = 8.675777522537765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.420214670317196e-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: 7:34:17 Bin 4 ray tracing: 13%|███▉ | ETA: 0:00:25 Bin 4 ray tracing: 26%|███████▊ | ETA: 0:00:14 Bin 4 ray tracing: 39%|███████████▋ | ETA: 0:00:09 Bin 4 ray tracing: 53%|███████████████▉ | ETA: 0:00:06 Bin 4 ray tracing: 67%|████████████████████ | ETA: 0:00:04 Bin 4 ray tracing: 79%|███████████████████████▊ | ETA: 0:00:02 Bin 4 ray tracing: 92%|███████████████████████████▊ | ETA: 0:00:01 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:10 Updating spectral results for spectral bin 4 Energy per ray: 0.0001853335835185918 Processing spectral bin 5/10 Bin 5 ray tracing: 13%|████ | ETA: 0:00:07 Bin 5 ray tracing: 27%|████████ | ETA: 0:00:06 Bin 5 ray tracing: 40%|████████████ | ETA: 0:00:05 Bin 5 ray tracing: 52%|███████████████▋ | ETA: 0:00:04 Bin 5 ray tracing: 66%|███████████████████▊ | ETA: 0:00:03 Bin 5 ray tracing: 79%|███████████████████████▊ | ETA: 0:00:02 Bin 5 ray tracing: 92%|███████████████████████████▊ | ETA: 0:00:01 Bin 5 ray tracing: 100%|██████████████████████████████| Time: 0:00:07 Updating spectral results for spectral bin 5 Energy per ray: 0.04303963948070305 Processing spectral bin 6/10 Bin 6 ray tracing: 13%|███▉ | ETA: 0:00:07 Bin 6 ray tracing: 26%|███████▉ | ETA: 0:00:06 Bin 6 ray tracing: 39%|███████████▉ | ETA: 0:00:05 Bin 6 ray tracing: 53%|███████████████▉ | ETA: 0:00:04 Bin 6 ray tracing: 66%|███████████████████▊ | ETA: 0:00:03 Bin 6 ray tracing: 78%|███████████████████████▌ | ETA: 0:00:02 Bin 6 ray tracing: 91%|███████████████████████████▎ | ETA: 0:00:01 Bin 6 ray tracing: 100%|██████████████████████████████| Time: 0:00:07 Updating spectral results for spectral bin 6 Energy per ray: 0.013246116789219256 Processing spectral bin 7/10 Bin 7 ray tracing: 13%|████ | ETA: 0:00:06 Bin 7 ray tracing: 27%|████████ | ETA: 0:00:05 Bin 7 ray tracing: 40%|████████████ | ETA: 0:00:05 Bin 7 ray tracing: 54%|████████████████▏ | ETA: 0:00:04 Bin 7 ray tracing: 67%|████████████████████▏ | ETA: 0:00:03 Bin 7 ray tracing: 80%|████████████████████████ | ETA: 0:00:02 Bin 7 ray tracing: 94%|████████████████████████████ | ETA: 0:00:00 Bin 7 ray tracing: 100%|██████████████████████████████| Time: 0:00:07 Updating spectral results for spectral bin 7 Energy per ray: 0.00021661482457376898 Processing spectral bin 8/10 Bin 8 ray tracing: 13%|████ | ETA: 0:00:07 Bin 8 ray tracing: 27%|████████ | ETA: 0:00:05 Bin 8 ray tracing: 40%|████████████▏ | ETA: 0:00:04 Bin 8 ray tracing: 54%|████████████████▎ | ETA: 0:00:03 Bin 8 ray tracing: 68%|████████████████████▎ | ETA: 0:00:02 Bin 8 ray tracing: 81%|████████████████████████▍ | ETA: 0:00:01 Bin 8 ray tracing: 95%|████████████████████████████▌ | ETA: 0:00:00 Bin 8 ray tracing: 100%|██████████████████████████████| Time: 0:00:07 Updating spectral results for spectral bin 8 Energy per ray: 1.0195075180910974e-6 Processing spectral bin 9/10 Bin 9 ray tracing: 14%|████ | ETA: 0:00:06 Bin 9 ray tracing: 27%|████████▏ | ETA: 0:00:05 Bin 9 ray tracing: 41%|████████████▎ | ETA: 0:00:04 Bin 9 ray tracing: 55%|████████████████▍ | ETA: 0:00:03 Bin 9 ray tracing: 68%|████████████████████▍ | ETA: 0:00:02 Bin 9 ray tracing: 81%|████████████████████████▌ | ETA: 0:00:01 Bin 9 ray tracing: 95%|████████████████████████████▌ | ETA: 0:00:00 Bin 9 ray tracing: 100%|██████████████████████████████| Time: 0:00:07 Updating spectral results for spectral bin 9 Energy per ray: 2.17242363082386e-9 Processing spectral bin 10/10 Bin 10 ray tracing: 14%|███▉ | ETA: 0:00:07 Bin 10 ray tracing: 27%|███████▊ | ETA: 0:00:06 Bin 10 ray tracing: 40%|███████████▌ | ETA: 0:00:05 Bin 10 ray tracing: 53%|███████████████▍ | ETA: 0:00:04 Bin 10 ray tracing: 66%|███████████████████▎ | ETA: 0:00:03 Bin 10 ray tracing: 79%|███████████████████████ | ETA: 0:00:02 Bin 10 ray tracing: 92%|██████████████████████████▋ | ETA: 0:00:01 Bin 10 ray tracing: 100%|█████████████████████████████| Time: 0:00:07 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: 27%|████████▊ | ETA: 0:00:03 Bin 1 progress: 60%|███████████████████▊ | ETA: 0:00:01 Bin 1 progress: 98%|████████████████████████████████▎| ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:03 Computing F matrix for spectral bin 2/10 Using 1 threads for spectral bin 2 Bin 2 progress: 36%|███████████▊ | ETA: 0:00:02 Bin 2 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 2 progress: 100%|█████████████████████████████████| Time: 0:00:02 Computing F matrix for spectral bin 3/10 Using 1 threads for spectral bin 3 Bin 3 progress: 36%|███████████▊ | ETA: 0:00:02 Bin 3 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 3 progress: 100%|█████████████████████████████████| Time: 0:00:02 Computing F matrix for spectral bin 4/10 Using 1 threads for spectral bin 4 Bin 4 progress: 31%|██████████▎ | ETA: 0:00:02 Bin 4 progress: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 4 progress: 100%|█████████████████████████████████| Time: 0:00:02 Computing F matrix for spectral bin 5/10 Using 1 threads for spectral bin 5 Bin 5 progress: 36%|███████████▊ | ETA: 0:00:02 Bin 5 progress: 71%|███████████████████████▌ | ETA: 0:00:01 Bin 5 progress: 100%|█████████████████████████████████| Time: 0:00:02 Computing F matrix for spectral bin 6/10 Using 1 threads for spectral bin 6 Bin 6 progress: 38%|████████████▌ | ETA: 0:00:02 Bin 6 progress: 76%|████████████████████████▉ | ETA: 0:00:01 Bin 6 progress: 100%|█████████████████████████████████| Time: 0:00:02 Computing F matrix for spectral bin 7/10 Using 1 threads for spectral bin 7 Bin 7 progress: 36%|███████████▊ | ETA: 0:00:02 Bin 7 progress: 73%|████████████████████████▎ | ETA: 0:00:01 Bin 7 progress: 100%|█████████████████████████████████| Time: 0:00:02 Computing F matrix for spectral bin 8/10 Using 1 threads for spectral bin 8 Bin 8 progress: 36%|███████████▊ | ETA: 0:00:02 Bin 8 progress: 69%|██████████████████████▊ | ETA: 0:00:01 Bin 8 progress: 100%|█████████████████████████████████| Time: 0:00:03 Computing F matrix for spectral bin 9/10 Using 1 threads for spectral bin 9 Bin 9 progress: 29%|█████████▌ | ETA: 0:00:03 Bin 9 progress: 67%|██████████████████████ | ETA: 0:00:01 Bin 9 progress: 100%|█████████████████████████████████| Time: 0:00:03 Computing F matrix for spectral bin 10/10 Using 1 threads for spectral bin 10 Bin 10 progress: 27%|████████▌ | ETA: 0:00:03 Bin 10 progress: 53%|█████████████████▏ | ETA: 0:00:02 Bin 10 progress: 78%|████████████████████████▉ | ETA: 0:00:01 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.010518564707839467 Iteration 10: norm(E*F-F'*E) = 5.683709284498492e-5 Iteration 20: norm(E*F-F'*E) = 6.759036373708808e-7 Iteration 30: norm(E*F-F'*E) = 8.675777522537765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.420214670317196e-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.011719251552722184 Iteration 10: norm(E*F-F'*E) = 3.824682577125913e-5 Iteration 20: norm(E*F-F'*E) = 3.3659513516737886e-7 Iteration 30: norm(E*F-F'*E) = 4.04324147682107e-9 Converged after 33 iterations. norm(E*F-F'*E) = 7.056419611966807e-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.010176576191344081 Iteration 10: norm(E*F-F'*E) = 3.9505367647341385e-5 Iteration 20: norm(E*F-F'*E) = 4.2072461440424105e-7 Iteration 30: norm(E*F-F'*E) = 5.30535390074006e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.936390825639801e-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.009887509081617509 Iteration 10: norm(E*F-F'*E) = 2.7241600036606982e-5 Iteration 20: norm(E*F-F'*E) = 2.414524497702094e-7 Iteration 30: norm(E*F-F'*E) = 2.927212639525717e-9 Converged after 33 iterations. norm(E*F-F'*E) = 5.167982034253975e-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.010438146478284277 Iteration 10: norm(E*F-F'*E) = 3.201560398755085e-5 Iteration 20: norm(E*F-F'*E) = 3.1817881215456386e-7 Iteration 30: norm(E*F-F'*E) = 4.170145928653666e-9 Converged after 33 iterations. norm(E*F-F'*E) = 7.501025301262153e-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.010453879173439278 Iteration 10: norm(E*F-F'*E) = 3.835021897153918e-5 Iteration 20: norm(E*F-F'*E) = 4.610269525866596e-7 Iteration 30: norm(E*F-F'*E) = 6.290500989342614e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.827244643771424e-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.011018515964291218 Iteration 10: norm(E*F-F'*E) = 4.355292371654655e-5 Iteration 20: norm(E*F-F'*E) = 4.744814828385083e-7 Iteration 30: norm(E*F-F'*E) = 6.259899138476903e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.746430695504521e-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.01006711484487524 Iteration 10: norm(E*F-F'*E) = 5.595070111487897e-5 Iteration 20: norm(E*F-F'*E) = 7.152026808510183e-7 Iteration 30: norm(E*F-F'*E) = 9.756650455498742e-9 Converged after 35 iterations. norm(E*F-F'*E) = 7.488926845402742e-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.010629161030126209 Iteration 10: norm(E*F-F'*E) = 4.0708508719348755e-5 Iteration 20: norm(E*F-F'*E) = 4.4886803797402584e-7 Iteration 30: norm(E*F-F'*E) = 6.080306982938471e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.696831139194944e-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.011034445435588038 Iteration 10: norm(E*F-F'*E) = 4.3578297535608444e-5 Iteration 20: norm(E*F-F'*E) = 4.693758115437744e-7 Iteration 30: norm(E*F-F'*E) = 5.943929657822851e-9 Converged after 35 iterations. norm(E*F-F'*E) = 4.4222431638866315e-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: 10%|███ | ETA: 0:00:13 Bin 4 ray tracing: 17%|█████▏ | ETA: 0:00:12 Bin 4 ray tracing: 24%|███████▏ | ETA: 0:00:11 Bin 4 ray tracing: 30%|█████████ | ETA: 0:00:11 Bin 4 ray tracing: 37%|███████████ | ETA: 0:00:10 Bin 4 ray tracing: 43%|█████████████ | ETA: 0:00:09 Bin 4 ray tracing: 50%|██████████████▉ | ETA: 0:00:08 Bin 4 ray tracing: 56%|████████████████▊ | ETA: 0:00:07 Bin 4 ray tracing: 63%|██████████████████▊ | ETA: 0:00:06 Bin 4 ray tracing: 70%|████████████████████▉ | ETA: 0:00:05 Bin 4 ray tracing: 77%|███████████████████████ | ETA: 0:00:04 Bin 4 ray tracing: 83%|█████████████████████████ | ETA: 0:00:03 Bin 4 ray tracing: 90%|██████████████████████████▉ | ETA: 0:00:02 Bin 4 ray tracing: 96%|████████████████████████████▊ | ETA: 0:00:01 Bin 4 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 4 Energy per ray: 0.0001853335835185918 Processing spectral bin 5/10 Bin 5 ray tracing: 7%|██ | ETA: 0:00:15 Bin 5 ray tracing: 13%|███▉ | ETA: 0:00:14 Bin 5 ray tracing: 19%|█████▊ | ETA: 0:00:13 Bin 5 ray tracing: 26%|███████▋ | ETA: 0:00:12 Bin 5 ray tracing: 32%|█████████▋ | ETA: 0:00:11 Bin 5 ray tracing: 39%|███████████▊ | ETA: 0:00:09 Bin 5 ray tracing: 46%|█████████████▊ | ETA: 0:00:08 Bin 5 ray tracing: 52%|███████████████▊ | ETA: 0:00:07 Bin 5 ray tracing: 59%|█████████████████▋ | ETA: 0:00:06 Bin 5 ray tracing: 65%|███████████████████▌ | ETA: 0:00:05 Bin 5 ray tracing: 72%|█████████████████████▋ | ETA: 0:00:04 Bin 5 ray tracing: 79%|███████████████████████▋ | ETA: 0:00:03 Bin 5 ray tracing: 86%|█████████████████████████▊ | ETA: 0:00:02 Bin 5 ray tracing: 93%|████████████████████████████ | ETA: 0:00:01 Bin 5 ray tracing: 99%|██████████████████████████████| ETA: 0:00:00 Bin 5 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 5 Energy per ray: 0.04303963948070305 Processing spectral bin 6/10 Bin 6 ray tracing: 6%|█▉ | ETA: 0:00:15 Bin 6 ray tracing: 13%|███▉ | ETA: 0:00:14 Bin 6 ray tracing: 19%|█████▊ | ETA: 0:00:13 Bin 6 ray tracing: 26%|███████▊ | ETA: 0:00:12 Bin 6 ray tracing: 32%|█████████▋ | ETA: 0:00:11 Bin 6 ray tracing: 39%|███████████▋ | ETA: 0:00:10 Bin 6 ray tracing: 45%|█████████████▌ | ETA: 0:00:09 Bin 6 ray tracing: 52%|███████████████▌ | ETA: 0:00:08 Bin 6 ray tracing: 58%|█████████████████▌ | ETA: 0:00:07 Bin 6 ray tracing: 65%|███████████████████▌ | ETA: 0:00:05 Bin 6 ray tracing: 72%|█████████████████████▋ | ETA: 0:00:04 Bin 6 ray tracing: 79%|███████████████████████▋ | ETA: 0:00:03 Bin 6 ray tracing: 85%|█████████████████████████▌ | ETA: 0:00:02 Bin 6 ray tracing: 92%|███████████████████████████▌ | ETA: 0:00:01 Bin 6 ray tracing: 98%|█████████████████████████████▌| ETA: 0:00:00 Bin 6 ray tracing: 100%|██████████████████████████████| Time: 0:00:15 Updating spectral results for spectral bin 6 Energy per ray: 0.013246116789219256 Processing spectral bin 7/10 Bin 7 ray tracing: 9%|██▋ | ETA: 0:00:12 Bin 7 ray tracing: 17%|█████ | ETA: 0:00:11 Bin 7 ray tracing: 25%|███████▌ | ETA: 0:00:10 Bin 7 ray tracing: 34%|██████████▏ | ETA: 0:00:08 Bin 7 ray tracing: 43%|████████████▉ | ETA: 0:00:07 Bin 7 ray tracing: 52%|███████████████▋ | ETA: 0:00:06 Bin 7 ray tracing: 62%|██████████████████▋ | ETA: 0:00:04 Bin 7 ray tracing: 72%|█████████████████████▊ | ETA: 0:00:03 Bin 7 ray tracing: 83%|████████████████████████▉ | ETA: 0:00:02 Bin 7 ray tracing: 93%|████████████████████████████ | ETA: 0:00:01 Bin 7 ray tracing: 100%|██████████████████████████████| Time: 0:00:10 Updating spectral results for spectral bin 7 Energy per ray: 0.00021661482457376898 Processing spectral bin 8/10 Bin 8 ray tracing: 9%|██▋ | ETA: 0:00:11 Bin 8 ray tracing: 18%|█████▍ | ETA: 0:00:09 Bin 8 ray tracing: 29%|████████▋ | ETA: 0:00:07 Bin 8 ray tracing: 39%|███████████▊ | ETA: 0:00:06 Bin 8 ray tracing: 49%|██████████████▊ | ETA: 0:00:05 Bin 8 ray tracing: 60%|██████████████████ | ETA: 0:00:04 Bin 8 ray tracing: 71%|█████████████████████▏ | ETA: 0:00:03 Bin 8 ray tracing: 81%|████████████████████████▍ | ETA: 0:00:02 Bin 8 ray tracing: 92%|███████████████████████████▌ | ETA: 0:00:01 Bin 8 ray tracing: 100%|██████████████████████████████| Time: 0:00:10 Updating spectral results for spectral bin 8 Energy per ray: 1.0195075180910974e-6 Processing spectral bin 9/10 Bin 9 ray tracing: 10%|███▏ | ETA: 0:00:09 Bin 9 ray tracing: 20%|██████▏ | ETA: 0:00:08 Bin 9 ray tracing: 30%|█████████▏ | ETA: 0:00:07 Bin 9 ray tracing: 40%|████████████▏ | ETA: 0:00:06 Bin 9 ray tracing: 51%|███████████████▎ | ETA: 0:00:05 Bin 9 ray tracing: 62%|██████████████████▌ | ETA: 0:00:04 Bin 9 ray tracing: 72%|█████████████████████▋ | ETA: 0:00:03 Bin 9 ray tracing: 83%|████████████████████████▊ | ETA: 0:00:02 Bin 9 ray tracing: 93%|███████████████████████████▉ | ETA: 0:00:01 Bin 9 ray tracing: 100%|██████████████████████████████| Time: 0:00:09 Updating spectral results for spectral bin 9 Energy per ray: 2.17242363082386e-9 Processing spectral bin 10/10 Bin 10 ray tracing: 11%|███ | ETA: 0:00:09 Bin 10 ray tracing: 21%|██████▏ | ETA: 0:00:08 Bin 10 ray tracing: 31%|█████████▏ | ETA: 0:00:07 Bin 10 ray tracing: 42%|████████████▎ | ETA: 0:00:06 Bin 10 ray tracing: 53%|███████████████▍ | ETA: 0:00:04 Bin 10 ray tracing: 64%|██████████████████▌ | ETA: 0:00:04 Bin 10 ray tracing: 73%|█████████████████████▎ | ETA: 0:00:03 Bin 10 ray tracing: 84%|████████████████████████▎ | ETA: 0:00:02 Bin 10 ray tracing: 94%|███████████████████████████▎ | ETA: 0:00:01 Bin 10 ray tracing: 100%|█████████████████████████████| Time: 0:00:09 Updating spectral results for spectral bin 10 Energy per ray: 1.5017824710260812e-5 Iter 1: T = 967.305417730812 K, relative_change = 0.03269458226918802 Iter 2: T = 936.6851767375069 K, relative_change = 0.03165519434920215 Iter 3: T = 908.107964181941 K, relative_change = 0.030508876691206905 Iter 20: T = 677.3854950701672 K, relative_change = 0.0061345815526851755 Iter 40: T = 654.5498334664484 K, relative_change = 0.00022294632810979255 Iter 60: T = 653.805168147979 K, relative_change = 6.859721834445629e-6 Iter 80: T = 653.7823407282485 K, relative_change = 2.09879284821979e-7 Converged in 96 iterations to T = 653.781664607872 K Iter 1: T = 970.2908177220922 K, relative_change = 0.02970918227790775 Iter 2: T = 942.744909332799 K, relative_change = 0.028389332235423476 Iter 3: T = 917.3178358626076 K, relative_change = 0.026971318771889914 Iter 20: T = 728.8025659939848 K, relative_change = 0.004113870304121787 Iter 40: T = 712.7524639682982 K, relative_change = 0.00014082408509574537 Iter 60: T = 712.2407713614839 K, relative_change = 4.323676239333337e-6 Iter 80: T = 712.2250977567156 K, relative_change = 1.322780061414647e-7 Converged in 93 iterations to T = 712.2246544182841 K Iter 1: T = 974.4162276910746 K, relative_change = 0.025583772308925373 Iter 2: T = 951.0216797778844 K, relative_change = 0.024008783154837877 Iter 3: T = 929.7416094701689 K, relative_change = 0.02237600967486416 Iter 20: T = 785.841123865936 K, relative_change = 0.0025842220723520863 Iter 40: T = 775.1766154664787 K, relative_change = 8.472227954183204e-5 Iter 60: T = 774.8420338989473 K, relative_change = 2.597404588315833e-6 Iter 80: T = 774.8317907669377 K, relative_change = 7.946110720214266e-8 Converged in 90 iterations to T = 774.8315240640918 K Iter 1: T = 970.3160345418481 K, relative_change = 0.02968396545815182 Iter 2: T = 942.79584114884 K, relative_change = 0.02836209277527006 Iter 3: T = 917.3948295687533 K, relative_change = 0.0269422185285995 Iter 20: T = 729.1933812889499 K, relative_change = 0.004101197964473577 Iter 40: T = 713.1867208241389 K, relative_change = 0.00014033905593309508 Iter 60: T = 712.6764828599453 K, relative_change = 4.308730107300002e-6 Iter 80: T = 712.6608538834979 K, relative_change = 1.3182069490730633e-7 Converged in 93 iterations to T = 712.6604118074702 K Iter 1: T = 969.3425863484777 K, relative_change = 0.030657413651522396 Iter 2: T = 940.826650304709 K, relative_change = 0.029417810014092546 Iter 3: T = 914.4130061736314 K, relative_change = 0.02807493189369469 Iter 20: T = 713.6350229639135 K, relative_change = 0.004633888802890332 Iter 40: T = 695.8142538077348 K, relative_change = 0.00016103629185737554 Iter 60: T = 695.242886990122 K, relative_change = 4.946849570024845e-6 Iter 80: T = 695.2253821815555 K, relative_change = 1.5134574124848712e-7 Converged in 94 iterations to T = 695.2248778810571 K Iter 1: T = 963.5301797037033 K, relative_change = 0.03646982029629668 Iter 2: T = 928.9359724429795 K, relative_change = 0.03590360529377701 Iter 3: T = 896.1837018035059 K, relative_change = 0.03525783435142411 Iter 20: T = 589.2915849351762 K, relative_change = 0.012088686273904546 Iter 40: T = 546.7331348980734 K, relative_change = 0.0005357858200266823 Iter 60: T = 545.232733965467 K, relative_change = 1.662098912906491e-5 Iter 80: T = 545.186603154664 K, relative_change = 5.086628465655474e-7 Iter 100: T = 545.185191773176 K, relative_change = 1.556042569861936e-8 Converged in 101 iterations to T = 545.185184647192 K Iter 1: T = 966.9355930501414 K, relative_change = 0.03306440694985864 Iter 2: T = 935.9303461980428 K, relative_change = 0.032065472690165804 Iter 3: T = 906.9537852223888 K, relative_change = 0.030960168236197606 Iter 20: T = 670.155041634734 K, relative_change = 0.006485114733781875 Iter 40: T = 646.1610260195697 K, relative_change = 0.00023822735599939458 Iter 60: T = 645.3753752358739 K, relative_change = 7.332820529414856e-6 Iter 80: T = 645.3512879506364 K, relative_change = 2.2435691571009126e-7 Converged in 96 iterations to T = 645.3505745114197 K Iter 1: T = 965.1854449863501 K, relative_change = 0.03481455501364984 Iter 2: T = 932.3456226801712 K, relative_change = 0.034024365448904324 Iter 3: T = 901.4510776550453 K, relative_change = 0.0331363651778777 Iter 20: T = 632.3763700897027 K, relative_change = 0.008666225732640429 Iter 40: T = 601.2103294483193 K, relative_change = 0.00034117672578605566 Iter 60: T = 600.1621486603333 K, relative_change = 1.0529982193316878e-5 Iter 80: T = 600.129981173477 K, relative_change = 3.222048663913805e-7 Converged in 98 iterations to T = 600.1290100158649 K Iter 1: T = 980.1422709917096 K, relative_change = 0.019857729008290408 Iter 2: T = 962.328239343813 K, relative_change = 0.01817494477599892 Iter 3: T = 946.437001439047 K, relative_change = 0.016513323889987795 Iter 20: T = 848.9673510224037 K, relative_change = 0.0014485531989214349 Iter 40: T = 842.5986761538657 K, relative_change = 4.604096957372609e-5 Iter 60: T = 842.4011296677166 K, relative_change = 1.4100990103753485e-6 Iter 80: T = 842.3950840257955 K, relative_change = 4.313712819437033e-8 Converged in 87 iterations to T = 842.3949495450864 K Iter 1: T = 976.4043392998155 K, relative_change = 0.023595660700184463 Iter 2: T = 954.9709913978613 K, relative_change = 0.02195130340912268 Iter 3: T = 935.6085911076311 K, relative_change = 0.02027538057662673 Iter 20: T = 809.4399567015928 K, relative_change = 0.0021045486664948657 Iter 40: T = 800.5469363197257 K, relative_change = 6.809190505760207e-5 Iter 60: T = 800.2692847362946 K, relative_change = 2.0866499656735114e-6 Iter 80: T = 800.2607858286174 K, relative_change = 6.383499845415941e-8 Converged in 89 iterations to T = 800.2605734763835 K Iter 1: T = 980.9181316271962 K, relative_change = 0.01908186837280377 Iter 2: T = 963.8443591545758 K, relative_change = 0.01740590975140556 Iter 3: T = 948.652450501386 K, relative_change = 0.015761786131648043 Iter 20: T = 856.5210443702374 K, relative_change = 0.001340824454901771 Iter 40: T = 850.5813321728194 K, relative_change = 4.2493757546154564e-5 Iter 60: T = 850.3972860145577 K, relative_change = 1.3013383963996326e-6 Iter 80: T = 850.3916537192019 K, relative_change = 3.980985894138925e-8 Converged in 86 iterations to T = 850.3915384231861 K Iter 1: T = 967.3070829693232 K, relative_change = 0.032692917030676796 Iter 2: T = 936.6885734830057 K, relative_change = 0.031653349825919275 Iter 3: T = 908.1131544771248 K, relative_change = 0.030506851278889156 Iter 20: T = 677.4175422594573 K, relative_change = 0.0061330697683919105 Iter 40: T = 654.586883930554 K, relative_change = 0.00022288112595665176 Iter 60: T = 653.8423948336789 K, relative_change = 6.8577039923622325e-6 Iter 80: T = 653.8195728299867 K, relative_change = 2.0981753615166725e-7 Converged in 96 iterations to T = 653.8188968700393 K Iter 1: T = 973.5114939853354 K, relative_change = 0.026488506014664556 Iter 2: T = 949.2160275319459 K, relative_change = 0.024956527584414374 Iter 3: T = 927.0462034324295 K, relative_change = 0.023355931059402854 Iter 20: T = 774.3682644300319 K, relative_change = 0.0028460030530301204 Iter 40: T = 762.757228908751 K, relative_change = 9.398548003301779e-5 Iter 60: T = 762.3919718323588 K, relative_change = 2.8820890931056225e-6 Iter 80: T = 762.3807886043888 K, relative_change = 8.817097040921177e-8 Converged in 91 iterations to T = 762.380487547512 K Iter 1: T = 970.0244000535008 K, relative_change = 0.029975599946499155 Iter 2: T = 942.2065532384693 K, relative_change = 0.028677471219793263 Iter 3: T = 916.5035837907315 K, relative_change = 0.02727954858665957 Iter 20: T = 724.6352192609937 K, relative_change = 0.0042512006495532195 Iter 40: T = 708.1152874052232 K, relative_change = 0.00014610306538452916 Iter 60: T = 707.5878340574847 K, relative_change = 4.4863719879748876e-6 Iter 80: T = 707.5716768846476 K, relative_change = 1.3725607916662175e-7 Converged in 93 iterations to T = 707.5712198674738 K Iter 1: T = 973.4804259831559 K, relative_change = 0.026519574016844082 Iter 2: T = 949.1539283111288 K, relative_change = 0.024989200627695007 Iter 3: T = 926.9533576210875 K, relative_change = 0.023389852823496836 Iter 20: T = 773.9652479944035 K, relative_change = 0.00285557061850998 Iter 40: T = 762.3198555193143 K, relative_change = 9.43265855883889e-5 Iter 60: T = 761.9534815164129 K, relative_change = 2.8925748568294616e-6 Iter 80: T = 761.9422640550989 K, relative_change = 8.849178265323779e-8 Converged in 91 iterations to T = 761.9419620766173 K Iter 1: T = 964.3098824446795 K, relative_change = 0.0356901175553205 Iter 2: T = 930.5444141030623 K, relative_change = 0.03501516364844916 Iter 3: T = 898.6725952585339 K, relative_change = 0.03425072286877265 Iter 20: T = 610.6933516439191 K, relative_change = 0.010241033775171094 Iter 40: T = 574.3274971835215 K, relative_change = 0.0004250831440971143 Iter 60: T = 573.078679367153 K, relative_change = 1.3148520460618425e-5 Iter 80: T = 573.0403240539592 K, relative_change = 4.0235636450845625e-7 Converged in 99 iterations to T = 573.0391576567524 K Iter 1: T = 963.5333598833362 K, relative_change = 0.03646664011666387 Iter 2: T = 928.9425413235738 K, relative_change = 0.03589996984012112 Iter 3: T = 896.1938814553703 K, relative_change = 0.03525369806139209 Iter 20: T = 589.3835759233385 K, relative_change = 0.01208003565575236 Iter 40: T = 546.8544411847148 K, relative_change = 0.0005352329708780664 Iter 60: T = 545.3552658026515 K, relative_change = 1.6603597409026694e-5 Iter 80: T = 545.3091729150349 K, relative_change = 5.081303670692392e-7 Iter 100: T = 545.3077626940441 K, relative_change = 1.5544136473682682e-8 Converged in 101 iterations to T = 545.3077555739194 K Iter 1: T = 969.34618187104 K, relative_change = 0.03065381812896 Iter 2: T = 940.8339352916695 K, relative_change = 0.02941389475980187 Iter 3: T = 914.4240563459998 K, relative_change = 0.02807071254023416 Iter 20: T = 713.6943584429931 K, relative_change = 0.00463174271436386 Iter 40: T = 695.8808526253866 K, relative_change = 0.00016095162019523983 Iter 60: T = 695.3097321448853 K, relative_change = 4.9442376436228745e-6 Iter 80: T = 695.292234897227 K, relative_change = 1.512658207346707e-7 Converged in 94 iterations to T = 695.2917308145651 K Iter 1: T = 966.4908513858966 K, relative_change = 0.03350914861410339 Iter 2: T = 935.0213809555866 K, relative_change = 0.032560546626162604 Iter 3: T = 905.561850245036 K, relative_change = 0.031506799000085996 Iter 20: T = 661.1444029065246 K, relative_change = 0.006949541611913754 Iter 40: T = 635.6199731034395 K, relative_change = 0.00025898493034821895 Iter 60: T = 634.7795913246919 K, relative_change = 7.976077288149065e-6 Iter 80: T = 634.7538209914359 K, relative_change = 2.440422455840552e-7 Converged in 97 iterations to T = 634.7530497025363 K Iter 1: T = 966.4899381344285 K, relative_change = 0.03351006186557149 Iter 2: T = 935.0195130686293 K, relative_change = 0.03256156512766719 Iter 3: T = 905.5589875323687 K, relative_change = 0.03150792590367911 Iter 20: T = 661.1255254657419 K, relative_change = 0.006950548260970537 Iter 40: T = 635.5977825975131 K, relative_change = 0.00025903057122518986 Iter 60: T = 634.7572816053788 K, relative_change = 7.97749242239554e-6 Iter 80: T = 634.7315076053321 K, relative_change = 2.440855531022904e-7 Converged in 97 iterations to T = 634.7307362066775 K Iter 1: T = 976.454468274302 K, relative_change = 0.02354553172569801 Iter 2: T = 955.0702436210435 K, relative_change = 0.021899868706680366 Iter 3: T = 935.7555410236266 K, relative_change = 0.0202233319762819 Iter 20: T = 810.0091433499473 K, relative_change = 0.0020938714943222743 Iter 40: T = 801.156180542379 K, relative_change = 6.772667499377508e-5 Iter 60: T = 800.879809251646 K, relative_change = 2.075437920572315e-6 Iter 80: T = 800.871349562692 K, relative_change = 6.34919800247056e-8 Converged in 89 iterations to T = 800.8711381903922 K Iter 1: T = 965.2189393534553 K, relative_change = 0.03478106064654467 Iter 2: T = 932.4144230543617 K, relative_change = 0.033986606521694905 Iter 3: T = 901.5570258061455 K, relative_change = 0.03309407971954679 Iter 20: T = 633.1637356396523 K, relative_change = 0.00861396133730306 Iter 40: T = 602.1695627940155 K, relative_change = 0.0003385389319908161 Iter 60: T = 601.1278593225836 K, relative_change = 1.044784872231263e-5 Iter 80: T = 601.0958914156791 K, relative_change = 3.1969099854870847e-7 Converged in 98 iterations to T = 601.0949262842594 K Iter 1: T = 964.603746543527 K, relative_change = 0.035396253456473015 Iter 2: T = 931.1495380033613 K, relative_change = 0.034681814848887375 Iter 3: T = 899.6070590476971 K, relative_change = 0.033874772706540636 Iter 20: T = 618.220351955627 K, relative_change = 0.009663375907166809 Iter 40: T = 583.767769472547 K, relative_change = 0.00039327650623198376 Iter 60: T = 582.5938490535949 K, relative_change = 1.215454911705353e-5 Iter 80: T = 582.5578049675025 K, relative_change = 3.7193036156875886e-7 Converged in 99 iterations to T = 582.5567088656464 K Iter 1: T = 964.2755303859041 K, relative_change = 0.0357244696140959 Iter 2: T = 930.4736377265918 K, relative_change = 0.03505418481975241 Iter 3: T = 898.5632309978596 K, relative_change = 0.03429479937411051 Iter 20: T = 609.7955872136649 K, relative_change = 0.010312284605497372 Iter 40: T = 573.1931273159597 K, relative_change = 0.0004290934008022317 Iter 60: T = 571.9349576267376 K, relative_change = 1.3273960428383871e-5 Iter 80: T = 571.8963136152067 K, relative_change = 4.061962616877537e-7 Converged in 99 iterations to T = 571.8951384372117 K Iter 1: T = 980.1821242233469 K, relative_change = 0.01981787577665311 Iter 2: T = 962.4062083467757 K, relative_change = 0.018135319383279076 Iter 3: T = 946.5510658873044 K, relative_change = 0.016474480652724843 Iter 20: T = 849.3602485391223 K, relative_change = 0.0014428235040311076 Iter 40: T = 843.0142709546582 K, relative_change = 4.5851789569320216e-5 Iter 60: T = 842.817439583912 K, relative_change = 1.4042980802500997e-6 Iter 80: T = 842.8114158378141 K, relative_change = 4.295966223544166e-8 Converged in 87 iterations to T = 842.8112818441682 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: 64%|█████████████████████▎ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010518564707839467 Iteration 10: norm(E*F-F'*E) = 5.683709284498492e-5 Iteration 20: norm(E*F-F'*E) = 6.759036373708808e-7 Iteration 30: norm(E*F-F'*E) = 8.675777522537765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.420214670317196e-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: 69%|██████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010261832974668806 Iteration 10: norm(E*F-F'*E) = 4.071631216626311e-5 Iteration 20: norm(E*F-F'*E) = 3.616098301207915e-7 Iteration 30: norm(E*F-F'*E) = 4.341811931869331e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.1781126012830594e-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: 72%|███████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010261832974668806 Iteration 10: norm(E*F-F'*E) = 4.071631216626311e-5 Iteration 20: norm(E*F-F'*E) = 3.616098301207915e-7 Iteration 30: norm(E*F-F'*E) = 4.341811931869331e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.1781126012830594e-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: 69%|██████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010261832974668806 Iteration 10: norm(E*F-F'*E) = 4.071631216626311e-5 Iteration 20: norm(E*F-F'*E) = 3.616098301207915e-7 Iteration 30: norm(E*F-F'*E) = 4.341811931869331e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.1781126012830594e-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: 66%|█████████████████████▋ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010261832974668806 Iteration 10: norm(E*F-F'*E) = 4.071631216626311e-5 Iteration 20: norm(E*F-F'*E) = 3.616098301207915e-7 Iteration 30: norm(E*F-F'*E) = 4.341811931869331e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.1781126012830594e-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: 69%|██████████████████████▊ | ETA: 0:00:00 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010261832974668806 Iteration 10: norm(E*F-F'*E) = 4.071631216626311e-5 Iteration 20: norm(E*F-F'*E) = 3.616098301207915e-7 Iteration 30: norm(E*F-F'*E) = 4.341811931869331e-9 Converged after 35 iterations. norm(E*F-F'*E) = 3.1781126012830594e-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: 67%|██████████████████████ | ETA: 0:00:01 Bin 1 progress: 100%|█████████████████████████████████| Time: 0:00:01 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.010518564707839467 Iteration 10: norm(E*F-F'*E) = 5.683709284498492e-5 Iteration 20: norm(E*F-F'*E) = 6.759036373708808e-7 Iteration 30: norm(E*F-F'*E) = 8.675777522537765e-9 Converged after 35 iterations. norm(E*F-F'*E) = 6.420214670317196e-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 7m12.4s Testing RayTraceHeatTransfer tests passed Testing completed after 446.81s PkgEval succeeded after 513.14s