This project is a GPU accelerated "vex" node in Blender's geometry nodes. Right now it works on arbitrary attributes even. I was frustrated with not having Houdini wrangle nodes, so I decided to take matters into my own hands.
Built using C++ to integrate directly into Blender's geometry node system, leveraging GPU acceleration for high-performance attribute manipulation. As you can see with the modifier profiler, it's quite fast.
Wrote a little solver for it to move the particles around and "collide" with the sphere. Substeps can go quite a bit high for better collisions.
There's a few buggy things I want to figure out here and there. This is 100k particles running in realtime (well, 24 fps - I want to see how much I can push this).
I don't think the bottleneck is my node but maybe also the simulation loop within geonodes, but I'm just guessing as I haven't explored that yet.