Matrix twotone

???

performance of incomplete LU preconditioners

Abbreviations for Preconditioners
ILUPACK	inverse-based Multilevel ILU. Uses the concept of ARMS and an inverse based partial ILUC
ILUPACK - MEDIUM VERSION	medium Schur-complement used, less sensitive w.r.t. drop tolerances
ILUPACK - SIMPLE VERSION	simple Schur-complement used, more sensitive w.r.t. drop tolerances, but often faster if it works
ILUTP	incomplete LU decomposition with pivoting from SPARSKIT, based on a recent binary search tree implementation
IB-ILUC	inverse-based incomplete LU decomposition WITHOUT pivoting, based on the Crout version. The norm of the inverse triangular factors is included in the dropping strategy
PARDISO	sparse direct solver. For details see the PARDISO home page

Abbreviations for Orderings
PQ	new PQ ordering for multilevel ILU
RCM	Reverse Cuthill-McKee
MMD	minimum degree
AMF	approximate inverse fill
ND	MeTiS multilevel (node) nested dissection

Performance - Best Computation Time

	Computation Time
time [sec]
	drop tolerance

	Memory Requirement
nz(L+U) ---------- nz(A)
	drop tolerance

Performance - Best Memory Requirement

	Computation Time
time [sec]
	drop tolerance

	Memory Requirement
nz(L+U) ---------- nz(A)
	drop tolerance

Performance - Most Robust Method

	Computation Time
time [sec]
	drop tolerance

	Memory Requirement
nz(L+U) ---------- nz(A)
	drop tolerance

ILUPACK - MEDIUM VERSION

ILUPACK - MEDIUM VERSION + PQ
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
2.5e-01	1.0e+01	1.7	3.6e+01	1.1e+01	4.70e+01	51
2.5e-01	2.5e+01	2.2	5.4e+01	1.1e+01	6.50e+01	44
2.5e-01	5.0e+01	2.3	8.6e+01	1.2e+01	9.80e+01	44
2.5e-01	1.0e+02	2.3	9.1e+01	1.2e+01	1.03e+02	51
1.0e-01	2.5e+01	2.4	7.2e+01	5.3e+00	7.73e+01	21
1.0e-01	5.0e+01	3.2	8.5e+01	6.0e+00	9.10e+01	20
1.0e-01	1.0e+02	3.1	9.0e+01	5.4e+00	9.54e+01	20
5.0e-02	1.0e+01	3.3	8.5e+01	2.1e+00	8.71e+01	8
1.0e-01	1.0e+01	2.4	5.6e+01	4.6e+00	6.06e+01	21

ILUPACK - MEDIUM VERSION + PQ + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-01	1.0e+01	2.5	5.3e+01	2.4e+01	7.70e+01	88
1.0e-01	2.5e+01	2.4	6.6e+01	2.3e+01	8.90e+01	89
1.0e-01	5.0e+01	2.3	8.8e+01	2.9e+01	1.17e+02	112
1.0e-01	1.0e+02	2.3	9.8e+01	2.7e+01	1.25e+02	111
5.0e-02	1.0e+01	2.9	8.9e+01	3.6e+00	9.26e+01	13
5.0e-02	2.5e+01	2.7	1.1e+02	3.7e+00	1.14e+02	14
5.0e-02	5.0e+01	2.9	1.1e+02	3.7e+00	1.14e+02	14
5.0e-02	1.0e+02	3.1	1.1e+02	3.9e+00	1.14e+02	13
2.5e-02	1.0e+01	2.9	9.9e+01	2.8e+00	1.02e+02	10
2.5e-02	2.5e+01	3.1	1.3e+02	3.3e+00	1.33e+02	10
2.5e-02	5.0e+01	3.2	1.4e+02	2.7e+00	1.43e+02	10
2.5e-02	1.0e+02	3.8	1.0e+02	3.4e+00	1.03e+02	10
1.0e-02	1.0e+01	3.2	1.2e+02	1.8e+00	1.22e+02	7

ILUPACK - MEDIUM VERSION + RCM
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-05	1.0e+02	11.1	4.2e+02	inf	Inf	inf

