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

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

with _exactly_ I mean to machine precision

Stefan
 
On May 21, 2010, at 6:50 AM, nek5000-users at lists.mcs.anl.gov wrote:

> 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
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users