Before revision:  The parallel implementation of the hybrid LBM code in a multi-CPU computer system with MPI is straightforward. For a domain decomposition algorithm on one million lattice nodes, the speedup is linear up to 168 MPI processes. While offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a sui generis scheme that can be easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.
After revision:  The parallel implementation of the hybrid LBM code in a multi-CPU computer system , as well as on GPUs, is straightforward and efficient. Along with offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a sui generis scheme that can be easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.