Electron beam dose distributions were calculated on a three-dimensional grid using three pencil-beam algorithms, each taking into account irregularities in field shape. The algorithms differ in that patient anatomy in either one, two, or three dimensions is used in the calculation of dose to a point. Algorithms were optimized for speed by such techniques as precalculation and storage of several quantities, reordering of pencil-beam and grid-point loops, selection of cut-off values for some calculated quantities, and invoking error function symmetries. Execution times for optimized versions of each of the algorithms as implemented on a three-dimensional treatment planning system were comparable for both the one- and two-dimensional heterogeneity correction requires an additional calculational loop over fan lines. Execution times for the three-dimensional heterogeneity correction were approximately a factor of four longer than those for the two-dimensional correction. For certain geometries, three-dimensional heterogeneity corrections were necessary to calculate dose distributions accurately, in spite of the additional cost in calculation times.