anchor.diff:  ----
   g++ -std=c++14 -fopenmp -mavx -mfma -DTHRUST_DEVICE_SYSTEM=THRUST_DEVICE_SYSTEM_OMP -Ipath/to/feltor/inc -Ipath/to/include test.cpp -o test
   export OMP_NUM_THREADS=4
  +./test
  +----
  +If you do not want any parallelization, you can use a single thread version
  +[source,sh]
  +----
  +g++ -std=c++14 -mavx -mfma -DTHRUST_DEVICE_SYSTEM=THRUST_DEVICE_SYSTEM_CPU -Ipath/to/feltor/inc -Ipath/to/include test.cpp -o test
   ./test
   ----
   
  
positive.code: ----
  cd path/to/feltor/inc/dg
  
  make blas_b device=omp #(for an OpenMP version)
  #or
  make blas_b device=gpu #(if you have a gpu and nvcc )
  ----
  
  Run the code with
  
negative.code: ----
  cd path/to/feltor/inc/dg
  
  make blas_mpib device=cpu  # (for MPI+CPU)
  # or
  make blas_mpib device=omp  # (for MPI+OpenMP)
  # or
  make blas_mpib device=gpu # (for MPI+GPU, requires CUDA-aware MPI installation)
  