init;

% Set grid
seti.dim = 3;           % two dimensional problem
seti.rCD = 0.2;         % size of computational domain [-rCD,rCD)^dim
seti.nCD = 256;         % number of discretization points for each dimension of CD
seti = setGrid(seti);   % grid on CD is stored in seti.grid; grid on ROI is stored in seti.gridROI
                        % discretization points in each dimension are stored in seti.nROI