Home / Store / Journals / Electronic Articles / Nuclear Science and Engineering / Volume 165 / Number 3 / Pages 283-291
Thomas E. Booth, James E. Gubernatis
Nuclear Science and Engineering / Volume 165 / Number 3 / Pages 283-291
Format:electronic copy (download)
Recently, we proposed a modified power iteration method that simultaneously determines the dominant and subdominant eigenvalues and eigenfunctions of a matrix or a continuous operator. One advantage of this method is the convergence rate to the dominant eigenfunction being [vertical bar]k3[vertical bar]/k1 instead of [vertical bar]k2[vertical bar]/k1, a potentially significant acceleration. One challenge for a Monte Carlo implementation of this method is that the second eigenfunction is represented by particles of both positive and negative weights that somehow must sum (cancel) to estimate the second eigenfunction faithfully. Our previous Monte Carlo work has demonstrated the improved convergence rate by using a point flux estimator method and a binning method to effect this cancellation. This paper presents an exact method that cancels over a region instead of at points or in small bins and has the potential of being significantly more efficient than the other two.
Your cart is empty.
Home|Invoice Payment|Nuclear Links