Merge branch 'master' into GPU

.github/workflows/CompatHelper.yml

@@ -0,0 +1,19 @@
+name: CompatHelper
+
+on:
+  schedule:
+    - cron: '00 00 * * *'
+
+jobs:
+  CompatHelper:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: julia-actions/setup-julia@latest
+        with:
+          version: 1.3
+      - name: Pkg.add("CompatHelper")
+        run: julia -e 'using Pkg; Pkg.add("CompatHelper")'
+      - name: CompatHelper.main()
+        env:
+          GITHUB_TOKEN: ${{ secrets.GITHUB_TOKEN }}
+        run: julia -e 'using CompatHelper; CompatHelper.main()'

.github/workflows/TagBot.yml

@@ -0,0 +1,11 @@
+name: TagBot
+on:
+  schedule:
+    - cron: 0 * * * *
+jobs:
+  TagBot:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: JuliaRegistries/TagBot@v1
+        with:
+          token: ${{ secrets.GITHUB_TOKEN }}

Project.toml

@@ -29,8 +29,8 @@ StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
-AWSCore = "0.6.8"
-AWSS3 = "0.6.8"
+AWSCore = "0.6"
+AWSS3 = "0.6"
 Adapt = "1.0.1"
 CSV = "0.5,0.6"
 CUDAnative = "3.0"

README.md

@@ -37,7 +37,7 @@ process_raw!(S1,fs)
 process_raw!(S2,fs)
 R = RawData.([S1,S2],cc_len,cc_step)
 detrend!.(R)
-taper!(R)
+taper!.(R)
 bandpass!.(R,freqmin,freqmax,zerophase=true)
 FFT = compute_fft.(R)
 whiten!.(FFT,freqmin,freqmax)

src/VelocityChange/Stretching.jl

@@ -71,6 +71,13 @@ function stretching(ref::AbstractArray,cur::AbstractArray,t::AbstractArray,
     dv = 100. * dtfiner[argmax(CCfiner)]
     cc = maximum(CCfiner) # Maximum correlation coefficient of the refined analysis
 
+    # due to CubicSplineInterpolation, cc can be more than 1, causing an error when computing sqrt(1-X^2).
+    if cc > 1.0
+        cc = 1.0
+    elseif cc < -1.0
+        cc = -1.0
+    end
+
     # Error computation based on Weaver, R., C. Hadziioannou, E. Larose, and M.
     # Campillo (2011), On the precision of noise-correlation interferometry,
     # Geophys. J. Int., 185(3), 1384?1392