[Nek5000-users] Seg fault (may be from Hpts routine)

[nek5000-users at lists.mcs.anl.gov]

Hi Stefan,

the latest revision (519) fixed the seg fault for me, but I keep getting error
messages about the history points not being within the mesh even though they
are.
WARNING: point not within mesh xy[z]: ! -0.2000000E+00  0.3000000E+00

Do you have the same issue?


Also, I was wondering if in the following from postpro.f
----
      ! locate points (iel,iproc,r,s,t)
      if(iffind) then
        call findpts(ipth,iTl(3,1),iTlS,
     &               iTl(1,1),iTlS,
     &               iTL(2,1),iTlS,
     &               rTl(nndim+2,1),rTlS,
     &               rTl(1,1),rTlS,
     &               rTl(2,1),rTlS,
     &               rTl(3,1),rTlS,
     &               rTl(4,1),rTlS,n)
        nh(ih) = n ! store number of points for a given handle
      endif
---
the 5th line should read
     &               iTl(2,1),iTlS,
(small "l" in iTl)?


Thanks,
Markus



Quoting nek5000-users at lists.mcs.anl.gov:

> Hello Stefan,
>
> I tried the updated hpts routine ( version 518 ) and was having seg fault
> issues. I tried under two different compilers ( intel and pgi ) and was
> getting the seg fault error in both cases. I tried the same case without
> calling the hpts routine and that works fine.
>
> Error -1 :
> <code>
> [i124-101:13913] *** Process received signal ***
> [i124-101:13913] Signal: Segmentation fault (11)
> [i124-101:13913] Signal code: Address not mapped (1)
> [i124-101:13913] Failing at address: 0x1e29b000
> </code>
>
> Is that an issue with compiler /mpi or something else ? Thanks for any help.
>
> Regards
> Shriram
>