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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML3.2 EN">
<HTML>
<HEAD> <link rel="canonical" href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MATSOLVERMKL_PARDISO.html" />
<META NAME="GENERATOR" CONTENT="DOCTEXT">
<TITLE>MATSOLVERMKL_PARDISO</TITLE>
</HEAD>
<BODY BGCOLOR="FFFFFF">
<div id="version" align=right><b>petsc-3.7.5 2017-01-01</b></div>
<div id="bugreport" align=right><a href="mailto:petsc-maint@mcs.anl.gov?subject=Typo or Error in Documentation &body=Please describe the typo or error in the documentation: petsc-3.7.5 v3.7.5 docs/manualpages/Mat/MATSOLVERMKL_PARDISO.html "><small>Report Typos and Errors</small></a></div>
<A NAME="MATSOLVERMKL_PARDISO"><H1>MATSOLVERMKL_PARDISO</H1></A>
A matrix type providing direct solvers (LU) for sequential matrices via the external package MKL_PARDISO. Works with <A HREF="../Mat/MATSEQAIJ.html#MATSEQAIJ">MATSEQAIJ</A> matrices
<P>
Use -pc_type lu -pc_factor_mat_solver_package mkl_pardiso to us this direct solver
<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_mkl_pardiso_65 </B></TD><TD>- Number of threads to use within MKL_PARDISO
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_66 </B></TD><TD>- Maximum number of factors with identical sparsity structure that must be kept in memory at the same time
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_67 </B></TD><TD>- Indicates the actual matrix for the solution phase
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_68 </B></TD><TD>- Message level information
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_69 </B></TD><TD>- Defines the matrix type. IMPORTANT: When you set this flag, iparm parameters are going to be set to the default ones for the matrix type
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_1 </B></TD><TD>- Use default values
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_2 </B></TD><TD>- Fill-in reducing ordering for the input matrix
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_4 </B></TD><TD>- Preconditioned CGS/CG
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_5 </B></TD><TD>- User permutation
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_6 </B></TD><TD>- Write solution on x
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_8 </B></TD><TD>- Iterative refinement step
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_10 </B></TD><TD>- Pivoting perturbation
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_11 </B></TD><TD>- Scaling vectors
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_12 </B></TD><TD>- Solve with transposed or conjugate transposed matrix A
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_13 </B></TD><TD>- Improved accuracy using (non-) symmetric weighted matching
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_18 </B></TD><TD>- Numbers of non-zero elements
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_19 </B></TD><TD>- Report number of floating point operations
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_21 </B></TD><TD>- Pivoting for symmetric indefinite matrices
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_24 </B></TD><TD>- Parallel factorization control
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_25 </B></TD><TD>- Parallel forward/backward solve control
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_27 </B></TD><TD>- Matrix checker
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_31 </B></TD><TD>- Partial solve and computing selected components of the solution vectors
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_34 </B></TD><TD>- Optimal number of threads for conditional numerical reproducibility (CNR) mode
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>-mat_mkl_pardiso_60 </B></TD><TD>- Intel MKL_PARDISO mode
</TD></TR></TABLE>
<P>
<P>
For more information please check mkl_pardiso manual
<P>
<H3><FONT COLOR="#CC3333">See Also</FONT></H3>
<A HREF="../PC/PCFactorSetMatSolverPackage.html#PCFactorSetMatSolverPackage">PCFactorSetMatSolverPackage</A>(), <A HREF="../Mat/MatSolverPackage.html#MatSolverPackage">MatSolverPackage</A>
<BR>
<P>
<P><B><P><B><FONT COLOR="#CC3333">Level:</FONT></B>beginner
<BR><FONT COLOR="#CC3333">Location:</FONT></B><A HREF="../../../src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.c.html#MATSOLVERMKL_PARDISO">src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.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>
|
