The algorithm for fast evaluation of the hysteresis loops of uniaxial or "textured" ferroelectric microcrystal or "grains" with long-range interactions is developed. Two types of ceramic microstructures are considered: 1. Random ceramics with complete isotropic distribution function of the crystallographic orientations of grains; and 2. Textured ceramics with anisotropic distribution function of the crystallographic orientations of grains. The qualitative analysis of the hysteretic behavior in terms of the grain distribution function is successfully demonstrated. Comparison of the calculated results with experimental data for Pb(Mg1/3Nb2/3)O3- PbTiO3 ceramics is presented.