ILUPACK - MEDIUM VERSION + RCM + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
2.5e-01	2.5e+01	3.3	2.0e+02	7.9e+01	2.79e+02	217
2.5e-01	1.0e+02	5.9	8.6e+02	1.5e+01	8.75e+02	35
1.0e-01	1.0e+01	4.7	1.9e+02	1.9e+01	2.09e+02	45
1.0e-01	2.5e+01	4.5	3.4e+02	2.4e+01	3.64e+02	62
1.0e-01	5.0e+01	8.0	7.7e+02	2.1e+01	7.91e+02	27
1.0e-01	1.0e+02	8.3	1.2e+03	1.0e+01	1.21e+03	14
5.0e-02	1.0e+01	5.0	2.3e+02	9.5e+00	2.40e+02	18
5.0e-02	2.5e+01	5.8	5.8e+02	5.7e+00	5.86e+02	14
5.0e-02	1.0e+02	9.1	1.4e+03	8.9e+00	1.41e+03	11
2.5e-02	1.0e+01	6.2	3.0e+02	4.6e+00	3.05e+02	10
2.5e-02	2.5e+01	10.4	8.4e+02	7.1e+00	8.47e+02	9

ILUPACK - MEDIUM VERSION + MMD
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-03	2.5e+01	3.1	8.8e+02	inf	Inf	inf

ILUPACK - MEDIUM VERSION + MMD + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
2.5e-01	2.5e+01	2.3	5.2e+02	7.5e+01	5.95e+02	243
2.5e-01	5.0e+01	3.0	7.3e+02	1.8e+01	7.48e+02	56
2.5e-01	1.0e+02	4.1	9.5e+02	1.4e+01	9.64e+02	31
1.0e-01	1.0e+01	2.3	4.0e+02	1.5e+01	4.15e+02	54
1.0e-01	2.5e+01	2.9	6.2e+02	8.9e+00	6.29e+02	29
1.0e-01	5.0e+01	4.2	9.1e+02	1.3e+01	9.23e+02	26
1.0e-01	1.0e+02	4.8	1.5e+03	7.5e+00	1.51e+03	17
5.0e-02	1.0e+01	2.7	4.7e+02	7.3e+00	4.77e+02	22
5.0e-02	2.5e+01	3.3	7.5e+02	5.5e+00	7.56e+02	17
5.0e-02	5.0e+01	4.6	1.0e+03	5.9e+00	1.01e+03	14
5.0e-02	1.0e+02	5.2	1.8e+03	5.0e+00	1.80e+03	12
2.5e-02	1.0e+01	3.3	5.3e+02	5.1e+00	5.35e+02	15
2.5e-02	2.5e+01	3.8	9.9e+02	4.0e+00	9.94e+02	11
1.0e-02	1.0e+01	4.0	6.9e+02	3.4e+00	6.93e+02	10
1.0e-02	2.5e+01	5.1	1.2e+03	3.9e+00	1.20e+03	9

ILUPACK - MEDIUM VERSION + AMF
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-05	1.0e+02	5.7	3.6e+02	inf	Inf	inf

ILUPACK - MEDIUM VERSION + AMF + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
5.0e-01	5.0e+01	2.5	2.8e+02	3.8e+01	3.18e+02	129
5.0e-01	1.0e+02	4.1	6.1e+02	1.1e+02	7.20e+02	269
2.5e-01	2.5e+01	2.5	1.8e+02	2.2e+01	2.02e+02	78
2.5e-01	5.0e+01	3.0	4.2e+02	1.5e+01	4.35e+02	54
2.5e-01	1.0e+02	4.6	8.4e+02	2.6e+01	8.66e+02	62
1.0e-01	1.0e+01	2.5	1.6e+02	1.3e+01	1.73e+02	47
1.0e-01	2.5e+01	3.2	3.5e+02	9.2e+00	3.59e+02	27
1.0e-01	5.0e+01	4.5	6.2e+02	8.2e+00	6.28e+02	21
5.0e-02	1.0e+01	3.0	2.8e+02	6.8e+00	2.87e+02	22
5.0e-02	2.5e+01	3.6	4.8e+02	5.2e+00	4.85e+02	16
2.5e-02	1.0e+01	4.2	3.8e+02	5.1e+00	3.85e+02	14
2.5e-02	2.5e+01	5.3	5.7e+02	5.0e+00	5.75e+02	11
1.0e-02	1.0e+01	5.0	5.3e+02	4.4e+00	5.34e+02	10
1.0e-02	2.5e+01	5.9	9.1e+02	4.3e+00	9.14e+02	8

