     container evaluate( BinaryOp f, UnaryOp g, unsigned p0, unsigned rounds);
     template< class BinaryOp, class UnaryOp>
     container evaluateAvg( BinaryOp f, UnaryOp g, unsigned p0, unsigned rounds);
+
     void eins(const Matrix& interp, const container& n, container& npe);
     void einsPlus( const container& n, container& npe);
     void einsMinus( const container& n, container& nme);
     void einsPlusT( const container& n, container& npe);
     void einsMinusT( const container& n, container& nme);
-    void centeredT( const container& f, container& dzf);
-    void forwardT( const container& f, container& dzf);
-    void backwardT( const container& f, container& dzf);
-
-    void symv( const container& f, container& dzTdzf);
+
+
     /**
      * @brief Returns the weights used to make the matrix symmetric 
      *

    */
    template <class Field, class Limiter>
    DZ(Field field, const dg::Grid3d<double>& grid, double hz, double eps = 1e-4, Limiter limit = DefaultLimiter(), dg::bc globalbcz = dg::DIR);
    /**
    * @brief Apply the derivative on a 3d vector
    *
    * @param f The vector to derive
    * @param dzf contains result on output (write only)
    */
    void operator()( const container& f, container& dzf);
    void forward( const container& f, container& dzf);
    void backward( const container& f, container& dzf);

    //void dz2d( const container& f, container& dzf);
    //void dzz2d( const container& f, container& dzzf);
    /**
    * @brief Set boundary conditions in the limiter region
    *

    is a limiter and 0 if there isn't
    * @param field The field to integrate
    * @param grid The grid on which to operate
    * @param hz 
    * @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
