[petsc-users] METIS 5 in PETSc 3, 3 [not directly related to PETSc]

[Tabrez Ali]

I was passing vwgt(:)=>null() instead of vwgt=>null() and so on. The 
following seems to work fine on all three compilers.

program test
  implicit none
  integer, parameter   :: nels=2, nnds=6, npel=4
  integer              :: eptr(nels+1), nodes(nels*npel), epart(nels), 
npart(nnds), n
  integer, pointer     :: vwgt=>null(), vsize=>null(), mopts=>null()
  real(8), pointer     :: tpwgts=>null()
  eptr=(/0,4,8/)
  nodes=(/0,1,2,3,1,4,5,2/)
  call 
METIS_PartMeshNodal(nels,nnds,eptr,nodes,vwgt,vsize,2,tpwgts,mopts,n,epart,npart)
  print*, npart; print*, epart
end program test

Thanks for all help/suggestions.

T

On 06/08/2012 11:04 AM, Sean Farley wrote:
>> Here it is for the 3 compilers.
>>
>> GNU FC http://bpaste.net/show/31029/
>>
>> Intel FC http://bpaste.net/show/QJ6k5jc52t7MGUlvrlEV/
>>
>> PGI FC http://bpaste.net/show/DZzY6Re5iXMgLSHydOgJ/
> You're going to need to build metis 5 with debugging options and step
> into the function to see why it's segfaulting. It might help (though
> it is a different function call) to see my mumps patch to interface
> with metis 5:
>
> http://petsc.cs.iit.edu/petsc/externalpackages/MUMPS_4.10.0/rev/907457aa0d14