%some grids
params.xnumintervals = 5;
params.ynumintervals = 5;
params.xrange = [0,1];
params.yrange = [0,1];

%grid = rectgrid(params);
%grid = onedgrid(params);
grid = triagrid(params);
grid = set_boundary_types(grid, params);

%discrete functions
dmodel.pdeg = 1;
dmodel.dimrange = 1;

df_info = feminfo(dmodel,grid);

uh = femdiscfunc([],df_info);
uh.dofs = ones(uh.ndofs,1);

uh.dofs = [1:uh.ndofs]';

params = [];
params.no_lines = 0;
params.subsampling_level = dmodel.pdeg - 1;
plot_discfunc(uh,params);
