Coordinate parametrization and spectral method optimisation for Beltrami field solver in stellarator geometry