FARGO3D: A new GPU-oriented MHD code

Pablo Benítez-Llambay, Frédéric S. Masset

Research output: Contribution to journalArticlepeer-review

150 Scopus citations


We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on the physics of protoplanetary disks and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite-difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on either graphical processing units (GPUs) or central processing units (CPUs), achieving large speed-up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for CPUs, and have them translated automatically for GPUs.

Original languageEnglish
Article number11
JournalAstrophysical Journal, Supplement Series
Issue number1
StatePublished - 2016
Externally publishedYes


  • Accretion, accretion disks
  • Hydrodynamics
  • Magnetohydrodynamics (MHD)
  • Methods: numerical
  • Protoplanetary disks


Dive into the research topics of 'FARGO3D: A new GPU-oriented MHD code'. Together they form a unique fingerprint.

Cite this