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

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

Markus,

Fortran is not case sensitive so that's not a problem.
The warning you get is a feature and not a bug.

It seems like this point is really outside the mesh (using a tolerance of 1e-13).
Are you sure you have a point GLL point at _exactly_ this location?

The reason I am asking is the following:
Mesh coordinates generated by e.g. prex, genbox are in single precision which can be tricky in this case.


Stefan


-----Original Message-----
From: nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov
Sent: Fri 5/21/2010 01:11
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Seg fault (may be from Hpts routine)
 
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
>


_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users