FINE_LEVEL parametry sterujące na siatce ostatecznej (najbardziej zaadaptowanej) 1 1 50 wybór_solwera liczba_poziomów liczba_wektorów_Kryłowa 4 1 1 wybór_preconditionera liczba_wewnętrznych_iteracji wybór_formatu_przechowania 100 0 1.e-6 maksymalna_liczba_iteracji wybór_rodzaju_zbieżności tolerancja_błędu 3 1 1 poziom_szczegółowości_wydruku liczba_iteracji_dla_wersji_multigrid: _pre _post Opcje do testowania MMNiT: 1. niezależnie od wyboru solwera: maksymalna_liczba_iteracji - można rozpocząć od 100, zwiększając w razie potrzeby, tak aby osiągnąć zbieżność (wartość 100 gwarantuje, że w przypadku braku zbieżności obliczenia szybko się zakończą) poziom_szczegółowości_wydruku - można zacząć od 3 (duża liczba wydruków informacyjnych), zmniejszając dla osiągnięcia większej czytelności (poziom 3 może być konieczny do uzyskania danych do sprawozdania) wybór_solwera: 10 - metoda Jacobiego (wybór_preconditionera 1) i Gaussa-Seidla (wybór_preconditionera 2) 1 - metoda GMRES wybór_preconditionera - 0 brak (wybór_formatu_przechowania -2) wybór_preconditionera - 1 metoda Jacobiego (wybór_formatu_przechowania -2) wybór_preconditionera - 2 metoda Gaussa-Seidla (wybór_formatu_przechowania -2)) wybór_preconditionera - 4 metoda niekompletnego rozkładu LU (wybór_formatu_przechowania 1) Pozostałe (do pozostawienia bez zmian) liczba_poziomów: 1 (wartości >1 tylko dla wersji multigrid) liczba_wektorów_Kryłowa: 50 (relatywnie duża wartość dla zapewnienia szybkiej zbieżności) liczba_wewnętrznych_iteracji: 1 (zazwyczaj optymalna dla preconditionerów BJ i BGS) wybór_rodzaju_zbieżności: 0 (względna redukcja normy residuum) tolerancja_błędu: 1.e-6 (rozsądne optimum dla zadañ liniowych, stacjonarnych) liczba_iteracji_dla_wersji_multigrid: _pre _post: 1 1 (znaczące tylko dla multigrid) ***************************************************************** Rozbudowana wersja pliku konfiguracyjnego z szeregiem opcji solwera (w tym dla wersji multigrid) FINE_LEVEL (the only level for single-level solvers) 1 1 50 single/multi-level GMRES, number of levels, number of Krylov vectors 1 1 -3 type of preconditioner, number of internal sweeps, block type 10 0 1e-3 max_iter, convergence type, convergence treshold 3 1 1 output level, number of pre-smooth and post-smooth sweeps COARSE_LEVEL (the coarsest level solver) 1 -1 100 single/multi-level GMRES, coarse pdeg, number of Krylov vectors 2 1 2 type of preconditioner, number of internal sweeps, block type 100 0 1e-12 max_iter, convergence type, convergence treshold 0 output level INTER_LEVELS (smoothers/preconditioners for intermediate levels) 1 1 -1 number of pre-smooth sweeps, number of post-smooth sweeps, coarse_pdeg 2 1 2 type of preconditioner, number of internal sweeps, block type 0 output level *************************************************************************** FOR MPI: number of internal sweeps = 1; block type = 1; 2 1 1 - with GS 1 1 1 - with BJ FOR GS AND PRISMATIC ELEMENTS TRY BLOCKS IN ORDER: 2=6, 5, 4 coarse_pdeg = -1 - use pdeg from fine level Options for GMRES solver (on fine level a,g,h,i,k can be overriden by fem code): a b c a-solver: 1-single level GMRES, 2-multi-level/multigrid GMRES b-type: 0-standard, 1-matrix free (not present in this version of code) c-number of Krylov vectors d e f d-preconditioner (only at the current level): 0-no, 1-block Jacobi, 2-block Gauss-Seidel (both as smoothers) 3-additive Schwarz (block Jacobi as preconditioner) 4-ILU(0) (as smoother), 5-ILU(0) (as preconditioner) e-number of internal sweeps for block Jacobi and block Gauss-Seidel f-type of blocks (most useful options): 1-one-element blocks 2,3,4,5,6,7-larger blocks with overlap (for block smoothers) g h i g-maximal number of single iterations (not restarts!) h-convergence criterion for the norm of residual: /* Types of convergence measures from lsd_mkb/lsd_mkb_core/lsh_mkb_core.h */ #define REL_RES_INI 0 #define ABS_RES 1 #define REL_RES_RHS 2 i-convergence treshold k k-output level 0-no information is printed to output 1-only error information printed 2-important (convergence) information printed 3-all available information printed Options for standard iterations solver (for fine level a,g,h,i,k as above): a b c a-solver: 10-standard iterations, 20-multigrid b-number of pre-smooth steps (only for multigrid on the finest level) c-number of post-smooth steps (only for multigrid on the finest level) d e f d-actual solver/smoother (only at the current level): 0-no, 1-block Jacobi, 2-block Gauss-Seidel (both as smoothers) 3-additive Schwarz (block Jacobi as preconditioner) 4-ILU(0) (as smoother), 5-ILU(0) (as preconditioner) e-number of internal sweeps for block Jacobi and block Gauss-Seidel f-type of blocks (most useful options): 0-one element blocks (blsiter) 1-one-element blocks (blliter) 2-larger blocks with no overlap (for ILU(0)) 4-larger blocks with overlap (for block smoothers) -1 -2 -3 -4 g h i g-maximal number of single iterations h-trash (max norm of update used as convergence criterion) i-convergence treshold k k-output level (the same as for the GMRES solver) Options for intermediate level solvers/smoothers: b c b-number of pre-smooth steps c-number of post-smooth steps d e f d-smoother at the current level: 0-no, 1-block Jacobi, 2-block Gauss-Seidel (both as smoothers) 3-additive Schwarz (block Jacobi as preconditioner) 4-ILU(0) (as smoother), 5-ILU(0) (as preconditioner) e-number of internal sweeps for block Jacobi and block Gauss-Seidel f-type of blocks (most useful options): 1-one-element blocks 2-larger blocks with no overlap (for ILU(0)) 4-larger blocks with overlap (for block smoothers) k k-output level (the same as for the GMRES solver)