<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Hong,</p>
<p>I checked out & compiled your new branch:
hzhang/fix-superlu_dist-reuse-factornumeric. Unfortunately it did
not solve the problem.<br>
</p>
Sorry.<br>
<br>
<div class="moz-cite-prefix">On 11/21/2016 04:43 AM, Hong wrote:<br>
</div>
<blockquote
cite="mid:CAGCphBu1UrhURvd0KW74jVxPX9D0ezEonPJ_Omp76XafW7RrHQ@mail.gmail.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<div dir="ltr">Anton,
<div>I pushed a fix </div>
<div><a moz-do-not-send="true"
href="https://bitbucket.org/petsc/petsc/commits/28865de08051eb99557d70672c208e14da23c8b1">https://bitbucket.org/petsc/petsc/commits/28865de08051eb99557d70672c208e14da23c8b1</a></div>
<div>in branch hzhang/fix-superlu_dist-reuse-factornumeric.</div>
<div>Can you give it a try to see if it works?</div>
<div>I do not have an example which produces your problem.</div>
<div><br>
</div>
<div><span style="font-size:12.8px">In your email, you asked
"Setting Options.Fact = DOFACT for all factorizations is
currently impossible via PETSc interface.</span><br
style="font-size:12.8px">
<span style="font-size:12.8px">The user is expected to choose
some kind of reuse model.</span><br style="font-size:12.8px">
<span style="font-size:12.8px">If you could add it, I (and
other users probably too) would really appreciate that."</span><br>
</div>
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">We do not allow user set
superlu' </span><span style="font-size:12.8px">Options.Fact
= DOFACT. If user changes matrix structure, then user must
call</span></div>
<div><span style="font-size:12.8px">KSPSetOperators() -> call
symbolic matrix factorization again, in which we set </span><span
style="font-size:12.8px">Options.Fact = DOFACT.</span></div>
</div>
</blockquote>
<br>
I have a conceptual question. How can sparsity (column) permutation
be reused if it's applied on top of the equilibration (row)
permutation? Symbolic factorization should be repeated anyway. Does
it run in some kind of faster update mode in this case? Please
correct me if I misunderstand something.<br>
<br>
I would still appreciate the full factorization even for the same
pattern without destroying KSP/PC object (just as a custom option).<br>
<br>
<blockquote
cite="mid:CAGCphBu1UrhURvd0KW74jVxPX9D0ezEonPJ_Omp76XafW7RrHQ@mail.gmail.com"
type="cite">
<div dir="ltr">
<div><span style="font-size:12.8px"><br>
</span></div>
<div><span style="font-size:12.8px">Hong</span></div>
</div>
</blockquote>
<br>
Thanks a lot,<br>
Anton<br>
</body>
</html>