Artifact 61c5d96355a05faf20f8fe315d14eceb839a70f1e90de59fa2d77e06f3f0bbee:

• Executable file mttroot/mtt/bin/trans/p/mtt_implicit.p — part of check-in [f8fd113135] at 1999-02-16 21:55:45 on branch origin/master — Initial revision (user: gawthrop@users.sourceforge.net, size: 510) [annotate] [blame] [check-ins using]

PROCEDURE mtt_implicit(VAR xnew,x,dx : StateVector;
		       VAR AA	     : StateMatrix;
		       VAR AAx	     : StateVector;
			   dt	     : REAL;
			   Nx	     : INTEGER);

VAR
   i : INTEGER;
   BB  : StateVector;

(*$I $MTTPATH/trans/p/mtt_ludcmp.p *)
(*$I $MTTPATH/trans/p/mtt_lubksb.p *)
(*$I $MTTPATH/trans/p/mtt_solve_lud.p *)

BEGIN{mtt_implicit}
   FOR i := 1 TO  Nx DO {BB is (1-A*dt)*x +dx*dt}
      BB[i] := AAx[i] + dt*dx[i];
   mtt_solve_lud(xnew,AA,BB,Nx); {Solve AAx=BB}
END;{mtt_implicit}