view +grid/evalOnStaggered.m @ 1336:0666629aa183 feature/D2_boundary_opt

Add methods for creating grids with different grid point distributions for each coordinate direction, and also supports constructing periodic grids
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 13 May 2022 13:26:16 +0200
parents 8aa0909125a4
children
line wrap: on
line source

function gf = evalOnStaggered(g, f, gridGroup, subGrids)

	default_arg('gridGroup', 1);

	gf = [];

	if isa(g, 'multiblock.Grid')

		nSubGrids = numel(g.grids{1}.gridGroups{gridGroup});
		default_arg('subGrids', 1:nSubGrids );

		for i = 1:g.nBlocks()
			for j = subGrids
				gf = [gf; grid.evalOn(g.grids{i}.gridGroups{gridGroup}{j}, f)];
			end
		end

	else

		nSubGrids = numel(g.gridGroups{gridGroup});
		default_arg('subGrids', 1:nSubGrids );

		for j = subGrids
			gf = [gf; grid.evalOn(g.grids{i}.gridGroups{gridGroup}{j}, f)];
		end

	end

end