Performance of the EXP drag solver as function of the number of fluids (gas and dust fluids) on one A100 GPU obtained with an implementation with a moderate level of optimisation in the Dyablo (left) and the Shamrock (right) codes. We provide the results obtained with the RK1 algorithm for comparison. The tests with Dyablo used static arrays and global memory. The tests with Shamrock used dynamical arrays and either the local (EXP-LOC) or the global (EXP-GLB) memory. For 8-16 dust fluids, the computational cost of the drag step is ∼ 10−100% of the hydro step, which makes it usable for astrophysical applications.