ILUPACK - MEDIUM VERSION + Multilevel ND
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter

ILUPACK - MEDIUM VERSION + Multilevel ND + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter

ILUPACK - SIMPLE VERSION

ILUPACK - SIMPLE VERSION + PQ
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-02	1.0e+01	3.1	6.6e+01	6.8e+01	1.34e+02	211
1.0e-02	2.5e+01	3.0	7.5e+01	6.1e+01	1.36e+02	208
1.0e-02	5.0e+01	2.9	6.9e+01	5.8e+01	1.27e+02	199
1.0e-02	1.0e+02	3.0	7.6e+01	7.3e+01	1.49e+02	206
5.0e-03	1.0e+01	3.3	6.6e+01	5.4e+01	1.20e+02	147
5.0e-03	2.5e+01	3.2	7.0e+01	5.4e+01	1.24e+02	155
5.0e-03	5.0e+01	3.2	6.6e+01	4.9e+01	1.15e+02	155
5.0e-03	1.0e+02	3.2	7.3e+01	5.1e+01	1.24e+02	153
2.5e-03	1.0e+01	3.1	7.1e+01	1.1e+01	8.20e+01	30
2.5e-03	2.5e+01	3.4	6.5e+01	1.1e+01	7.60e+01	30
2.5e-03	5.0e+01	3.5	8.1e+01	1.1e+01	9.20e+01	30
2.5e-03	1.0e+02	3.5	7.2e+01	9.7e+00	8.17e+01	30
1.0e-03	1.0e+01	3.8	8.2e+01	2.0e+00	8.40e+01	6

ILUPACK - SIMPLE VERSION + PQ + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-02	1.0e+01	2.8	5.0e+01	1.5e+02	2.00e+02	496
1.0e-02	5.0e+01	2.7	5.5e+01	5.1e+01	1.06e+02	170
1.0e-02	1.0e+02	2.8	6.5e+01	5.5e+01	1.20e+02	183
5.0e-03	1.0e+01	2.9	5.7e+01	5.1e+01	1.08e+02	163
5.0e-03	2.5e+01	3.4	7.3e+01	5.2e+01	1.25e+02	141
5.0e-03	5.0e+01	3.4	7.7e+01	6.4e+01	1.41e+02	184
5.0e-03	1.0e+02	3.4	7.6e+01	5.8e+01	1.34e+02	193
2.5e-03	1.0e+01	3.1	6.3e+01	7.8e+00	7.08e+01	22
2.5e-03	2.5e+01	3.3	6.8e+01	6.8e+00	7.48e+01	24
2.5e-03	5.0e+01	3.7	9.3e+01	9.4e+00	1.02e+02	24
2.5e-03	1.0e+02	3.7	7.5e+01	9.3e+00	8.43e+01	24
1.0e-03	1.0e+01	3.3	6.5e+01	3.2e+00	6.82e+01	12
1.0e-03	2.5e+01	3.2	8.7e+01	3.8e+00	9.08e+01	12
1.0e-03	5.0e+01	3.3	7.2e+01	3.2e+00	7.52e+01	11
1.0e-03	1.0e+02	3.3	8.2e+01	3.5e+00	8.55e+01	11
1.0e-04	1.0e+01	3.7	8.7e+01	1.3e+00	8.83e+01	4

ILUPACK - SIMPLE VERSION + RCM
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-05	1.0e+02	6.8	1.3e+02	inf	Inf	inf

ILUPACK - SIMPLE VERSION + RCM + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
5.0e-02	1.0e+02	9.0	5.2e+02	1.4e+01	5.34e+02	16
2.5e-02	2.5e+01	5.1	2.0e+02	2.2e+02	4.20e+02	477
2.5e-02	5.0e+01	6.3	3.1e+02	9.2e+00	3.19e+02	17
2.5e-02	1.0e+02	9.1	5.6e+02	9.2e+00	5.69e+02	12
1.0e-02	2.5e+01	5.5	2.2e+02	7.4e+00	2.27e+02	17
1.0e-02	5.0e+01	9.5	5.5e+02	9.0e+00	5.59e+02	10
1.0e-02	1.0e+02	9.6	5.4e+02	7.2e+00	5.47e+02	8
5.0e-03	1.0e+01	5.9	2.0e+02	4.2e+01	2.42e+02	86
5.0e-03	2.5e+01	6.3	2.9e+02	6.0e+00	2.96e+02	13
2.5e-03	1.0e+01	6.2	2.1e+02	9.8e+00	2.20e+02	20

