A GPU-based acceleration for the direct boundary integral equation method (GPU-DBIE) to extrapolate solar coronal magnetic fields is developed, which is about 1000 times faster than the original DBIE. The 3-d coronal magnetic field is reconstructed for NOAA 11158 on 14-Feb-2011 with the SDO/HMI vector magnetogram as bottom boundary condition. The extrapolated results agree well with the projected SDO/AIA, EUV loops and the STEREO EUV sideviews, which verifies the correctness of our GPU-DBIE method. It is also found that the group of bright EUV loops along magnetic neutral lines agree well with current lines, which may have played an important role in the flare process of the active region.