causticsEngineering % julia create_mesh.jl cat_posing.jpg
-0.7084825184623098
0.9797413140103279
Loss:
Minimum: -0.7084825184623098
Maximum: 0.9797413140103279
(512, 512)
(512, 512)
(512, 512)
6.591187573040844
1.192422233106962
0.09291080010935611
0.0076503943138856835
0.0006274343139272808
5.145908379627362e-5
Convergence reached at step 3328.0 with max_update of 9.976433165049815e-6
Overall min_t:0.9882087086735224
-0.8705387614429418
1.0434134720282262
Loss:
Minimum: -0.8705387614429418
Maximum: 1.0434134720282262
(512, 512)
(512, 512)
(512, 512)
3.3736189059429162
0.6537082038600845
0.09809572805409107
0.0547325150128704
0.051162139277199
0.05086932351347559
0.05084530815605421
0.05084333856013063
0.050843177003521756
0.05084316379427653
Overall min_t:0.4750608801904638
-1.1313385534689124
1.1098316875462115
Loss:
Minimum: -1.1313385534689124
Maximum: 1.1098316875462115
(512, 512)
(512, 512)
(512, 512)
2.5819520332699244
0.5188760642923861
0.09672104006757247
0.06366445053157743
0.06094363912439287
0.060720498158988744
0.06070219718314486
0.06070069625486606
0.06070057315941933
0.06070056305410539
Overall min_t:0.22785528518448867
0.00029296875
div: -4.312328351034622e-5
Have all the divergences
0.3894483656565549
0.20231833167764193
0.124059012168825
0.0790930799778991
0.052047210613660494
0.036055963315667405
0.03205074059894223
0.03046537914099536
0.029400282745299523
0.028763086266614675
0.02837094564435181
0.02813328101987218
0.027989636398846145
0.027902794036071408
0.027850261189980855
0.02781841869539476
0.027799155643885923
0.02778746079170332
0.027780368819212083
0.027776066796648107
0.027773458094809737
0.027771875952155054
0.02777091663889291
0.02777033489090996
0.02776998210226679
Example heights: -47.027921842063016 and -47.027921842063016 and 2.972078157936984
Specs: 513 513 526338 2048 524288 1052672
We now have 526338 valid nodes
We've filled up 524288 triangles
We've filled up 1048576 triangles
But every once in a while, changing absolutely nothing about the code or source image, I see output like this:
julia create_mesh.jl cat_posing.jpg
-0.7084825184623098
0.9797413140103279
Loss:
Minimum: -0.7084825184623098
Maximum: 0.9797413140103279
(512, 512)
(512, 512)
(512, 512)
Convergence reached at step 1.0 with max_update of 0
Overall min_t:10000
NaN
NaN
Loss:
Minimum: NaN
Maximum: NaN
(512, 512)
(512, 512)
(512, 512)
Convergence reached at step 1.0 with max_update of 0
Overall min_t:10000
NaN
NaN
Loss:
Minimum: NaN
Maximum: NaN
(512, 512)
(512, 512)
(512, 512)
Convergence reached at step 1.0 with max_update of 0
Overall min_t:10000
0.00029296875
div: NaN
Have all the divergences
Convergence reached at step 1.0 with max_update of 0
Example heights: NaN and NaN and NaN
Specs: 513 513 526338 2048 524288 1052672
We now have 526338 valid nodes
We've filled up 524288 triangles
We've filled up 1048576 triangles
So the code is definitely not running deterministically. Is it possibly due to some uninitialized values somewhere?