Codes used in the scientific publication: M. Mangeat, S. Chatterjee, J. D. Noh, and H. Rieger, Emergent complex phases in a discrete flocking model with reciprocal and non-reciprocal interactions, submitted (2024). A preprint is available on arXiv.
For each model, a C++ code to compute the numerical simulations of the microscopic model and a C++ code to compute the numerical solutions of the hydrodynamic equations are available in this repository. Some additional Python codes are also available to generate movies of the system dynamics.
Exportations: density snapshots and profiles shown in the different figures of the paper.
Compile: g++ filename.cpp -fopenmp -lgsl -lgslcblas -lm -O3 -s -o filename.out.
Run: ./filename.out -parameter=value.
Generate the movie: python filename.py -parameter=value.
Codes for the reciprocal two-species active Ising model without species flip.
List of parameters for the numerical simulations (TSAIM_omp.cpp): beta, D, v, theta, rho0, mag0, LX, LY, init, tmax, ran, threads (details as comments in the code).
Codes to generate the movies for the numerical simulations: figure_TSAIM_dynamics1d.py and figure_TSAIM_dynamics2d_rect.py.
List of parameters for the numerical solutions of hydrodynamic equations (TSAIM_hydro_omp.cpp): beta, epsilon, rho0, mag0, LX, init, dt, tmax, dx, threads (details as comments in the code).
Code to generate the movies for numerical solutions of hydrodynamic equations: figure_TSAIM_hydro_dynamics1d.py.
Codes for the reciprocal two-species active Ising model with species flip.
List of parameters for the numerical simulations (TSAIM_species_omp.cpp): beta1, beta2, D, v, theta, rho0, mag0, gamma, LX, LY, init, tmax, ran, threads (details as comments in the code).
Codes to generate the movies for the numerical simulations: figure_TSAIM_species_dynamics1d.py and figure_TSAIM_species_dynamics2d_rect.py.
List of parameters for the numerical solutions of hydrodynamic equations (TSAIM_species_hydro_omp.cpp): beta1, beta2, epsilon, rho0, mag0, gamma, LX, init, dt, tmax, dx, threads (details as comments in the code).
Code to generate the movies for numerical solutions of hydrodynamic equations: figure_TSAIM_species_hydro_dynamics1d.py.
Codes for the non-reciprocal two-species active Ising model.
List of parameters for the numerical simulations (NRTSAIM_omp.cpp): beta, D, v, theta, rho0, mag0, JAB, JBA, LX, LY, init, tmax, ran, threads (details as comments in the code).
Codes to generate the movies for the numerical simulations: figure_NRTSAIM_dynamics1d.py and figure_NRTSAIM_dynamics2d_rect.py.
List of parameters for the numerical solutions of hydrodynamic equations (NRTSAIM_hydro_omp.cpp): beta, epsilon, rho0, mag0, JAB, JBA, LX, init, dt, tmax, dx, threads (details as comments in the code).
Code to generate the movies for numerical solutions of hydrodynamic equations: figure_NRTSAIM_hydro_dynamics1d.py.