     pol(     g), 
     invert_pol( omega, omega.size(), p.eps_pol), 
     w3d( dg::create::w3d(g)), v3d( dg::create::v3d(g)), one( g.size(), 1.),
-#if def APAR
-    maxwell(A,w3d,v3d),
-    invert_maxwell(omega, omega.size(), p.eps_pol), //chaneg here to eps_maxwell in input file
-#endif
     p(p),
     gp(gp)
 {
+
 }
 template< class container>
 const container& Feltor<container>::compute_vesqr( container& potential)

            dg::blas2::gemv( LaplacianM_perp, x[i], temp);
            dg::blas2::gemv( LaplacianM_perp, temp, y[i]);
            dg::blas1::axpby( -p.nu_perp, y[i], 0., y[i]); // - nu_perp lapl_RZ (lapl_RZ (lnN,U)) //factor MISSING!?!
            
            //additional heavy damping zone
//             dg::blas1::pointwiseDot( lapiris_, temp, temp); //N_i U_i
//             dg::blas1::axpby( -0.1, temp, 1., y[i]); // - nu_perp lapl_RZ (lapl_RZ (lnN,U)) //factor MISSING!?!

//             dg::blas2::gemv( LaplacianM_para, x[i], temp);
//             dg::blas1::axpby(  -p.nu_parallel, temp, 1., y[i]);
        }
        //add parallel resistivity
        std::vector<container>  expy(2);
        expy[0].resize( x[0].size()), expy[1].resize( x[1].size());

    const container& compute_vesqr( container& potential);
    //extrapolates and solves for phi[1], then adds square velocity ( omega)
    const container& compute_psi( container& potential);
    const container& polarisation( const std::vector<container>& y); //solves polarisation equation

    container chi, omega;

    const container binv, curvR, curvZ, gradlnB;
    const container iris, source, damping;
    const container w3d, v3d, one;
    std::vector<container> phi, curvphi, dzphi, expy;
    std::vector<container> dzy, curvy; 

    //matrices and solvers
    Matrix A; 
    dg::DZ<container> dz;
    dg::ArakawaX< dg::DMatrix, container>    arakawa; 
    //dg::Polarisation2dX< thrust::host_vector<value_type> > pol; //note the host vector
    dg::Polarisation< dg::DMatrix, container, container > pol; //note the host vector
    dg::Invert<container> invert_pol;

    const Parameters p;
    const solovev::GeomParameters gp;

