ZINZINBIN
diff --git a/Collapse file
‎example.ipynb‎
Copy file name to clipboardExpand all lines: example.ipynb
+11-11Lines changed: 11 additions & 11 deletions b/Collapse file
‎example.ipynb‎
Copy file name to clipboardExpand all lines: example.ipynb
+11-11Lines changed: 11 additions & 11 deletions
diff --git a/Collapse file
‎pic_1d_case_symp.py‎ ‎pic_1d_bump_on_tail.py‎pic_1d_case_symp.py renamed to pic_1d_bump_on_tail.py
Copy file name to clipboardExpand all lines: pic_1d_bump_on_tail.py
+5-5Lines changed: 5 additions & 5 deletions b/Collapse file
‎pic_1d_case_symp.py‎ ‎pic_1d_bump_on_tail.py‎pic_1d_case_symp.py renamed to pic_1d_bump_on_tail.py
Copy file name to clipboardExpand all lines: pic_1d_bump_on_tail.py
+5-5Lines changed: 5 additions & 5 deletions
diff --git a/Collapse file
‎pic_2d_case.py‎
Copy file name to clipboardExpand all lines: pic_2d_case.py b/Collapse file
‎pic_2d_case.py‎
Copy file name to clipboardExpand all lines: pic_2d_case.py
diff --git a/Collapse file
‎result/PIC_dist.png‎
Copy file name to clipboard
46.2 KB
Display the source diff
Display the rich diff b/Collapse file
‎result/PIC_dist.png‎
Copy file name to clipboard
46.2 KB
Display the source diff
Display the rich diff
diff --git a/Collapse file
‎result/bumpontail.gif‎
Copy file name to clipboard
18.5 MB
Display the source diff
Display the rich diff b/Collapse file
‎result/bumpontail.gif‎
Copy file name to clipboard
18.5 MB
Display the source diff
Display the rich diff
diff --git a/Collapse file
‎result/hamiltonian_dist.png‎
Copy file name to clipboard
25.3 KB
Display the source diff
Display the rich diff b/Collapse file
‎result/hamiltonian_dist.png‎
Copy file name to clipboard
25.3 KB
Display the source diff
Display the rich diff
diff --git a/Collapse file
‎result/hamiltonian_leapfrog.png‎
Copy file name to clipboard
4.6 KB
Display the source diff
Display the rich diff b/Collapse file
‎result/hamiltonian_leapfrog.png‎
Copy file name to clipboard
4.6 KB
Display the source diff
Display the rich diff
diff --git a/Collapse file
‎result/hamiltonian_midpoint.png‎
Copy file name to clipboard
4.52 KB
Display the source diff
Display the rich diff b/Collapse file
‎result/hamiltonian_midpoint.png‎
Copy file name to clipboard
4.52 KB
Display the source diff
Display the rich diff
diff --git a/Collapse file
‎result/simulation_dist.gif‎
Copy file name to clipboard
1.17 MB
Display the source diff
Display the rich diff b/Collapse file
‎result/simulation_dist.gif‎
Copy file name to clipboard
1.17 MB
Display the source diff
Display the rich diff
diff --git a/Collapse file
‎src/PIC1D.py‎
Copy file name to clipboardExpand all lines: src/PIC1D.py
+1-1Lines changed: 1 addition & 1 deletion b/Collapse file
‎src/PIC1D.py‎
Copy file name to clipboardExpand all lines: src/PIC1D.py
+1-1Lines changed: 1 addition & 1 deletion
@@ -1,5 +1,5 @@
 import argparse
-from src.SymPIC1D import PICsolver
+from src.PIC1D_dist import PICsolver
 from src.dist import BumpOnTail1D
 
 def parsing():
@@ -8,8 +8,8 @@ def parsing():
     parser.add_argument("--num_mesh", type = int, default = 400)
     parser.add_argument("--t_min", type = float, default = 0)
     parser.add_argument("--t_max", type = float, default = 50)
-    parser.add_argument("--dt", type = float, default = 1.0)
-    parser.add_argument("--L", type = float, default = 50)
+    parser.add_argument("--dt", type = float, default = 1)
+    parser.add_argument("--L", type = float, default = 10)
     parser.add_argument("--n0", type = float, default = 1.0)
     parser.add_argument("--eta", type = float, default = 10.0)
     parser.add_argument("--a", type = float, default = 0.3)
@@ -19,7 +19,7 @@ def parsing():
     parser.add_argument("--gamma", type = float, default = 1.0)
     parser.add_argument("--use_animation", type = bool, default = True)
     parser.add_argument("--plot_freq", type = int, default = 10)
-    parser.add_argument("--save_dir", type=str, default="./result/simulation_symplectic.gif")
+    parser.add_argument("--save_dir", type=str, default="./result/")
     args = vars(parser.parse_args())
     return args
 
@@ -28,7 +28,7 @@ def parsing():
     args = parsing()
 
     # Initial distribution: Bump-On-Tail distribution   
-    dist = BumpOnTail1D(eta = args['eta'], a = args['a'], v0 = args['v0'], sigma = args['sigma'], beta = args['beta'], n_samples=args['num_particle'])
+    dist = BumpOnTail1D(eta = args['eta'], a = args['a'], v0 = args['v0'], sigma = args['sigma'], beta = args['beta'], n_samples=args['num_particle'], L = args['L'])
 
     # PIC solver
     solver = PICsolver(
 
@@ -227,7 +227,7 @@ def solve(self):
         plt.figure(figsize = (6,4))
         plt.plot(np.arange(len(E_list)), E_list, 'b')
         plt.xlabel("Time step")
-        plt.ylabel("Hamiltonian")
+        plt.ylabel("Relative error ($(H(t)-H(t=0))/H(t=0)$)")
         plt.tight_layout()
 
         if self.save_dir is not None: