Risinātājs un priekšnosacīšana¶

FrontISTR lineāros risinātājus konfigurē, izvēloties tiešo vai iteratīvo metodi; iteratīvās metodes lineāro vienādojumu risināšanai tiek kombinētas ar priekšnosacītāju. MPC apstrādes metode un kontakta BP eliminācija ir papildu izvēles, kas nosaka, kā tiek apstrādāta lineāra sistēma ar daudzpunktu vai kontakta ierobežojumiem.

Funkciju pārskats¶

Risinātāja un priekšnosacīšanas iestatījumus veido šādas izvēles. Vispirms kā lineāro risinātāju izvēlieties iteratīvu vai tiešu metodi. Priekšnosacītāju kombinē tikai ar iteratīvu metodi.

Paralēlai izpildei MPI paralēlismu, kura pamatā ir apgabala dekompozīcija, var kombinēt ar OpenMP pavedienu paralēlismu. Iteratīvās metodes ir standarta izvēle liela mēroga paralēlām analīzēm, savukārt tiešās metodes izmanto, ja ir piesaistītas nepieciešamās bibliotēkas un vajadzīgs robustāks risinājums.

Risinātāja un priekšnosacītāja izvēle¶

Vispirms izlemiet, vai izmantot iteratīvu vai tiešu metodi. Ja izvēlēta iteratīva metode, pēc tam izvēlieties priekšnosacītāju.

Standarta analīzēm parasti pietiek ar noklusējuma konverģences toleranci, maksimālo iterāciju skaitu un priekšnosacītāja diagonālās korekcijas koeficientu. Tikai tad, ja konverģence ir lēna vai risinājums diverģē, secīgi pārskatiet iteratīvo metodi, priekšnosacītāju, MPC apstrādes metodi un kontakta BP elimināciju. Parametru vērtības un sintaksi skatiet !SOLVER atslēgvārda rokasgrāmatā.

Iteratīvās metodes¶

Iteratīvās metodes risina lielas retinātas sistēmas, atkārtoti atjauninot risinājumu ar matricas–vektora reizinājumiem un priekšnosacīšanu. FrontISTR nodrošina šādas iteratīvās metodes.

CG pieņem, ka matrica ir simetriska pozitīvi definita. Ja berzes kontakts, nesimetriska ierobežojumu apstrāde, saistītie efekti vai līdzīgi faktori rada izteiktu nesimetriju, apsveriet BiCGSTAB, GMRES, GPBiCG vai citu nesimetrisku risinātāju. Izmantojot GMRES, Krilova apakštelpas palielināšana var uzlabot konverģenci, bet palielina arī atmiņas patēriņu.

Ja ir norādīts METHOD2, FrontISTR var pārslēgties uz alternatīvu iteratīvo metodi, ja CG diverģē vai pārtrauc darboties. Šī funkcija ir rezerves variants gadījumiem, kad CG ir primārais risinātājs; alternatīvā metode izmanto tos pašus pārējos iestatījumus un datu rindas kā primārais risinātājs.

Priekšnosacīšana¶

Priekšnosacīšana pārveido koeficientu matricu, lai uzlabotu iteratīvās metodes konverģenci, un tiek piemērota tikai tad, ja izvēlēta iteratīva metode. Tiešās metodes priekšnosacīšanu neizmanto. Pat vienai un tai pašai iteratīvajai metodei priekšnosacītāja izvēle var būtiski mainīt konverģencei nepieciešamo iterāciju skaitu un aprēķina laiku.

SSOR ir standarta izvēle, un OpenMP paralēlisma gadījumā tiek izmantota versija ar daudzkrāsu kārtošanu. Diagonālā mērogošana ir viegla, lai gan tās konverģences uzlabojums ir atkarīgs no uzdevuma. BILU ir realizēta kā bloku nepilnā LU faktorizācija un atbalsta arī matricas ar vispārīgu brīvības pakāpju skaitu. BILU saimei ir ceļš, kas iteratīvās metodes diverģences gadījumā automātiski palielina diagonālās korekcijas koeficientu SIGMA_DIAG un mēģina vēlreiz; vērtības norādīšanu skatiet atslēgvārdu rokasgrāmatā.

AMG izmanto Trilinos-ML, tādēļ būvēšanas laikā jāiespējo ML. Var norādīt izlīdzinātājus, daudzrežģu ciklus, rupjināšanas metodes un citus iestatījumus, bet detalizētas vērtības skatiet atslēgvārdu rokasgrāmatā. Ja SAINV izmanto paralēlā analīzē ar kontakta vai MPC ierobežojumiem, konverģenci var ietekmēt apgabala dekompozīcijā izmantotais pārklājuma dziļums.

Tiešās metodes¶

Tiešās metodes risina lineārus vienādojumus, faktorizējot koeficientu matricu. Tā kā tās nav atkarīgas no iterāciju skaita, tās ir robustas un var būt stabila izvēle analīzēm ar kontakta vai ierobežojumu nosacījumiem. Tomēr to atmiņas patēriņš palielinās līdz ar uzdevuma izmēru.

Tiešās metodes var izmantot gan simetriskām, gan nesimetriskām matricām. MUMPS un MKL nepieciešams būvēšanas laikā sasaistīt attiecīgās bibliotēkas. Vidēs bez šīm bibliotēkām tās nevar izmantot, tādēļ nepieciešamās atkarības un CMake opcijas skatiet būvēšanas rokasgrāmatā.

FrontISTR satur arī iebūvētu tiešā risinātāja ceļu, kas neizmanto ārēju bibliotēku, taču parastām analīzēm, kurām vajadzīga tiešā metode, parasti vispirms jāapsver MUMPS vai MKL.

Tiešajām metodēm priekšnosacīšana netiek piemērota. Pat ja !SOLVER ir norādīts priekšnosacītājs, tiešā risinātāja ceļā tas netiek izmantots.

MPC apstrādes metodes¶

MPC apstrādes metode nosaka, kā lineārajos vienādojumos tiek iekļautas ar daudzpunktu ierobežojumiem saistītās brīvības pakāpes. Pat ar vienādiem ierobežojumu vienādojumiem ievades datos apstrādes metode maina koeficientu matricas īpašības un konverģences uzvedību.

Soda metode ir noklusējums tiešajām metodēm, bet tieša BP eliminācija — iteratīvajām metodēm. MPC-CG metode ir mantota savietojamības opcija un jaunām analīzēm nav ieteicama; parasti izmantojiet tiešu BP elimināciju vai izvēlētās tiešās metodes noklusējuma apstrādi.

MPC apstrādes metodi var skaidri norādīt ar MPCMETHOD elementā !SOLVER. Katras vērtības nozīmi un sintaksi skatiet atslēgvārdu rokasgrāmatā.

Kontakta BP eliminācija¶

Kontakta BP eliminācija samazina SLAGRANGE kontakta ieviestās papildu kontakta brīvības pakāpes pirms sistēmas nodošanas lineārajam risinātājam. Šo kontakta brīvības pakāpju eliminācija var uzlabot iteratīvās metodes apstrādātās lineārās sistēmas īpašības.

Noklusējuma automātiskajā režīmā kontakta BP eliminācija ir iespējota iteratīvajām metodēm un atspējota tiešajām metodēm. Tiešajām metodēm ir paredzēts ceļš, kas tieši risina paplašināto sistēmu ar kontakta brīvības pakāpēm.

Ja CONTACT_ELIM ir norādīts skaidri, kontakta BP elimināciju var veikt arī tad, ja izmanto tiešo metodi. Kontakta tipa, kontakta pāra un kontakta algoritma ievades izvēles skatiet lapā Kontakts un iegulšana.

Palīgfunkcijas¶

Nosacījuma skaitļa novērtēšanu, matricas izgāznes un žurnalēšanu var izmantot, lai pētītu risinātāja konverģenci un veiktspēju. Parastām analīzēm šīs funkcijas nav nepieciešamas un tās galvenokārt iespējo atkļūdošanai vai veiktspējas novērtēšanai.

Matricas izgāšana ļauj ar ārējiem rīkiem pārbaudīt analīzē samontēto koeficientu matricu. Ja ir iespējots DUMPEXIT, analīze beidzas pēc matricas un labās puses saglabāšanas. Nosacījuma skaitļa novērtēšana ir eksperimentāla funkcija, kas pieejama ar CG un GMRES, un to var izmantot konverģences uzvedības diagnosticēšanai.

Saistītās tēmas¶

• Lineārais risinātājs (teorija) — iteratīvo metožu un priekšnosacīšanas matemātiskie principi.
• Nelineārā iterācija un laika integrēšana — Ņūtona iterāciju, kontakta iterāciju un laika integrēšanas ietvars, kas izsauc lineāro risinātāju.
• Apgabala dekompozīcija — MPI paralēlā analīze un pārklājuma dziļuma iestatījumi.
• Kontakts un iegulšana — kontakta tipu un kontakta algoritmu izvēle.
• !SOLVER — risinātāja un priekšnosacītāja parametru vērtības un sintakse.
• Obligātās un neobligātās atkarību bibliotēkas — atkarības, piemēram, MUMPS, MKL (Intel PARDISO), Trilinos-ML un LAPACK.