@article { author = {Goli, Mehdi}, title = {Comparison of row-wise and column-wise methods in ultra-high degree/order spherical harmonic transforms}, journal = {Iranian Journal of Geophysics}, volume = {15}, number = {1}, pages = {55-65}, year = {2021}, publisher = {Iranian Geophysical Society}, issn = {2008-0336}, eissn = {2783-168X}, doi = {10.30499/ijg.2020.241659.1283}, abstract = {Spherical harmonic transform (SHT) provides a standard tool for spectral analysis of data measured on a spherical manifold particularly in geodesy and geophysics, such as Earth’s gravity/magnetic field, topography, atmosphere, etc. Thanks to recent satellite missions such as SRTM, the spatial resolution of geoscience data decreases, leading to an increase in the degree/order of spherical harmonic. The accurate and fast computation of fully normalized associated Legendre functions (fnALF) is the main part of spherical harmonic computations.In polar region, for high degrees (n>2000), the arithmetic underflow may occur and Extended Range Arithmetic (ERA) can be applied to fix underflow problem. In Belikov method the four terms row-wise reccurrence formula is used for fnALF computation without any numerical issues. Previous studies have only addressed the efficiency of the ERA and Belikov method without mentioning to their performnace in spherical harmonic synthesis (SHS) and analysis (SHA). This study aims at investigating the effect of row-wise and colunm-wise fnALF on SHS/SHA for high harmonic degree/order. For this purpose, we developed parallel FORTRAN95 software that utilizes Belikov and ERA to perform fnALF in the row-wise and column-wise SHT, respectively. For SHA, exact Gauss-Legendre quadrature was used that allows one to determine the level error of SHT. In the first part of numerical results, the computation time and error of Belikov and ERA methods for calculating fnALF are compared in serial and parallel processing. Based on the results of the developed software, Belikov method is approximately 2 times faster than ERA for degrees 2160 to 64800 in both serial and parallel mode. Also for mid-latitudes, the errors for both methods are nearly the same. In polar region the error in ERA grows, while Belikov method has not latitudinal dependencies. To evaluate the performance of fnALF method on error in SHS/SHA, calculations were performed using random spherical harmonic coefficients (Cnm , Snm) between [-1, +1]. Using SHS and SHA in a closed cycle, the computed coefficients (C׳nm , S׳nm) were compared with original coefficients (Cnm , Snm). The RMS and maximum error for different degrees from N=2160 to 21600 were computed. The level error of two methods in terms of RMS is in the same level. Numerical results show that in serial mode, the row-wise SHS/SHA is approximately 1.8 times faster than column-wise method. However, in parallel computing with 10 threads, the speed up factor decreases to about 1.2.  }, keywords = {associated Legendre function,spherical harmonic transforms,Belikov method,Extended Range Arithmetic (ERA)}, title_fa = {مقایسه روش‌های سطری و ستونی تبدیل هماهنگ کروی از درجه/مرتبه خیلی‌بالا}, abstract_fa = {تبدیل هماهنگ کروی، ابزار استانداردی را برای تحلیل/ساخت طیفی داده­های ژئودتیکی و ژئوفیزیکی نظیر میدان ثقل و میدان مغناطیسی زمین، توپوگرافی و جوّ فراهم می­آورد. به لطف مأموریت­های اخیر فضایی نظیر SRTM، تفکیک مکانی داده­های زمین مرجع نظیر توپوگرافی، پیوسته رو به کاهش است که این موضوع، افزایش درجه/مرتبه مدل هماهنگ را در­پی­دارد. محاسبه توابع وابسته لژاندر در درجات بالا (2000n>) با مشکل عددی خطای پاریز روبه‌رو می­شود. برای رفع این مشکل، روش­های مختلف سطری و ستونی محاسبه توابع وابسته لژاندر در مطالعات پیشین پیشنهاد شده است. هدف از مطالعه حاضر، بررسی کارایی این روش­ها در تبدیل تحلیل و ساخت هماهنگ کروی است. الگوریتم سطری و ستونی تحلیل و ساخت هماهنگ کروی در یک نرم­افزار به زبان برنامه­نویسی فرترن و با توان پردازش موازی اجرا شده است. برای محاسبه توابع وابسته لژاندر در حالت سطری از روش بلیکوف (Belikov) و در حالت ستونی از روش حساب برد مبسوط (ERA) استفاده شد. نتایج عددی مبتنی بر نرم­افزار توسعه­داده­شده نشان می­دهد هر دو روش سطری و ستونی در بازه درجه/مرتبه 2160 تا 21600، دقت یکسانی در تبدیل هماهنگ کروی دارند. همچنین سرعت روش سطری نسبت به ستونی در پردازش سری و موازی به‌ترتیب 8/1 و 2/1 برابر سریع­تر است.}, keywords_fa = {توابع وابسته لژاندر,تبدیل هماهنگ کروی,روش بلیکوف,روش ERA}, url = {https://www.ijgeophysics.ir/article_120827.html}, eprint = {https://www.ijgeophysics.ir/article_120827_2978838de12f888be052410eb367e1c3.pdf} }