blob: e2a8d191989cefaef9453836db24ffe49b9f2154 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
|
--- spack-src/src/io_tools.F90 2021-10-08 13:30:58.000000000 +0900
+++ spack-src/src/io_tools.F90.new 2021-10-08 13:30:30.000000000 +0900
@@ -12,6 +12,9 @@
#if defined (__INTEL_COMPILER)
USE IFPORT
#endif
+#if defined(__FUJITSU)
+use service_routines
+#endif
IMPLICIT NONE
PRIVATE
--- spack-src/src/excor.F90 2021-10-08 16:04:49.000000000 +0900
+++ spack-src/src/excor.F90.new 2021-10-08 16:12:28.000000000 +0900
@@ -11,6 +11,9 @@
MODULE excor
+#if defined(__FUJITSU)
+ USE, intrinsic :: ieee_arithmetic
+#endif
USE io_tools
USE Tools
USE globalmath
@@ -589,8 +592,15 @@
deallocate(grad,gradmag,gxc,dgxcdr,dfxcdgbg)
! not sure why/if this is needed
do i=1,n
+#if defined(__FUJITSU)
+ if (ieee_support_nan(tmpv(i))) then
+ if (ieee_is_nan(tmpv(i))) tmpv(i)=0.d0
+ if (ieee_is_nan(exci(i))) exci(i)=0.d0
+ endif
+#else
if (isnan(tmpv(i))) tmpv(i)=0.d0
if (isnan(exci(i))) exci(i)=0.d0
+#endif
enddo
elseif (libxc_ismgga()) then
write(std_out,*) ' atompaw not yet available for mgga -- stop '
|