ILUPACK - SIMPLE VERSION + MMD
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-05	1.0e+02	6.6	1.0e+03	inf	Inf	inf

ILUPACK - SIMPLE VERSION + MMD + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-01	5.0e+01	4.1	4.0e+02	1.5e+02	5.50e+02	443
1.0e-01	1.0e+02	4.7	6.3e+02	2.8e+01	6.58e+02	78
5.0e-02	2.5e+01	3.2	3.6e+02	8.6e+01	4.46e+02	302
5.0e-02	5.0e+01	4.4	4.8e+02	2.6e+01	5.06e+02	74
5.0e-02	1.0e+02	5.1	7.5e+02	1.0e+01	7.60e+02	27
2.5e-02	2.5e+01	3.8	4.9e+02	1.6e+01	5.06e+02	54
2.5e-02	5.0e+01	5.0	5.6e+02	9.0e+00	5.69e+02	24
2.5e-02	1.0e+02	5.9	8.6e+02	4.6e+00	8.65e+02	12
1.0e-02	1.0e+01	3.8	3.9e+02	1.3e+01	4.03e+02	42
1.0e-02	2.5e+01	5.0	4.8e+02	4.9e+00	4.85e+02	13
1.0e-02	5.0e+01	5.6	8.1e+02	4.6e+00	8.15e+02	12
1.0e-02	1.0e+02	6.7	1.2e+03	3.3e+00	1.20e+03	8
5.0e-03	1.0e+01	4.2	5.4e+02	5.0e+00	5.45e+02	15
5.0e-03	2.5e+01	5.6	6.1e+02	4.0e+00	6.14e+02	10
5.0e-03	5.0e+01	6.0	8.5e+02	4.0e+00	8.54e+02	10
2.5e-03	1.0e+01	4.7	7.5e+02	3.7e+00	7.54e+02	11

ILUPACK - SIMPLE VERSION + AMF
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-05	1.0e+02	5.6	1.6e+02	inf	Inf	inf

ILUPACK - SIMPLE VERSION + AMF + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
2.5e-01	1.0e+02	4.4	1.8e+02	1.6e+02	3.40e+02	453
1.0e-01	5.0e+01	4.2	1.4e+02	2.1e+01	1.61e+02	50
1.0e-01	1.0e+02	5.1	2.4e+02	1.2e+01	2.52e+02	31
5.0e-02	2.5e+01	3.5	1.0e+02	1.3e+01	1.13e+02	47
5.0e-02	5.0e+01	4.7	1.8e+02	9.5e+00	1.90e+02	25
5.0e-02	1.0e+02	5.9	2.9e+02	6.2e+00	2.96e+02	16
2.5e-02	1.0e+01	4.0	1.1e+02	1.9e+01	1.29e+02	59
2.5e-02	2.5e+01	5.0	1.8e+02	7.1e+00	1.87e+02	16
2.5e-02	5.0e+01	5.2	2.5e+02	5.9e+00	2.56e+02	13
2.5e-02	1.0e+02	6.4	3.1e+02	4.3e+00	3.14e+02	11
1.0e-02	1.0e+01	4.6	1.4e+02	7.9e+00	1.48e+02	24
1.0e-02	2.5e+01	5.6	2.2e+02	4.0e+00	2.24e+02	10
1.0e-02	5.0e+01	5.9	2.9e+02	4.0e+00	2.94e+02	9
5.0e-03	1.0e+01	5.3	1.9e+02	4.2e+00	1.94e+02	12

ILUPACK - SIMPLE VERSION + Multilevel ND
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter

