def plotGlass(c): xradii = numpy.linspace(0, properties.glassHeight, num=properties.interPoints) fi_height_volume, fi_radius_height, fi_side_height, fi_volume_height = integrator.interStates(c) xvolume = integrator.interFunc( fi_height_volume, integrator.totalVolume, fi_volume_height[-1] / properties.interPoints, "fi_height_volume" ) yvolume = integrator.interFunc(fi_radius_height, xvolume, properties.interDHeight, "fi_radius_height") plt.plot(xradii, fi_radius_height, "k", [xvolume], [yvolume], "ro") plt.xlabel("height / m") plt.ylabel("radius / m") plt.show()
def plotGlass(c): xradii = numpy.linspace(0, properties.glassHeight, num=properties.interPoints) fi_height_volume, fi_radius_height, fi_side_height, fi_volume_height = integrator.interStates( c) xvolume = integrator.interFunc( fi_height_volume, integrator.totalVolume, fi_volume_height[-1] / properties.interPoints, "fi_height_volume") yvolume = integrator.interFunc(fi_radius_height, xvolume, properties.interDHeight, "fi_radius_height") plt.plot(xradii, fi_radius_height, 'k', [xvolume], [yvolume], 'ro') plt.xlabel("height / m") plt.ylabel("radius / m") plt.show()
def getVolume(self): _, _, _, fi_volume_height = integrator.interStates(self) return fi_volume_height[-1]