1 /* ---------------------------------------------------------------- */
2 /* berechnet die lauflaenge von punkt zu punkt und summiert */
3 /* diese auf (s). */
4 /* ---------------------------------------------------------------- */
5 #include <extUtil.h>
6
7
8
linelength(double * x,double * y,double * z,int n,double * s)9 void linelength(double *x, double *y, double *z, int n, double *s )
10 {
11 register int i, j;
12 double p0[3], p1[3], v01[3], ds;
13
14 j=0;
15 s[j]=0.;
16 for (i=0; i<n-1; i++)
17 {
18 j++;
19 p0[0]=x[i];
20 p0[1]=y[i];
21 p0[2]=z[i];
22 p1[0]=x[j];
23 p1[1]=y[j];
24 p1[2]=z[j];
25 v_result( p0, p1, v01 );
26 ds=v_betrag( v01 );
27 s[j]=s[i]+ds;
28 }
29 }
30
31