Performance optimization on lots of small linear algebra operations

Thank you for the suggestion, perhaps I will integrate my approach wih the CoordinateTransformations. In reality I need non-linear transformations, such as a free-form non-linear (already implemented) where transformations are represented as a field of 3d displacement vectors : https://github.com/vfonov/Minc2.jl/blob/main/src/geo_tools.jl#L147
Perhaps you could recommend some package for that too?