23
24
25
26
27
28
29
30
31 newton_log_file1="/home/karol/mofem_install/users_modules_debug/multifield_plasticity/test_vertical_flow_form_3D/log_classical_bc"
32 newton_log_file2="/home/karol/mofem_install/users_modules_debug/multifield_plasticity/test_vertical_flow_form_3D/log_constraint_bc"
33
34 newton_data1=pd.read_fwf(newton_log_file1, header=None)
35 newton_data2=pd.read_fwf(newton_log_file2, header=None)
36
37 plt.plot(newton_data1[8].to_numpy(), ls='solid',lw=line_width, marker='o', markersize=markersize, c='blue', label="classical bc")
38 plt.plot(newton_data2[8].to_numpy(), ls='solid', lw=line_width, marker='x', markersize=markersize, c='orange', label="constraint bc")
39
40 plt.legend(loc='upper right')
41
42
43
44
45
46 plt.title('Newton method convergence')
47 plt.ylabel('absolute residual')
48 plt.xlabel('accumulated iterations')
49 plt.yscale('log')
50 plt.grid(True)
51 plt.show()
52 plt.savefig("my_log.png")
53
54
55