gradparx2.hpp

Go to the documentation of this file.

#ifndef GRADPARX2_HPP
#define GRADPARX2_HPP

#include "globals.hpp"
#include "grid.hpp"
#include "boundary.hpp"

struct GradParXTmp{
    double sgn;
    Projection<DeviceType> lplane_to_neighbors;
    Projection<DeviceType> rplane_to_neighbors;
    View<int*,CLayout,DeviceType> is_inside;

    GradParXTmp(){}

    GradParXTmp(const Grid<DeviceType>& grid, const Boundary& boundary, double sml_bt_sign)
      : sgn(-(sml_bt_sign<0.0 ? -1.0 : 1.0)), // Reverse sign since negative gradient is desired
        lplane_to_neighbors(grid.ff_lplane_to_neighbors.template mirror<DeviceType>()),
        rplane_to_neighbors(grid.ff_rplane_to_neighbors.template mirror<DeviceType>()),
        is_inside(NoInit("is_inside"),grid.nnode)

    {
        // get is_inside for every node
        is_inside = boundary.get_view_int_not_in_list<DeviceType>(grid.nnode);
    }
};

void negative_gradparx2(const GradParXTmp& gptx, const View<double**,CLayout,DeviceType>& input, const View<double**,CLayout,DeviceType>& output, Order derivative_order);

#endif