In our group, 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. We have successfully used an implementation of this method, in Julia, to study a) phase crystals in d-wave superconductors and b) transport in open superconducting hybrid structures. Our aim with this follow up project is to develop a package for this finite element method that is entirely written in C++, highly parallelised and more user-friendly. Tests on a local workstation indicate a speedup of at least a factor of 10, and we hope to increase this further through better parallelisation support. The aim of this small project is to test the scaling of this new code on a HPC cluster environment and the study of example systems such as superconductor-magnet hybrid structures or d-wave film under the influence of external bias.