* installing to library ‘/home/hornik/tmp/R.check/r-devel-gcc/Work/build/Packages’
* installing *source* package ‘UBL’ ...
** package ‘UBL’ successfully unpacked and MD5 sums checked
** using staged installation
** libs
using C compiler: ‘gcc-13 (Debian 13.2.0-23) 13.2.0’
using Fortran compiler: ‘GNU Fortran (Debian 13.2.0-23) 13.2.0’
make[1]: Entering directory '/home/hornik/tmp/scratch/Rtmp5Lmp0U/R.INSTALL2123f120e73681/UBL/src'
gcc-13 -I"/home/hornik/tmp/R.check/r-devel-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c UBL_init.c -o UBL_init.o
gfortran-13 -fpic -g -O2 -Wall -mtune=native -c neighbours.f90 -o neighbours.o
neighbours.f90:8:2:
8 | nnom, Cl, distm, numD, res)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:52:2:
52 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:54:2:
54 | integer, parameter :: dp = kind(1.0d0)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:57:2:
57 | integer , intent(in) :: p, k, n, nnum, nnom, Cl
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:58:2:
58 | real (dp), intent(in) :: tgtData(1,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:59:2:
59 | real (dp), intent(in) :: numData(nnum,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:60:2:
60 | integer , intent(in) :: nomData(nnom,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:61:2:
61 | integer, intent(inout) :: res(k,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:62:2:
62 | real (dp), intent(inout) :: distm(n,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:63:2:
63 | real (dp), intent(inout) :: numD(nnum,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:66:2:
66 | integer :: i,j,l,bestIndex
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:67:2:
67 | integer :: used(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:68:2:
68 | real (dp) :: bestDist
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:69:2:
69 | real (dp) :: ranges(nnum,2)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:70:2:
70 | real (dp) :: mean(nnum), sd(nnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:72:2:
72 | distm=0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:73:2:
73 | numD=numData
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:76:2:
76 | if (p >= 1) then ! p-norm
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:77:2:
77 | call callpNorm(p, distm, numD, nnum, n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:78:2:
78 | else if (p == 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:79:2:
79 | call callChebyshev(distm, numD, nnum, n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:80:2:
80 | else if (p == -1) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:81:2:
81 | call callCanberra(distm, numD, nnum, n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:82:2:
82 | else if (p == -2) then ! overlap metric: only for nominal attributes
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:83:2:
83 | call callOverlap(distm, nomData, nnom, n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:84:2:
84 | else if (p == -3) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:85:2:
85 | call callHEOM(distm, numD, nnum, nomData, nnom, n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:86:2:
86 | else if (p == -4) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:87:2:
87 | call callHVDM(distm, numD, nomData, nnum, nnom, tgtData, n, Cl)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:97:2:
97 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:102:2:
102 | do i=1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:103:2:
103 | used = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:104:2:
104 | used(i) = 1 ! the example is not a neighbour of himself
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:105:2:
105 | do j=1,k
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:106:2:
106 | bestDist = huge(bestDist)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:107:2:
107 | bestIndex = -1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:108:2:
108 | do l=1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:109:2:
109 | if (distm(l,i) < bestDist .and. used(l) == 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:110:2:
110 | bestDist = distm(l,i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:111:2:
111 | bestIndex = l
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:112:2:
112 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:113:2:
113 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:114:2:
114 | used(bestIndex) = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:115:2:
115 | res(j,i) = bestIndex
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:116:2:
116 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:117:2:
117 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:120:2:
120 | contains
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:123:2:
123 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:125:2:
125 | integer , intent(in) :: n, nnum, p
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:126:2:
126 | real (dp), intent(in) :: numD(nnum,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:127:2:
127 | real (dp), intent(inout) :: distm(n,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:130:2:
130 | integer :: i,j
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:132:2:
132 | do j=1,n-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:133:2:
133 | do i=j+1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:134:2:
134 | distm(i,j) = distm(i,j) + pNorm(p, numD(:,j), numD(:,i), nnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:135:2:
135 | distm(j,i) = distm(i,j)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:136:2:
136 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:137:2:
137 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:141:2:
141 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:143:2:
143 | integer , intent(in) :: n, nnum
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:144:2:
144 | real (dp), intent(in) :: numD(nnum,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:145:2:
145 | real (dp), intent(inout) :: distm(n,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:148:2:
148 | integer :: i, j
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:150:2:
150 | do j=1,n-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:151:2:
151 | do i=j+1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:152:2:
152 | distm(i,j) = distm(i,j) + chebyshev(numD(:,j), numD(:,i), nnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:153:2:
153 | distm(j,i) = distm(i,j)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:154:2:
154 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:155:2:
155 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:159:2:
159 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:161:2:
161 | integer, intent(in) :: nnum, nnom, n, Cl
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:162:2:
162 | real(dp), intent(inout) :: distm(n,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:163:2:
163 | real (dp), intent(in) :: numD(nnum,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:164:2:
164 | integer, intent(in) :: nomData(nnom,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:165:2:
165 | real (dp), intent(in) :: tgtData(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:168:2:
168 | integer :: i,j
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:169:2:
169 | real (dp):: sd(nnum), mean(nnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:171:2:
171 | sd = 0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:172:2:
172 | mean = 0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:173:2:
173 | if (nnum /= 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:174:2:
174 | do i=1,nnum
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:175:2:
175 | mean(i) = sum(numD(i,:))/n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:176:2:
176 | sd(i) = sqrt(sum((numD(i,:)-mean(i))**2) / (n-1))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:177:2:
177 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:178:2:
178 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:180:2:
180 | do j=1,n-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:181:2:
181 | do i=j,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:182:2:
182 | distm(i,j) = distm(i,j) + &
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:183:2:
183 | HVDM(numD(:,j), numD(:,i), sd, nomData(:,:),&
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:184:2:
184 | nnum, nnom, tgtData, n, i, j, Cl)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:185:2:
185 | distm(j,i) = distm(i,j)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:186:2:
186 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:187:2:
187 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:192:2:
192 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:194:2:
194 | integer, intent(in) :: n, nnom
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:195:2:
195 | integer, intent(in) :: nomData(nnom,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:196:2:
196 | real(dp), intent(inout) :: distm(n,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:199:2:
199 | integer :: i, j
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:200:2:
200 | do j=1,n-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:201:2:
201 | do i=j+1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:202:2:
202 | distm(i,j) = distm(i,j) +&
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:203:2:
203 | overlap(nomData(:,j), nomData(:,i), nnom)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:204:2:
204 | distm(j,i)=distm(i,j)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:205:2:
205 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:206:2:
206 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:210:2:
210 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:212:2:
212 | integer, intent(in) :: n, nnum
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:213:2:
213 | real (dp), intent(in) :: numD(nnum,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:214:2:
214 | real(dp), intent(inout) :: distm(n,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:217:2:
217 | integer :: i, j
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:218:2:
218 | do j=1,n-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:219:2:
219 | do i=j+1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:220:2:
220 | distm(i,j) = distm(i,j) + canberra(numD(:,j), numD(:,i), nnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:221:2:
221 | distm(j,i) = distm(i,j)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:222:2:
222 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:223:2:
223 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:237:2:
237 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:239:2:
239 | integer, intent(in) :: n, nnum, nnom
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:240:2:
240 | real (dp), intent(in) :: numD(nnum,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:241:2:
241 | integer, intent(in) :: nomData(nnom,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:242:2:
242 | real(dp), intent(inout) :: distm(n,n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:245:2:
245 | integer :: i, j
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:246:2:
246 | real (dp) :: ranges(nnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:248:2:
248 | do i=1,nnum
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:249:2:
249 | ranges(i) = maxval(numD(i,:)) - minval(numD(i,:))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:250:2:
250 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:252:2:
252 | do j=1,n-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:253:2:
253 | do i=j+1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:254:2:
254 | distm(i,j) = distm(i,j) + sqrt(HEOMnum(numD(:,j), &
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:255:2:
255 | numD(:,i), nnum, ranges) + &
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:256:2:
256 | overlap(nomData(:,j), nomData(:,i), nnom))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:257:2:
257 | distm(j,i) = distm(i,j)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:258:2:
258 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:259:2:
259 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:265:2:
265 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:266:2:
266 | integer, intent(in) :: p, d
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:267:2:
267 | real (dp), intent(in) :: a(d), b(d)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:269:2:
269 | pNorm = (sum((abs(a-b))**p))**(1.0/dble(p))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:274:2:
274 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:275:2:
275 | integer, intent(in) :: d
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:276:2:
276 | real (dp), intent(in) :: a(d), b(d), ranges(d)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:277:2:
277 | integer :: i
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:279:2:
279 | real (dp), PARAMETER :: eps = 1d-30
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:281:2:
281 | HEOMnum =0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:282:2:
282 | do i=1, d
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:284:2:
284 | if (ranges(i) > eps) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:285:2:
285 | HEOMnum = HEOMnum +(abs(a(i)-b(i))/ranges(i))**2
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:286:2:
286 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:287:2:
287 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:292:2:
292 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:293:2:
293 | integer, intent(in) :: d
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:294:2:
294 | real (dp), intent(in) :: a(d), b(d)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:297:2:
297 | chebyshev = maxval(abs(a-b))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:302:2:
302 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:303:2:
303 | integer, intent(in) :: d
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:304:2:
304 | integer, intent(in) :: a(d), b(d)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:305:2:
305 | integer :: i,s
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:307:2:
307 | s = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:308:2:
308 | do i=1,d
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:309:2:
309 | if(a(i) /= b(i))then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:310:2:
310 | s = s + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:311:2:
311 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:312:2:
312 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:314:2:
314 | overlap = s
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:319:2:
319 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:320:2:
320 | integer, intent(in) :: d
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:321:2:
321 | real (dp), intent(in) :: a(d), b(d)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:323:2:
323 | canberra = sum(abs(a-b)/(abs(a)+abs(b)))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:329:2:
329 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:330:2:
330 | integer, intent(in) :: dimnum, dimnom, n, i, j, Cl
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:331:2:
331 | real (dp), intent(in) :: numa(dimnum), numb(dimnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:332:2:
332 | integer, intent(in) :: nom(dimnom, n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:333:2:
333 | real (dp), intent(in) :: tgtData(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:334:2:
334 | real (dp), intent(in) :: sd(dimnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:336:2:
336 | real (dp) :: resnum, resnom
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:339:2:
339 | if (dimnom == 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:340:2:
340 | resnum = HVDMnum(numa, numb, dimnum, sd)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:341:2:
341 | resnom = 0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:342:2:
342 | else if (dimnum == 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:343:2:
343 | resnum = 0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:344:2:
344 | resnom = HVDMnom(nom, dimnom, tgtData, n, i, j, Cl)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:345:2:
345 | else
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:346:2:
346 | resnum = HVDMnum(numa, numb, dimnum, sd)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:347:2:
347 | resnom = HVDMnom(nom, dimnom, tgtData, n, i, j, Cl)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:348:2:
348 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:350:2:
350 | HVDM = sqrt(resnum + resnom)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:355:2:
355 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:356:2:
356 | integer, intent(in) :: dimnum
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:357:2:
357 | real (dp), intent(in) :: numa(dimnum), numb(dimnum), sd(dimnum)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:360:2:
360 | integer :: i
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:362:2:
362 | real (dp), PARAMETER :: eps = 1d-30
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:364:2:
364 | HVDMnum = 0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:365:2:
365 | do i=1, dimnum
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:366:2:
366 | if (sd(i) > eps) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:367:2:
367 | HVDMnum = HVDMnum + ((abs(numa(i)-numb(i))/(4*sd(i)))**2)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:368:2:
368 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:369:2:
369 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:374:2:
374 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:375:2:
375 | integer, intent(in) :: dimnom, n, i, j, Cl
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:376:2:
376 | integer, intent(in) :: nomdata(dimnom, n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:377:2:
377 | real (dp), intent(in) :: tgtData(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:379:2:
379 | real (dp) :: res
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:380:2:
380 | integer :: att, k, l, p
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:381:2:
381 | real (dp) :: ci, cj, NCli, NClj, resAt, finalresAt
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:383:2:
383 | finalresAt = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:385:2:
385 | do att=1,dimnom
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:386:2:
386 | resAt = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:387:2:
387 | ci = 0.0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:388:2:
388 | cj = 0.0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:389:2:
389 | do p=1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:390:2:
390 | if (nomdata(att,p)==nomdata(att,i)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:391:2:
391 | ci = ci + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:392:2:
392 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:393:2:
393 | if (nomdata(att,p)==nomdata(att,j)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:394:2:
394 | cj = cj + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:395:2:
395 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:396:2:
396 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:397:2:
397 | do k=1,Cl
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:398:2:
398 | NCli = 0.0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:399:2:
399 | NClj = 0.0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:400:2:
400 | do l=1,n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:401:2:
401 | if (nomdata(att,l)==nomdata(att,i) .and. tgtData(l)==k) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:402:2:
402 | NCli = NCli + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:403:2:
403 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:404:2:
404 | if (nomdata(att,l)==nomdata(att,j) .and. tgtData(l)==k) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:405:2:
405 | NClj = NClj + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:406:2:
406 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:407:2:
407 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:408:2:
408 | resAt = resAt + (abs((NCli/ci) - (NClj/cj)))**2
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:409:2:
409 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:410:2:
410 | finalresAt = finalresAt + resAt
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:411:2:
411 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:413:2:
413 | HVDMnom = finalresAt
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
neighbours.f90:379:31:
379 | real (dp) :: res
| 1
Warning: Unused variable ‘res’ declared at (1) [-Wunused-variable]
neighbours.f90:70:26:
70 | real (dp) :: mean(nnum), sd(nnum)
| 1
Warning: Unused variable ‘mean’ declared at (1) [-Wunused-variable]
neighbours.f90:69:30:
69 | real (dp) :: ranges(nnum,2)
| 1
Warning: Unused variable ‘ranges’ declared at (1) [-Wunused-variable]
neighbours.f90:70:36:
70 | real (dp) :: mean(nnum), sd(nnum)
| 1
Warning: Unused variable ‘sd’ declared at (1) [-Wunused-variable]
gfortran-13 -fpic -g -O2 -Wall -mtune=native -c phi.f90 -o phi.o
phi.f90:34:2:
34 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:36:2:
36 | integer, parameter :: dp = kind(1.0d0)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:39:2:
39 | integer , intent(in) :: n, lparms, method, npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:40:2:
40 | real (dp), intent (in) :: y(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:41:2:
41 | real (dp), intent (in) :: phiParms(lparms)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:42:2:
42 | real (dp), intent (inout) :: yPhi(n), ydPhi(n), yddPhi(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:45:2:
45 | real (dp) :: x(npts), a(npts), b(npts), c(npts), d(npts) !,&
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:47:2:
47 | integer :: i
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:49:2:
49 | x=0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:50:2:
50 | a=0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:51:2:
51 | b=0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:52:2:
52 | c=0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:53:2:
53 | d=0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:54:2:
54 | i=1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:55:2:
55 | call phiSplInit(lparms, phiParms, npts, x, a, b, c, d)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:57:2:
57 | call phiEval(n, y, method, npts, x, a, b, c, d, yPhi, ydPhi, yddPhi)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:63:2:
63 | integer, intent(in) :: n, method, npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:64:2:
64 | real (dp), intent(in) :: y(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:65:2:
65 | real(dp), intent(inout) :: yPhi(n), ydPhi(n), yddPhi(n)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:66:2:
66 | real (dp), intent(in) :: x(npts), a(npts), b(npts), c(npts), d(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:69:2:
69 | integer :: i
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:70:2:
70 | real (dp) :: yval, yvald, yvaldd
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:71:2:
71 | do i=1, n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:72:2:
72 | call phiSplValue(y(i), method, n, npts, x, a, b, c, d, yval, yvald, yvaldd)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:73:2:
73 | yPhi(i) = yval
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:74:2:
74 | ydPhi(i) = yvald
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:75:2:
75 | yddPhi(i) = yvaldd
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:76:2:
76 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:83:2:
83 | integer , intent(in) :: lparms, npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:84:2:
84 | real (dp), intent(in) :: phiParms(lparms)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:85:2:
85 | real (dp), intent(inout) :: x(npts), a(npts), b(npts), c(npts), d(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:88:2:
88 | integer :: i, n
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:89:2:
89 | real (dp) :: newm(npts), y(npts), m(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:91:2:
91 | do i=1, npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:92:2:
92 | x(i) = phiParms(3 * i - 2)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:93:2:
93 | y(i) = phiParms(3 * i - 1)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:94:2:
94 | m(i) = phiParms(3 * i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:95:2:
95 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:97:2:
97 | newm=m
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:98:2:
98 | call pchipSet(npts, x, y, m, newm, c, d)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:100:2:
100 | a=y
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:101:2:
101 | b=newm
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:108:2:
108 | integer, intent(in) :: method, n, npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:109:2:
109 | real (dp), intent(in) :: y, x(npts), a(npts), b(npts), c(npts), d(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:110:2:
110 | real(dp), intent(inout) :: yval, yvald, yvaldd
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:113:2:
113 | integer :: extrapol ! //linear
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:115:2:
115 | extrapol = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:116:2:
116 | call pchipVal(npts, x, a, b, c, d, y, extrapol, yval, yvald, yvaldd)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:160:2:
160 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:162:2:
162 | integer , intent(in) :: npts, extrapol
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:163:2:
163 | real (dp), intent(in) :: x(npts), a(npts), b(npts), c(npts), d(npts), xval
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:164:2:
164 | real (dp), intent(inout) :: yval, yvald, yvaldd
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:167:2:
167 | integer :: i, rightmostClosed, allInside, mfl, left
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:168:2:
168 | real (dp) :: s
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:170:2:
170 | i = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:171:2:
171 | left = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:172:2:
172 | rightmostClosed = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:173:2:
173 | allInside = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:174:2:
174 | mfl = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:175:2:
175 | s=0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:177:2:
177 | call findInterval(x, npts, xval, rightmostClosed, allInside, i,&
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:178:2:
178 | mfl, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:179:2:
179 | i = left
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:182:2:
182 | if (extrapol == 0 .and. (i == 0 .or. i == npts))then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:183:2:
183 | if (i == 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:184:2:
184 | i = i + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:185:2:
185 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:186:2:
186 | yval = a(i) + b(i) * (xval - x(i))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:187:2:
187 | yvald = b(i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:188:2:
188 | yvaldd = 0 ! // This is the case where ydd does not exist
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:189:2:
189 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:190:2:
190 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:192:2:
192 | s = xval -x(i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:193:2:
193 | yval = a(i) + s * (b(i) + s * (c(i) + s * d(i)))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:194:2:
194 | yvald = b(i) + s * (2 * c(i) + s * (3 * d(i)))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:195:2:
195 | yvaldd = 2 * c(i) + s * 6 * d(i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:227:2:
227 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:229:2:
229 | integer, intent(in) :: npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:230:2:
230 | real(dp), intent(inout) :: x(npts), y(npts), newm(npts), c(npts), d(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:231:2:
231 | real (dp), intent(in) :: m(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:234:2:
234 | integer :: i
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:235:2:
235 | real (dp) :: h(npts), delta(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:238:2:
238 | do i=1, npts-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:239:2:
239 | h(i) = x(i+1) - x(i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:240:2:
240 | delta(i) = (y(i+1) - y(i))/ h(i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:241:2:
241 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:243:2:
243 | call pchipSlopeMonoFC(npts, newm, delta) ! adjusts newm values
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:245:2:
245 | do i=1, npts-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:246:2:
246 | c(i) = (3*delta(i) - 2 * newm(i) - newm(i+1)) / h(i)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:247:2:
247 | d(i) = (newm(i) - 2 * delta(i) + newm(i+1)) / (h(i) * h(i))
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:248:2:
248 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:279:2:
279 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:281:2:
281 | integer, intent(in) :: npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:282:2:
282 | real (dp), intent(in) :: delta(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:283:2:
283 | real(dp), intent(inout) :: newm(npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:286:2:
286 | integer :: k, k1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:287:2:
287 | real (dp) :: Sk, alpha, beta, a2b3, ab23, tauS
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:289:2:
289 | do k=1,npts-1 ! modify both newm(i) and newm(i+1) if needed
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:290:2:
290 | Sk = delta(k)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:291:2:
291 | k1 = k + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:293:2:
293 | if (abs(Sk) == 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:294:2:
294 | newm(k) = 0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:295:2:
295 | newm(k1) = 0.0d0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:296:2:
296 | else
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:297:2:
297 | alpha = newm(k) / Sk
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:298:2:
298 | beta = newm(k1) / Sk
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:299:2:
299 | if (abs(newm(k)) /= 0 .and. alpha < 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:300:2:
300 | newm(k) = - newm(k)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:301:2:
301 | alpha = newm(k) / Sk
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:302:2:
302 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:304:2:
304 | if (abs(newm(k1)) /= 0 .and. beta < 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:305:2:
305 | newm(k1) = -newm(k1)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:306:2:
306 | beta = newm(k1) / Sk
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:307:2:
307 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:309:2:
309 | a2b3 = 2 * alpha + beta - 3
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:310:2:
310 | ab23 = alpha + 2 * beta - 3
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:312:2:
312 | if (a2b3 > 0 .and. ab23 > 0 .and. &
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:313:2:
313 | alpha * (a2b3 + ab23) < a2b3 * a2b3) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:315:2:
315 | tauS = 3 * Sk / sqrt(alpha * alpha + beta * beta)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:316:2:
316 | newm(k) = tauS * alpha
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:317:2:
317 | newm(k1) = tauS * beta
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:318:2:
318 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:319:2:
319 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:320:2:
320 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:326:2:
326 | mflag, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:353:2:
353 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:355:2:
355 | integer , intent(in) :: npts, rightmostClosed, allInside
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:356:2:
356 | integer , intent(inout) :: mflag, left, ilo
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:357:2:
357 | real (dp), intent(in) :: xt(npts), x
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:360:2:
360 | integer :: istep, ihi
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:362:2:
362 | if (ilo <= 0) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:363:2:
363 | if (x < xt(1)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:364:2:
364 | call leftBoundary(mflag, allInside, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:365:2:
365 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:366:2:
366 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:367:2:
367 | ilo = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:368:2:
368 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:370:2:
370 | ihi = ilo + 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:372:2:
372 | if (ihi >= npts) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:373:2:
373 | if (x >= xt(npts)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:374:2:
374 | call rightBoundary(xt, npts, x, mflag, allInside, &
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:375:2:
375 | rightmostClosed, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:376:2:
376 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:377:2:
377 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:378:2:
378 | if (npts <= 1) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:379:2:
379 | call leftBoundary(mflag, allInside, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:380:2:
380 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:381:2:
381 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:382:2:
382 | ilo = npts - 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:383:2:
383 | ihi = npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:384:2:
384 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:386:2:
386 | if (x < xt(ihi)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:387:2:
387 | if (x >= xt(ilo)) then ! /* `lucky': same interval as last time */
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:388:2:
388 | mflag = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:389:2:
389 | left = ilo
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:390:2:
390 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:391:2:
391 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:393:2:
393 | istep = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:394:2:
394 | do while (ilo > 1)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:395:2:
395 | ihi = ilo
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:396:2:
396 | ilo = ihi - istep
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:397:2:
397 | if (x >= xt(ilo)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:398:2:
398 | call gotoL50(xt, npts, x, ilo, ihi, mflag, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:399:2:
399 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:400:2:
400 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:401:2:
401 | istep = istep * 2
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:402:2:
402 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:403:2:
403 | ilo = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:404:2:
404 | if (x < xt(1)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:405:2:
405 | call leftBoundary(mflag, allInside, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:406:2:
406 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:407:2:
407 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:409:2:
409 | else
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:411:2:
411 | istep = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:412:2:
412 | do while (ihi < npts)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:413:2:
413 | ilo = ihi
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:414:2:
414 | ihi = ilo + istep
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:415:2:
415 | if (ihi >= npts) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:416:2:
416 | exit
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:417:2:
417 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:418:2:
418 | if (x < xt(ihi)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:419:2:
419 | call gotoL50(xt, npts, x, ilo, ihi, mflag, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:420:2:
420 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:421:2:
421 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:422:2:
422 | istep = istep * 2
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:423:2:
423 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:425:2:
425 | if (x >= xt(npts)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:426:2:
426 | call rightBoundary(xt, npts, x, mflag, allInside, &
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:427:2:
427 | rightmostClosed, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:428:2:
428 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:429:2:
429 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:430:2:
430 | ihi = npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:431:2:
431 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:433:2:
433 | call gotoL50(xt, npts, x, ilo, ihi, mflag, left)
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:438:2:
438 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:440:2:
440 | integer, intent(inout) :: mflag, left
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:441:2:
441 | integer , intent(in) :: allInside
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:443:2:
443 | mflag = -1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:444:2:
444 | if (allInside == 1) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:445:2:
445 | left = 1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:446:2:
446 | else
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:447:2:
447 | left = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:448:2:
448 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:449:2:
449 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:454:2:
454 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:456:2:
456 | integer , intent(in) :: npts, allInside, rightmostClosed
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:457:2:
457 | integer, intent(inout) :: mflag, left
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:458:2:
458 | real (dp), intent(in) :: xt(npts), x
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:460:2:
460 | mflag = +1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:461:2:
461 | if (allInside == 1 .or. (rightmostClosed == 1 .and. x == xt(npts))) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:462:2:
462 | left = npts-1
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:463:2:
463 | else
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:464:2:
464 | left = npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:465:2:
465 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:466:2:
466 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:470:2:
470 | implicit none
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:472:2:
472 | integer , intent(in) :: npts
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:473:2:
473 | integer, intent(inout) :: mflag, left, ilo, ihi
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:474:2:
474 | real (dp), intent(in) :: xt(npts), x
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:477:2:
477 | integer :: middle
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:480:2:
480 | do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:481:2:
481 | middle = (ilo + ihi) / 2
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:482:2:
482 | if (middle == ilo) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:483:2:
483 | mflag = 0
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:484:2:
484 | left = ilo
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:485:2:
485 | return
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:486:2:
486 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:488:2:
488 | if (x >= xt(middle)) then
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:489:2:
489 | ilo = middle
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:490:2:
490 | else
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:491:2:
491 | ihi = middle
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:492:2:
492 | end if
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:493:2:
493 | end do
| 1
Warning: Nonconforming tab character at (1) [-Wtabs]
phi.f90:205:33:
205 | subroutine pchipSet(npts, x, y, m, newm, c, d)
| 1
Warning: Unused dummy argument ‘m’ at (1) [-Wunused-dummy-argument]
phi.f90:106:33:
106 | subroutine phiSplValue( y, method, n, npts, x, a, b, c, d, yval, yvald, yvaldd)
| 1
Warning: Unused dummy argument ‘method’ at (1) [-Wunused-dummy-argument]
phi.f90:106:36:
106 | subroutine phiSplValue( y, method, n, npts, x, a, b, c, d, yval, yvald, yvaldd)
| 1
Warning: Unused dummy argument ‘n’ at (1) [-Wunused-dummy-argument]
phi.f90:88:17:
88 | integer :: i, n
| 1
Warning: Unused variable ‘n’ declared at (1) [-Wunused-variable]
gcc-13 -shared -L/home/hornik/tmp/R.check/r-devel-gcc/Work/build/lib -Wl,-O1 -o UBL.so UBL_init.o neighbours.o phi.o -lgfortran -lm -lquadmath -L/home/hornik/tmp/R.check/r-devel-gcc/Work/build/lib -lR
make[1]: Leaving directory '/home/hornik/tmp/scratch/Rtmp5Lmp0U/R.INSTALL2123f120e73681/UBL/src'
make[1]: Entering directory '/home/hornik/tmp/scratch/Rtmp5Lmp0U/R.INSTALL2123f120e73681/UBL/src'
make[1]: Leaving directory '/home/hornik/tmp/scratch/Rtmp5Lmp0U/R.INSTALL2123f120e73681/UBL/src'
installing to /home/hornik/tmp/R.check/r-devel-gcc/Work/build/Packages/00LOCK-UBL/00new/UBL/libs
** R
** data
*** moving datasets to lazyload DB
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** testing if installed package can be loaded from temporary location
** checking absolute paths in shared objects and dynamic libraries
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (UBL)