| 18 |
|
set mat=mark |
| 19 |
|
if ($#argv > 2) set mat=$3 |
| 20 |
|
set npts=`wc -l < $pts` |
| 21 |
< |
(head -1 $pts ; cat $pts) | lam - $pts \ |
| 21 |
> |
(head -1 $pts ; cat $pts) | rlam - $pts \ |
| 22 |
|
| rcalc -e '$1=d($1,$2,$3,$4,$5,$6)' -e "cond=$npts+.5-recno" \ |
| 23 |
|
-e 's(x):x*x;d(x0,y0,z0,x1,y1,z1):sqrt(s(x1-x0)+s(y1-y0)+s(z1-z0))' \ |
| 24 |
|
| total -1 -r > /tmp/run$$.dat |
| 25 |
< |
lam /tmp/run$$.dat $pts | tabfunc -i xp yp zp > /tmp/path$$.cal |
| 25 |
> |
rlam /tmp/run$$.dat $pts | tabfunc -i xp yp zp > /tmp/path$$.cal |
| 26 |
|
set tmax=`tail -1 /tmp/run$$.dat` |
| 27 |
|
set nsteps=`ev "floor($tmax/$step)"` |
| 28 |
|
echo $mat > /tmp/tri$$.fmt |
