From fb15c36443b0eb0d9c0d3f9bf1833b9f7e459013 Mon Sep 17 00:00:00 2001
From: David A Minton <daminton@purdue.edu>
Date: Mon, 15 May 2023 12:14:00 -0400
Subject: [PATCH] Added Swifter vs Swiftest comparison for RMVS

---
 examples/.gitignore                           |  1 +
 examples/rmvs_swifter_comparison/.gitignore   |  3 ++
 .../1pl_1tp_encounter/.gitignore              |  3 ++
 .../1pl_1tp_encounter/init_cond.py            | 38 +++++++++++++++++++
 4 files changed, 45 insertions(+)
 create mode 100644 examples/rmvs_swifter_comparison/.gitignore
 create mode 100644 examples/rmvs_swifter_comparison/1pl_1tp_encounter/.gitignore
 create mode 100755 examples/rmvs_swifter_comparison/1pl_1tp_encounter/init_cond.py

diff --git a/examples/.gitignore b/examples/.gitignore
index 90ecccf87..dbdf17577 100644
--- a/examples/.gitignore
+++ b/examples/.gitignore
@@ -9,3 +9,4 @@
 !solar_impact
 !Swifter_Swiftest
 !swiftest_performance
+!rmvs_swifter_comparison
diff --git a/examples/rmvs_swifter_comparison/.gitignore b/examples/rmvs_swifter_comparison/.gitignore
new file mode 100644
index 000000000..ecf4c57c8
--- /dev/null
+++ b/examples/rmvs_swifter_comparison/.gitignore
@@ -0,0 +1,3 @@
+*
+!.gitignore
+!1pl_1tp_encounter
diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/.gitignore b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/.gitignore
new file mode 100644
index 000000000..80c7e7bd6
--- /dev/null
+++ b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/.gitignore
@@ -0,0 +1,3 @@
+*
+!.gitignore
+!init_cond.py
diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/init_cond.py b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/init_cond.py
new file mode 100755
index 000000000..9d997f720
--- /dev/null
+++ b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/init_cond.py
@@ -0,0 +1,38 @@
+#!/usr/bin/env python3
+import numpy as np
+import swiftest 
+
+# Simple initial conditions of a circular planet with one test particle in a close encounter state
+# Simulation start, stop, and output cadence times
+t_0	  = 0 # simulation start time
+deltaT	= 0.25 * swiftest.JD2S / swiftest.YR2S   # simulation step size
+end_sim = 0.15
+t_print = deltaT  #output interval to print results
+
+iout = int(np.ceil(t_print / deltaT))
+
+radius = np.double(4.25875607065041e-05)
+Gmass = np.double(0.00012002693582795244940133) 
+apl = np.longdouble(1.0)
+atp = np.longdouble(1.01)
+vpl = np.longdouble(2 * np.pi)
+vtp = np.longdouble(2 * np.pi / np.sqrt(atp))
+
+p_pl = np.array([apl, 0.0, 0.0], dtype=np.double)
+v_pl = np.array([0.0, vpl, 0.0], dtype=np.double)
+
+p_tp = np.array([atp, 0.0, 0.0], dtype=np.double)
+v_tp = np.array([0.0, vtp, 0.0], dtype=np.double)
+
+rhill = np.double(apl * 0.0100447248332378922085)
+
+sim = swiftest.Simulation(simdir="swiftest_sim",init_cond_format="XV", general_relativity=False, integrator="RMVS")
+sim.clean()
+sim.add_solar_system_body(["Sun"])
+sim.add_body(name=["Planet"], rh=[p_pl], vh=[v_pl], Gmass=[Gmass], radius=[radius], rhill=[rhill])
+sim.add_body(name=["TestParticle"], rh=[p_tp], vh=[v_tp])
+sim.set_parameter(tstart=t_0, tstop=end_sim, dt=deltaT, istep_out=iout, dump_cadence=0)
+sim.save()
+
+sim.set_parameter(simdir="swifter_sim",codename="Swifter",init_cond_file_type="ASCII",output_format="XV",output_file_type="REAL8")
+sim.save()
\ No newline at end of file