[petsc-users] MatMatSolve in sequential call

[Hong Zhang]

Alexander:

>
>
> Sorry, but I forgot to mention that this is double complex symmetric
> matrix. Did you load and treat it like that?
>
I'm not aware it, thus I used real precision :-(

>
>
>  I'm convinced that matrix A is corrupted or obtained from an incorrect
> numerical model.
> Also note, your matrix A has size of approx. 60k, which is too large for
> sequential direct solver.
>
>
> MUMPS takes only several minutes and 6 GB of memory to factorize it.
> This factorization gives residual on the order of 10e-12 and solution is
> indeed correct.
>
> Nevertheless, you're right, there is numerical null-space in this matrix
> since it comes
> from the discretization of equation that contains curl curl operator, but
> practically this
> case is not really the worst one.
>
>
>  The error on MatMatSolve() when calling
> ' -pc_type cholesky -pc_factor_mat_solver_package mumps'
> is a bug in our petsc-mumps interface. I've pushed a fix to petsc-3.3
> http://petsc.cs.iit.edu/petsc/releases/petsc-3.3/rev/8badc49a596e
>
>  Thanks for your report.
>
>
> Thank you for the patch.
>
> --
> Regards,
> Alexander
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120821/2a7976e7/attachment.html>