def __init__(self, function, x0, y0, X, n):
self.x = []
self.y = []
self.n = n
self.x0 = x0
self.X = X
self.c1 = ((1 / math_exp(x0) * y0) - 1) / math_exp(x0)
self.function = function
def __init__(self, alpha, mu):
self.f = lambda x_, y_: x_ ** alpha * math_exp(mu * (1 - y_))
self.g = lambda x_, y_: x_
self.f_symbolic = x ** alpha * sympy_exp(mu * (1 - y))
self.g_symbolic = x
self.alpha = alpha
self.mu = mu
def exp(self):
v = Quaternion(0, self.b, self.c, self.d)
norm_v = v.norm()
unit_v = v.unit()
rhs = unit_v * sin(norm_v) + cos(norm_v)
return rhs * math_exp(self.a)
def probability_from_log_odds(log_odds: float) -> float:
"""
Returns a probability from a log odds.
Args:
log_odds: ln(o); natural log of o
Returns:
float: p
"""
return probability_from_odds(math_exp(log_odds))
def probability_from_log_prob(log_p: float) -> float:
"""
Returns a probability from a log probability.
Args:
log_p: ln(p); natural log of p
Returns:
float: p
"""
return math_exp(log_p)
def exp(other):
if hasattr(other, 'shape'):
new = [0] * other.shape[0]
for i, line in enumerate(other):
new[i] = map(math_exp, line)
return Matrix(new)
if is_math(other):
return math_exp(other)
if is_iter(other):
new_ = list()
for item in other:
new_.append(exp(item))
return new_
raise TypeError('expects an iterable or numeric for exp(), got %s' %
type(other))
def exp(v):
if v == None:
return Null
return math_exp(v)
def qexp(a):
number = math_exp(a.number)
return Q(number, "exp0(%s)", a.units, abs(a.uncert * number), a.prefu, (a,))
def _update_callback(self):
# print("-***********************************")
self.function_graph.clear()
self.fig_errors.clear()
self.controller.set_parameters(self.initial_value_x.get(),
self.initial_value_y.get(),
self.max_X.get(), self.steps_n.get())
obj = self.controller.compute()
axis_function = self.function_graph.add_subplot(111)
axis_error = self.fig_errors.add_subplot(111)
plots = [[None, None, None, None],
["Exact", "Euler", "Improved Euler", "Runge Kutte"]]
# print(len(obj["Exact"][0]), len(obj["Exact"][1]))
if (self.exact_checkbox_val.get()):
plots[0][0] = axis_function.plot(obj["Function"]["Exact"][0],
obj["Function"]["Exact"][1],
label="Exact",
color='green')
if (self.euler_checkbox_val.get()):
plots[0][1] = axis_function.plot(obj["Function"]["Euler"][0],
obj["Function"]["Euler"][1],
label="Euler",
color='red')
# print(len(obj["LocalError"]["Euler"][0]),len(obj["LocalError"]["Euler"][1]))
if (self.local_global_error_checkbox_val.get()):
axis_error.plot(obj["GlobalError"]["Euler"][0],
obj["GlobalError"]["Euler"][1],
color='red')
else:
axis_error.plot(obj["LocalError"]["Euler"][0],
obj["LocalError"]["Euler"][1],
color='red')
if (self.improved_euler_checkbox_val.get()):
plots[0][2] = axis_function.plot(
obj["Function"]["ImprovedEuler"][0],
obj["Function"]["ImprovedEuler"][1],
label="Improved Euler",
color='blue')
if (self.local_global_error_checkbox_val.get()):
axis_error.plot(obj["GlobalError"]["ImprovedEuler"][0],
obj["GlobalError"]["ImprovedEuler"][1],
color='blue')
else:
axis_error.plot(obj["LocalError"]["ImprovedEuler"][0],
obj["LocalError"]["ImprovedEuler"][1],
color='blue')
if (self.runge_kutte_checkbox_val.get()):
plots[0][3] = axis_function.plot(obj["Function"]["RungeKutte"][0],
obj["Function"]["RungeKutte"][1],
label="Runge Kutte",
color='yellow')
if (self.local_global_error_checkbox_val.get()):
axis_error.plot(obj["GlobalError"]["RungeKutte"][0],
obj["GlobalError"]["RungeKutte"][1],
color='yellow')
else:
axis_error.plot(obj["LocalError"]["RungeKutte"][0],
obj["LocalError"]["RungeKutte"][1],
color='yellow')
# self.fig_errors.legend(handless=[])
available_plots = [[], []]
for i in range(4):
if (plots[0][i] is not None):
available_plots[0].append(plots[0][i][0])
available_plots[1].append(plots[1][i])
# print(available_plots[0])
# print(available_plots[1])
if (len(available_plots[0]) > 0):
axis_function.legend(handles=available_plots[0],
labels=available_plots[1],
loc='upper right')
axis_error.legend(handles=available_plots[0],
labels=available_plots[1],
loc='upper right')
self.canvas_function_graph.draw()
self.canvas_error_graph.draw()
self.label.config(
text="Where c1 = {:}".format((
(1 / math_exp(self.initial_value_x.get()) *
self.initial_value_y.get()) - 1) /
math_exp(self.initial_value_x.get())))
def __init__(self,
root,
controller,
x0=-4,
y0=1,
X=4,
n=100,
gui_app_title="ODE_GUI"):
# internal parameters
self.root = root
self.update_flag = 0
self.gui_app_geometry = "1000x700"
self.gui_app_title = gui_app_title
self.initial_value_x = DoubleVar(value=x0)
self.initial_value_y = DoubleVar(value=y0)
self.steps_n = DoubleVar(value=n)
self.max_X = DoubleVar(value=X)
self.exact_checkbox_val = IntVar()
self.euler_checkbox_val = IntVar()
self.improved_euler_checkbox_val = IntVar()
self.runge_kutte_checkbox_val = IntVar()
self.local_global_error_checkbox_val = IntVar()
self.controller = controller
# GUI construction
self.root.geometry(self.gui_app_geometry)
self.root.title(self.gui_app_title)
self.function_graph_label = Label(self.root,
text="Function graph Figure")
self.function_graph_label.place(x=170, y=30)
self.function_graph = Figure(figsize=(4.75, 5), dpi=100)
# t = np.arange(0, 3, .01)
# self.function_graph.add_subplot(111).plot(t, 2 * np.sin(2 * np.pi * t))
self.canvas_function_graph = FigureCanvasTkAgg(
self.function_graph, master=self.root) # A tk.DrawingArea.
self.canvas_function_graph.draw()
self.canvas_function_graph.get_tk_widget().place(x=10, y=50)
self.toolbar2 = NavigationToolbar2Tk(self.canvas_function_graph,
self.root)
self.toolbar2.place(x=5, y=660)
self.canvas_function_graph.get_tk_widget().place(x=10, y=50)
self.canvas_function_graph.mpl_connect("key_press_event",
self.on_key_press)
self.errors_label = Label(self.root, text="Errors Figure")
self.errors_label.place(x=700, y=30)
self.fig_errors = Figure(figsize=(4.75, 5), dpi=100)
# t = np.arange(0, 3, .01)
# self.fig_errors.add_subplot(111).plot(t, 2 * np.sin(2 * np.pi * t))
self.canvas_error_graph = FigureCanvasTkAgg(
self.fig_errors, master=self.root) # A tk.DrawingArea.
self.canvas_error_graph.draw()
self.canvas_error_graph.get_tk_widget().place(x=515, y=50)
self.toolbar = NavigationToolbar2Tk(self.canvas_error_graph, self.root)
self.toolbar.update()
self.toolbar.place(x=750, y=660)
self.canvas_error_graph.get_tk_widget().place(x=515, y=50)
self.canvas_error_graph.mpl_connect("key_press_event",
self.on_key_press)
self.label = Label(self.root, text="Initial Value (x0)", wraplength=50)
self.label.place(x=15, y=560)
self.initial_value_entry = Entry(self.root,
textvariable=self.initial_value_x,
width=5)
self.initial_value_entry.place(x=20, y=620)
self.label = Label(self.root, text="Initial Value (y0)", wraplength=50)
self.label.place(x=75, y=560)
self.initial_value_entry = Entry(self.root,
textvariable=self.initial_value_y,
width=5)
self.initial_value_entry.place(x=70, y=620)
self.label = Label(self.root, text="Steps (n)", wraplength=50)
self.label.place(x=135, y=560)
self.initial_value_entry = Entry(self.root,
textvariable=self.steps_n,
width=5)
self.initial_value_entry.place(x=130, y=620)
self.label = Label(self.root, text="Max (X)", wraplength=50)
self.label.place(x=195, y=560)
self.initial_value_entry = Entry(self.root,
textvariable=self.max_X,
width=5)
self.initial_value_entry.place(x=190, y=620)
# self.label = Label(self.root, text="Exact", wraplength=50)
# self.label.place(x=700, y=560)
self.exact_checkbox = Checkbutton(self.root,
variable=self.exact_checkbox_val,
text="Exact",
command=self._update_callback)
self.exact_checkbox.place(x=700, y=580)
# self.label = Label(self.root, text="Euler", wraplength=50)
# self.label.place(x=760, y=560)
self.euler_checkbox = Checkbutton(self.root,
text="Euler",
variable=self.euler_checkbox_val,
command=self._update_callback)
self.euler_checkbox.place(x=760, y=580)
# self.label = Label(self.root, text="Improved Euler", wraplength=70)
# self.label.place(x=840, y=560)
self.improved_euler_checkbox = Checkbutton(
self.root,
text="Improved Euler",
variable=self.improved_euler_checkbox_val,
command=self._update_callback)
self.improved_euler_checkbox.place(x=830, y=580)
# self.label = Label(self.root, text="Runge Kutte", wraplength=50)
# self.label.place(x=920, y=560)
self.runge_kutte_checkbox = Checkbutton(
self.root,
text="Runge Kutte",
variable=self.runge_kutte_checkbox_val,
command=self._update_callback)
self.runge_kutte_checkbox.place(x=890, y=580)
self.local_global_error_checkbox = Checkbutton(
self.root,
text="Global",
variable=self.local_global_error_checkbox_val,
command=self._update_callback)
self.local_global_error_checkbox.place(x=930, y=580)
self.update_button = Button(self.root,
text="Update",
command=self._update_callback)
self.update_button.place(x=820, y=620)
self.close_button = Button(self.root, text="Close", command=self._quit)
self.close_button.place(x=920, y=620)
self.label = Label(self.root, text="Equation")
self.label.place(x=590, y=580)
img = ImageTk.PhotoImage(Image.open("assets/equation.png"))
panel = Label(root, image=img)
panel.image = img
panel.place(x=550, y=610)
self.label = Label(self.root, text="Exact Solution")
self.label.place(x=380, y=580)
img = ImageTk.PhotoImage(Image.open("assets/solution.png"))
panel = Label(root, image=img)
panel.image = img
panel.place(x=370, y=600)
self.label = Label(self.root,
text="Where c1 = {:}".format(
((1 / math_exp(x0) * y0) - 1) / math_exp(x0)))
self.label.place(x=350, y=670)
self._update_callback()
def qexp(a):
number = math_exp(a.number)
return Q(number, "exp(%s)", a.units, abs(a.uncert * number), a.prefu, (a,))
def test_prim_exp(x):
return math_exp(x)
def exp(x):
try:
return math_exp(x)
except OverflowError:
return INF
def exp(x):
try:
return math_exp(x)
except OverflowError:
return inf
