Using a bihomogeneous resultant to find the singularities of rational space curves

We provide a new technique to detect the singularities of rational space curves. Given a rational parametrization of a space curve, we first compute a μ-basis for the parametrization. From this μ-basis we generate three planar algebraic curves of different bidegrees whose intersection points correspond to the parameters of the singularities. To find these intersection points, we construct a new sparse resultant matrix for these three bivariate polynomials. We then compute the parameter values corresponding to the singularities by applying Gaussian elimination to this resultant matrix. Let ν_Q denote the multiplicity of the singular point Q, and let n be the degree of the curve. We find that when ∑ν_Q≤2n-3, the last nonzero row after Gaussian elimination represents a univariate polynomial whose roots are exactly the parameter values of the singularities with the correct multiplicity. Otherwise the last two nonzero rows represent two bivariate polynomials whose common roots provide the parameter values of the singularities. We also show that if R is this resultant matrix, then size(R)-rank(R) gives the total multiplicity ∑ν_Q(ν_Q-1) of all the singular points including the infinitely near singular points of a rational space curve and we provide bounds on the expression ∑ν_Q(ν_Q-1) for the total multiplicity of all the singular points of a rational space curve. To verify our results, we present several examples to illustrate our methods.