It's hard to tell from the info you provided, but you seem to be playing fast and loose with your type declarations. METIS is expecting real_t, which is a 32 bit real if you haven't changed the definition in metis.h. I know this has caused me problems in the past.
On Fri, Jun 8, 2012 at 7:38 AM, Tabrez Ali <stali at geology.wisc.edu> wrote: > Sorry about a question not directly related to PETSc but has anyone here > been able to use the METIS 5.0 (that PETSc 3.3/dev downloads/builds) with > Fortran? There has been an API change from 4 to 5 but I am having some > trouble and METIS manual/forums havent been useful. > > For example consider the simple code (below) that partitions a two element > mesh made of linear quads into two. The elements are numbered 0 1 2 3 and 1 > 4 5 2.It works fine with GNU FC (no valgrind errors). With Intel FC it > works fines (though valgrind throws a bunch of errors). However with PGI > compilers I get a segfault. > > 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,7/) > 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 > > According to the manual moving from METIS 4 to 5 only involves passing > some additional nulls. I am not sure what I missed. > > http://glaros.dtc.umn.edu/**gkhome/fetch/sw/metis/manual.**pdf<http://glaros.dtc.umn.edu/gkhome/fetch/sw/metis/manual.pdf> > > Thanks in advance. > > Tabrez > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120608/ad6dcd5f/attachment.html>
