Index: mttroot/mtt/bin/trans/p/mtt_implicit.p ================================================================== --- mttroot/mtt/bin/trans/p/mtt_implicit.p +++ mttroot/mtt/bin/trans/p/mtt_implicit.p @@ -8,13 +8,11 @@ VAR i,ii,j,jj : INTEGER; BB,xsub : StateVector; AAsub : StateMatrix; -(*$I $MTTPATH/trans/p/mtt_ludcmp.p *) -(*$I $MTTPATH/trans/p/mtt_lubksb.p *) -(*$I $MTTPATH/trans/p/mtt_solve_lud.p *) +(*$I $MTTPATH/trans/p/mtt_solve.p *) BEGIN{mtt_implicit} ii := 0; FOR i := 1 TO Nx DO {BB is (1-A*dt)*x +dx*dt} BEGIN @@ -30,11 +28,11 @@ END; END; END; Nx := ii; {New size} - mtt_solve_lud(xsub,AAsub,BB,Nx); {Solve AAx=BB} + mtt_solve(xsub,AAsub,BB,Nx); {Solve AAx=BB} ii := 0; FOR i := 1 TO Nx DO IF open[i]<0.5 THEN BEGIN