Algorithmen zur Parallelisierung von Stoffsimulationssystemen?

Eine einfache Stoffsimulation kann mit dem folgenden Algorithmus durchgeführt werden:

Als Test habe ich es in JavaScript implementiert . Leider skaliert das schlecht. Die Leistung sinkt sehr schnell mit der Anzahl der Partikel und die Grenze ist sehr niedrig. Gibt es eine Möglichkeit, diesen Algorithmus zu parallelisieren? Könnten Sie es als einen einfachen Pseudocode beschreiben?