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 = 0;
  plasma.for_all_nonadiabatic_species([&](Species<DeviceType>& species){
      total_ptl += species.get_total_n_ptl();
    }, Plasma::NoDevicePtl);

  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