!SOLVER¶
Həlledicinin idarə edilməsi
Məcburi idarəetmə verilənləridir.
Parametrlər¶
METHOD = Həll üsulu (CG, BiCGSTAB, GMRES, GPBiCG, CR, DIRECT, DIRECTmkl, MUMPS)
DIRECT : Kontakt analizindən başqa hallarda birbaşa üsul (yalnız ardıcıl emal) (hazırda istifadə edilə bilməz)
DIRECTmkl : Intel MKL ilə birbaşa üsul
MUMPS : MUMPS ilə birbaşa üsul
Birbaşa üsul seçildikdə verilənlər sətirləri nəzərə alınmır.
3 sərbəstlik dərəcəli məsələlər üçün iterativ üsullarda OpenMP ilə axın paralelliyi istifadə edilə bilər
PRECOND = İterativ üsul üçün önşərtləndirmə üsulu (1､2､3､5､10､11､12､22)
1, 2 : (Blok) SSOR (yalnız 3 sərbəstlik dərəcəli məsələlər üçün çoxrəngli sıralama ilə)
3 : (Blok) diaqonal miqyaslama
5 : ML çoxsəviyyəli önşərtləndirmə paketi ilə AMG
10 : Blok ILU(0)
11 : Blok ILU(1)
12 : Blok ILU(2)
22 : HEC-MW-yə daxil edilmiş SA-AMG (hamarlanmış aqreqasiyalı AMG) ilə AMG
10､11､12 yalnız 3 sərbəstlik dərəcəli məsələlərdə istifadə oluna bilər
22 yalnız LAPACK aktiv edilmiş yığımda istifadə oluna bilər
OpenMP axın paralelliyi zamanı SSOR, diaqonal miqyaslama, ML və ya SA-AMG tövsiyə olunur
ITERLOG = İterativ həlledicinin yığılma tarixçəsinin çıxarılıb-çıxarılmaması (YES/NO)(standart: NO)
TIMELOG = Həlledicinin hesablama vaxtının çıxarılıb-çıxarılmaması (YES/NO/VERBOSE)(standart: NO)
LOGLEVEL = Önşərtləndiricinin diaqnostik məlumat çıxış səviyyəsi (tam ədəd)(standart: çıxış yoxdur)
ITERLOG və TIMELOG-dan asılı olmayaraq göstərilir
1 və yuxarı qiymətdə AMG tipli önşərtləndiricilərin (PRECOND=5､22) iyerarxiyasının qurulması və digər məlumatlar çıxarılır
2 və yuxarı qiymətdə daha ətraflı məlumat (məsələn, yaddaş istifadəsi) çıxarılır
USEJAD = Vektor prosessorları üçün sıralamadan istifadə edilib-edilməməsi (YES/NO)(standart: NO)
Yalnız 3 sərbəstlik dərəcəli məsələdə iterativ üsul istifadə edildikdə etibarlıdır
SCALING = Matrisin diaqonal elementlərini 1 edən miqyaslamanın tətbiq edilib-edilməməsi (YES/NO)(standart: NO)
Yalnız 3 sərbəstlik dərəcəli məsələdə iterativ üsul istifadə edildikdə etibarlıdır
DUMPTYPE = Matris dump formatı (NONE, MM, CSR, BSR) (əsasən sazlama üçün)
NONE : dump etmə (standart)
MM : Matrix Market formatı
CSR : Compressed Sparse Row (CSR) formatı
BSR : Blocked CSR formatı
DUMPEXIT = Matris dump-dan dərhal sonra proqramın dayandırılması (YES/NO)(standart: NO)
MPCMETHOD = Çoxnöqtəli məhdudiyyət şərtlərinin emal üsulu (1､3)
1: Cərimə üsulu (birbaşa üsul istifadə edildikdə standart)
3: Aşkar sərbəstlik dərəcəsi aradan qaldırma üsulu (iterativ üsul istifadə edildikdə standart)
Əvvəllər seçilə bilən MPC-CG üsulu silinib; 2 göstərildikdə xəbərdarlıq verilir və müvafiq susmaya görə üsul istifadə olunur.
ESTCOND = Şərt ədədinin qiymətləndirilməsi tezliyi (sınaq xarakterli)
Göstərilən hər iterasiyada və iterasiyalar bitdikdə şərt ədədi qiymətləndirilir
0 olduqda qiymətləndirmə aparılmır
RECOMPUTE_RESIDUAL = Qalıq vektorunun yenidən hesablanma tezliyi
Yalnız METHOD üçün CG, BiCGSTAB, GPBiCG və ya PipeCG göstərildikdə qüvvədədir
İterativ üsullar qalıq vektorunu rekurent düsturla yenilədiyinə görə yuvarlaqlaşdırma xətaları toplanır
Göstərilən hər iterasiya intervalında qalıq vektoru {b}-[A]{x} ifadəsindən yenidən hesablanır
Dəyər 0 olduqda hər üsul üçün standart dəyər istifadə olunur (CG: 50, BiCGSTAB: 100, GPBiCG: 20, PipeCG: 50)
Mənfi dəyər göstərildikdə dövri yenidən hesablama aparılmır
Yaxınsama meyarı ödənildikdə yoxlama üçün yenidən hesablama bu parametrdən asılı olmayaraq həmişə aparılır
Periodu həddən artıq qısaltmaq matris-vektor hasillərinin sayını artırmaqla yanaşı
yaxınsamanın özünü pisləşdirə və ya yaxınsamanın baş verməməsinə səbəb ola bilər
Standart dəyər dəyişdirilərsə, yığılma tarixçəsini və son nisbi qalığı yoxlayın
METHOD2 = İkinci həll üsulu (BiCGSTAB, GMRES, GPBiCG) (sınaq xarakterli)
Yalnız METHOD üçün CG göstərildikdə etibarlıdır
CG üsulu dağıldıqda avtomatik olaraq keçid edilir və həll davam etdirilir
Digər parametrlər və verilənlər sətirlərindəki eyni məlumatlar istifadə olunur
CONTACT_ELIM = Kontakt analizində sərbəstlik dərəcələrinin aradan qaldırılıb-qaldırılmamasını göstərir (0,1)
0: Sərbəstlik dərəcələrinin aradan qaldırılması yalnız iterativ üsul istifadə edildikdə aparılır (standart)
1: Həmişə (birbaşa üsul istifadə edildikdə də) sərbəstlik dərəcələrinin aradan qaldırılması aparılır
** 2-ci sətirdən etibarən **
| Dəyişən adı | Atribut | Məzmun |
|---|---|---|
| NITER | I | İterasiya sayı (standart: 100) |
| iterPREmax | I | Additive Schwarz ilə önşərtləndirmənin təkrar sayı (standart: 1) (Tövsiyə olunan qiymət 1-dir (paralel hesablamada bəzən 2 də səmərəli ola bilər)) |
| NREST | I | Krılov altfəzalarının sayı (standart: 10) (Yalnız həll üsulu kimi GMRES seçildikdə etibarlıdır) |
| NCOLOR_IN | I | Çoxrəngli sıralamada hədəf rəng sayı (standart: 10) (Yalnız OpenMP axınlarının sayı 2 və ya daha çox olduqda etibarlıdır) |
| RECYCLEPRE | I | Önşərtləndiricinin qurulma məlumatının təkrar istifadə sayı (standart: 3) (Yalnız qeyri-xətti analizdə etibarlıdır) |
| Dəyişən adı | Atribut | Məzmun |
|---|---|---|
| RESID | R | Kəsmə xətası (standart qiymət: 1.0e-8) |
| SIGMA_DIAG | R | Önşərtləndirmə matrisi hesablanarkən diaqonal elementlərə tətbiq olunan əmsal (standart qiymət: 1.0) (ILU önşərtləndirməsində sıfıra bölmə və ya dağılma baş verərsə, 1.0-dan böyük qiymət verməklə həll mümkün ola bilər) |
| SIGMA | R | İstifadə edilmir (standart qiymət: 0.0) |
PRECOND=5 olduqda (buraxıla bilər)¶
PRECOND üçün 5 və ya 22-dən başqa qiymət göstərildikdə 4-cü sətir nəzərə alınmır.
(4-cü sətir) ML_CoarseSolver, ML_Smoother, ML_MGCycle, ML_MaxLevels, ML_CoarseningScheme, ML_NumSweep
| Dəyişən adı | Atribut | Məzmun |
|---|---|---|
| ML_CoarseSolver | I | ML-in kobud tor həlledicisi (1: hamarlaşdırıcı ilə əvəz et, 2: KLU (ardıcıl birbaşa üsul), 3: MUMPS (paralel birbaşa üsul)) (standart qiymət: 1) (Tövsiyə olunan qiymət çətin məsələlər üçün 3 və ya 2, adi məsələlər üçün 1-dir) |
| ML_Smoother | I | ML hamarlaşdırıcısı (1: Chebyshev, 2: SymBlockGaussSeidel, 3: Jacobi) (standart qiymət: 1) (Tövsiyə olunan qiymət 1-dir) |
| ML_MGCycle | I | ML çoxsəviyyəli dövrü (1: V-cycle, 2: W-cycle, 3: Full-V-cycle) (standart qiymət: 1) (Tövsiyə olunan qiymət çətin məsələlər üçün 2, adi məsələlər üçün 1-dir) |
| ML_MaxLevels | I | ML-in maksimum səviyyə sayı (standart qiymət: 10) (Xüsusilə çətin məsələlər üçün kobud tor həlledicisini birbaşa üsul seçərək 2 (yaddaş çatmırsa 3), adi məsələlər üçün 10 tövsiyə olunur) |
| ML_CoarseningScheme | I | ML kobudlaşdırma sxemi (1: UncoupledMIS, 2: METIS, 4: Zoltan, 5: DD) (standart qiymət: 1) (Tövsiyə olunan qiymət 1 və ya 5-dir) |
| ML_NumSweep | I | ML hamarlaşdırıcısının sweep sayı (Chebyshev olduqda polinom dərəcəsi) (standart qiymət: 2) (Tövsiyə olunan qiymət Chebyshev üçün 2, SymBlockGaussSeidel üçün 1-dir) |
PRECOND=22 olduqda (buraxıla bilər)¶
PRECOND üçün 22-dən başqa qiymət göstərildikdə 5-ci sətir nəzərə alınmır. 4-cü və 5-ci sətirlərin hər biri bütövlükdə buraxıla bilər, həmçinin hər sətirin sonundakı elementlər də buraxıla bilər. 0 göstərilən element üçün standart qiymət istifadə olunur.
4-cü sətirdəki 1-dən 7-yə qədər elementlər PRECOND=5 (ML) seçim sətiri ilə eyni ardıcıllıqdadır; buna görə ML seçim sətiri dəyişdirilmədən təkrar istifadə oluna bilər. 8-dən 10-a qədər elementlər SA-AMG-yə məxsusdur.
(4-cü sətir) SA_CoarseSolver, SA_Smoother, SA_MGCycle, SA_MaxLevels, (istifadə edilmir), SA_ChebyDegree, SA_CoarseSize, SA_MaxAggSize, SA_GalerkinLowMem, (istifadə edilmir)
| Dəyişən adı | Atribut | Məzmun |
|---|---|---|
| SA_CoarseSolver | I | Kobud tor həlledicisi (0: avtomatik, 1: hamarlaşdırıcı ilə əvəz et, 2: sıx matrisli birbaşa üsul, 3: MUMPS (paralel birbaşa üsul)) (standart qiymət: 0) (Avtomatik rejimdə MUMPS aktiv edilmiş yığımda MUMPS, digər hallarda sıx matrisli birbaşa üsul istifadə olunur) |
| SA_Smoother | I | Hamarlaşdırıcı (0, 1: Chebyshev) (standart qiymət: 0) (Yalnız Chebyshev dəstəklənir. Başqa qiymət göstərilərsə, xəbərdarlıq çıxarılır və Chebyshev istifadə olunur) |
| SA_MGCycle | I | Çoxsəviyyəli dövr (0: standart (W-cycle), 1: V-cycle, 2: W-cycle) (standart qiymət: 0) (Full-V-cycle dəstəklənmir) |
| SA_MaxLevels | I | Maksimum səviyyə sayı (standart qiymət: 20) |
| (istifadə edilmir) | I | ML kobudlaşdırma sxeminə uyğun mövqe. SA-AMG həmişə əlaqəsiz aqreqasiyadan istifadə etdiyinə görə 0-dan fərqli qiymət xəbərdarlıq yaradır və nəzərə alınmır |
| SA_ChebyDegree | I | Chebyshev hamarlaşdırıcısının polinom dərəcəsi (standart qiymət: 2) |
| SA_CoarseSize | I | Kobud torun dayandırma ölçüsü (standart qiymət: kobud tor həlledicisi kimi MUMPS istifadə edilirsə 50000, əks halda 100) |
| SA_MaxAggSize | I | Aqreqatın maksimum ölçüsü (standart qiymət: 96) |
| SA_GalerkinLowMem | I | Kobud tor matrisinin hesablama üsulu (0: adi (sürətə üstünlük verir), müsbət qiymət: yaddaşa qənaət) (standart qiymət: 0) (Qurulma zamanı yaddaş çatmırsa, müsbət qiymət göstərin) |
| (istifadə edilmir) | I | Gələcək genişləndirmələr üçün ayrılıb |
5-ci sətir yalnız SA-AMG üçündür və ML tərəfindən oxunmur.
(5-ci sətir) SA_Theta, SA_ChebyAlpha, SA_Safety, SA_TaperK, SA_AggOrder, SA_MinAggSize, SA_Verify, SA_DumpVTK
| Dəyişən adı | Atribut | Məzmun |
|---|---|---|
| SA_Theta | R | Güclü əlaqələrin müəyyən edilməsi üçün hədd (standart qiymət: 0.0) (0.0 olduqda matrisin sıfırdan fərqli strukturu birbaşa qraf kimi istifadə olunur) |
| SA_ChebyAlpha | R | Chebyshev hamarlaşdırıcısının hədəflədiyi məxsusi qiymətlər intervalının aşağı həddini təyin edən nisbət (standart qiymət: 20.0) |
| SA_Safety | R | Maksimum məxsusi qiymətin təxmininə tətbiq olunan təhlükəsizlik əmsalı (standart qiymət: 1.1) |
| SA_TaperK | R | Kobudlaşdırmanın tədricilik əmsalı K (standart qiymət: 100) (2-ci səviyyədən etibarən ən azı təxminən K aqreqat qalması üçün kobudlaşdırma daha tədricən aparılır. Mənfi qiymət tədriciliyi söndürür) |
| SA_AggOrder | R | Aqreqasiyanın baxış ardıcıllığı (0: standart (BFS ardıcıllığı), 1: BFS ardıcıllığı, 2: düyün nömrələrinin hash ardıcıllığı (sınaq xarakterli), 3: dərəcənin artan ardıcıllığı (sınaq xarakterli), 4: dərəcənin azalan ardıcıllığı (sınaq xarakterli), mənfi qiymət: düyün nömrələri ardıcıllığı) (standart qiymət: 0) |
| SA_MinAggSize | R | Aqreqatın minimum ölçüsü (standart qiymət: 3) |
| SA_Verify | R | Qurulma zamanı özünüyoxlamanın aparılması (0: aparılmır, müsbət qiymət: aparılır) (standart qiymət: 0) (əsasən sazlama üçün) |
| SA_DumpVTK | R | Ən yuxarı səviyyəli aqreqasiyanın VTK formatında çıxarılması (0: çıxarılmır, müsbət qiymət: çıxarılır) (standart qiymət: 0) (əsasən sazlama üçün) |
İstifadə nümunəsi¶
SSOR önşərtləndirməli CG üsulundan istifadə edin, maksimum iterasiya sayını 10000, kəsmə xətasını 1.0e-8 təyin edin
SSOR önşərtləndirməli GMRES üsulundan istifadə edin, Krılov altfəzası sayını 40, çoxrəngli sıralamanın hədəf rəng sayını 100 təyin edin
ILU(0) önşərtləndirməli CG üsulundan istifadə edin və önşərtləndirmə matrisi hesablanarkən diaqonal elementlərə tətbiq olunan əmsalı 1.1 təyin edin
CG üsulunun önşərtləndiricisini ML ilə AMG üsulu edin
CG üsulunun önşərtləndiricisini ML ilə AMG üsulu, kobud tor həlledicisini isə MUMPS edin (çətin məsələlər üçün)
CG üsulunun önşərtləndiricisini ML ilə AMG üsulu, çoxsəviyyəli dövrü isə W-cycle edin (çətin məsələlər üçün)
CG üsulunun önşərtləndiricisini ML ilə AMG üsulu, kobud tor həlledicisini MUMPS və maksimum səviyyə sayını 2 edin (çox çətin məsələlər üçün)
CG üsulunun önşərtləndiricisi kimi HEC-MW-yə daxil edilmiş SA-AMG istifadə edin
CG üsulunun önşərtləndiricisi kimi SA-AMG, kobud tor həlledicisi kimi MUMPS istifadə edin
Sürtünməli kontakt analizində BiCGSTAB üsulunu SA-AMG önşərtləndirməsi ilə istifadə edin və AMG önşərtləndiricisinin iyerarxiyasını çıxarın