Insight Toolkit (ITK) -- Official Repository. ITK builds on a proven, spatially-oriented architecture for processing, segmentation, and registration of scientific images in two, three, or more dimensions.
Thanks @dzenanz I think this PR might even improve the performance of the FEM module, but I did not have time to do bench-marking. It's a style improvement anyway 😇
Note that the type of these local matrices is
MatrixType = vnl_matrix<Float>
, so such a matrix is just empty, when default-constructed.Following C++ Core Guidelines, February 15, 2024, "Prefer scoped objects, don’t heap-allocate unnecessarily"