diff --git a/B2R_Berggren_poles.jl b/calculations/B2R_Berggren_poles.jl
similarity index 96%
rename from B2R_Berggren_poles.jl
rename to calculations/B2R_Berggren_poles.jl
index 8757ebb..0757c69 100644
--- a/B2R_Berggren_poles.jl
+++ b/calculations/B2R_Berggren_poles.jl
@@ -1,5 +1,5 @@
 using Plots
-include("p_space.jl")
+include("../p_space.jl")
 
 # contour
 p, w = get_mesh([0, 0.4 - 0.15im, 0.8, 6], [128, 128, 128])
@@ -57,4 +57,4 @@ scatter!(real.(extrapolate_E), imag.(extrapolate_E), label="extrapolated")
 scatter!(real.(basis_E), imag.(basis_E), m=:x, label="Berggren basis")
 xlims!(-0.3,0.3)
 ylims!(-0.120,0.020)
-savefig("temp/" * string(rand(UInt16)) * ".pdf")
+savefig("../temp/" * string(rand(UInt16)) * ".pdf")
diff --git a/B2R_comparison.jl b/calculations/B2R_comparison.jl
similarity index 95%
rename from B2R_comparison.jl
rename to calculations/B2R_comparison.jl
index b73ae30..5b3f228 100644
--- a/B2R_comparison.jl
+++ b/calculations/B2R_comparison.jl
@@ -1,5 +1,5 @@
 using Plots
-include("p_space.jl")
+include("../p_space.jl")
 
 berggren_mesh = get_mesh([0, 0.4 - 0.15im, 0.8, 6], [128, 128, 128])
 csm_mesh = get_mesh([0, 8 - 3im], [512])
@@ -47,5 +47,5 @@ for mesh in (berggren_mesh, csm_mesh)
     plot!(real.(mesh_E), imag.(mesh_E), label="contour")
     xlims!(-0.3,0.3)
     ylims!(-0.120,0.020)
-    savefig("temp/" * string(rand(UInt16)) * ".pdf")
+    savefig("../temp/" * string(rand(UInt16)) * ".pdf")
 end
\ No newline at end of file
diff --git a/XZ_technique.jl b/calculations/XZ_technique.jl
similarity index 96%
rename from XZ_technique.jl
rename to calculations/XZ_technique.jl
index c739da6..7ae8efe 100644
--- a/XZ_technique.jl
+++ b/calculations/XZ_technique.jl
@@ -1,6 +1,6 @@
 using LinearAlgebra, Plots
-include("ho_basis.jl")
-include("p_space.jl")
+include("../ho_basis.jl")
+include("../p_space.jl")
 
 μω_gen = 0.5 * exp(-1im * 0.47 * pi)
 μ = 0.5