view +multiblock/setAllInterfaceTypes.m @ 1083:a22b23098021 feature/laplace1dUpdate

Divide interface penalty strength by 4, to make it sharp. Stability verified by eigenvalue computations.
author Martin Almquist <malmquist@stanford.edu>
date Tue, 05 Mar 2019 12:41:50 -0800
parents dd95470d4baf
children
line wrap: on
line source

% Create interface configuration with a single type for all interfaces
%    g -- multiblock grid
% type -- type for all interfaces
function intfTypes = setAllInterfaceTypes(g, type)
	intfTypes = cell(g.nBlocks(), g.nBlocks());
	for i = 1:g.nBlocks()
		for j = 1:g.nBlocks()
			if ~isempty(g.connections{i,j})
				intfTypes{i,j} = type;
			end
		end
	end
end