<< spilu_ilutplot Graphics spilu_plotparameter >>

Spilu >> Spilu > Graphics > spilu_ilutpplot

spilu_ilutpplot

Plots the sensitivity of ILUT for A.

Calling Sequence

spilu_ilutpplot(A)
data = spilu_ilutpplot(A)
data = spilu_ilutpplot(A,N)
[data,hLfil] = spilu_ilutpplot(...)
[data,hLfil,hDrop] = spilu_ilutpplot(...)
[data,hAlpha,hDrop,hPermtol] = spilu_ilutpplot(...)
[data,hAlpha,hDrop,hPermtol,hBloc] = spilu_ilutpplot(...)

Parameters

A:

a n-by-n sparse matrix

N:

a 1-by-1 matrix of doubles, integer value, the number of points in the plot (default N = 100).

data:

a N-by-6 matrix of doubles, the values of the parameters

hLfil:

a graphics handle, the lfil plot

hDrop:

a graphics handle, the drop plot

hPermtol:

a graphics handle, the permtol plot

hBloc:

a graphics handle, the bloc plot

Description

For one particular matrix, plot norm(A-L*U)/norm(A) and nnz(L)+nnz(U), depending on lfil and drop for the ILUT method.

The columns of the array data are:

If the decomposition does not work for one value of the parameter, then norm(A-L*U)/norm(A) and nnz(L)+nnz(U) are set to Infinity.

Examples

A = [
10.     -6.    -1.     0.     3.
-6.     34.     0.    10.     0.
-1.      0.     1.     0.     0.
0.     10.     0.    10.     0.
3.      0.     0.     0.     3.
];
A = sparse(A);
spilu_ilutpplot(A);

// See on a 225-by-225 sparse matrix
path = spilu_getpath (  );
filename = fullfile(path,"tests","matrices","pde225.mtx");
A=mmread(filename);
[data,hLfil,hDrop,hPermtol,hBloc] = spilu_ilutpplot(A);

Authors


Report an issue
<< spilu_ilutplot Graphics spilu_plotparameter >>