OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
diagainv_example.m File Reference

Go to the source code of this file.

Typedefs

using complex_arithmetic = false

Functions

 if (complex_arithmetic) id
size(mat, 1), size(mat, 1) speye ()
30 ICNTL ()
end diagonal values have been computed in the (sparse) matrix id.SOL
end diagonal values have been computed in which has the same structure as id RHS Compare diagonal of inverse computed by Mumps and by matlab disp (' ')
 disp ('Computing 2-norm of error on diagonal of inverse:')
 norm (diag(diag(diag(inv(mat))) -id.SOL), 2) % destroy mumps instance id.JOB

Variables

initialization of a matlab MUMPS structure id = initmumps
end load lhr01
 mat = Problem.A
end JOB

Typedef Documentation

◆ complex_arithmetic

using complex_arithmetic = false

Definition at line 4 of file diagainv_example.m.

Function Documentation

◆ disp() [1/2]

end diagonal values have been computed in which has the same structure as id RHS Compare diagonal of inverse computed by Mumps and by matlab disp ( ' ' )

◆ disp() [2/2]

disp ( 'Computing 2-norm of error on diagonal of inverse:' )

◆ ICNTL()

30 ICNTL ( )
virtual

Definition at line 51 of file dmumps.m.

◆ if()

◆ norm()

norm ( diag(diag(diag(inv(mat))) -id.SOL) ,
2  )

◆ speye()

size(mat, 1), size(mat, 1) speye ( )
virtual

◆ the()

end diagonal values have been computed in the ( sparse )

Variable Documentation

◆ id

else id = initmumps

Definition at line 7 of file diagainv_example.m.

◆ JOB

end JOB
Initial value:
= 4 means analysis+factorization
id.JOB = 4
initmumps id
#define JOB
Definition mumpsmex.c:33

Definition at line 21 of file diagainv_example.m.

◆ lhr01

end load lhr01

Definition at line 13 of file diagainv_example.m.

◆ mat

if(complex_arithmetic) To test complex version mat = Problem.A

Definition at line 14 of file diagainv_example.m.