     dg::blas1::axpby( -2., temp0, +1., tempM, dzzf); 
 }
 
-/*
-template<class M, class container>
-void DZ<M,container>::dz2d( const container& f, container& dzf)
-{
-    assert( &f != &dzf);
-    View ghostPV( ghostP.begin(), ghostP.end());
-    View ghostMV( ghostM.begin(), ghostM.end());
-    cView fp( f.cbegin(), f.cend());
-    cView fm( f.cbegin(), f.cend());
-
-    cusp::multiply( plus, fp, ghostPV);
-    cusp::multiply( minus, fm, ghostMV );
-    dg::blas1::axpby( 1., ghostP, -1., ghostM);
-    dg::blas1::pointwiseDivide( ghostM, hz_plane, dzf);
-}
-template< class M, class container >
-void DZ<M,container>::dzz2d( const container& f, container& dzzf)
-{
-    assert( &f != &dzzf);
-
-    View ghostPV( ghostP.begin(), ghostP.end());
-    View ghostMV( ghostM.begin(), ghostM.end());
-    cView fp( f.cbegin() , f.cend());
-    cView fm( f.cbegin() , f.cend());
-
-    cusp::multiply( plus, fp, ghostPV);
-    cusp::multiply( minus, fm, ghostMV );
-    dg::blas1::pointwiseDivide( ghostP, hp_plane, ghostP);
-    dg::blas1::pointwiseDivide( ghostP, hz_plane, ghostP);
-    dg::blas1::pointwiseDivide( f,      hp_plane, dzzf);
-    dg::blas1::pointwiseDivide( dzzf,   hm_plane, dzzf);
-    dg::blas1::pointwiseDivide( ghostM, hm_plane, ghostM);
-    dg::blas1::pointwiseDivide( ghostM, hz_plane, ghostM);
-    dg::blas1::axpby( 2., ghostP, +2., ghostM); //fp+fm
-    dg::blas1::axpby( -2., dzzf, +1., ghostM, dzzf);
-}
-*/
 
 template< class M, class container>
 template< class BinaryOp>

    is a limiter and 0 if there isn't
    * @param field The field to integrate
    * @param grid The grid on which to operate
    * @param eps Desired accuracy of runge kutta
    * @param limit Instance of the limiter class (Default is a limiter everywhere)
    * @param globalbcz Choose NEU or DIR. Defines BC in parallel on box

*
* @ingroup dz
* @tparam Matrix The matrix class of the interpolation matrix
* @tparam container The container-class on which the interpolation matrix operates on (does not need to be dg::HVec)
*/
template< class Matrix = dg::DMatrix, class container=thrust::device_vector<double> >
struct DZ
