|
|
1.1 root 1: real u(3), d(3,3), xyz(3,5), work(25)
2: data u / .5, .7, 0 /
3: data d / 0., 0., 0.,
4: $ 1., 0., 1.,
5: $ 0., 1., 10. /
6: data xyz / 0., 0., 0.,
7: $ 1., 0., 0.,
8: $ 0., 1., 10.,
9: $ 0., 2., 20.5,
10: $ 0., 3., 29.4 /
11: call interp(u,d(1,1), d(1,2),d(1,3))
12: write(6,1001)u(3)
13: 1001 format(' interp:',f8.3)
14: c call l2fit(5,xyz,u,work)
15: c write(6,1002)u
16: 1002 format(' l2fit:',3f8.3)
17: end
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.