2020-06-28 23:07:29 +08:00
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://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" / >
2021-10-26 03:34:27 +08:00
< meta http-equiv = "X-UA-Compatible" content = "IE=11" / >
2022-06-08 02:13:03 +08:00
< meta name = "generator" content = "Doxygen 1.9.4" / >
2020-06-28 23:07:29 +08:00
< meta name = "viewport" content = "width=device-width, initial-scale=1" / >
< title > Algorithms_in_C: numerical_methods/qr_decomposition.c File Reference< / 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 = "../../navtree.css" rel = "stylesheet" type = "text/css" / >
< script type = "text/javascript" src = "../../resize.js" > < / script >
< script type = "text/javascript" src = "../../navtreedata.js" > < / script >
< script type = "text/javascript" src = "../../navtree.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/x-mathjax-config" >
2021-10-26 03:34:27 +08:00
MathJax.Hub.Config({
extensions: ["tex2jax.js", "TeX/AMSmath.js", "TeX/AMSsymbols.js"],
jax: ["input/TeX","output/HTML-CSS"],
2020-06-28 23:07:29 +08:00
});
< / script >
< script type = "text/javascript" async = "async" src = "https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.7/MathJax.js?config=TeX-MML-AM_CHTML/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 >
2022-01-15 04:11:27 +08:00
< tr id = "projectrow" >
< td id = "projectalign" >
2021-10-26 03:34:27 +08:00
< div id = "projectname" > Algorithms_in_C< span id = "projectnumber" >   1.0.0< / span >
2020-06-28 23:07:29 +08:00
< / div >
< div id = "projectbrief" > Set of algorithms implemented in C.< / div >
< / td >
< / tr >
< / tbody >
< / table >
< / div >
<!-- end header part -->
2022-06-08 02:13:03 +08:00
<!-- Generated by Doxygen 1.9.4 -->
2020-06-28 23:07:29 +08:00
< script type = "text/javascript" >
2021-10-26 03:34:27 +08:00
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699& dn=expat.txt MIT */
var searchBox = new SearchBox("searchBox", "../../search",'Search','.html');
2020-06-28 23:07:29 +08:00
/* @license-end */
< / script >
< script type = "text/javascript" src = "../../menudata.js" > < / script >
< script type = "text/javascript" src = "../../menu.js" > < / script >
< script type = "text/javascript" >
2021-10-26 03:34:27 +08:00
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699& dn=expat.txt MIT */
2020-06-28 23:07:29 +08:00
$(function() {
initMenu('../../',true,false,'search.php','Search');
$(document).ready(function() { init_search(); });
});
2021-10-26 03:34:27 +08:00
/* @license-end */
< / script >
2020-06-28 23:07:29 +08:00
< div id = "main-nav" > < / div >
< / div > <!-- top -->
< div id = "side-nav" class = "ui-resizable side-nav-resizable" >
< div id = "nav-tree" >
< div id = "nav-tree-contents" >
< div id = "nav-sync" class = "sync" > < / div >
< / div >
< / div >
< div id = "splitbar" style = "-moz-user-select:none;"
class="ui-resizable-handle">
< / div >
< / div >
< script type = "text/javascript" >
2021-10-26 03:34:27 +08:00
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699& dn=expat.txt MIT */
2020-06-28 23:07:29 +08:00
$(document).ready(function(){initNavTree('d5/d23/qr__decomposition_8c.html','../../'); initResizable(); });
/* @license-end */
< / script >
< div id = "doc-content" >
<!-- 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 class = "header" >
< div class = "summary" >
< a href = "#func-members" > Functions< / a > < / div >
2021-10-26 03:34:27 +08:00
< div class = "headertitle" > < div class = "title" > qr_decomposition.c File Reference< / div > < / div >
2020-06-28 23:07:29 +08:00
< / div > <!-- header -->
< div class = "contents" >
2021-10-26 03:34:27 +08:00
< p > Program to compute the < a href = "https://en.wikipedia.org/wiki/QR_decomposition" target = "_blank" > QR decomposition< / a > of a given matrix.
2020-06-28 23:07:29 +08:00
< a href = "#details" > More...< / a > < / p >
2020-06-28 23:30:37 +08:00
< div class = "textblock" > < code > #include < math.h> < / code > < br / >
2020-06-28 23:07:29 +08:00
< code > #include < stdio.h> < / code > < br / >
< code > #include < stdlib.h> < / code > < br / >
< code > #include < time.h> < / code > < br / >
2020-06-28 23:30:37 +08:00
< code > #include " < a class = "el" href = "../../d4/d68/qr__decompose_8h_source.html" > qr_decompose.h< / a > " < / code > < br / >
2020-06-28 23:07:29 +08:00
< / div > < div class = "textblock" > < div class = "dynheader" >
Include dependency graph for qr_decomposition.c:< / div >
< div class = "dyncontent" >
2022-06-08 02:13:03 +08:00
< div class = "center" > < iframe scrolling = "no" frameborder = "0" src = "../../db/d45/qr__decomposition_8c__incl.svg" width = "395" height = "202" > < p > < b > This browser is not able to show SVG: try Firefox, Chrome, Safari, or Opera instead.< / b > < / p > < / iframe >
2020-06-28 23:07:29 +08:00
< / div >
< / div >
< / div > < table class = "memberdecls" >
2021-10-26 03:34:27 +08:00
< tr class = "heading" > < td colspan = "2" > < h2 class = "groupheader" > < a id = "func-members" name = "func-members" > < / a >
2020-06-28 23:07:29 +08:00
Functions< / h2 > < / td > < / tr >
2021-10-26 03:34:27 +08:00
< tr class = "memitem:a840291bc02cba5474a4cb46a9b9566fe" > < td class = "memItemLeft" align = "right" valign = "top" > int  < / td > < td class = "memItemRight" valign = "bottom" > < a class = "el" href = "../../d5/d23/qr__decomposition_8c.html#a840291bc02cba5474a4cb46a9b9566fe" > main< / a > (void)< / td > < / tr >
< tr class = "memdesc:a840291bc02cba5474a4cb46a9b9566fe" > < td class = "mdescLeft" >   < / td > < td class = "mdescRight" > main function < a href = "../../d5/d23/qr__decomposition_8c.html#a840291bc02cba5474a4cb46a9b9566fe" > More...< / a > < br / > < / td > < / tr >
2020-06-28 23:07:29 +08:00
< tr class = "separator:a840291bc02cba5474a4cb46a9b9566fe" > < td class = "memSeparator" colspan = "2" >   < / td > < / tr >
< / table >
< a name = "details" id = "details" > < / a > < h2 class = "groupheader" > Detailed Description< / h2 >
2021-10-26 03:34:27 +08:00
< div class = "textblock" > < p > Program to compute the < a href = "https://en.wikipedia.org/wiki/QR_decomposition" target = "_blank" > QR decomposition< / a > of a given matrix. < / p >
< dl class = "section author" > < dt > Author< / dt > < dd > < a href = "https://github.com/kvedala" target = "_blank" > Krishna Vedala< / a > < / dd > < / dl >
< / div > < h2 class = "groupheader" > Function Documentation< / h2 >
< a id = "a840291bc02cba5474a4cb46a9b9566fe" name = "a840291bc02cba5474a4cb46a9b9566fe" > < / a >
< h2 class = "memtitle" > < span class = "permalink" > < a href = "#a840291bc02cba5474a4cb46a9b9566fe" > ◆ < / a > < / span > main()< / h2 >
< div class = "memitem" >
< div class = "memproto" >
< table class = "memname" >
< tr >
< td class = "memname" > int main < / td >
< td > (< / td >
< td class = "paramtype" > void  < / td >
< td class = "paramname" > < / td > < td > )< / td >
< td > < / td >
< / tr >
< / table >
< / div > < div class = "memdoc" >
< p > main function < / p >
2022-01-15 04:11:27 +08:00
< div class = "fragment" > < div class = "line" > < span class = "lineno" > 19< / span > {< / div >
< div class = "line" > < span class = "lineno" > 20< / span > < span class = "keywordtype" > double< / span > **A;< / div >
< div class = "line" > < span class = "lineno" > 21< / span > < span class = "keywordtype" > unsigned< / span > < span class = "keywordtype" > int< / span > ROWS, COLUMNS;< / div >
< div class = "line" > < span class = "lineno" > 22< / span > < / div >
< div class = "line" > < span class = "lineno" > 23< / span > printf(< span class = "stringliteral" > " Enter the number of rows and columns: " < / span > );< / div >
< div class = "line" > < span class = "lineno" > 24< / span > scanf(< span class = "stringliteral" > " %u %u" < / span > , & ROWS, & COLUMNS);< / div >
< div class = "line" > < span class = "lineno" > 25< / span > < span class = "keywordflow" > if< / span > (ROWS < COLUMNS)< / div >
< div class = "line" > < span class = "lineno" > 26< / span > {< / div >
< div class = "line" > < span class = "lineno" > 27< / span > fprintf(stderr,< / div >
< div class = "line" > < span class = "lineno" > 28< / span > < span class = "stringliteral" > " Number of rows must be greater than or equal to " < / span > < / div >
< div class = "line" > < span class = "lineno" > 29< / span > < span class = "stringliteral" > " number of columns.\n" < / span > );< / div >
< div class = "line" > < span class = "lineno" > 30< / span > < span class = "keywordflow" > return< / span > -1;< / div >
< div class = "line" > < span class = "lineno" > 31< / span > }< / div >
< div class = "line" > < span class = "lineno" > 32< / span > < / div >
< div class = "line" > < span class = "lineno" > 33< / span > printf(< span class = "stringliteral" > " Enter matrix elements row-wise:\n" < / span > );< / div >
< div class = "line" > < span class = "lineno" > 34< / span > < / div >
< div class = "line" > < span class = "lineno" > 35< / span > A = (< span class = "keywordtype" > double< / span > **)< a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a725f50ecaf1959d96de79b36b4788fee" > malloc< / a > (ROWS * < span class = "keyword" > sizeof< / span > (< span class = "keywordtype" > double< / span > *));< / div >
< div class = "line" > < span class = "lineno" > 36< / span > < span class = "keywordflow" > for< / span > (< span class = "keywordtype" > int< / span > i = 0; i < ROWS; i++)< / div >
< div class = "line" > < span class = "lineno" > 37< / span > A[i] = (< span class = "keywordtype" > double< / span > *)< a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a725f50ecaf1959d96de79b36b4788fee" > malloc< / a > (COLUMNS * < span class = "keyword" > sizeof< / span > (< span class = "keywordtype" > double< / span > ));< / div >
< div class = "line" > < span class = "lineno" > 38< / span > < / div >
< div class = "line" > < span class = "lineno" > 39< / span > < span class = "keywordflow" > for< / span > (< span class = "keywordtype" > int< / span > i = 0; i < ROWS; i++)< / div >
< div class = "line" > < span class = "lineno" > 40< / span > < span class = "keywordflow" > for< / span > (< span class = "keywordtype" > int< / span > j = 0; j < COLUMNS; j++) scanf(< span class = "stringliteral" > " %lf" < / span > , & A[i][j]);< / div >
< div class = "line" > < span class = "lineno" > 41< / span > < / div >
< div class = "line" > < span class = "lineno" > 42< / span > < a class = "code hl_function" href = "../../d4/d68/qr__decompose_8h.html#a90562ce8c3707401e9c5809dece68d6a" > print_matrix< / a > (A, ROWS, COLUMNS);< / div >
< div class = "line" > < span class = "lineno" > 43< / span > < / div >
< div class = "line" > < span class = "lineno" > 44< / span > < span class = "keywordtype" > double< / span > **R = (< span class = "keywordtype" > double< / span > **)< a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a725f50ecaf1959d96de79b36b4788fee" > malloc< / a > (< span class = "keyword" > sizeof< / span > (< span class = "keywordtype" > double< / span > *) * ROWS);< / div >
< div class = "line" > < span class = "lineno" > 45< / span > < span class = "keywordtype" > double< / span > **Q = (< span class = "keywordtype" > double< / span > **)< a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a725f50ecaf1959d96de79b36b4788fee" > malloc< / a > (< span class = "keyword" > sizeof< / span > (< span class = "keywordtype" > double< / span > *) * ROWS);< / div >
< div class = "line" > < span class = "lineno" > 46< / span > < span class = "keywordflow" > if< / span > (!Q || !R)< / div >
< div class = "line" > < span class = "lineno" > 47< / span > {< / div >
< div class = "line" > < span class = "lineno" > 48< / span > perror(< span class = "stringliteral" > " Unable to allocate memory for Q & R!" < / span > );< / div >
< div class = "line" > < span class = "lineno" > 49< / span > < span class = "keywordflow" > return< / span > -1;< / div >
< div class = "line" > < span class = "lineno" > 50< / span > }< / div >
< div class = "line" > < span class = "lineno" > 51< / span > < span class = "keywordflow" > for< / span > (< span class = "keywordtype" > int< / span > i = 0; i < ROWS; i++)< / div >
< div class = "line" > < span class = "lineno" > 52< / span > {< / div >
< div class = "line" > < span class = "lineno" > 53< / span > R[i] = (< span class = "keywordtype" > double< / span > *)< a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a725f50ecaf1959d96de79b36b4788fee" > malloc< / a > (< span class = "keyword" > sizeof< / span > (< span class = "keywordtype" > double< / span > ) * COLUMNS);< / div >
< div class = "line" > < span class = "lineno" > 54< / span > Q[i] = (< span class = "keywordtype" > double< / span > *)< a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a725f50ecaf1959d96de79b36b4788fee" > malloc< / a > (< span class = "keyword" > sizeof< / span > (< span class = "keywordtype" > double< / span > ) * ROWS);< / div >
< div class = "line" > < span class = "lineno" > 55< / span > < span class = "keywordflow" > if< / span > (!Q[i] || !R[i])< / div >
< div class = "line" > < span class = "lineno" > 56< / span > {< / div >
< div class = "line" > < span class = "lineno" > 57< / span > perror(< span class = "stringliteral" > " Unable to allocate memory for Q & R." < / span > );< / div >
< div class = "line" > < span class = "lineno" > 58< / span > < span class = "keywordflow" > return< / span > -1;< / div >
< div class = "line" > < span class = "lineno" > 59< / span > }< / div >
< div class = "line" > < span class = "lineno" > 60< / span > }< / div >
< div class = "line" > < span class = "lineno" > 61< / span > < / div >
< div class = "line" > < span class = "lineno" > 62< / span > clock_t t1 = clock();< / div >
< div class = "line" > < span class = "lineno" > 63< / span > < a class = "code hl_function" href = "../../d4/d68/qr__decompose_8h.html#a45c7640d9d22c89c11beb1f567843c56" > qr_decompose< / a > (A, Q, R, ROWS, COLUMNS);< / div >
< div class = "line" > < span class = "lineno" > 64< / span > < span class = "keywordtype" > double< / span > dtime = (double)(clock() - t1) / CLOCKS_PER_SEC;< / div >
< div class = "line" > < span class = "lineno" > 65< / span > < / div >
< div class = "line" > < span class = "lineno" > 66< / span > < a class = "code hl_function" href = "../../d4/d68/qr__decompose_8h.html#a90562ce8c3707401e9c5809dece68d6a" > print_matrix< / a > (R, ROWS, COLUMNS);< / div >
< div class = "line" > < span class = "lineno" > 67< / span > < a class = "code hl_function" href = "../../d4/d68/qr__decompose_8h.html#a90562ce8c3707401e9c5809dece68d6a" > print_matrix< / a > (Q, ROWS, COLUMNS);< / div >
< div class = "line" > < span class = "lineno" > 68< / span > printf(< span class = "stringliteral" > " Time taken to compute: %.4g sec\n" < / span > , dtime);< / div >
< div class = "line" > < span class = "lineno" > 69< / span > < / div >
< div class = "line" > < span class = "lineno" > 70< / span > < span class = "keywordflow" > for< / span > (< span class = "keywordtype" > int< / span > i = 0; i < ROWS; i++)< / div >
< div class = "line" > < span class = "lineno" > 71< / span > {< / div >
< div class = "line" > < span class = "lineno" > 72< / span > < a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a9cc854374299a1dd933bf62029761768" > free< / a > (A[i]);< / div >
< div class = "line" > < span class = "lineno" > 73< / span > < a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a9cc854374299a1dd933bf62029761768" > free< / a > (R[i]);< / div >
< div class = "line" > < span class = "lineno" > 74< / span > < a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a9cc854374299a1dd933bf62029761768" > free< / a > (Q[i]);< / div >
< div class = "line" > < span class = "lineno" > 75< / span > }< / div >
< div class = "line" > < span class = "lineno" > 76< / span > < a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a9cc854374299a1dd933bf62029761768" > free< / a > (A);< / div >
< div class = "line" > < span class = "lineno" > 77< / span > < a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a9cc854374299a1dd933bf62029761768" > free< / a > (R);< / div >
< div class = "line" > < span class = "lineno" > 78< / span > < a class = "code hl_define" href = "../../d2/ddd/malloc__dbg_8h.html#a9cc854374299a1dd933bf62029761768" > free< / a > (Q);< / div >
< div class = "line" > < span class = "lineno" > 79< / span > < span class = "keywordflow" > return< / span > 0;< / div >
< div class = "line" > < span class = "lineno" > 80< / span > }< / div >
2021-10-26 03:34:27 +08:00
< div class = "ttc" id = "amalloc__dbg_8h_html_a725f50ecaf1959d96de79b36b4788fee" > < div class = "ttname" > < a href = "../../d2/ddd/malloc__dbg_8h.html#a725f50ecaf1959d96de79b36b4788fee" > malloc< / a > < / div > < div class = "ttdeci" > #define malloc(bytes)< / div > < div class = "ttdoc" > This macro replace the standard malloc function with malloc_dbg.< / div > < div class = "ttdef" > < b > Definition:< / b > malloc_dbg.h:18< / div > < / div >
< div class = "ttc" id = "amalloc__dbg_8h_html_a9cc854374299a1dd933bf62029761768" > < div class = "ttname" > < a href = "../../d2/ddd/malloc__dbg_8h.html#a9cc854374299a1dd933bf62029761768" > free< / a > < / div > < div class = "ttdeci" > #define free(ptr)< / div > < div class = "ttdoc" > This macro replace the standard free function with free_dbg.< / div > < div class = "ttdef" > < b > Definition:< / b > malloc_dbg.h:26< / div > < / div >
< div class = "ttc" id = "aqr__decompose_8h_html_a45c7640d9d22c89c11beb1f567843c56" > < div class = "ttname" > < a href = "../../d4/d68/qr__decompose_8h.html#a45c7640d9d22c89c11beb1f567843c56" > qr_decompose< / a > < / div > < div class = "ttdeci" > void qr_decompose(double **A, double **Q, double **R, int M, int N)< / div > < div class = "ttdoc" > Decompose matrix using Gram-Schmidt process.< / div > < div class = "ttdef" > < b > Definition:< / b > qr_decompose.h:142< / div > < / div >
< div class = "ttc" id = "aqr__decompose_8h_html_a90562ce8c3707401e9c5809dece68d6a" > < div class = "ttname" > < a href = "../../d4/d68/qr__decompose_8h.html#a90562ce8c3707401e9c5809dece68d6a" > print_matrix< / a > < / div > < div class = "ttdeci" > void print_matrix(double **A, int M, int N)< / div > < div class = "ttdoc" > function to display matrix on stdout< / div > < div class = "ttdef" > < b > Definition:< / b > qr_decompose.h:22< / div > < / div >
< / div > <!-- fragment --> < div class = "dynheader" >
Here is the call graph for this function:< / div >
< div class = "dyncontent" >
2021-12-10 01:57:52 +08:00
< div class = "center" > < iframe scrolling = "no" frameborder = "0" src = "../../d5/d23/qr__decomposition_8c_a840291bc02cba5474a4cb46a9b9566fe_cgraph.svg" width = "486" height = "164" > < p > < b > This browser is not able to show SVG: try Firefox, Chrome, Safari, or Opera instead.< / b > < / p > < / iframe >
2021-10-26 03:34:27 +08:00
< / div >
< / div >
< / div >
< / div >
< / div > <!-- contents -->
2020-06-28 23:07:29 +08:00
< / div > <!-- doc - content -->
<!-- start footer part -->
< div id = "nav-path" class = "navpath" > <!-- id is needed for treeview function! -->
< ul >
< li class = "navelem" > < a class = "el" href = "../../dir_9c6faab82c22511b50177aa2e38e2780.html" > numerical_methods< / a > < / li > < li class = "navelem" > < a class = "el" href = "../../d5/d23/qr__decomposition_8c.html" > qr_decomposition.c< / a > < / li >
2022-06-08 02:13:03 +08:00
< li class = "footer" > Generated by < a href = "https://www.doxygen.org/index.html" > < img class = "footer" src = "../../doxygen.svg" width = "104" height = "31" alt = "doxygen" / > < / a > 1.9.4 < / li >
2020-06-28 23:07:29 +08:00
< / ul >
< / div >
< / body >
< / html >
