Newer
Older
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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.11"/>
<title>LASs - Linear Algebra Routines on OmpSs: src/ddss_dnpgetrf.c Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/searchdata.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
$(document).ready(function() { init_search(); });
</script>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({
extensions: ["tex2jax.js"],
jax: ["input/TeX","output/HTML-CSS"],
});
</script><script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td id="projectalign" style="padding-left: 0.5em;">
<div id="projectname">LASs - Linear Algebra Routines on OmpSs
 <span id="projectnumber">1.0.0</span>
</div>
<div id="projectbrief">LASs</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.11 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
<li>
<div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.png"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
</span>
</div>
</li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>File Members</span></a></li>
</ul>
</div>
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_68267d1309a1af8e8297ef4c3efbcdba.html">src</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">ddss_dnpgetrf.c</div> </div>
</div><!--header-->
<div class="contents">
<a href="ddss__dnpgetrf_8c.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="preprocessor">#include "../include/lass.h"</span></div><div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"></span></div><div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">/**</span></div><div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment"> *</span></div><div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> * @file ddss_dnpgetrf.c</span></div><div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"> *</span></div><div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment"> * @brief LASs-DDSs ddss_dnpgetrf routine.</span></div><div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment"> *</span></div><div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment"> * LASs-DDSs is a software package provided by:</span></div><div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment"> * Barcelona Supercomputing Center - Centro Nacional de Supercomputacion</span></div><div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment"> *</span></div><div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment"> * @author Pedro Valero-Lara pedro.valero@bsc.es</span></div><div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment"> * @author Boro Sofranac boro.sofranac@bsc.es</span></div><div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment"> * @date 2018-04-06</span></div><div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment"> * @reviewer </span></div><div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment"> * @modified </span></div><div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment"> *</span></div><div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment"> **/</span></div><div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment"></span></div><div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment">/**</span></div><div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment"> * </span></div><div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="comment"> * @ingroup DDSS</span></div><div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="comment"> * </span></div><div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="comment"> * Performs the LU factorization without pivoting of a general M-by-N matrix A:</span></div><div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="comment"> *</span></div><div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="comment"> * A = L * U</span></div><div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="comment"> *</span></div><div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="comment"> * where L is a lower triangular ( lower trapezoidal if M > N ) matrix with</span></div><div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="comment"> * unit diagonal elements and U is an upper triangular ( upper trapezoidal </span></div><div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="comment"> * if M < N ) matrix.</span></div><div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="comment"> *</span></div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="comment">**/</span></div><div class="line"><a name="l00033"></a><span class="lineno"> 33</span> <span class="comment"></span></div><div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment">/**</span></div><div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="comment"> * </span></div><div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="comment"> * @param[in]</span></div><div class="line"><a name="l00037"></a><span class="lineno"> 37</span> <span class="comment"> * M int.</span></div><div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="comment"> * M specifies the number of rows of the matrix A. M >= 0.</span></div><div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment"> * </span></div><div class="line"><a name="l00040"></a><span class="lineno"> 40</span> <span class="comment"> * @param[in]</span></div><div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="comment"> * N int.</span></div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span> <span class="comment"> * N specifies the number of columns of the matrix A. N >= 0.</span></div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="comment"> *</span></div><div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="comment"> * @param[in,out]</span></div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span> <span class="comment"> * A double *.</span></div><div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="comment"> * A is a pointer to a regular matrix of dimension M-by-N.</span></div><div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="comment"> * On exit, if return value is Success, the matrix A is overwriten by </span></div><div class="line"><a name="l00048"></a><span class="lineno"> 48</span> <span class="comment"> * the factors L and U. The unit diagonal elements of L are not stored.</span></div><div class="line"><a name="l00049"></a><span class="lineno"> 49</span> <span class="comment"> * </span></div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span> <span class="comment"> * @param[in]</span></div><div class="line"><a name="l00051"></a><span class="lineno"> 51</span> <span class="comment"> * LDA int.</span></div><div class="line"><a name="l00052"></a><span class="lineno"> 52</span> <span class="comment"> * LDA specifies the number of columns of A ( row-major order ).</span></div><div class="line"><a name="l00053"></a><span class="lineno"> 53</span> <span class="comment"> * LDA must be at least max( 1, N ).</span></div><div class="line"><a name="l00054"></a><span class="lineno"> 54</span> <span class="comment"> * </span></div><div class="line"><a name="l00055"></a><span class="lineno"> 55</span> <span class="comment"> **/</span></div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span> <span class="comment"></span></div><div class="line"><a name="l00057"></a><span class="lineno"> 57</span> <span class="comment">/**</span></div><div class="line"><a name="l00058"></a><span class="lineno"> 58</span> <span class="comment"> *</span></div><div class="line"><a name="l00059"></a><span class="lineno"> 59</span> <span class="comment"> * @retval Success successful exit</span></div><div class="line"><a name="l00060"></a><span class="lineno"> 60</span> <span class="comment"> * @retval NoSuccess unsuccessful exit</span></div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span> <span class="comment"> *</span></div><div class="line"><a name="l00062"></a><span class="lineno"> 62</span> <span class="comment"> **/</span></div><div class="line"><a name="l00063"></a><span class="lineno"> 63</span> <span class="comment"></span></div><div class="line"><a name="l00064"></a><span class="lineno"> 64</span> <span class="comment">/**</span></div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span> <span class="comment"> *</span></div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span> <span class="comment"> * @sa knpdgetrf</span></div><div class="line"><a name="l00067"></a><span class="lineno"> 67</span> <span class="comment"> *</span></div><div class="line"><a name="l00068"></a><span class="lineno"> 68</span> <span class="comment"> **/</span></div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span> </div><div class="line"><a name="l00070"></a><span class="lineno"><a class="line" href="lass_8h.html#a665c631275c07b1dc9554d982e493a5e"> 70</a></span> <span class="keywordtype">int</span> <a class="code" href="ddss__dnpgetrf_8c.html#a665c631275c07b1dc9554d982e493a5e">ddss_dnpgetrf</a>( <span class="keywordtype">int</span> M, <span class="keywordtype">int</span> N, <span class="keywordtype">double</span> *A, <span class="keywordtype">int</span> LDA )</div><div class="line"><a name="l00071"></a><span class="lineno"> 71</span> {</div><div class="line"><a name="l00072"></a><span class="lineno"> 72</span> </div><div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  <span class="comment">// Argument checking</span></div><div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  <span class="keywordflow">if</span> ( M < 0 )</div><div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  {</div><div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  fprintf( stderr, <span class="stringliteral">"Illegal value of M, in ddss_dnpgetrf code\n"</span> );</div><div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  <span class="keywordflow">return</span> NoSuccess;</div><div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  }</div><div class="line"><a name="l00079"></a><span class="lineno"> 79</span> </div><div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  <span class="keywordflow">if</span> ( N < 0 )</div><div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  {</div><div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  fprintf( stderr, <span class="stringliteral">"Illegal value of N, in ddss_dnpgetrf code\n"</span> );</div><div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <span class="keywordflow">return</span> NoSuccess;</div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  }</div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span> </div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <span class="keywordflow">if</span> ( LDA < MAX( 1, N ) )</div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  {</div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  fprintf( stderr, <span class="stringliteral">"Illegal value of LDA, in ddss_dnpgetrf code\n"</span> );</div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  <span class="keywordflow">return</span> NoSuccess;</div><div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  }</div><div class="line"><a name="l00091"></a><span class="lineno"> 91</span> </div><div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  <span class="comment">// Qick return</span></div><div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <span class="keywordflow">if</span> ( MIN( M, N ) == 0 )</div><div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  {</div><div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  <span class="keywordflow">return</span> Success;</div><div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  }</div><div class="line"><a name="l00097"></a><span class="lineno"> 97</span> </div><div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  <span class="keywordflow">return</span> <a class="code" href="kdnpgetrf_8c.html#aad8c6b9542e74c481faff1ff5e03850f">kdnpgetrf</a>( M, N, A, LDA );</div><div class="line"><a name="l00099"></a><span class="lineno"> 99</span> </div><div class="line"><a name="l00100"></a><span class="lineno"> 100</span> } </div><div class="ttc" id="ddss__dnpgetrf_8c_html_a665c631275c07b1dc9554d982e493a5e"><div class="ttname"><a href="ddss__dnpgetrf_8c.html#a665c631275c07b1dc9554d982e493a5e">ddss_dnpgetrf</a></div><div class="ttdeci">int ddss_dnpgetrf(int M, int N, double *A, int LDA)</div><div class="ttdef"><b>Definition:</b> <a href="ddss__dnpgetrf_8c_source.html#l00070">ddss_dnpgetrf.c:70</a></div></div>
<div class="ttc" id="kdnpgetrf_8c_html_aad8c6b9542e74c481faff1ff5e03850f"><div class="ttname"><a href="kdnpgetrf_8c.html#aad8c6b9542e74c481faff1ff5e03850f">kdnpgetrf</a></div><div class="ttdeci">enum LASS_RETURN kdnpgetrf(int M, int N, double *A, int LDA)</div><div class="ttdef"><b>Definition:</b> <a href="kdnpgetrf_8c_source.html#l00074">kdnpgetrf.c:74</a></div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated by  <a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/>
</a> 1.8.11
</small></address>
</body>
</html>