torus torus intersection论文笔记

第一部分是圆环面交线的拓扑分析
圆环面上重要的几个圆：内圆、外圆、大圆、小圆。
1.圆环面参数域中交线的原象：经过变换后T2转移到xoy平面上，T2隐式方程带入T1参数方程中，并通过正切半角公式转换参数，得F(u,v)。
2.交线原象的特征点：左右上下边界点（成对出现）（记为1和2）、转折点（对于某一参数的偏导为0，另一偏导不为0）（记为3）、奇异点（两个偏导为0）（判别式判定）（isolated point，cusp，branch point（记为4））。（这些特征点不是Torus的特性）
给出一个定理：交线原象的显式分段，端点为特征点，分段无交。
特征点的拓扑特征：SSP SEP DSP DEP SSEP DSEP。
3.21个类别：1、2（斜率正负）、3（二阶导数正负）、4、123（正负）、124（同正、同负、异号）、12（正负）、13（正负）、23（正负）、14、24（同正、同负、异号）。
4.用扰动方法处理退化情况：
斜率为0：求解F($u_0$+$\delta$,v) = 0 和 F($u_0$-$\delta$,v) = 0，v2-v1的符号是斜率符号。
二阶导数为0或特征点是cusp：求解F($u_0$+$\delta$,v) = 0 和 F($u_0$-$\delta$,v) = 0，若F($u_0$-$\delta$,v) = 0有两根而F($u_0$+$\delta$,v) = 0无根，二阶导数为负；若F($u_0$+$\delta$,v) = 0有两根而F($u_0$-$\delta$,v) = 0无根，二阶导数为正；若两个方程都有根，选择忽略该点。
两斜率：$F_{uu}=0$则k1为0，求解F($u_0$+$\delta$,v) = 0 和 F($u_0$-$\delta$,v) = 0，v2-v1的符号是k1符号。$F_{vv}=0$则k2为无穷，求解F(u,$v_0$+$\delta$) = 0 和 F(u,$v_0$-$\delta$) = 0，u2-u1的符号是k2符号。
5.构造原象曲线：有一套算法，定义数据结构存储特征点，记录每个特征点沿v轴往下的分支条数。

第二部分是细化曲线：
这是个自适应的方法。
首先估计误差，AB两点的T1对应参数分别为(u1,v1)、(u2,v2)，求解F($\frac{u_1+u_2}{2}$,v) = 0，得到值$v_0$，此点为C，把C到线段AB的距离记为误差。
之后一套算法，大概是不断细分曲线的方法来逼近误差。

第三部分是计算交线里的圆