1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML3.2 EN">
<HTML>
<HEAD>
<META NAME="GENERATOR" CONTENT="DOCTEXT">
<TITLE>MATAIJMUMPS</TITLE>
</HEAD>
<BODY BGCOLOR="FFFFFF">
<A NAME="MATAIJMUMPS"><H1>MATAIJMUMPS</H1></A>
<A HREF="../Mat/MATAIJMUMPS.html#MATAIJMUMPS">MATAIJMUMPS</A> = "aijmumps" - A matrix type providing direct solvers (LU) for distributed and sequential matrices via the external package MUMPS. If MUMPS is installed (see the manual for instructions
on how to declare the existence of external packages),
a matrix type can be constructed which invokes MUMPS solvers.
After calling <A HREF="../Mat/MatCreate.html#MatCreate">MatCreate</A>(...,A), simply call <A HREF="../Mat/MatSetType.html#MatSetType">MatSetType</A>(A,<A HREF="../Mat/MATAIJMUMPS.html#MATAIJMUMPS">MATAIJMUMPS</A>).
<P>
If created with a single process communicator, this matrix type inherits from <A HREF="../Mat/MATSEQAIJ.html#MATSEQAIJ">MATSEQAIJ</A>.
Otherwise, this matrix type inherits from <A HREF="../Mat/MATMPIAIJ.html#MATMPIAIJ">MATMPIAIJ</A>. Hence for single process communicators,
<A HREF="../Mat/MatSeqAIJSetPreallocation.html#MatSeqAIJSetPreallocation">MatSeqAIJSetPreallocation</A> is supported, and similarly <A HREF="../Mat/MatMPIAIJSetPreallocation.html#MatMPIAIJSetPreallocation">MatMPIAIJSetPreallocation</A> is supported
for communicators controlling multiple processes. It is recommended that you call both of
the above preallocation routines for simplicity. One can also call <A HREF="../Mat/MatConvert.html#MatConvert">MatConvert</A> for an inplace
conversion to or from the <A HREF="../Mat/MATSEQAIJ.html#MATSEQAIJ">MATSEQAIJ</A> or <A HREF="../Mat/MATMPIAIJ.html#MATMPIAIJ">MATMPIAIJ</A> type (depending on the communicator <A HREF="../Sys/size.html#size">size</A>)
without data copy.
<P>
<H3><FONT COLOR="#CC3333">Options Database Keys</FONT></H3>
<TABLE border="0" cellpadding="0" cellspacing="0">
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_type aijmumps </B></TD><TD>- sets the matrix type to "aijmumps" during a call to <A HREF="../Mat/MatSetFromOptions.html#MatSetFromOptions">MatSetFromOptions</A>()
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_sym <0,1,2> </B></TD><TD>- 0 the matrix is unsymmetric, 1 symmetric positive definite, 2 symmetric
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_4 <0,1,2,3,4> </B></TD><TD>- print level
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_6 <0,...,7> </B></TD><TD>- matrix prescaling options (see MUMPS User's Guide)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_7 <0,...,7> </B></TD><TD>- matrix orderings (see MUMPS User's Guide)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_9 <1,2> </B></TD><TD>- A or A^T x=b to be solved: 1 denotes A, 2 denotes A^T
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_10 <n> </B></TD><TD>- maximum number of iterative refinements
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_11 <n> </B></TD><TD>- error analysis, a positive value returns statistics during -ksp_view
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_12 <n> </B></TD><TD>- efficiency control (see MUMPS User's Guide)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_13 <n> </B></TD><TD>- efficiency control (see MUMPS User's Guide)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_14 <n> </B></TD><TD>- efficiency control (see MUMPS User's Guide)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_icntl_15 <n> </B></TD><TD>- efficiency control (see MUMPS User's Guide)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_cntl_1 <delta> </B></TD><TD>- relative pivoting threshold
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_cntl_2 <tol> </B></TD><TD>- stopping criterion for refinement
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mumps_cntl_3 <adelta> </B></TD><TD>- absolute pivoting threshold
</TD></TR></TABLE>
<P>
<P>
<H3><FONT COLOR="#CC3333">See Also</FONT></H3>
<A HREF="../Mat/MATSBAIJMUMPS.html#MATSBAIJMUMPS">MATSBAIJMUMPS</A>
<BR><P><B><P><B><FONT COLOR="#CC3333">Level:</FONT></B>beginner
<BR><FONT COLOR="#CC3333">Location:</FONT></B><A HREF="../../../src/mat/impls/aij/mpi/mumps/mumps.c.html#MATAIJMUMPS">src/mat/impls/aij/mpi/mumps/mumps.c</A>
<BR><A HREF="./index.html">Index of all Mat routines</A>
<BR><A HREF="../../index.html">Table of Contents for all manual pages</A>
<BR><A HREF="../singleindex.html">Index of all manual pages</A>
</BODY></HTML>
|
