derivativesT.h

Go to the documentation of this file.

#pragma once
// must be included after both derivatives.h, derivativesX.h and mpi_derivatives.h are included
// else there is a problem with ADL
// https://stackoverflow.com/questions/40875453/declare-function-after-template-defined
namespace dg{
namespace create{
template<class Topology>
auto dx( const Topology& g, dg::bc bc, dg::direction dir = centered)
{
    static_assert( Topology::ndim() > 0 );
    return derivative( 0, g, bc, dir);
}
template<class Topology>
auto dy( const Topology& g, dg::bc bc, dg::direction dir = centered)
{
    static_assert( Topology::ndim() > 1 );
    return derivative( 1, g, bc, dir);
}
template<class Topology>
auto dz( const Topology& g, dg::bc bc, dg::direction dir = centered)
{
    static_assert( Topology::ndim() > 2 );
    return derivative( 2, g, bc, dir);
}
 
template<class Topology>
auto jumpX( const Topology& g, bc bc)
{
    static_assert( Topology::ndim() > 0 );
    return jump( 0, g, bc);
}
template<class Topology>
auto jumpY( const Topology& g, bc bc)
{
    static_assert( Topology::ndim() > 1 );
    return jump( 1, g, bc);
}
template<class Topology>
auto jumpZ( const Topology& g, bc bc)
{
    static_assert( Topology::ndim() > 2 );
    return jump( 2, g, bc);
}
template<class Topology>
auto derivative( unsigned coord, const Topology& g, direction dir = centered)
{
    return derivative( coord, g, g.bcx(), dir);
}
template<class Topology>
auto jump( unsigned coord, const Topology& g)
{
    return jump( coord, g, g.bcx());
}
template<class Topology>
auto dx( const Topology& g, dg::direction dir = centered)
{
    static_assert( Topology::ndim() > 0 );
    return derivative( 0, g, g.bcx(), dir);
}
template<class Topology>
auto dy( const Topology& g, dg::direction dir = centered)
{
    static_assert( Topology::ndim() > 1 );
    return derivative( 1, g, g.bcy(), dir);
}
template<class Topology>
auto dz( const Topology& g, dg::direction dir = centered)
{
    static_assert( Topology::ndim() > 2 );
    return derivative( 2, g, g.bcz(), dir);
}
 
template<class Topology>
auto jumpX( const Topology& g)
{
    static_assert( Topology::ndim() > 0 );
    return jump( 0, g, g.bcx());
}
template<class Topology>
auto jumpY( const Topology& g)
{
    static_assert( Topology::ndim() > 1 );
    return jump( 1, g, g.bcy());
}
template<class Topology>
auto jumpZ( const Topology& g)
{
    static_assert( Topology::ndim() > 2 );
    return jump( 2, g, g.bcz());
}
} //namespace create
}//namespace dg