[petsc-users] SuperLU_dist issue in 3.7.4

[Anton Popov]

Hong,

I checked out & compiled your new branch: 
hzhang/fix-superlu_dist-reuse-factornumeric. Unfortunately it did not 
solve the problem.

Sorry.

On 11/21/2016 04:43 AM, Hong wrote:
> Anton,
> I pushed a fix
> https://bitbucket.org/petsc/petsc/commits/28865de08051eb99557d70672c208e14da23c8b1
> in branch hzhang/fix-superlu_dist-reuse-factornumeric.
> Can you give it a try to see if it works?
> I do not have an example which produces your problem.
>
> In your email, you asked "Setting Options.Fact = DOFACT for all 
> factorizations is currently impossible via PETSc interface.
> The user is expected to choose some kind of reuse model.
> If you could add it, I (and other users probably too) would really 
> appreciate that."
>
> We do not allow user set superlu' Options.Fact = DOFACT. If user 
> changes matrix structure, then user must call
> KSPSetOperators() -> call symbolic matrix factorization again, in 
> which we set Options.Fact = DOFACT.

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.

I would still appreciate the full factorization even for the same 
pattern without destroying KSP/PC object (just as a custom option).

>
> Hong

Thanks a lot,
Anton
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161126/b6d79f2b/attachment.html>