
ORIGINAL ARTICLE 



Year : 2010  Volume
: 35
 Issue : 2  Page : 104112 

Performance evaluation of an algorithm for fast optimization of beam weights in anatomybased intensity modulated radiotherapy
Vaitheeswaran Ranganathan^{1}, VK Sathiya Narayanan^{2}, Janhavi R Bhangle^{2}, Kamlesh K Gupta^{2}, Sumit Basu^{2}, Vikram Maiya^{2}, Jolly Joseph^{2}, Amit Nirhali^{2}
^{1} Healthcare Sector, Siemens Ltd, 403A, Senapati Bapat Road, Pune 411 016, India ^{2} Department of Radiation Oncology, Ruby Hall Clinic, 40Sasoon Road, Pune411 001, India
Date of Submission  23Sep2009 
Date of Decision  03Dec2009 
Date of Acceptance  09Dec2009 
Date of Web Publication  3Apr2010 
Correspondence Address: Vaitheeswaran Ranganathan Healthcare Sector, Siemens Ltd, 403A, ICC Trade Tower, Tower B, Senapati Bapat Road, Pune 411016 India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.62203
Abstract   
This study aims to evaluate the performance of a new algorithm for optimization of beam weights in anatomybased intensity modulated radiotherapy (IMRT). The algorithm uses a numerical technique called GaussianElimination that derives the optimum beam weights in an exact or noniterative way. The distinct feature of the algorithm is that it takes only fraction of a second to optimize the beam weights, irrespective of the complexity of the given case. The algorithm has been implemented using MATLAB with a Graphical User Interface (GUI) option for convenient specification of dose constraints and penalties to different structures. We have tested the numerical and clinical capabilities of the proposed algorithm in several patient cases in comparison with KonRad^{®} inverse planning system. The comparative analysis shows that the algorithm can generate anatomybased IMRT plans with about 50% reduction in number of MUs and 60% reduction in number of apertures, while producing dose distribution comparable to that of beamletbased IMRT plans. Hence, it is clearly evident from the study that the proposed algorithm can be effectively used for clinical applications.
Keywords: Anatomybased MLC fields, aperturebased IMRT, Gaussian Elimination method, IMRT
How to cite this article: Ranganathan V, Sathiya Narayanan V K, Bhangle JR, Gupta KK, Basu S, Maiya V, Joseph J, Nirhali A. Performance evaluation of an algorithm for fast optimization of beam weights in anatomybased intensity modulated radiotherapy. J Med Phys 2010;35:10412 
How to cite this URL: Ranganathan V, Sathiya Narayanan V K, Bhangle JR, Gupta KK, Basu S, Maiya V, Joseph J, Nirhali A. Performance evaluation of an algorithm for fast optimization of beam weights in anatomybased intensity modulated radiotherapy. J Med Phys [serial online] 2010 [cited 2019 Apr 25];35:10412. Available from: http://www.jmp.org.in/text.asp?2010/35/2/104/62203 
Introduction   
In inverse planning for IMRT, the conventional practice is to define the constraints for tumor and normal tissues based on which the optimal fluence is calculated.^{ [1],[2] } The optimized fluence is delivered using an ordered MLC shapesequence, made by a "leaf sequencer program".^{ [3],[4],[5]} This method is known as fluencebased inverse planning or Beamletbased Inverse Planning (BBIP). The main limitation of this method is that the optimization step does not take into account the delivery constraints imposed by MLC. Therefore the delivery fluence maps may differ from the optimized fluence maps due to physical and mechanical constraints of the MLC. This discrepancy degrades the treatment delivery efficiency and also increases the treatment verification load. It may also produce a leaf sequence that generally requires a large number of apertures and monitor units (MUs). Increase in number of MUs increases the whole body dose due to leakage radiation. Various investigators^{ [6],[7],[8],[9] } and recently Hal^{ [10]} have expressed concern over the the greater probability of cancer induction due to increase in the whole body scatter with IMRT.
Several techniques have been proposed to simplify IMRT plans. An alternative approach to BBIP is Aperturebased Inverse planning (ABIP), where only required numbers of deliverable apertures are generated and their weights get optimized. Also, plans with predefined anatomybased MLC fields (anatomybased ABIP), as proposed by many authors,^{ [11],[12],[13],[14],[15]} could be considered to reduce both the treatment complexity and verification burden. The optimization of the beam weights in anatomybased ABIP was addressed by many investigators by using some form of Monte Carlo technique^{ [16],[17],[18]} simultaneous iterative back projection method^{ [19],[20],[21]} or other algebraic techniques.^{ [22],[23] } The efficiency of such approaches or algorithms depends upon the ability to match the BBIP dose distributions along with significant reduction in no. of beam segments and MUs.
In this communication, we present an algorithm to optimize beam weights in anatomybased IMRT. The distinct feature of the algorithm is that it takes only fraction of a second for optimizing the beam weights, irrespective of the complexity of the given case. The algorithm uses a numerical technique called GaussianElimination that derives the optimum beam weights in an exact or noniterative way. The algorithm has been implemented using MATLAB with a Graphical User Interface (GUI) option for convenient specification of dose constraints and penalties to different structures. We initially tested the performance of the algorithm in different patient cases in terms of convergence, consistency and optimization speed. We also tested the clinical capabilities of the proposed algorithm in several patient cases. To date, we have treated three patients planned using the algorithm in our clinic. In this article, we limited the clinical performance analysis to two patient cases, performed in comparison with beamletbased IMRT plans generated using KonRad^{® } inverse planning system.
Materials and Methods   
Design of anatomybased MLC shapes
We use a simple anatomybased segmentation method ^{[11],[12],[13],[14],[15]} to manually generate the anatomybased MLC fields. In this method, the first field is adjusted to the projection of the target with an appropriate margin using the standard beam'seyeview (BEV) display, blocking the organsatrisk (OAR) present within the BEV of the target results in the subsequent fields per gantry/couch angle combination, depending on the anatomy. In this approach, the intensity of each aperture is allowed to vary continuously. In combination, these apertures produce complex intensity maps.^{ [24]}
Optimization description
To optimize beam weights in IMRT, the objective is to minimize the difference between the prescribed and calculated dose distributions.^{ [25]} The optimization is based on the following quadratic dosebased cost function.
Here F(d (x)) is the cost function, where d(x) is a dose distribution among the voxels of the patient model and x is the vector of parameters to be selected (beam weights) and N is the no. of voxels. D_{ i} = D_{ i} (x) and ^{ p} D_{ i} are the calculated and prescribed maximum dose in voxel i, respectively. Here the array P_{ i} contains the penalties for different structures for violation of dose limits. For instance, a voxel in the target is penalized if it is overdosed (>D_{ max} ) as well as underdosed ( min ). Likewise, a voxel in the normal structure get penalized if it receives a dose more than its D_{}max .
In order to minimize the cost function, a set of optimum beam weights is to be derived. In our approach, instead of attempting to solve the problem in an iterative manner, an exact analytical model is created that describes the whole optimization problem with a set of linear equations solving which will minimize the cost function for the given dose constraints and penalties. A set of simultaneous linear equations that describes the optimization problem is shown in equation 2:
Where the general term D_{ mn} denotes the dose to an m^{ th} voxel from a n^{ th} aperturesegment before optimization and W_{ 1} , W_{ 2}…W_{ n} are the optimized aperture weights that will produce a prescribed dose of d_{ m} at the m^{ th} voxel.
In order to solve Equation 2, a GaussianElimination code is written along with additional modules to smoothening of the beam weights. In linear algebra, Gaussian elimination is a simple yet powerful algorithm that can be used to determine the exact solutions of a system of linear equations. Gaussian Elimination derives the solution through a series of parameter elimination and backsubstitution processes that derives the result without taking any iteration.^{ [26],[27],[28],[29] } Our beam weight optimization scheme can be expressed using Equation 3.
where W is a 1D array of optimum aperture weights, D_{ 0} is a 2D array of dose values at each voxel due to each aperture before optimization, and D_{ p} is a 1D array of dose constraints to each voxel.
Dose calculation
Our algorithm has been implemented in MATLAB. Patient contours are first generated in the CMSXiO^{®} (4.3.1) treatment planning system. In this planning system the dose is calculated using a fast convolution superposition algorithm.^{ [30]} Also, the user defined anatomy fields are generated and doses are estimated by the planning system initially for a set of equal weightings of the fields. Sample points are then manually distributed in each structure based on the anatomy. Then the dose values obtained at those user defined sample points are loaded into MATLAB for optimization. The main advantage of this approach is that there is no need for separate dose calculation software for optimization. Because, once the optimum beam weights are calculated for a given set of dose constraints and penalties, the corresponding dose values at each sample point can be derived just by reversing the operation of the algorithm. Once the doses at each sample point are obtained, a subroutine incorporated with the algorithm immediately calculates the current cost function value. After the final optimization, the resulting beam weights are again exported to CMSXiO planning system and the dose is recalculated in order to confirm that the dose values derived by our algorithm are accurate. The DVHs, dose distributions and MU calculations presented in the paper were calculated using CMS XiO and KonRad planning systems.
Graphical user interface
As part of the work, a Graphical user Interface (GUI) was developed. The GUI was written in MATAB with DICOM compatibility for enabling communication between the planning system and the algorithm. The GUI contains a series of push buttons, checkboxes and popup menus to implement the following functions: (1) import: importing dose values obtained at sample points; (2) Optimization parameters: inputting dose constrains and penalties for different structures; (3) Optimization process: performing optimization and reoptimization; and (4) Output: display of the final beam weights and cost values.
Numerical performance analysis
By far, applying of Gaussian Elimination in optimization of IMRT has not been attempted by any researchers, because the performance of such noniterative algorithms may degrade for large scale optimizations problems (E.g. 10^{ 4}  10^{ 6} simultaneous equations). However, the speed and accuracy of Gaussian Elimination method for relatively small (E.g. 10^{ 1}  10^{ 2 } equations) and medium (E.g. 10^{ 2}  10^{ 4} equations) scale optimization problems can be better than any iterative method. In order to verify the actual numerical capability of the algorithm in the context of beam weight optimization, an analysis was initially performed for a variety of patient cases, in which we evaluated the characteristics of the algorithm in terms of convergence, consistency and speed. To perform the analysis, different data sets (A, B, C, D, E, F, G and H) are generated that belong to different patient cases (H&N 1, H&N 2, Brain 1, Brain 2, Brain 3, Abdomen 1, Abdomen 2 and Abdomen 3 respectively). Each data set is represented using a cost function, whose basic structure resembles the one shown in equation 1. The data sets comprise the user specified treatment parameters (no. of gantry angles, no. of apertures, etc) and optimization parameters (dose constraints and penalties). The number of gantry angles, number of apertures and their shapes were adapted to the anatomy of the given case as described in earlier Section. Then the dose values obtained at the user defined points were loaded into MATLAB for optimization. A Javabased Data Base Management System (DBMS) was separately used for maintaining the data sets with appropriate connectivity to the optimization algorithm. Then the whole problem was described using a set of linear simultaneous equations (as shown in equation 2). Then we attempted to minimize those cost functions by solving the linear equations (optimizing beam weights) for each data set.
Clinical performance analysis
We have studied the clinical performance of the proposed algorithm in eight patient cases (2 H&N, 3 Brain and 3 Abdomen) out of which three cases planed using the algorithm have been treated in our clinic. In this communication, we have presented a detailed account on the clinical performance analysis done for two patient cases (H&N and brain) planned using the algorithm. We used a sampling density (ñ) of 1 point/cm^{ 3} (approximately) for which the cost function has roughly converged in most of the patient cases. For the sake of comparison, beamletbased IMRT plans were also generated for the same cases presented, using the KonRad inverse planning system, which uses a pencil beam algorithm for dose calculation. We used a grid size of 3mm for the dose calculation in KonRad. Also, we used 3 mm CT slices for the study cases. The gantry, collimator and couch angles are kept same for both plans to make the comparison more realistic. We used an intensity level of 7 in the BBIP plans presented in this study. Moreover, the parameters of the objective function are also kept same for both plans, as the clinical objectives for both plans are obviously the same.
Case A
The case is a typical concave H and N lesion of volume 379cc encircling spinal cord and surrounded by parotids on both sides. The difficulty in this case is to cover the target with a prescription dose of 50.4 Gy as homogeneously as possible and at the same time maintain the spinal cord D_{ max} Strictly within 45 Gy. The treatment goals for this plan are summarized in [Table 1]. We developed an anatomybased ABIP plan with seven 6 MV coplanar beams of angles 0^{ o} , 45^{ o} , 100^{ o} , 150^{ o} , 210^{ o} , 260^{ o} , 315^{ o} . The number of apertures/gantry angle was set to be 3 to 5 and hence totally 30 apertures were included in the optimization with appropriate collimator angles.
Case B
The case is a temporal brain lesion of volume 273cc surrounded by critical organs such as brainstem, chiasm, and optic nerves. The prescription for the PTV was 50 Gy along with tight constraints to the nearby OARs. The difficulty in this case is because of the close proximity of the target to optic chiasm, left optic nerve and brainstem. Our intension was to give uniform target coverage and reduce D_{ max} of chiasm and brainstem to less than 50 Gy. The treatment goals for this plan are summarized in [Table 1]. We used a noncoplanar beam arrangement with seven 6 MV beams in ABIP and BBIP plans. We set three to five apertures per gantry angle in ABIP plan that leads a total of 28 apertures.
Results   
Numerical performance analysis
[Table 2] shows the result of the numerical performance analysis in detail, performed in different data sets/patient cases. The reduction in the cost function in all data sets is illustrated in [Figure 1]. The data sets A and E actually represent the Patient Cases A and B taken for the clinical analysis described earlier section. In patient cases A and B, the percent reduction in the cost function is 78 and 80 respectively. Such a reduction in the cost function in both cases resulted in a nearly acceptable dose distribution. It is admitted that the percent reduction in the cost function value does not necessarily mean an equal amount of reduction in the required outcomes. However, in a single criteria optimization, the reduction in the cost function can be used as an approximate indication of the corresponding improvement in the dose distribution.
In our approach, further improvement in dose distribution is derived in successive trials by modifying the penalties of different subobjectives according to the situation as usually performed in most of the optimization systems. Moreover, we also investigated the consistency of the Gaussian Elimination method by repeating the optimizations for multiple times. Our observation is that the variation in the results is less than 0.3% in almost all the data sets presented. The time taken for each optimization was of fractions of CPU seconds in a 1.8 GHz Pentium microprocessor. Therefore, it is evident from this analysis that the characteristics of the algorithm are numerically acceptable for a broad category of cases in the context of beam weight optimization in anatomybased IMRT.
Clinical performance analysis
Case A
[Figure 2] presents the dose distribution on a transversal slice for the two plans, namely ABIP and BBIP. [Figure 3] shows the dosevolume histograms for the structures of interest, allowing comparison between the two plans. The isodose lines of 50.4 Gy (prescription dose), 40 Gy and 30 Gy are shown in the [Figure 2] indicated that the required curvatures in the dose distributions are obtained in ABIP plan as well. [Table 3] summarizes the results in terms of dosevolume indices for both plans. The total MU in ABIP and BBIP plans are 360 and 768 respectively. The total beam segments in ABIP and BBIP are 30 and 84 respectively. [Figure 4] shows the DVHs obtained for PTV and spinal cord at different optimizations trials. In each trial, we changed the penalties of the structures to improve the results. It took 6 optimization trials with different combination of penalties to get an acceptable DVH and dose distribution in ABIP plan.
Case B
[Figure 5] presents the dose distribution on a transversal slice for the two plans and [Figure 6] shows the dosevolume histograms for the structures of interest. The isodose lines of 50 Gy (prescription), 40 Gy, 30 Gy and 20 Gy are shown for both plans in [Figure 7] indicate that ABIP plan is comparable to BBIP plan in terms of target coverage, OAR sparing and spillage control. [Table 3] summarizes the results in terms of dosevolume indices for both plans. The total MU in ABIP and BBIP is 342 and 652 respectively. The total beam segments in ABIP and BBIP are 28 and 66 respectively. [Figure 8] illustrates how the difference in the no. of apertures/gantry angle is reflected in the final cost value. Acceptable dose distribution was obtained only above 3 apertures/gantry angle. Also, it took nine optimization trials in the ABIP plan to get a clinically acceptable dose distribution.
Discussion   
The method described here looks attractive in many aspects. Since the algorithm optimizes beam weights in fractions of a second, the user is given an opportunity to repeat the optimization multiple times with different combinations of penalties by analyzing a plan in multiple clinical, physical and technical viewpoints in a shorter time. The method offers a GUI option for conveniently specifying different dose limits for the tumor and normal tissues as well as assigning individual penalties for reaching the specified objectives. Another advantage in our optimization method is that we did not employ a separate dose calculation algorithm for optimization. We use an existing treatment planning system (CMS XiO) to calculate dose at the user defined points for a set of uniform beam weights and use those dose values as a standard reference. Afterwards, the dose values resulting from each optimization are derived by just performing an inverse Gaussian Elimination operation over those pregenerated reference dose values. This approach reduces a significant amount of time during optimization and allows to directly utilizing the dose calculation accuracy of the existing planning system for the optimization.
It has been clearly demonstrated in the numerical analysis that the numerical performance of the algorithm in terms of convergence, consistency and speed are suitable to perform beam weight optimization in anatomybased IMRT for a wide range of patient cases. The clinical performance analysis shows that the required concavities are produced in the isodose lines of interest so that it conforms to the tumor volume while sparing OARs effectively. Also it is evident from the study that the ABIP plans generated using the algorithm are comparable to BBIP plans in terms of target coverage, OAR sparing and spillage. However, it can be observed from the clinical results that in some parameters ABIP plans are ranking behind BBIP plans. For instance, in Case A, the PTV D_{ mean} is about 4% higher in ABIP plan as compared to BBIP plan and the spinal card D_{ max} in ABIP is 5% higher in ABIP plan than BBIP plan. Likewise, in Case B, the D_{ max} of chiasm, brainstem and Lt. Optic nerve was about 4% higher in ABIP plan as compared to BBIP plan. However, all the parameters were stringently kept within the desired limits in ABIP plans for both the cases. As suggested in recent investigations,^{ [31],[32] } it will be necessary to optimize the beam orientations, couch and wedge angles to reach more competitive plans.
The number of MUs and apertures are significantly lesser in the ABIP plans as compared to their BBIP counter plans. For instance, in Case A, the reduction in the number of MUs and segments in ABIP plan was 53% (BBIP:768/ABIP:360) and 64% (BBIP:84/ABIP:30) respectively as compared to BBIP plan. Likewise, in Case B, the reduction in MU and no of segments in ABIP plan was 48% (BBIP:652/ ABIP:342) and 57% (BBIP:66/ABIP:28) respectively as compared to BBIP plan. The observed reduction in the number of MUs and apertures in the ABIP plans can significantly improve the treatment delivery efficiency and reduce the treatment verification burden.
It is accepted that the optimization methodology mentioned in this study is designed only for quickly finding a local optimum instead of locating a global optimum. But, in the cases studied, we could not observe any particular problem in producing desired results by not getting the global optimum in the optimization. This demonstrates that locating a global optimum for the underlying nonconvex optimization problem may not be necessarily of any significant importance, at least from a clinical point of view.
Conclusion   
An algorithm for fast optimization of beam weights in anatomybased IMRT has been proposed. The numerical as well as clinical performance of the algorithm has been investigated in different patient cases. Our results show that one is able to generate anatomybased IMRT plans using the proposed algorithm that are comparable to BBIP plans in terms of dose distribution and superior to the same in terms of monitor units and number of apertures. It is evident from the study that the proposed algorithm could effectively produce satisfactory plans meeting the clinical objectives, while the verification could remain simple.
Acknowledgment   
The authors would like to thank Health Care Sector, Siemens Ltd, India for providing access to MATLAB software for this research work. The authors would also like to thank Dr. Kumara Swamy for his initial discussions about the research.
References   
1.  Spirou SV, Chui CS. A gradient inverse planning algorithm with dosevolume constraints. Med Phys 1998;25:32133. [PUBMED] 
2.  Wu Q, Mohan R. Algorithm and functionality of an intensity modulated radiotherapy optimization system. Med Phys 2000;27:70111. [PUBMED] 
3.  Spirou SV, Chui CS. Generation of arbitrary intensity profiles by dynamic jaws or multileaf collimators. Med Phys 1994;21:103141. [PUBMED] 
4.  Xia P, Verhey LJ. Multileaf collimator leaf sequencing algorithm for intensity modulated beams with multiple static segments. Med Phys 1998;25:142434. [PUBMED] 
5.  Otto K, Clark BG. Enhancement of IMRT delivery through MLC rotation. Phys Med Biol 2002;47:39974017. [PUBMED] [FULLTEXT] 
6.  Kry SF, Salehpour M, Followil DD, Stovall M, Kuban DA, White RA, et al. The calculated risk of fatal secondary malignancies from intensitymodulated radiation therapy. Int J Radiat Oncol Biol Phys 2005;62:1195203. 
7.  Hall EJ, Wuu CS. Radiationinduced second cancers: The impact of 3DCRT and 1MRT. Int J Radiat Oncol Biol Phys 2003;56:838. [PUBMED] [FULLTEXT] 
8.  Verellen D, Vanhavere F. Risk assessment of radiationinduced malignancies based on wholebody equivalent dose estimates for IMRT treatment in the head and neck region. Radiother Oncol 1999;53:199203. [PUBMED] [FULLTEXT] 
9.  Followill D, Geis P, Boyer A. Estimates of wholebody dose equivalent produced by beam intensity modulated conformal therapy. Int J Radiat Oncol Biol Phys 1997;38:66772. [PUBMED] [FULLTEXT] 
10.  Hall EJ. Intensitymodulated radiation therapy, protons, and the risk of second cancers. Int J Radiat Oncol Biol Phys 2006;65:17. [PUBMED] [FULLTEXT] 
11.  Webb S. Optimization by simulated annealing of threedimensional conformal treatment planning for radiation fields defined by a multileaf collimator. Phys Med Biol 1991;36:12016. [PUBMED] [FULLTEXT] 
12.  Mageras G, Mohan R. Application of fast simulated annealing to optimization of conformal radiation treatments. Med Phys 1993;20:63947. 
13.  Xiao Y, Galvin J, Hossain M, Valicenti R. An optimized forward planning technique for intensity modulated radiation therapy. Med Phys 2000;27:20939. [PUBMED] 
14.  De Gersem W, Claus F, De Wagter C, De Neve W. An anatomy based segmentation tool for intensitymodulated radiation therapy and its application to headandneck cancer. Int J Radiat Oncol Biol Phys 2001;51:84959. [PUBMED] [FULLTEXT] 
15.  Bednarz G, Michalski D, Houser C, Saiful Huq M, Xiao Y, Rani Anne P, et al. The use of mixedinteger programming for inverse treatment planning with predefined field segments. Phys Med Biol 2002;47:223545. 
16.  Webb S. Optimization of conformal radiotherapy dose distributions by simulated annealing. Phys Med Biol 1989;34:134979. [PUBMED] [FULLTEXT] 
17.  Morrill SM, Lane RG, Rosen II. Treatment planning optimization using constrained simulated annealing. Phys Med Biol 1991;36:134161. [PUBMED] [FULLTEXT] 
18.  Mohan R, Mageras GS, Baldwin B, Brewster LJ, Kutcher GJ, Leibel S, et al. Clinically relevant optimization of 3D conformal treatments. Med Phys 1992;19:93344. [PUBMED] 
19.  Bortfeld T, Urkelbach J, Boesecke R, Schlegel W. Methods of image reconstruction from projections applied to conformation therapy. Phys Med Biol 1990;35:142334. 
20.  Powlis WD, Altscheuler MD, Censor Y, Buhe EL. Semiautomated radiotherapy treatment planning with a mathematical model to satisfy treatment goals. Int J Radiat Oncol Biol Phys 1989;16:2716. 
21.  Xing L, Hamilton RJ, Pelizzri C, Chen GT. A three dimensional algorithm for optimizing beam weights and wedge filters. Med Phys 1998;25:185865. 
22.  Starkschall G. A constrained leastsquares optimization method for external beam radiation therapy treatment planning. Med Phys 1984;11:65965. [PUBMED] 
23.  Bahr GK, Kereiaks J, Horwitz H, Finney R, Galvin J, Goode K. The method of linear programming applied to radiation treatment planning. Radiology 1968;91:68693. 
24.  Shepard DM, Earl MA, Li XA, Naqvi S, Yu C. Direct aperture optimization: A turnkey solution for stepandshoot IMRT. Med Phys 2002;29:6. 
25.  Lu R, Radke RJ, Happersett L, Yang J, Chui CS, Yorke E, et al. Reducedorder parameter optimization for simplifying prostate IMRT planning. Phys Med Biol 2007;52:84970. [PUBMED] 
26.  Isaacson E, Keller HB. Analysis of Numerical Methods. New York: Willey; 1966. 
27.  Westlake JR. A handbook of numerical matrix inversion and solution of linear equations. New York: Willey; 1968. 
28.  Chawla MM. A parallel Gaussian elimination method for general linear systems. Int J Comput Math 1992;42:1029265. 
29.  Grimes RG. Solving systems of large dense linear equations. J Supercomput 1988;1:2919. 
30.  Bedford JL. Speed versus accuracy in a fast convolution photon dose calculation for conformal radiotherapy. Phys Med Biol 2002;47:347584. [PUBMED] [FULLTEXT] 
31.  Beaulieu F, Tremblay D, Beaulieu L. Optimizing beam orientations and virtual wedges in multifield conformal radiotherapy. Med Phys 2002;29:1337. 
32.  Beaulieu F, Tremblay D, Roy R, Beaulieu L. Simultaneous optimization of beam orientations, wedge filters and field weights for inverse planning with anatomybased MLC fields. Med Phys 2004;31:154657. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8]
[Table 1], [Table 2], [Table 3]
This article has been cited by  1 
A hybrid algorithm for instant optimization of beam weights in anatomybased intensity modulated radiotherapy: A performance evaluation study 

 Vaitheeswaran, R., Sathiya Narayanan, V.K., Bhangle, J.R., Nirhali, A., Kumar, N., Basu, S., Maiya, V.   Journal of Medical Physics. 2011; 36(2): 8594   [Pubmed]  


