A method to refine crystal structural parameters and charge density distribution using convergent-beam electron diffraction is described, which is based on the least-squares fitting between full dynamical calculations and energy-filtered intensities of two-dimensional CBED patterns. A function of parallel computation using a computer cluster was implemented to extend the applications of the method. Crystal structural parameters of the intermediate phase of hexagonal barium titanate (h-BaTiO3) was refined by the method.