hyper_viscosity.hpp

Go to the documentation of this file.

#ifndef HYPER_VISCOSITY_HPP
#define HYPER_VISCOSITY_HPP
 
#include "grid.hpp"
#include "domain_decomposition.hpp"
#include "boundary.hpp"
#include "field.hpp"
#include "plasma.hpp"
 
class HyperViscosity{
    public:
 
    bool rad_on;  
    bool parallel_on;  
 
    double parallel;  
 
    Matrix<DeviceType> rad_mat;
 
    HyperViscosity() : rad_on(false), parallel_on(false) {}
 
    HyperViscosity(NLReader::NamelistReader& nlr, const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid, Plasma& plasma);
 
    void apply_radial(const View<double*, CLayout, DeviceType>& view) const;
    void apply_radial(const View<double*, CLayout, DeviceType>& input, const View<double*, CLayout, DeviceType>& output) const;
#ifdef USE_GPU
    void apply_radial(const View<double*, CLayout, HostType>& view) const;
    void apply_radial(const View<double*, CLayout, HostType>& input, const View<double*, CLayout, HostType>& output) const;
#endif
    void apply_radial(const View<ScalarField*, CLayout, HostType>& input, const View<ScalarField*, CLayout, HostType>& output) const;
 
    void apply(const MagneticField<DeviceType>& magnetic_field, const Grid<DeviceType>& grid, const DomainDecomposition<DeviceType>& pol_decomp, const Boundary& boundary, const View<double*, CLayout, HostType>& field);
 
};
 
#endif