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"

class HyperViscosity{
    public:

    bool rad_on;  
    bool parallel_on;  

    double rad;  
    double parallel;  

    HyperViscosity() : rad_on(false), parallel_on(false) {}

    HyperViscosity(NLReader::NamelistReader& nlr);

    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