We have developed a finite-element method to simulate equilibrium properties and transport in two-dimensional (2D) superconducting systems by solving the quasiclassical equations of motion. The code is in C++, highly parallelized, and can now treat systems with non-trivial spin degree of freedom. We want to use this to study the interplay of conventional superconductors and magnets hosting skyrmions. The latter are complicated spatial variations or patterns of the magnetic field in a material. These patterns can be described in terms of topological excitations, and are typically rotating or winding around a "core". Recent experimental work suggests that the combination of these two classes of materials leads to the creation of vortex-antivortex pairs in a superconductor. Our goal is to investigate these pairs in more detail using quasiclassical theory, which allows for e.g. lower temperatures than studied in the available literature, and understand the energetics of vortex-antivortex creation in such materials.