sparserhs__example_8m_source.html

%Example of using MUMPS in matlab with sparse right hansd side


% initialization of a matlab MUMPS structure

id = initmumps;

id = dmumps(id);

load lhr01;

mat = Problem.A;

% JOB = 6 means analysis+facto+solve

id.JOB = 6;

% we set the rigth hand side

id.RHS = ones(size(mat,1),2);

id.RHS(:,2) = 2*id.RHS(:,2);

id.RHS = sparse(id.RHS);

%call to mumps

id = dmumps(id,mat);

if(norm(mat*id.SOL - id.RHS,'inf') > sqrt(eps))

    disp('WARNING : precision may not be OK');

else

    disp('SOLUTION OK');

end

norm(mat*id.SOL - id.RHS,'inf')

% solution OK

% destroy mumps instance

id.JOB = -2;

id = dmumps(id)


MUMPS
!This file is part of MUMPS
Definition cmumps_root.h:2

lhr01
end load lhr01
Definition diagainv_example.m:13

mat
mat
Definition diagainv_example.m:14

the
end diagonal values have been computed in the(sparse) matrix id.SOL

end
end[inform, rinform, sol, inst, schur, redrhs, pivnul_list, sym_perm, uns_perm, icntl, cntl, colsca_out, rowsca_out, keep_out, dkeep_out]
Definition dmumps.m:40

SOL
id SOL
Definition dmumps.m:43

id
initmumps id
Definition multiplerhs_example.m:4

RHS
we set the rigth hand side id RHS
Definition multiplerhs_example.m:11

dmumps
****************************************This help menu gives details about the use of dmumps
Definition mumps_help.m:2

JOB
#define JOB
Definition mumpsmex.c:33

norm
end norm(mat *id.SOL - id.RHS, 'inf') % solution OK % destroy mumps instance id.JOB

disp
if(norm(mat *id.SOL - id.RHS, 'inf') > sqrt(eps)) disp('WARNING els disp)('SOLUTION OK')
Definition sparserhs_example.m:19