Changeset 14387 for vendors/PPR/src/p5e.h90
- Timestamp:
- 2021-02-03T16:51:56+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
vendors/PPR/src/p5e.h90
r14212 r14387 82 82 83 83 integer, parameter :: NSIZ = +6 84 real*8 , parameter :: ZERO = 1. d-1484 real*8 , parameter :: ZERO = 1.e-14 85 85 86 86 head = +4 ; tail = npos - 3 … … 121 121 122 122 edge(ivar,ipos) = & 123 & + ( 1. d0 / 60.d0) * &123 & + ( 1.e0 / 60.e0) * & 124 124 & fdat(1,ivar,ipos-3) & 125 & - ( 8. d0 / 60.d0) * &125 & - ( 8.e0 / 60.e0) * & 126 126 & fdat(1,ivar,ipos-2) & 127 & + (37. d0 / 60.d0) * &127 & + (37.e0 / 60.e0) * & 128 128 & fdat(1,ivar,ipos-1) & 129 & + (37. d0 / 60.d0) * &129 & + (37.e0 / 60.e0) * & 130 130 & fdat(1,ivar,ipos+0) & 131 & - ( 8. d0 / 60.d0) * &131 & - ( 8.e0 / 60.e0) * & 132 132 & fdat(1,ivar,ipos+1) & 133 & + ( 1. d0 / 60.d0) * &133 & + ( 1.e0 / 60.e0) * & 134 134 & fdat(1,ivar,ipos+2) 135 135 136 136 dfdx(ivar,ipos) = & 137 & - ( 1. d0 / 90.d0) * &137 & - ( 1.e0 / 90.e0) * & 138 138 & fdat(1,ivar,ipos-3) & 139 & + ( 5. d0 / 36.d0) * &139 & + ( 5.e0 / 36.e0) * & 140 140 & fdat(1,ivar,ipos-2) & 141 & - (49. d0 / 36.d0) * &141 & - (49.e0 / 36.e0) * & 142 142 & fdat(1,ivar,ipos-1) & 143 & + (49. d0 / 36.d0) * &143 & + (49.e0 / 36.e0) * & 144 144 & fdat(1,ivar,ipos+0) & 145 & - ( 5. d0 / 36.d0) * &145 & - ( 5.e0 / 36.e0) * & 146 146 & fdat(1,ivar,ipos+1) & 147 & + ( 1. d0 / 90.d0) * &147 & + ( 1.e0 / 90.e0) * & 148 148 & fdat(1,ivar,ipos+2) 149 149 … … 185 185 & + delh(-1) ) / xhat 186 186 xmap(-1) = - delh(-1) / xhat 187 xmap(+0) = + 0. d0187 xmap(+0) = + 0.e0 188 188 xmap(+1) = + delh(+0) / xhat 189 189 xmap(+2) = +( delh(+0) &
Note: See TracChangeset
for help on using the changeset viewer.