<< spilu_iludpplot Graphics spilu_ilutplot >>

Spilu >> Spilu > Graphics > spilu_ilukplot

spilu_ilukplot

Plots the sensitivity of ILUK for A.

Calling Sequence

spilu_ilukplot(A)
data = spilu_ilukplot(A)
data = spilu_ilukplot(A,N)
[data,hLfil] = spilu_ilukplot(...)

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

Description

For one particular matrix, plot norm(A-L*U)/norm(A) and nnz(L)+nnz(U), depending on lfil for the ILUK 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 = [
30.    8.     7.     0.    0.
8.     18.    0.     0.   -1.
7.     0.     14.    0.    0.
0.     0.     0.     1.    0.
0.    -1.     0.     0.    2.
];
A = sparse(A);
spilu_ilukplot(A);

// See on a 225-by-225 sparse matrix
path = spilu_getpath (  );
filename = fullfile(path,"tests","matrices","pde225.mtx");
A=mmread(filename);
[data,hLfil] = spilu_ilukplot(A);
hLfil.children(1).data_bounds(2,1) = 20;
hLfil.children(2).data_bounds(2,1) = 20;

Authors


Report an issue
<< spilu_iludpplot Graphics spilu_ilutplot >>