Abstract
Polynomial switching systems such as multivariate splines provide accurate fitting while retaining an algebraic representation and offering arbitrary degrees of smoothness; yet, application of sum-of-squares techniques for local stability analysis is computationally demanding for a large number of subdomains. This communiqué presents an algorithm for region of attraction estimation that is confined to those subdomains actually covered by the estimate, thereby significantly reducing computation time. Correctness of the results is subsequently proven and the run time is approximated in terms of the number of total and covered subdomains. Application to longitudinal aircraft motion concludes the study.