PSYDAC is a high-performance Python library for isogeometric analysis [1]. It is an academic open-source project [https://github.com/pyccel/psydac] created by numerical mathematicians for computational plasma physics applications [2, 3, 4]. It can solve a general system of partial differential equations in weak form, which is specified by the user through the domain-specific language provided by the SymPDE library [https://github.com/pyccel/sympde]. It supports finite element exterior calculus (FEEC) with tensor-product spline spaces, and handles multi-patch geometries in a variety of ways. PSYDAC automatically generates Python code for the assembly of user-defined functionals, as well as linear and bilinear forms contained in the weak formulation of the problem. This Python code is then accelerated to C/Fortran speed using Pyccel [5]. The library can perform large parallel computations on distributed-memory supercomputers using MPI and OpenMP. In this contribution, we present strong and weak parallel scaling of PSYDAC on some typical 3D simulations of solid and structural mechanics.