1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
/*
**************************************************************************
description
--------------------
copyright : (C) 2002 by Andreas Zehender
email : [email protected]
**************************************************************************
**************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
* *
**************************************************************************/
#include "pmsplinesegment.h"
PMVector PMSplineSegment::point( double t ) const
{
double t2 = t * t;
double t3 = t2 * t;
return PMVector( m_a[0]*t3 + m_b[0]*t2 + m_c[0]*t + m_d[0],
m_a[1]*t3 + m_b[1]*t2 + m_c[1]*t + m_d[1] );
}
PMVector PMSplineSegment::gradient( double t ) const
{
double t2 = t * t;
return PMVector( 3*m_a[0]*t2 + 2*m_b[0]*t + m_c[0],
3*m_a[1]*t2 + 2*m_b[1]*t + m_c[1] );
}
void PMSplineSegment::calculateLinear( const PMVector& p1, const PMVector& p2 )
{
m_a[0] = 0.0;
m_b[0] = 0.0;
m_c[0] = -1.0 * p1[0] + 1.0 * p2[0];
m_d[0] = 1.0 * p1[0];
m_a[1] = 0.0;
m_b[1] = 0.0;
m_c[1] = -1.0 * p1[1] + 1.0 * p2[1];
m_d[1] = 1.0 * p1[1];
}
void PMSplineSegment::calculateQuadratic( const PMVector& p1,
const PMVector& p2,
const PMVector& p3 )
{
m_a[0] = 0.0;
m_b[0] = 0.5 * p1[0] - 1.0 * p2[0] + 0.5 * p3[0];
m_c[0] = -0.5 * p1[0] + 0.5 * p3[0];
m_d[0] = 1.0 * p2[0];
m_a[1] = 0.0;
m_b[1] = 0.5 * p1[1] - 1.0 * p2[1] + 0.5 * p3[1];
m_c[1] = -0.5 * p1[1] + 0.5 * p3[1];
m_d[1] = 1.0 * p2[1];
}
void PMSplineSegment::calculateCubic( const PMVector& p1, const PMVector& p2,
const PMVector& p3, const PMVector& p4 )
{
m_a[0] = -0.5 * p1[0] + 1.5 * p2[0] - 1.5 * p3[0] + 0.5 * p4[0];
m_b[0] = p1[0] - 2.5 * p2[0] + 2.0 * p3[0] - 0.5 * p4[0];
m_c[0] = -0.5 * p1[0] + 0.5 * p3[0];
m_d[0] = p2[0];
m_a[1] = -0.5 * p1[1] + 1.5 * p2[1] - 1.5 * p3[1] + 0.5 * p4[1];
m_b[1] = p1[1] - 2.5 * p2[1] + 2.0 * p3[1] - 0.5 * p4[1];
m_c[1] = -0.5 * p1[1] + 0.5 * p3[1];
m_d[1] = p2[1];
}
void PMSplineSegment::calculateBezier( const PMVector& p1, const PMVector& p2,
const PMVector& p3, const PMVector& p4 )
{
m_a[0] = - p1[0] + 3.0 * p2[0] - 3.0 * p3[0] + p4[0];
m_b[0] = 3.0 * p1[0] - 6.0 * p2[0] + 3.0 * p3[0];
m_c[0] = -3.0 * p1[0] + 3.0 * p2[0];
m_d[0] = p1[0];
m_a[1] = - p1[1] + 3.0 * p2[1] - 3.0 * p3[1] + p4[1];
m_b[1] = 3.0 * p1[1] - 6.0 * p2[1] + 3.0 * p3[1];
m_c[1] = -3.0 * p1[1] + 3.0 * p2[1];
m_d[1] = p1[1];
}
void PMSplineSegment::calculateQuadricBezier( const PMVector& p1,
const PMVector& p2,
const PMVector& p3 )
{
m_a[0] = 0;
m_b[0] = p1[0] - 2.0 * p2[0] + p3[0];
m_c[0] = -2.0 * p1[0] + 2.0 * p2[0];
m_d[0] = p1[0];
m_a[1] = 0;
m_b[1] = p1[1] - 2.0 * p2[1] + p3[1];
m_c[1] = -2.0 * p1[1] + 2.0 * p2[1];
m_d[1] = p1[1];
}
|
