Skip to content

Commit

Permalink
Using ?PossibleZeroQ and closing secitons #50
Browse files Browse the repository at this point in the history
  • Loading branch information
Philip-Lynch committed Aug 2, 2023
1 parent 1f6492a commit 29ddd77
Show file tree
Hide file tree
Showing 4 changed files with 17 additions and 17 deletions.
18 changes: 9 additions & 9 deletions Kernel/FourVelocity.m
Original file line number Diff line number Diff line change
Expand Up @@ -17,22 +17,22 @@
Begin["`Private`"];


(* ::Subsection:: *)
(* ::Subsection::Closed:: *)
(*Error messages*)


KerrGeoFourVelocity::parametrization = "Parameterization error: `1`"


(* ::Section:: *)
(* ::Section::Closed:: *)
(*Schwarzschild*)


(* ::Subsection::Closed:: *)
(*Circular, Equatorial*)


KerrGeoVelocityMino[(0|0.),p_,(0|0.),x_,initPhases_ ]:= Module[{En,L,Q,r,z,r1,r2,r3,r4,kr,zp,zm,kz, \[CapitalUpsilon]r, \[CapitalUpsilon]z,
KerrGeoVelocityMino[a_?PossibleZeroQ,p_,e_?PossibleZeroQ,x_,initPhases_ ]:= Module[{En,L,Q,r,z,r1,r2,r3,r4,kr,zp,zm,kz, \[CapitalUpsilon]r, \[CapitalUpsilon]z,
qr, qz, \[Lambda]local ,qr0, qz0, rprime, zprime, \[CapitalDelta], \[CapitalSigma], \[Omega], utContra,urContra,u\[Theta]Contra,uzContra,u\[Phi]Contra, utCo, urCo, u\[Theta]Co, u\[Phi]Co},

(*Constants of Motion*)
Expand Down Expand Up @@ -64,7 +64,7 @@
]


(* ::Section:: *)
(* ::Section::Closed:: *)
(*Kerr*)


Expand Down Expand Up @@ -130,18 +130,18 @@



(* ::Subsection:: *)
(* ::Subsection::Closed:: *)
(*Equatorial (Darwin)*)


(* ::Subsubsection:: *)
(*Circular Case*)


KerrGeoVelocityDarwin[a_,p_,(0|0.),x_,initPhases_]:= Module[{ut,ur,u\[Theta],u\[Phi], MinoVelocities,utContra,urContra,u\[Theta]Contra,u\[Phi]Contra,
KerrGeoVelocityDarwin[a_,p_,e_?PossibleZeroQ,x_,initPhases_]:= Module[{ut,ur,u\[Theta],u\[Phi], MinoVelocities,utContra,urContra,u\[Theta]Contra,u\[Phi]Contra,
utCo,urCo,u\[Theta]Co,u\[Phi]Co,utUp,urUp,u\[Theta]Up,u\[Phi]Up, utDown,urDown,u\[Theta]Down,u\[Phi]Down},

MinoVelocities = KerrGeoVelocityMino[a,p,0,x,{0,0}];
MinoVelocities = KerrGeoVelocityMino[a,p,e,x,{0,0}];

utUp="\!\(\*SuperscriptBox[\(u\), \(t\)]\)"; urUp="\!\(\*SuperscriptBox[\(u\), \(r\)]\)";
u\[Theta]Up="\!\(\*SuperscriptBox[\(u\), \(\[Theta]\)]\)"; u\[Phi]Up="\!\(\*SuperscriptBox[\(u\), \(\[Phi]\)]\)";
Expand All @@ -167,7 +167,7 @@
]


(* ::Subsubsection::Closed:: *)
(* ::Subsubsection:: *)
(*Eccentric Case*)


Expand Down Expand Up @@ -227,7 +227,7 @@
]


(* ::Section:: *)
(* ::Section::Closed:: *)
(*KerrGeoFourVelocity Wrapper*)


Expand Down
2 changes: 1 addition & 1 deletion Kernel/KerrGeoOrbit.m
Original file line number Diff line number Diff line change
Expand Up @@ -136,7 +136,7 @@
]


(* ::Section:: *)
(* ::Section::Closed:: *)
(*Kerr*)


Expand Down
10 changes: 5 additions & 5 deletions Kernel/OrbitalFrequencies.m
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@
Begin["`Private`"];


(* ::Section:: *)
(* ::Section::Closed:: *)
(*Roots of the radial and polar equations*)


Expand Down Expand Up @@ -70,7 +70,7 @@
]


(* ::Section:: *)
(* ::Section::Closed:: *)
(*Orbital Frequencies*)


Expand Down Expand Up @@ -121,11 +121,11 @@
KerrGeoProperFrequencyFactor[a_?PossibleZeroQ ,p_,e_,x_]:=(p^2 ((1+e) (28+4 e^2+(-12+p) p)-((1+e) (-4+p) (-6+2 e+p) EllipticE[(4 e)/(-6+2 e+p)]+2 (6+2 e-p) (3+e^2-p) EllipticPi[(2 e (-4+p))/((1+e) (-6+2 e+p)),(4 e)/(-6+2 e+p)])/EllipticK[(4 e)/(-6+2 e+p)]))/(2 (-1+e) (1+e)^2 (-4+p)^2)


(* ::Subsection:: *)
(* ::Subsection::Closed:: *)
(*Kerr*)


(* ::Subsubsection:: *)
(* ::Subsubsection::Closed:: *)
(*KerrGeoMinoFrequencyr*)


Expand All @@ -147,7 +147,7 @@



(* ::Subsubsection:: *)
(* ::Subsubsection::Closed:: *)
(*KerrGeoMinoFrequency\[Theta]*)


Expand Down
4 changes: 2 additions & 2 deletions Kernel/SpecialOrbits.m
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
(*Define usage for public functions*)


(* ::Section:: *)
(* ::Section::Closed:: *)
(*Create Package*)


Expand Down Expand Up @@ -602,7 +602,7 @@
]


(* ::Section:: *)
(* ::Section::Closed:: *)
(*Close the package*)


Expand Down

0 comments on commit 29ddd77

Please sign in to comment.