mc_tm_1d_c.m

来自「Matlab下实现的CFD算法」· M 代码 · 共 25 行

function time=MC_TM_1D_C(N,dt,t,pN,Cx)
% MACCORMACK`S TECHNIQUE 
% Conservation Form
% Time Marching

tic;
initial_data(t,N);
for h=1:t
    preCal(h,N,Cx);
    predictor_step(N,h,dt);
    BoundaryConditions(h+1,N,pN);
    preCal(h+1,N,Cx);
    corrector_step(N,h,dt);
    BoundaryConditions(h+1,N,pN);
    Cal_rVT(h);
end
time=toc;