 #include "functions.h"
 #include "interpolation.cuh"
 #include "draw/host_window.h"
+
 
 int main()
 {

        std::cout << "Norm diff "<<sqrt( normdiff2)<<"\n"; //=sqrt(gradlnB - dz(ln(B)))
        double reldiff2 =sqrt( dg::blas2::dot( diff2, w3d, diff2)/normsol2 );
        std::cout << "Relative Difference Is "<<reldiff2 <<"\n";
        dzerrfile1 << reldiff  << std::endl;
        dzerrfile2 << reldiff2 << std::endl;
    }
    dzerrfile1.close();
    dzerrfile2.close();
    }

    return 0;

#define DG_DEBUG

#include <iostream>
#include <iomanip>
#include <sstream>
