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

536 lines
16 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
class color_types(IntEnum):
default = 0
black = auto()
red = auto()
green = auto()
yellow = auto()
blue = auto()
cyan = auto()
light_gray = auto()
dark_gray = auto()
light_red = auto()
light_green = auto()
light_yellow = auto()
light_blue = 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 = True, axislabel: 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()
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
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=width // 2))
2022-02-01 23:20:23 +08:00
strm = ""
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 - 4)
2022-02-01 23:20:23 +08:00
ylabelstrm = ""
ylabellength = 0
if border and axislabel 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:
2022-02-01 23:20:23 +08:00
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
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 border:
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 axisunitslabel:
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:
2022-02-01 23:20:23 +08:00
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 axisunitslabel:
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:
2022-02-01 23:20:23 +08:00
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 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 axisunitslabel:
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:
2022-02-01 23:20:23 +08:00
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):
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 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 i < (height - 4):
strm += "\n"
2022-02-01 23:20:23 +08:00
i += 4
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 = True, axislabel: 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,
axislabel=axislabel, 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 = True, axislabel: 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, axislabel=axislabel, 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 = True, axislabel: 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
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,
axislabel=axislabel, 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 = True, axislabel: 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, axislabel=axislabel, axisunitslabel=axisunitslabel, style=style, color=color, title=title)