Mode solver based on Gegenbauer polynomial expansion for cylindrical structures with arbitrary cross sections.