[Nek5000-users] Seg fault (may be from Hpts routine)
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Thu May 20 23:50:54 CDT 2010
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
More information about the Nek5000-users
mailing list