view +multiblock/setAllInterfaceTypes.m @ 1317:34997aced843 feature/poroelastic

Add some interface forcing penalties in ElasticCurvilinearAnisotropic
author Martin Almquist <malmquist@stanford.edu>
date Sun, 26 Jul 2020 20:06:06 -0700
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