The LEMan code is a linear global full wave solver that computes perturbed electromagnetic fields in 3D equilibria with nested flux surfaces using the potential formulation. Fourier decomposition in the Boozer angular coordinates is employed and cubic Hermite finite elements are used for the radial discretisation. In stellarator systems, LEMan requires large shared memory capability that can approach 64Gb because the number of Fourier modes required to converge ICRH wave structures exceeds 1300. The matrix inversion in both TERPSICHORE and LEMan is accomplished with the “Burn-At-Both-Ends” (BABE) algorithm which entail an OpenMP parallelisation scheme limited to 2 nodes