Hi, thanks for the advice. The template idea is great. I will implement it in the future.

I was also thinking about the vector arrays. Although I still don't know how to make it work. I'm interested in implementing something like this Julia package: https://juliaarrays.github.io/StructArrays.jl/stable/ but for C++. The problem is that C++ does not support reflection, therefore is a bit complicated. But it's on my to-do list. For now, I guess I can do a specific SOA implementation.

I know the angle calculation is not fast. But I don't want to sacrifice precision.