Numerical solution to the Ginzburg-Landau (GL) equation becomes infeas
ible as the GL parameter kappa and the number of GL vortices increase
to a physically interesting regime. It is in this regime that we focus
our attention to design a simulated annealing approach for minimizing
the associated renormalized energy to study the GL vortex configurati
ons. This approach applied to a commonly used model produces vortex co
nfigurations which are of interest in recent studies.