Table-and-Graph-Libs/python/graphs.py

560 lines
17 KiB
Python
Raw Normal View History

#!/usr/bin/env python3
2022-02-01 23:20:23 +08:00
# -*- coding: utf-8 -*-
# Teal Dulcet, CS546
from __future__ import division, print_function, unicode_literals
import sys
import math
import shutil
from fractions import Fraction
2022-02-01 23:20:23 +08:00
import textwrap
from enum import Enum, IntEnum, auto
2022-02-01 23:20:23 +08:00
from wcwidth import wcswidth
from typing import List, Tuple, Optional, Sequence, Callable
import locale
locale.setlocale(locale.LC_ALL, '')
class style_types(IntEnum):
ASCII = 0
basic = auto()
light = auto()
heavy = auto()
double = auto()
light_dashed = auto()
heavy_dashed = auto()
2022-02-01 23:20:23 +08:00
styles = [
["-", "|", "+", "+", "+", "+", "+", "+", "+", "+", "+"], # ASCII
["", "|", "+", "+", "+", "+", "+", "+", "+", "+", "+"], # Basic
["", "", "", "", "", "", "", "", "", "", ""], # Light
["", "", "", "", "", "", "", "", "", "", ""], # Heavy
["", "", "", "", "", "", "", "", "", "", ""], # Double
["", "", "", "", "", "", "", "", "", "", ""], # Light Dashed
["", "", "", "", "", "", "", "", "", "", ""] # Heavy Dashed
# [" ", " ", " ", " ", " ", " ", " ", " ", " ", " ", " "]] #No border
2022-02-01 23:20:23 +08:00
]
class color_types(IntEnum):
default = 0
black = auto()
red = auto()
green = auto()
yellow = auto()
blue = auto()
magenta = auto()
cyan = auto()
light_gray = auto()
dark_gray = auto()
light_red = auto()
light_green = auto()
light_yellow = auto()
light_blue = auto()
light_magenta = auto()
light_cyan = auto()
white = auto()
colors = ["\033[39m", "\033[30m", "\033[31m", "\033[32m", "\033[33m",
"\033[34m", "\033[35m", "\033[36m", "\033[37m", "\033[90m",
"\033[91m", "\033[92m", "\033[93m", "\033[94m", "\033[95m",
"\033[96m", "\033[97m"]
dots = [
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
"", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
""]
2022-02-01 23:20:23 +08:00
values = [[0x1, 0x2, 0x4, 0x40], [0x8, 0x10, 0x20, 0x80]]
fractions = {
"¼": Fraction(1, 4),
"½": Fraction(1, 2),
"¾": Fraction(3, 4),
"": Fraction(1, 7),
"": Fraction(1, 9),
"": Fraction(1, 10),
"": Fraction(1, 3),
"": Fraction(2, 3),
"": Fraction(1, 5),
"": Fraction(2, 5),
"": Fraction(3, 5),
"": Fraction(4, 5),
"": Fraction(1, 6),
"": Fraction(5, 6),
"": Fraction(1, 8),
"": Fraction(3, 8),
"": Fraction(5, 8),
"": Fraction(7, 8)
2022-02-01 23:20:23 +08:00
}
constants = {
"π": math.pi,
"e": math.e
}
MAX = sys.float_info.radix ** sys.float_info.mant_dig - 1
2022-02-01 23:20:23 +08:00
def strcol(astr: str) -> int:
2022-02-01 23:20:23 +08:00
"""Returns the number of columns that the given string would take up if printed."""
width = wcswidth(astr)
2022-02-01 23:20:23 +08:00
if width == -1:
print(
"\nError! wcswidth failed. Nonprintable wide character.",
file=sys.stderr)
2022-02-01 23:20:23 +08:00
sys.exit(1)
return width
# return len(astr)
2022-02-01 23:20:23 +08:00
def outputfraction(number: float) -> Tuple[int, str]:
2022-02-01 23:20:23 +08:00
"""Convert fractions and constants to Unicode characters"""
output = False
strm = ""
n = abs(number)
2022-02-01 23:20:23 +08:00
if n <= MAX:
# fractionpart, intpart = math.modf(number)
# fractionpart = abs(fractionpart)
intpart, fractionpart = divmod(Fraction(number).limit_denominator(), 1)
2022-02-01 23:20:23 +08:00
for fraction, value in fractions.items():
if abs(fractionpart - value) <= sys.float_info.epsilon * n:
if intpart == 0 and number < 0:
2022-02-01 23:20:23 +08:00
strm += "-"
elif intpart != 0:
strm += "{0:n}".format(intpart)
2022-02-01 23:20:23 +08:00
strm += fraction
2022-02-01 23:20:23 +08:00
output = True
break
if n > sys.float_info.epsilon and not output:
for constant, value in constants.items():
if not output and number % value <= sys.float_info.epsilon * n:
intpart = number / value
if intpart == -1:
strm += "-"
elif intpart != 1:
strm += "{0:.{prec}n}".format(intpart,
prec=sys.float_info.dig)
strm += constant
output = True
break
2022-02-01 23:20:23 +08:00
if not output:
strm += "{0:n}".format(number)
2022-02-01 23:20:23 +08:00
length = strcol(strm)
return length, strm
def graph(height: int, width: int, xmin: float, xmax: float, ymin: float, ymax: float, array: List[List[int]], border: bool = False, axis: bool = True, axislabel: bool = True, axistick: bool = True,
axisunitslabel: bool = True, style: style_types = style_types.light, title: Optional[str] = None, check: bool = True) -> int:
2022-02-01 23:20:23 +08:00
"""Output graph"""
if not array:
return 1
if height == 0:
return 1
if width == 0:
return 1
w = shutil.get_terminal_size()
aheight = height // 4
awidth = width // 2
2022-02-01 23:20:23 +08:00
if check:
if aheight > w.lines:
print("The height of the graph ({0}) is greater then the height of the terminal ({1}).".format(
aheight, w.lines), file=sys.stderr)
return 1
2022-02-01 23:20:23 +08:00
if awidth > w.columns:
print("The width of the graph ({0}) is greater then the width of the terminal ({1}).".format(
awidth, w.columns), file=sys.stderr)
return 1
2022-02-01 23:20:23 +08:00
if xmin >= xmax:
print("xmin must be less than xmax.", file=sys.stderr)
return 1
if ymin >= ymax:
print("ymin must be less than ymax.", file=sys.stderr)
return 1
xstep = (xmax - xmin) / width
ystep = (ymax - ymin) / height
xaxis = 0 if xmin > 0 else width if xmax < 0 else width - xmax / xstep
yaxis = height if ymin > 0 else 0 if ymax < 0 else ymax / ystep
xdivisor = 2 * 4 * (width // 160 + 2)
ydivisor = 2 * 4 * (height // 160 + 2)
2022-02-01 23:20:23 +08:00
if title:
print(textwrap.fill(title, width=awidth))
2022-02-01 23:20:23 +08:00
strm = ""
if border:
strm += styles[style][2]
strm += styles[style][0] * awidth
strm += styles[style][4] + "\n"
2022-02-01 23:20:23 +08:00
i = 0
while i < height:
ayaxis = i <= yaxis and i + 4 > yaxis if yaxis <= height - \
4 else i < yaxis and i + 4 >= yaxis
yaxislabel = i <= yaxis + 4 and i + 4 > yaxis + \
4 if yaxis <= height - 4 else i < yaxis - 4 and i + 4 >= yaxis
2022-02-01 23:20:23 +08:00
ylabelstrm = ""
ylabellength = 0
if axis and axislabel and axistick and axisunitslabel and yaxis >= 0 and yaxis <= height:
2022-02-01 23:20:23 +08:00
output = False
label = 0.0
adivisor = -ydivisor if i < yaxis else ydivisor
2022-02-01 23:20:23 +08:00
k = yaxis + adivisor
while (k >= i if i < yaxis else k < i + 4) and i >= 4 and not output:
if i <= k and i + 4 > k:
label = ymax - (height if k > height else k) * ystep
2022-02-01 23:20:23 +08:00
output = True
k += adivisor
if output:
ylabellength, ylabelstrm = outputfraction(label)
2022-02-01 23:20:23 +08:00
ylabellength *= 2
if border:
strm += styles[style][1]
2022-02-01 23:20:23 +08:00
j = 0
while j < width:
axaxis = j < xaxis and j + 2 >= xaxis if xaxis >= 2 else j <= xaxis and j + 2 > xaxis
xaxislabel = j < xaxis - 2 and j + 2 >= xaxis - \
2 if xaxis >= 2 else j <= xaxis + 2 and j + 2 > xaxis + 2
2022-02-01 23:20:23 +08:00
output = False
if axis:
2022-02-01 23:20:23 +08:00
if axaxis and ayaxis:
strm += styles[style][6]
output = True
elif axaxis:
if i == 0:
strm += styles[style][4]
output = True
elif i >= height - 4:
strm += styles[style][10]
output = True
elif axislabel and axistick:
adivisor = -ydivisor if i < yaxis else ydivisor
2022-02-01 23:20:23 +08:00
k = yaxis + adivisor
while (k >= i if i < yaxis else k < i + 4) and i >= 4 and not output:
if i <= k and i + 4 > k:
strm += styles[style][7 if xaxis >= 2 else 5]
2022-02-01 23:20:23 +08:00
output = True
k += adivisor
if not output:
strm += styles[style][1]
2022-02-01 23:20:23 +08:00
output = True
elif ayaxis:
if j == 0:
strm += styles[style][2]
output = True
elif j >= width - 2:
strm += styles[style][4]
output = True
elif axislabel and axistick:
adivisor = -xdivisor if j < xaxis else xdivisor
2022-02-01 23:20:23 +08:00
k = xaxis + adivisor
while (k >= j if j < xaxis else k < j + 2) and j < width - 4 and not output:
if j <= k and j + 2 > k:
strm += styles[style][3 if yaxis <=
height - 4 else 9]
2022-02-01 23:20:23 +08:00
output = True
k += adivisor
if not output:
strm += styles[style][0]
2022-02-01 23:20:23 +08:00
output = True
elif yaxislabel and xaxislabel and axislabel and axistick and axisunitslabel and ymin <= 0 and ymax >= 0 and xmin <= 0 and xmax >= 0:
2022-02-01 23:20:23 +08:00
strm += "0"
output = True
elif (j >= width - 2 if xaxis <= width - 2 else j == 0) and yaxislabel and axislabel:
2022-02-01 23:20:23 +08:00
strm += "x"
output = True
elif yaxislabel and axislabel and axistick and axisunitslabel:
2022-02-01 23:20:23 +08:00
label = 0.0
adivisor = -xdivisor if j < xaxis else xdivisor
2022-02-01 23:20:23 +08:00
if j < xaxis:
j += 2
k = xaxis + adivisor
while (k >= j if j < xaxis else k < j + 2) and j < width - 2 and not output:
if j <= k and j + 2 > k:
label = (width if k > width else k) * xstep + xmin
2022-02-01 23:20:23 +08:00
output = True
k += adivisor
if adivisor < 0:
j -= 2
if output:
output = False
length, astrm = outputfraction(label)
2022-02-01 23:20:23 +08:00
length *= 2
if (j >= xaxis or j + length < xaxis - 4) and j + length < width - 2:
2022-02-01 23:20:23 +08:00
strm += astrm
if length > 2:
j += length - 2
if adivisor < 0:
output = True
else:
j += 2
elif (i == 0 if yaxis >= 4 else i >= height - 4) and xaxislabel and axislabel:
2022-02-01 23:20:23 +08:00
strm += "y"
output = True
elif ylabellength and (xaxislabel if xaxis < 2 else j < xaxis - ylabellength and j + 2 >= xaxis - ylabellength) and (yaxis >= 4 or i < height - 4) and axislabel and axistick and axisunitslabel:
2022-02-01 23:20:23 +08:00
strm += ylabelstrm
output = True
if ylabellength > 2:
j += ylabellength - 2
if not output:
dot = 0
color = 0
for k in range(min(2, width - j)):
for l in range(min(4, height - i)):
dot += (1 if array[j + k][i + l] else 0) * values[k][l]
if color:
if array[j + k][i + l] and color != array[j + k][i + l]:
color = 1
else:
color = array[j + k][i + l]
if color:
color -= 1
if color:
strm += colors[color]
strm += dots[dot]
2022-02-01 23:20:23 +08:00
if color:
strm += colors[0]
j += 2
if border:
strm += styles[style][1]
if i < height - 4 or border:
strm += "\n"
2022-02-01 23:20:23 +08:00
i += 4
if border:
strm += styles[style][8]
strm += styles[style][0] * awidth
strm += styles[style][10]
2022-02-01 23:20:23 +08:00
print(strm)
return 0
def arrays(height: int, width: int, xmin: float, xmax: float, ymin: float, ymax: float, aarrays: Sequence[Sequence[Sequence[float]]], border: bool = False, axis: bool = True, axislabel: bool = True, axistick: bool = True, axisunitslabel: bool = True,
style: style_types = style_types.light, color: color_types = color_types.red, title: Optional[str] = None, check: bool = True) -> int:
2022-02-01 23:20:23 +08:00
"""Convert one or more arrays to graph and output"""
if not aarrays:
return 1
if not all(len(x) == 2 for aarray in aarrays for x in aarray):
print("Error: The arrays must have two columns.")
return 1
w = shutil.get_terminal_size()
if height == 0:
height = w.lines * 4
if width == 0:
width = w.columns * 2
if check:
aheight = height // 4
awidth = width // 2
2022-02-01 23:20:23 +08:00
if aheight > w.lines:
print("The height of the graph ({0}) is greater then the height of the terminal ({1}).".format(
aheight, w.lines), file=sys.stderr)
return 1
2022-02-01 23:20:23 +08:00
if awidth > w.columns:
print("The width of the graph ({0}) is greater then the width of the terminal ({1}).".format(
awidth, w.columns), file=sys.stderr)
return 1
2022-02-01 23:20:23 +08:00
if xmin == 0 and xmax == 0:
xmin = min(i[0] for aarray in aarrays for i in aarray)
xmax = max(i[0] for aarray in aarrays for i in aarray)
if ymin == 0 and ymax == 0:
ymin = min(i[1] for aarray in aarrays for i in aarray)
ymax = max(i[1] for aarray in aarrays for i in aarray)
if xmin >= xmax:
print("xmin must be less than xmax.", file=sys.stderr)
return 1
if ymin >= ymax:
print("ymin must be less than ymax.", file=sys.stderr)
return 1
xstep = (xmax - xmin) / width
ystep = (ymax - ymin) / height
xaxis = width - xmax / xstep
yaxis = ymax / ystep
2022-02-01 23:20:23 +08:00
aaarray = [[0 for j in range(height)] for i in range(width)]
for j, aarray in enumerate(aarrays):
acolor = color + 1 if len(aarrays) == 1 else j % (len(colors) - 2) + 3
2022-02-01 23:20:23 +08:00
for i in aarray:
if i[0] >= xmin and i[0] < xmax and i[1] >= ymin and i[1] < ymax:
x = int(i[0] / xstep + xaxis)
y = int(yaxis - i[1] / ystep - 1)
2022-02-01 23:20:23 +08:00
if aaarray[x][y]:
if aaarray[x][y] != acolor:
aaarray[x][y] = 1
else:
aaarray[x][y] = acolor
return graph(height, width, xmin, xmax, ymin, ymax, aaarray, border=border, axis=axis, axislabel=axislabel,
axistick=axistick, axisunitslabel=axisunitslabel, style=style, title=title)
2022-02-01 23:20:23 +08:00
def array(height: int, width: int, xmin: float, xmax: float, ymin: float, ymax: float, aarray: Sequence[Sequence[float]], border: bool = False, axis: bool = True, axislabel: bool = True, axistick: bool = True,
axisunitslabel: bool = True, style: style_types = style_types.light, color: color_types = color_types.red, title: Optional[str] = None) -> int:
2022-02-01 23:20:23 +08:00
"""Convert single array to graph and output"""
return arrays(height, width, xmin, xmax, ymin, ymax, [aarray], border=border, axis=axis, axislabel=axislabel,
axistick=axistick, axisunitslabel=axisunitslabel, style=style, color=color, title=title)
2022-02-01 23:20:23 +08:00
def functions(height: int, width: int, xmin: float, xmax: float, ymin: float, ymax: float, afunctions: Sequence[Callable[[float], float]], border: bool = False, axis: bool = True, axislabel: bool = True, axistick: bool = True,
axisunitslabel: bool = True, style: style_types = style_types.light, color: color_types = color_types.red, title: Optional[str] = None, check: bool = True) -> int:
2022-02-01 23:20:23 +08:00
"""Convert one or more functions to graph and output"""
if not afunctions:
2022-02-01 23:20:23 +08:00
return 1
w = shutil.get_terminal_size()
if height == 0:
height = w.lines * 4
if width == 0:
width = w.columns * 2
if check:
aheight = height // 4
awidth = width // 2
2022-02-01 23:20:23 +08:00
if aheight > w.lines:
print("The height of the graph ({0}) is greater then the height of the terminal ({1}).".format(
aheight, w.lines), file=sys.stderr)
return 1
2022-02-01 23:20:23 +08:00
if awidth > w.columns:
print("The width of the graph ({0}) is greater then the width of the terminal ({1}).".format(
awidth, w.columns), file=sys.stderr)
return 1
2022-02-01 23:20:23 +08:00
if xmin >= xmax:
print("xmin must be less than xmax.", file=sys.stderr)
return 1
if ymin >= ymax:
print("ymin must be less than ymax.", file=sys.stderr)
return 1
rows = width
xstep = (xmax - xmin) / width
ystep = (ymax - ymin) / height
xaxis = width - xmax / xstep
yaxis = ymax / ystep
xres = 2
2022-02-01 23:20:23 +08:00
array = [[0 for j in range(height)] for i in range(width)]
for j, function in enumerate(afunctions):
acolor = color + \
1 if len(afunctions) == 1 else j % (len(colors) - 2) + 3
2022-02-01 23:20:23 +08:00
for i in (x / xres for x in range(rows * xres)):
x = i * xstep + xmin
2022-02-01 23:20:23 +08:00
y = function(x)
if x >= xmin and x < xmax and y >= ymin and y < ymax:
ax = int(x / xstep + xaxis)
ay = int(yaxis - y / ystep - 1)
2022-02-01 23:20:23 +08:00
if array[ax][ay]:
if array[ax][ay] != acolor:
array[ax][ay] = 1
else:
array[ax][ay] = acolor
return graph(height, width, xmin, xmax, ymin, ymax, array, border=border, axis=axis, axislabel=axislabel,
axistick=axistick, axisunitslabel=axisunitslabel, style=style, title=title)
2022-02-01 23:20:23 +08:00
def function(height, width, xmin: float, xmax: float, ymin: float, ymax: float, afunction: Callable[[float], float], border: bool = False, axis: bool = True, axislabel: bool = True, axistick: bool = True,
axisunitslabel: bool = True, style: style_types = style_types.light, color: color_types = color_types.red, title: Optional[str] = None) -> int:
2022-02-01 23:20:23 +08:00
"""Convert single function to function array and output"""
return functions(height, width, xmin, xmax, ymin, ymax, [afunction], border=border, axis=axis, axislabel=axislabel,
axistick=axistick, axisunitslabel=axisunitslabel, style=style, color=color, title=title)