istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

added test in pol_limit for anpl=0 to avoid 0/0 expressions

Browse Source
This commit is contained in:
Daniele Micheletti 2015-10-14 10:26:05 +00:00
parent af5fb208b2
commit e1ba175efb
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/gray.f
View File

@ -5083,8 +5083,14 @@ c
exom=(ffo*csgam-ui*anpl*sngam)/sqrt(deno) exom=(ffo*csgam-ui*anpl*sngam)/sqrt(deno)
eyom=(-ffo*sngam-ui*anpl*csgam)/sqrt(deno) eyom=(-ffo*sngam-ui*anpl*csgam)/sqrt(deno)
c c
exxm=(ffx*csgam-ui*anpl*sngam)/sqrt(denx) c if anpl=0 the expressions for exxm and eyxm are 0/0
eyxm=(-ffx*sngam-ui*anpl*csgam)/sqrt(denx) if (anpl.eq.0.0_wp_) then
exxm=-ui*sngam
eyxm=-ui*csgam
else
exxm=(ffx*csgam-ui*anpl*sngam)/sqrt(denx)
eyxm=(-ffx*sngam-ui*anpl*csgam)/sqrt(denx)
end if
c c
if (sox.lt.0.0_wp_) then if (sox.lt.0.0_wp_) then
ext=exom ext=exom