1 function diffusivity = diffusivity_tensor_zero(glob,U,params)
2 %
function diffusivity = diffusivity_tensor_zero(glob,U,params)
4 %
function computing the diffusivity pointwise evaluation in the point
5 % sequences indicated by global coordinates in the columns of the matrix glob.
7 % generated fields of diffusivity:
8 % epsilon: upper bound on diffusivity value
9 % K : diffusivity tensor, i.e. a sparse square matrix of size
14 if ~isempty(glob) && size(glob,1) < size(glob,2)
15 warning('coordinates in variable glob are given row-wise, but expected them to be column-wise');
23 diffusivity.epsilon = 0;