Index: mttroot/mtt/bin/trans/m/mtt_implicit.m ================================================================== --- mttroot/mtt/bin/trans/m/mtt_implicit.m +++ mttroot/mtt/bin/trans/m/mtt_implicit.m @@ -1,21 +1,26 @@ -function [x] = mtt_implicit(x,dx,AA,AAx,dt,Nx,open); # Implicit update - -## ############################################################### -## ## Version control history -## ############################################################### -## ## $Id$ -## ## $Log$ -## ## Revision 1.2 1999/04/20 00:58:22 peterg -## ## Set open-switch states to zero -## ## -## ## Revision 1.1 1999/04/02 06:23:20 peterg -## ## Initial revision -## ## -## ############################################################### - -## Copyright (C) 1999 by P.J. Gawthrop +function [x] = mtt_implicit(x,dx,AA,AAx,dt,Nx,open); + ## x = mtt_implicit(x,dx,AA,AAx,dt,Nx,open); + ## Implicit integration update + + ## ############################################################### + ## ## Version control history + ## ############################################################### + ## ## $Id$ + ## ## $Log$ + ## ## Revision 1.3 1999/04/20 06:14:53 peterg + ## ## Reorder to make equivalent to .p version + ## ## + ## ## Revision 1.2 1999/04/20 00:58:22 peterg + ## ## Set open-switch states to zero + ## ## + ## ## Revision 1.1 1999/04/02 06:23:20 peterg + ## ## Initial revision + ## ## + ## ############################################################### + + ## Copyright (C) 1999 by P.J. Gawthrop I = nozeros(!open.*[1:Nx]'); # Indices of states to update x(I) = AA(I,I)\(AAx(I) + dx(I)*dt); # Implicit update (except open switches); I_open = nozeros(open.*[1:Nx]'); # Indices of open switches x(I_open) = 0.0; # Open switches have zero state