In this paper, we report our efforts to model split-ring resonators (SRRs) and their large arrays accurately and efficiently in a sophisticated simulation environment based on recent advances in the computational electromagnetics. The resulting linear system obtained from the simultaneous discretization of the geometry and Maxwell's equations is solved iteratively with the multilevel fast multipole algorithm. As an example, we present an array of 125 SRRs showing a negative effective permeability about 92 GHz.