On-Line Geometric Modeling Notes

Eigenvalues and Eigenvectors
of the Refinement Matrix


Overview

The basic operation in the development of subdivision curves is the refinement procedure. In the cubic case, the points of the refinement can be specified as vertex and edge points and the refinement procedure can be specified by a matrix operation. In this case the refinement matrix A is defined to be
displaymath37
In these notes, we develop the eigenvalues and eigenvectors of the refinement matrix which play an important role in the analysis of subdivision curves. This matrix is also diagonalizable and we use the eigenvalues and eigenvectors to calculate this diagonal decomposition. Finally, the diagonal decomposition allows an easy calculation of the inverse of the matrix.

For a postscript version of these notes look here


The Eigenvalues of the Refinement Matrix

One of the eigenvalues is easy to calculate. Since the rows of A all sum to one, we have
displaymath51
and so 1 is an eigenvalue of A, with corresponding eigenvector tex2html_wrap_inline223. The other eigenvalues can be calculated from the characteristic polynomial and turn out to be tex2html_wrap_inline225 and tex2html_wrap_inline227.

The (right) eigenvectors for the eigenvalues tex2html_wrap_inline229 can be calculated to be
displaymath78
respectively, and in a similar fashion the left eigenvectors turn out to be
eqnarray91


Diagonalization of the Matrix

We will let L be the tex2html_wrap_inline233 matrix whose rows are the left eigenvectors of A,
displaymath121
and R be the tex2html_wrap_inline239 matrix whose columns are the right eigenvectors of A,
displaymath144
noting that tex2html_wrap_inline243 Then since the tex2html_wrap_inline245 matrix A has 3 distinct eigenvalues, it is diagonalizable [1] and can be written as
displaymath153
where tex2html_wrap_inline249 is the diagonal matrix whose diagonal elements are the eigenvalues of A.
displaymath155


The Inverse of the Refinement Matrix

Since the refinement matrix can be written as tex2html_wrap_inline253, it is easy to generate the inverse of A.
eqnarray165
since tex2html_wrap_inline257. The matrix tex2html_wrap_inline259 is just
displaymath174
The inverse of A then works out to be
displaymath180


Summary

The eigenvalues and eigenvectors of the refinement matrix are straightforward to calculate. Since this matrix is well conditioned it can be diagonalized and written in form tex2html_wrap_inline263 which will be very useful when analyzing subdivision curves.


References

1
GOLUB, G. H., AND VAN LOAN, C. F. Matrix Computations, 2nd ed. Johns Hopkins University Press, 1989.

2
HALSTEAD, M., KASS, M., AND DEROSE, T. Efficient, fair interpolation using Catmull-Clark surfaces. In Computer Graphics (SIGGRAPH '93 Proceedings) (Aug. 1993), J. T. Kajiya, Ed., vol. 27, pp. 35-44.


This document maintained by Ken Joy

Comments to the Author

All contents copyright (c) 1996, 1997
Computer Science Department,
University of California, Davis
All rights reserved.



Ken Joy
Tue May 19 07:12:52 PDT 1998