simpleFOM.hpp

Go to the documentation of this file.

#ifndef SIMPLEFOM_HPP
#define SIMPLEFOM_HPP

#include <iostream>
#include "plasma.hpp"
#include "globals.hpp"
#include "timer_macro.hpp"

extern "C" int get_sml_mype();
extern "C" int get_sml_totalpe();

// Simple performance figure-of-merit (FOM)

template <class DeviceType>
void simpleFOM(Plasma& plasma, DomainDecomposition<DeviceType>& pol_decomp, Simulation<DeviceType>& sml)
{
  double loop_time = 0;
  GPTLget_wallclock("MAIN_LOOP", -1, &loop_time);
  double total_ptl = plasma.all_species[0].n_ptl * pol_decomp.mpi.nranks; // 0 = electrons; assume all species have same particle count
  if (is_rank_zero()) {
    std::cout << "SimpleFOM = 1/(loop_time/(total_ptl * sml.mstep)) =  1/(" << loop_time << "/(" << total_ptl << " * " << sml.mstep << ")) = "
              << 1.0/(loop_time/(total_ptl * sml.mstep)) << std::endl;
  }
  
  return;
}

#endif