Skip to content
Snippets Groups Projects
Commit 61b1df1e authored by Timo Koch's avatar Timo Koch
Browse files

Merge branch 'fix/network1d3d_convergencetest' into 'master'

[test][convergence] Change condition of checking error i twice to checking...

See merge request !1634
parents a3fdcd62 059f2e2e
No related branches found
No related tags found
1 merge request!1634[test][convergence] Change condition of checking error i twice to checking...
...@@ -87,7 +87,7 @@ for permIndex in range(0, len(k)): ...@@ -87,7 +87,7 @@ for permIndex in range(0, len(k)):
if isnan(l2_matrix[i]) or isinf(l2_matrix[i]): if isnan(l2_matrix[i]) or isinf(l2_matrix[i]):
sys.stderr.write("l2 error norm is not a number!\n") sys.stderr.write("l2 error norm is not a number!\n")
sys.exit(1) sys.exit(1)
if not (l2_matrix[i] < 1e-12 or l2_matrix[i] < 1e-12): if not (l2_matrix[i] < 1e-12 or l2_matrix[i+1] < 1e-12):
rateP = (log(l2_matrix[i])-log(l2_matrix[i+1]))/(log(eps[i])-log(eps[i+1])) rateP = (log(l2_matrix[i])-log(l2_matrix[i+1]))/(log(eps[i])-log(eps[i+1]))
message = "{}\t{:0.4e}\t{:0.4e}\t{:0.4e}\n".format(i, eps[i], l2_matrix[i], rateP) message = "{}\t{:0.4e}\t{:0.4e}\t{:0.4e}\n".format(i, eps[i], l2_matrix[i], rateP)
logfile.write(message) logfile.write(message)
......
...@@ -39,7 +39,7 @@ logfile.write("-"*50 + "\n") ...@@ -39,7 +39,7 @@ logfile.write("-"*50 + "\n")
for i in range(len(error)-1): for i in range(len(error)-1):
if isnan(error[i]) or isinf(error[i]): if isnan(error[i]) or isinf(error[i]):
continue continue
if not (error[i] < 1e-12 or error[i] < 1e-12): if not (error[i] < 1e-12 or error[i+1] < 1e-12):
rate = (log(error[i])-log(error[i+1]))/(log(hmax[i])-log(hmax[i+1])) rate = (log(error[i])-log(error[i+1]))/(log(hmax[i])-log(hmax[i+1]))
message = "{}\t{:0.4e}\t{:0.4e}\t{:0.4e}\n".format(i, hmax[i], error[i], rate) message = "{}\t{:0.4e}\t{:0.4e}\t{:0.4e}\n".format(i, hmax[i], error[i], rate)
logfile.write(message) logfile.write(message)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment