Make some changes
This commit is contained in:
25
src/plot_general_alg_error.py
Normal file
25
src/plot_general_alg_error.py
Normal file
@@ -0,0 +1,25 @@
|
||||
import matplotlib.pyplot as plt
|
||||
|
||||
def main():
|
||||
_, axs = plt.subplots(2)
|
||||
for i in range(6):
|
||||
x = []
|
||||
abs_err = []
|
||||
rel_err = []
|
||||
with open(f"output/error/out_{10**(i+1)}.txt", "r") as f:
|
||||
lines = f.readlines()
|
||||
for line in lines:
|
||||
x_i, abs_err_i, rel_err_i = line.strip().split(",")
|
||||
x.append(float(x_i))
|
||||
abs_err.append(float(abs_err_i))
|
||||
rel_err.append(float(rel_err_i))
|
||||
|
||||
axs[0].plot(x, abs_err, label=f"abs_err {10**(i+1)} steps")
|
||||
axs[1].plot(x, rel_err, label=f"rel_err {10**(i+1)} steps")
|
||||
|
||||
axs[0].legend(bbox_to_anchor=(1.05, 1), loc=2, borderaxespad=0.)
|
||||
axs[1].legend(bbox_to_anchor=(1.05, 1), loc=2, borderaxespad=0.)
|
||||
plt.savefig("../latex/images/problem8.pdf", bbox_inches="tight")
|
||||
|
||||
if __name__ == "__main__":
|
||||
main()
|
||||
Reference in New Issue
Block a user