1 function [DT,N_T,dx,dt,LAMBDA]=HMM_CFL_BL(DX,T,epsilon,M)
3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 %%% Bestimmung der Ortsschrittweiten %%%
5 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7 %
function [DT,N_T,dx,dt,LAMBDA]=HMM_CFL_BL(DX,T,epsilon,M)
9 % Input: DX makroskopische Ortsschrittweite
10 % T Endzeit (fuer Makrorechnung)
11 % epsilon Regularisierungsparameter
12 % Output: DT makroskopische Zeitschrittweite
13 % N_T Anzahl makroskopische Ortsschritte
14 % dx, dt mikroskopische Orts- bzw. Zeitschrittweite
15 % LAMBDA Abkuerzungen fuer Verhaeltnisse
18 % mikroskopische Ortsschrittweiten
21 h = -(2*s*M.*(s-1))./(s.^2+M-2*M*s+M*s.^2).^2;
24 dt = 0.95*min(dx^2/(2*epsilon),dx/w);
26 DT=0.95*(DX/w); % makroskopische Zeitschrittweite mit CFL-Bedingung fuer grobes Gitter