1 |
/* Copyright (c) 1986 Regents of the University of California */ |
2 |
|
3 |
#ifndef lint |
4 |
static char SCCSid[] = "$SunId$ LBL"; |
5 |
#endif |
6 |
|
7 |
/* |
8 |
* hermite.c - routines for 3D hermite curves. |
9 |
* |
10 |
* 10/29/85 |
11 |
*/ |
12 |
|
13 |
#include <stdio.h> |
14 |
|
15 |
|
16 |
hermite3(hp, p0, p1, r0, r1, t) /* compute point on hermite curve */ |
17 |
double hp[3]; /* returned hermite point */ |
18 |
double p0[3]; /* first endpoint */ |
19 |
double p1[3]; /* second endpoint */ |
20 |
double r0[3]; /* tangent at p0 */ |
21 |
double r1[3]; /* tangent at p1 */ |
22 |
double t; /* position parameter */ |
23 |
{ |
24 |
register int i; |
25 |
double tmh[4]; |
26 |
|
27 |
tmh[0] = (2.0*t - 3.0)*t*t + 1.0; |
28 |
tmh[1] = (-2.0*t + 3.0)*t*t; |
29 |
tmh[2] = ((t - 2.0)*t + 1.0)*t; |
30 |
tmh[3] = (t - 1.0)*t*t; |
31 |
|
32 |
for (i = 0; i < 3; i++) |
33 |
hp[i] = p0[i]*tmh[0] + p1[i]*tmh[1] + |
34 |
r0[i]*tmh[2] + r1[i]*tmh[3]; |
35 |
} |
36 |
|
37 |
|
38 |
htan3(ht, p0, p1, r0, r1, t) /* compute tangent on hermite curve */ |
39 |
double ht[3]; /* returned hermite tangent */ |
40 |
double p0[3]; /* first endpoint */ |
41 |
double p1[3]; /* second endpoint */ |
42 |
double r0[3]; /* tangent at p0 */ |
43 |
double r1[3]; /* tangent at p1 */ |
44 |
double t; /* position parameter */ |
45 |
{ |
46 |
register int i; |
47 |
double tpmh[4]; |
48 |
|
49 |
tpmh[0] = (6.0*t - 6.0)*t; |
50 |
tpmh[1] = (-6.0*t + 6.0)*t; |
51 |
tpmh[2] = (3.0*t - 4.0)*t + 1.0; |
52 |
tpmh[3] = (3.0*t - 2.0)*t; |
53 |
|
54 |
for (i = 0; i < 3; i++) |
55 |
ht[i] = p0[i]*tpmh[0] + p1[i]*tpmh[1] + |
56 |
r0[i]*tpmh[2] + r1[i]*tpmh[3]; |
57 |
} |