for i in xrange(len(xs)):
xy = np.array((xs[i], ys[i]), dtype=np.double)
uv = timoBeam.displacement(xy[0], xy[1], 0, None)
if is_vs:
exp_uvs.append(uv[1])
else:
exp_uvs.append(uv[0])
if None is axis:
axis = xs
ax.plot(axis, act_uvs, "r", label="Numerical")
ax.plot(axis, exp_uvs, "b", label="Precise")
if __name__ == "__main__":
WeakformProcessor2D = JClass("net.epsilony.tsmf.process.WeakformProcessor2D")
processor = WeakformProcessor2D.genTimoshenkoProjectProcess()
processor.process()
processor.solve()
pp = processor.postProcessor()
fig = plt.figure()
ax = fig.add_subplot(221)
ax.set_title("displacement $v$ along axis $x$ ($y=0$)")
ax.set_xlabel("$x$")
ax.set_ylabel("$v(x,0)$")
y = 0
plot_vs_by_y(pp, ax, y)
ax2 = fig.add_subplot(222)
ax2.set_title("displacement $u$\n along the left edge ($x=0$)")
ax2.set_xlabel("$y$")