ILUPACK - SIMPLE VERSION + Multilevel ND + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-01	5.0e+01	4.7	2.7e+02	5.5e+01	3.25e+02	114
1.0e-01	1.0e+02	4.9	2.7e+02	2.0e+01	2.90e+02	30
5.0e-02	2.5e+01	3.7	1.3e+02	1.1e+02	2.40e+02	238
5.0e-02	5.0e+01	5.4	3.1e+02	1.2e+01	3.22e+02	27
5.0e-02	1.0e+02	5.4	3.2e+02	6.8e+00	3.27e+02	16
2.5e-02	2.5e+01	4.3	1.8e+02	1.3e+01	1.93e+02	30
2.5e-02	5.0e+01	5.8	3.5e+02	1.2e+01	3.62e+02	16
2.5e-02	1.0e+02	6.2	4.3e+02	8.7e+00	4.39e+02	12
1.0e-02	1.0e+01	4.3	1.5e+02	1.8e+01	1.68e+02	34
1.0e-02	2.5e+01	6.1	3.5e+02	9.4e+00	3.59e+02	16
1.0e-02	5.0e+01	6.9	4.2e+02	7.2e+00	4.27e+02	10
1.0e-02	1.0e+02	7.1	4.5e+02	3.0e+00	4.53e+02	7
5.0e-03	1.0e+01	4.8	1.8e+02	6.6e+00	1.87e+02	17
2.5e-03	1.0e+01	5.4	1.7e+02	3.8e+00	1.74e+02	12

ILUTP

ILUTP + RCM
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-04		18.4	6.7e+02	7.1e+00	6.77e+02	11
1.0e-05		19.7	7.5e+02	2.9e+00	7.53e+02	5

ILUTP + RCM + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-03		-0.0	2.9e+03	-	-	failure
1.0e-05		35.6	3.2e+03	7.8e+00	3.21e+03	8

ILUTP + MMD
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-04		30.0	2.7e+03	2.2e+01	2.72e+03	24
1.0e-05		33.6	3.0e+03	5.9e+00	3.01e+03	6

ILUTP + MMD + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-04		7.2	2.7e+02	3.5e+01	3.05e+02	115
1.0e-05		7.7	3.0e+02	1.8e+00	3.02e+02	6

ILUTP + AMF
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-01		-0.0	2.0e+03	-	-	failure
5.0e-01		4.8	8.5e+01	inf	Inf	inf

ILUTP + AMF + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-01		-0.0	1.8e+03	-	-	failure
5.0e-01		1.6	1.7e+01	inf	Inf	inf

ILUTP + Multilevel ND
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-03		0.0	1.3e+03	-	-	failure
1.0e-02		18.7	7.9e+02	inf	Inf	inf

ILUTP + Mutilevel ND + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-04		15.2	4.7e+02	1.3e+02	6.00e+02	185
1.0e-05		16.3	5.0e+02	3.5e+00	5.04e+02	7

ILUC

ILUC + RCM
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
2.5e-01		-0.0	3.5e+02	-	-	failure
5.0e-01		14.7	2.1e+02	NaN	NaN	NaN

ILUC + RCM + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
2.5e-02		-0.0	3.8e+02	-	-	failure
1.0e-03		19.5	3.4e+02	5.5e+01	3.95e+02	83
1.0e-04		21.7	3.6e+02	1.3e+01	3.73e+02	17
1.0e-05		24.7	4.2e+02	5.7e+00	4.26e+02	7

ILUC + MMD
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
5.0e-01		-0.0	6.9e+02	-	-	failure

ILUC + MMD + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-03		6.4	1.8e+02	2.0e+01	2.00e+02	73
1.0e-04		7.0	1.9e+02	1.8e+00	1.92e+02	7

ILUC + AMF
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-03		-0.0	1.0e+03	-	-	failure
1.0e-02		12.1	4.9e+02	NaN	NaN	NaN

ILUC + AMF + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
5.0e-01		-0.0	1.0e+03	-	-	failure

ILUC + Multilevel ND
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-07		21.8	3.4e+02	NaN	NaN	NaN

ILUC + Multilevel ND + MC64
parameters			time [sec]
drop tol	condest	fill-in	ILU	GMRES(30)	total	iter
1.0e-03		8.3	1.4e+02	4.5e+01	1.85e+02	123
1.0e-04		9.0	1.4e+02	4.8e+00	1.45e+02	9

PARDISO

PARDISO
parameters			time [sec]
		fill-in	LU fac.	LU solve	total	iter
		20.6	4.1e+02	1.3e+00	4.2e+02	0

PARDISO + MC64
parameters			time [sec]
		fill-in	LU fac.	LU solve	total	iter
		22.1	3.5e+02	1.3e+00	3.5e+02	0

m.bollhoefer@tu-bs.de

Last modified: October 11, 2011