A high-performance method has been proposed for calculating the shielding current density in a hightemperature superconducting thin film. After spatially discretized, the initial-boundary-value problem of the shielding current density is reduced to a system of first-order ordinary differential equations that has a strong nonlinearity. However, the system cannot be always solved by means of the Runge-Kutta method even when an adaptive step-size control algorithm is incorporated to the method. In order to suppress an overflow in the algorithm, the following method is proposed: the J-E constitutive relation is modified so that its solution may satisfy the original constitutive relation. A numerical code for analyzing the shielding current density has been developed on the basis of the proposed method and the inductive method has been investigated by use of the code.