function v = inflow_masked(v, p_inflow, vessel, mask)
% Inflow condizionato a regioni definite da 'mask'
% mask vale 1 dove l'inflow Ã¨ permesso, 0 altrove
% vessel vale 1 nei punti dei vasi

a = v(:);
b = vessel(:) .* mask(:);  % solo dove entrambe sono 1

a(b == 1) = a(b == 1) + (rand(sum(b), 1) < p_inflow);

v(:) = a;
