Merge pull request #4 from AshtonSBradley/add-current-W>0-in-2D

+ rotating frame current
AshtonSBradley · Apr 3, 2023 · 20a1d9a · 20a1d9a · AshtonSBradley · Apr 3, 2023
2 parents c0e243f + 5bbe329
commit 20a1d9a
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 8 deletions.
diff --git a/.github/workflows/CI.yml b/.github/workflows/CI.yml
@@ -10,12 +10,13 @@ jobs:
       fail-fast: false
       matrix:
         version:
-          - '1.6'
-          - 'nightly'
+          # - '1.6'
+          - '1.8'
+          # - 'nightly'
         os:
           - ubuntu-latest
         arch:
-          - x64
+          # - x64
           - x86
     steps:
       - uses: actions/checkout@v2

diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "QuantumFluidSpectra"
 uuid = "8ddfec37-9569-48bd-b1f3-d8731767ee3a"
 authors = ["Ashton Bradley <[email protected]> and contributors"]
-version = "0.1.4"
+version = "0.1.5"
 
 [deps]
 FFTW = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"

diff --git a/src/analysis.jl b/src/analysis.jl
@@ -41,11 +41,12 @@ function current(psi::Psi{1})
     return jx
 end
 
-function current(psi::Psi{2})
-	@unpack ψ = psi 
+function current(psi::Psi{2},Ω = 0)
+	@unpack ψ,X = psi 
+    x,y = X
     ψx,ψy = gradient(psi)
-	jx = @. imag(conj(ψ)*ψx) # TODO add rot frame -n(Wxr) term, and in 3D
-	jy = @. imag(conj(ψ)*ψy)
+	jx = @. imag(conj(ψ)*ψx) + Ω*y'  
+	jy = @. imag(conj(ψ)*ψy) - Ω*x
 	return jx,jy
 end