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

[Tabrez Ali]

Thanks for your answer. Unfortunately its not that.

Based on past experience I am sure I am overlooking something very 
simple but I cant seem to find out what. Btw METIS 4 worked fine before.

On 06/08/2012 08:39 AM, John Mousel wrote:
> 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 
> <mailto: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
>
>     Thanks in advance.
>
>     Tabrez
>

-- 
No one trusts a model except the one who wrote it; Everyone trusts an observation except the one who made it- Harlow Shapley

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120608/360c0bac/attachment.html>