* installing *source* package 'UBL' ... ** using staged installation ** libs using C compiler: 'gcc.exe (GCC) 12.2.0' using Fortran compiler: 'GNU Fortran (GCC) 12.2.0' gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools43/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c UBL_init.c -o UBL_init.o gfortran -pedantic -O2 -mfpmath=sse -msse2 -mstackrealign -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] gfortran -pedantic -O2 -mfpmath=sse -msse2 -mstackrealign -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] gcc -shared -s -static-libgcc -o UBL.dll tmp.def UBL_init.o neighbours.o phi.o -Ld:/rtools43/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools43/x86_64-w64-mingw32.static.posix/lib -lgfortran -lm -lquadmath -LD:/RCompile/recent/R/bin/x64 -lR installing to d:/RCompile/CRANincoming/R-devel/lib/00LOCK-UBL/00new/UBL/libs/x64 ** 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 ** testing if installed package can be loaded from final location ** testing if installed package keeps a record of temporary installation path * DONE (UBL)