Overview
Comment: | Deleted the commented out stuff |
---|---|
Downloads: | Tarball | ZIP archive | SQL archive |
Timelines: | family | ancestors | descendants | both | origin/master | trunk |
Files: | files | file ages | folders |
SHA3-256: |
6fe802cf1f351ae7e72655519c249901 |
User & Date: | gawthrop@users.sourceforge.net on 1998-08-15 09:33:25 |
Other Links: | branch diff | manifest | tags |
Context
1998-08-15
| ||
13:40:47 | Initial revision check-in: b93f24a29b user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
09:33:25 | Deleted the commented out stuff check-in: 6fe802cf1f user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
09:30:05 | Commented out the cariabel iteration stuff check-in: f4415c4d81 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
Changes
Modified mttroot/mtt/bin/trans/p/mtt_sparse.p from [f61ecb354a] to [4534575ce1].
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | PROCEDURE mtt_sparse( b : glnarray; n : integer; VAR x : glnarray; VAR rsq : real); {* ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ ## Revision 1.2 1998/08/15 08:26:30 peterg ## This is variable interations version - now going on to fixed ## iterations ## ## Revision 1.1 1998/08/13 14:58:35 peterg ## Initial revision ## | > > > | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | PROCEDURE mtt_sparse( b : glnarray; n : integer; VAR x : glnarray; VAR rsq : real); {* ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ ## Revision 1.3 1998/08/15 09:30:05 peterg ## Commented out the cariabel iteration stuff ## ## Revision 1.2 1998/08/15 08:26:30 peterg ## This is variable interations version - now going on to fixed ## iterations ## ## Revision 1.1 1998/08/13 14:58:35 peterg ## Initial revision ## |
︙ | ︙ | |||
27 28 29 30 31 32 33 | in the main routine. They must also provide two routines, PROCEDURE asub(x: glnarray; VAR y: glnarray; n: integer); and PROCEDURE atsub(x: glnarray; VAR z: glnarray; n: integer); which calculate A*x and (A transpose)*x *) LABEL 1,99; CONST | < < < < < < < < < < < < < < < < < < < < < < < < < < < < < < < < | 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 | in the main routine. They must also provide two routines, PROCEDURE asub(x: glnarray; VAR y: glnarray; n: integer); and PROCEDURE atsub(x: glnarray; VAR z: glnarray; n: integer); which calculate A*x and (A transpose)*x *) LABEL 1,99; CONST iters = 500; VAR j,iter,irst: integer; rp,gg,gam,eps2,dgg,bsq,anum,aden: real; g,h,xi,xj: glnarray; BEGIN mtt_asub(x,xi,n); rp := 0.0; bsq := 0.0; FOR j := 1 TO n DO BEGIN xi[j] := xi[j]-b[j]; END; mtt_atsub(xi,g,n); FOR j := 1 TO n DO BEGIN g[j] := -g[j]; h[j] := g[j] END; FOR iter := 1 TO iters DO BEGIN mtt_asub(h,xi,n); anum := 0.0; aden := 0.0; FOR j := 1 TO n DO BEGIN anum := anum+g[j]*h[j]; aden := aden+sqr(xi[j]) END; anum := anum/aden; FOR j := 1 TO n DO BEGIN xi[j] := x[j]; x[j] := x[j]+anum*h[j] END; mtt_asub(x,xj,n); FOR j := 1 TO n DO BEGIN xj[j] := xj[j]-b[j]; END; mtt_atsub(xj,xi,n); gg := 0.0; dgg := 0.0; FOR j := 1 TO n DO BEGIN gg := gg+sqr(g[j]); dgg := dgg+(xi[j]+g[j])*xi[j] END; gam := dgg/gg; FOR j := 1 TO n DO BEGIN g[j] := -xi[j]; h[j] := g[j]+gam*h[j] END END; END; |