1 <?xml version="1.0" encoding="UTF-8"?>
2 <refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" xmlns:scilab="http://www.scilab.org" xml:lang="en" xml:id="umf_lusolve">
4 <refname>umf_lusolve</refname>
5 <refpurpose>solve a linear sparse system given the LU factors </refpurpose>
9 <synopsis>[x] = umf_lusolve(LU_ptr, b [, st, A])</synopsis>
12 <title>Arguments</title>
17 <para>a pointer to umf lu factors (L,U,p,q,R) </para>
23 <para>a real or complex column vector or a matrix (multiple rhs)</para>
29 <para>(optional) a string "Ax=b" (default) or "Ax'=b"
30 (to be written "Ax''=b" in scilab language: a quote in a string
38 <para>(optional) the sparse square matrix corresponding to the LU factors (LU_ptr must be got with LU_ptr = umf_lufact(A))</para>
44 <para>a column vector or a matrix in case of multiple rhs ( x(:,i) is solution of A x(:,i) = b(:,i) or A'x(:,i) = b(:,i) )</para>
50 <title>Description</title>
52 This function must be used in conjunction with <link linkend="umf_lufact">umf_lufact</link> which
53 computes the LU factors of a sparse matrix. The optional
54 <literal>st</literal> argument lets us choose between the solving of Ax=b (general case)
55 or of A'x=b (sometimes useful). If you give the 4th argument then
56 iterative refinement will be also proceeded (as in umfpack) to
57 give a better numerical solution.
61 <title>Examples</title>
63 see the example section of <link linkend="umf_lufact">umf_lufact</link>
66 <refsection role="see also">
67 <title>See Also</title>
68 <simplelist type="inline">
70 <link linkend="umfpack">umfpack</link>
73 <link linkend="umf_lufact">umf_lufact</link>
76 <link linkend="umf_luget">umf_luget</link>
79 <link linkend="umf_ludel">umf_ludel</link>
82 <link linkend="umf_luinfo">umf_luinfo</link>