Commit 75be7679 authored by Taddeüs Kroes's avatar Taddeüs Kroes

improc ass4: Updated plot results.

parent 97a4fa53
...@@ -16,7 +16,7 @@ def exit_with_usage(): ...@@ -16,7 +16,7 @@ def exit_with_usage():
def Gauss(s): def Gauss(s):
"""Sample a two-dimensional Gaussian function of scale s.""" """Sample a two-dimensional Gaussian function of scale s."""
size = int(ceil(3 * s)) size = int(ceil(3 * s))
r = 2 * size r = 2 * size + 1
W = zeros((r, r)) W = zeros((r, r))
t = float(s) ** 2 t = float(s) ** 2
a = 1 / (2 * pi * t) a = 1 / (2 * pi * t)
...@@ -32,7 +32,7 @@ def Gauss(s): ...@@ -32,7 +32,7 @@ def Gauss(s):
def Gauss1(s): def Gauss1(s):
"""Sample a one-dimensional Gaussian function of scale s.""" """Sample a one-dimensional Gaussian function of scale s."""
size = int(ceil(3 * s)) size = int(ceil(3 * s))
r = 2 * size r = 2 * size + 1
W = zeros((r,)) W = zeros((r,))
t = float(s) ** 2 t = float(s) ** 2
a = 1 / (2 * pi * t) a = 1 / (2 * pi * t)
...@@ -60,7 +60,7 @@ def gD(F, s, iorder, jorder): ...@@ -60,7 +60,7 @@ def gD(F, s, iorder, jorder):
lambda x: -(x ** 2 - s ** 2) * e ** -(x ** 2 / (2 * s ** 2)) \ lambda x: -(x ** 2 - s ** 2) * e ** -(x ** 2 / (2 * s ** 2)) \
/ (2 * pi * s ** 6)] / (2 * pi * s ** 6)]
size = int(ceil(3 * s)) size = int(ceil(3 * s))
r = 2 * size r = 2 * size + 1
iterator = map(float, range(r)) iterator = map(float, range(r))
W = zeros((r, r)) W = zeros((r, r))
Fx = funcs[iorder] Fx = funcs[iorder]
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment