macroscopic.hpp

Go to the documentation of this file.

#ifndef VIENNASHE_SHE_POSTPROC_MACROSCOPIC_HPP
#define VIENNASHE_SHE_POSTPROC_MACROSCOPIC_HPP
 
/* ============================================================================
   Copyright (c) 2011-2022, Institute for Microelectronics,
                            Institute for Analysis and Scientific Computing,
                            TU Wien.
 
                            -----------------
     ViennaSHE - The Vienna Spherical Harmonics Expansion Boltzmann Solver
                            -----------------
 
                    http://viennashe.sourceforge.net/
 
   License:         MIT (X11), see file LICENSE in the base directory
=============================================================================== */
 
 
 
#include <math.h>
#include <fstream>
#include <iostream>
 
#include "viennagrid/mesh/mesh.hpp"
 
#include "viennashe/forwards.h"
#include "viennashe/physics/constants.hpp"
 
namespace viennashe
{
  template <typename DeviceType,
            typename MacroscopicQuantityAccessor,
            typename ContainerType>
  void write_macroscopic_quantity_to_container(DeviceType const & device,
                                               MacroscopicQuantityAccessor const & quantity,
                                               ContainerType & quantity_container)
  {
    typedef typename DeviceType::mesh_type              MeshType;
 
    typedef typename viennagrid::result_of::const_cell_range<MeshType>::type    CellContainer;
    typedef typename viennagrid::result_of::iterator<CellContainer>::type       CellIterator;
 
    MeshType const & mesh = device.mesh();
 
    CellContainer cells(mesh);
    for (CellIterator cit = cells.begin();
        cit != cells.end();
        ++cit)
    {
      quantity_container.at(static_cast<std::size_t>(cit->id().get())) = quantity(*cit);
    }
  }
 
} //namespace viennashe
 
#endif