Improving the numerical technique for computing the accumulated distribution of a quadratic form in normal variables

This paper is concerned with the technique of numerically evaluating the cumulative distribution function of a quadratic form in normal variables. The efficiency of two new truncation bounds and all existing truncation bounds are investigated. We also find that the suggestion in the literature for further splitting truncation errors might reduce computational efficiency, and the optimum splitting rate could be different in different situations. A practical solution is provided. The paper also discusses a modified secant algorithm for finding the critical value of the distribution at any given significance level.