Overview
Comment: | Initial revision |
---|---|
Downloads: | Tarball | ZIP archive | SQL archive |
Timelines: | family | ancestors | descendants | both | origin/master | trunk |
Files: | files | file ages | folders |
SHA3-256: |
451efc65f57637c34c5247992a7ca191 |
User & Date: | gawthrop@users.sourceforge.net on 1998-08-13 08:51:57 |
Other Links: | branch diff | manifest | tags |
Context
1998-08-13
| ||
08:52:33 | Added implicit method check-in: c97c09e8d0 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
08:51:57 | Initial revision check-in: 451efc65f5 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
08:40:40 | Initial revision check-in: b2a1f79f26 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
Changes
Added mttroot/mtt/bin/trans/p/mtt_solve.p version [7fcb2d509d].
> > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 | PROCEDURE mtt_solve(VAR x : StateVector; VAR A : StateMatrix; VAR B : StateVector; n : integer; Small : real); { % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Id$ % % $Log$ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% } VAR i : integer; wmax,wmin : real; w : StateVector ; v : StateMatrix; (*$I svdcmp.p *) (*$I svbksb.p *) BEGIN{mtt_solve} (* decompose matrix A using SVD *) svdcmp(A,n,n,w,v); (* find maximum singular value *) wmax := 0.0; FOR i := 1 to n DO BEGIN IF (w[i] > wmax) THEN wmax := w[i] END; (* define "small" *) wmin := wmax*Small; (* zero the "small" singular values *) FOR i := 1 to n DO BEGIN IF (w[i] < wmin) THEN w[i] := 0.0 END; (* backsubstitute for B *) svbksb(A,w,v,n,n,B,x); END{mtt_solve}; |