numlib/src/eig.pas

99 var            i, sr, nsr : ArbInt;
107   sr:=sizeof(ArbFloat); nsr:=n*sr;
108   getmem(d, nsr); getmem(cd, nsr); getmem(dh, nsr); getmem(cdh, nsr);
109   getmem(u, n*nsr);
112   move(d^[1], dh^[1], nsr); move(cd^[1], cdh^[1], nsr);
115   freemem(d, nsr); freemem(cd, nsr); freemem(dh, nsr); freemem(cdh, nsr);
116   freemem(u, n*nsr);
121 var          i, sr, nsr : ArbInt;
129   sr:=sizeof(ArbFloat); nsr:=n*sr;
130   getmem(d, nsr); getmem(cd, nsr); getmem(u, n*nsr);
134   freemem(d, nsr); freemem(cd, nsr); freemem(u, n*nsr);
139 var   nsr : ArbInt;
146   nsr:=n*sizeof(ArbFloat);
147   getmem(d, nsr); getmem(cd, nsr);
150   freemem(d, nsr); freemem(cd, nsr)
155 var      i, sr, nsr : ArbInt;
163   sr:=sizeof(ArbFloat); nsr:=n*sr;
164   getmem(d, nsr); getmem(cd, nsr); getmem(u, n*nsr);
169   freemem(d, nsr); freemem(cd, nsr); freemem(u, n*nsr) { toegevoegd 3 apr 92 }
174 var               sr, nsr : ArbInt;
181   sr:=sizeof(ArbFloat); nsr:=n*sr;
182   pd:=@d; pcd:=@cd; getmem(dh, nsr); getmem(cdh, nsr);
183   move(pd^[1], dh^[1], nsr); move(pcd^[1], cdh^[2], (n-1)*sr);
187       move(pd^[1], dh^[1], nsr); move(pcd^[1], cdh^[2], (n-1)*sr);
190   freemem(dh, nsr); freemem(cdh, nsr);
195 var               sr, nsr : ArbInt;
203   term:=1; sr:=sizeof(ArbFloat); nsr:=n*sr; getmem(cdh, nsr);
206   freemem(cdh, nsr)
211 var             i, sr, nsr : ArbInt;
219   sr:=sizeof(ArbFloat); nsr:=n*sr;
220   getmem(cdh, nsr);
222   for i:=1 to n do fillchar(px^[(i-1)*rwidth+1], nsr, 0);
225   freemem(cdh, nsr);
248       uwidth, sr, nsr : ArbInt;
254   sr:=sizeof(ArbFloat); nsr:=n*sr; uwidth:=l+1;
255   getmem(u, uwidth*nsr); getmem(d, nsr); getmem(cd, nsr);
259   freemem(u, uwidth*nsr); freemem(d, nsr); freemem(cd, nsr);
265             sr, nsr, uwidth : ArbInt;
271   sr:=sizeof(ArbFloat); nsr:=n*sr; uwidth:=l+1;
272   getmem(u, uwidth*nsr); getmem(d, nsr); getmem(cd, nsr);
276   freemem(u, uwidth*nsr); freemem(d, nsr); freemem(cd, nsr)
282            sr, nsr, uwidth : ArbInt;
288   sr:=sizeof(ArbFloat); nsr:=n*sr; uwidth:=l+1;
289   getmem(u, uwidth*nsr); getmem(d, nsr); getmem(cd, nsr);
293   freemem(u, uwidth*nsr); freemem(d, nsr); freemem(cd, nsr)