Stable implicit numerical algorithm of time-dependent Ginzburg-Landau theory coupled with thermal effect for vortex behaviors in hybrid superconductor systems | Synapse