[petsc-dev] Problem with MUMPS+metis

[John Fettig]

On Mon, Feb 20, 2012 at 3:21 PM, Sean Farley <sean at mcs.anl.gov> wrote:

> This hangs for me in the METIS ordering stage
>> (src/ksp/ksp/examples/tutorials/ex2):
>>
>> mpirun -np 1 ./ex2 -n 100 -m 100 -ksp_type preonly -pc_type lu
>> -pc_factor_mat_solver_package mumps -mat_mumps_icntl_7 5 -info
>
>
> Thanks, I'll look into this.
>

It seems to be a problem with metis' interpretation of the vertex weights.
If you set the vwgt=0x0 in a debugger when stopped at metis_nodend_ it
completes successfully.  Also, if you set the weights to 1 inside
dmumps_part2.F like:

         DO I=1,NCMP
           FILS(I) = 1
         ENDDO

it also runs fine.  I think having zero vertex weights is probably wrong
anyways.

Regards,
John
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120220/aa7441c2/attachment.html>