Skip to content

Latest commit

 

History

History
1291 lines (836 loc) · 43.6 KB

PYTHON-README.md

File metadata and controls

1291 lines (836 loc) · 43.6 KB

math-as-code (Python version)

Chinese translation of javascript version (中文版)
JavaScript version (English)

This is a reference to ease developers into mathematical notation by showing comparisons with Python code.

Motivation: Academic papers can be intimidating for self-taught graphics programmers and data wranglers :)

This guide is not yet finished. If you see errors or want to contribute, please open a ticket or send a PR.

Note: There will be some discussion of vectorization with numpy.

Note: We're only looking at python3, if there are significant differences with python2 please bring them up in issues.

foreword

Mathematical symbols can mean different things depending on the author, context and the field of study (linear algebra, set theory, etc). This guide may not cover all uses of a symbol. In some cases, real-world references (blog posts, publications, etc) will be cited to demonstrate how a symbol might appear in the wild.

For a more complete list, refer to Wikipedia - List of Mathematical Symbols.

For simplicity, many of the code examples here operate on floating point values and are not numerically robust. For more details on why this may be a problem, see Robust Arithmetic Notes by Mikola Lysenko.

contents

variable name conventions

There are a variety of naming conventions depending on the context and field of study, and they are not always consistent. However, in some of the literature you may find variable names to follow a pattern like so:

  • s - italic lowercase letters for scalars (e.g. a number)
  • x - bold lowercase letters for vectors (e.g. a 2D point)
  • A - bold uppercase letters for matrices (e.g. a 3D transformation)
  • θ - italic lowercase Greek letters for constants and special variables (e.g. polar angle θ, theta)

This will also be the format of this guide.

Numpy

Numpy is a powerful array programming library, which in python can be interpreted as a domain specific language (DSL). Sometimes it's helpful to think of math in python as two languages sharing a namespace, with special syntax sugar to access one or the other. This will be important in our vectors and matrices section, because slightly different python syntax means different speeds at large input. The convention is import numpy as np, so when you see np.something you know that we're working in numpy.

equals symbols

There are a number of symbols resembling the equals sign =. Here are a few common examples:

  • = is for equality (values are the same)
  • is for inequality (value are not the same)
  • is for approximately equal to (π ≈ 3.14159)
  • := is for definition (A is defined as B)

In Python:

## equality
2 == 3

## inequality
2 != 3

## approximately equal
import math
math.isclose(math.pi, 3.14159) # math.isclose doesn't have a third argument for tolerance, so this is false

from numpy.testing import assert_almost_equal
assert_almost_equal(math.pi, 3.14159, 1e-5) # we gave it a the tolerance we want, 5 decimal places. 
# This is actually a unit test, equivalent to "assert isclose(x,y)", read on for more. 

def almost_equal(x, y, epsilon=7): 
  ''' you can make your own! 
  in numpy, 1e-7 is the default epsilon
  '''
  return abs(x - y) < 10 ** -epsilon

Read more: programmers got this idea from the epsilon-delta definition of limit

Note: subclasses of unittest.TestCase come with their own assertAlmostEqual.

Warning: please don't use exact == equality on floats!

In mathematical notation, you might see the :=, =: and = symbols being used for definition.2

For example, the following defines x to be another name for 2kj.

equals1

In python, we define our variables and provide aliases with =.

x = 2 * k * j

Assignment in python is generally mutable besides special cases like Tuple.

Note: Some languages have pre-processor #define statements, which are closer to a mathematical define.

Notice that def is a form of := as well.

def plus(x, y): 
  return x + y

The following, on the other hand, represents equality:

equals2

Important: the difference between = and == can be more obvious in code than it is in math literature! In python, a = is an instruction. You're telling the machine to interact with the namespace, add something to it or change something in it. In python, when you write == you're asking the machine "may I have a bool?". In math, the former case is either covered by := or =, while the latter case is usually =, and you might have to do some disambiguating in your reading.

In math, when I write 1 + 1 = 2 I'm making a judgment. It's not that i'm asking the world (or the chalkboard) for a bool, it's that I'm keeping track of my beliefs. This distinction is the foundation of unit tests or assertions.

# assert in python takes an expression that lands in bool and a string to be printed if it turns out false. 
assert plus(1, 1) == 2, "DANGER: PHYSICS IS BROKEN. PLEASE STAY INSIDE. "

It's important to know when a falsehood ought to crash a program vs. when you just want a boolean value. To understand this better, read this.

square root and complex numbers

A square root operation is of the form:

squareroot

In programming we use a sqrt function, like so:

import math
print(math.sqrt(2))
# Out: 1.4142135623730951

import numpy as np
print(np.sqrt(2))
# Out: 1.4142135623730951

Complex numbers are expressions of the form complex, where a is the real part and b is the imaginary part. The imaginary number i is defined as:

imaginary.

Vanilla python has a complex constructor, and a standard module cmath for working with them.

complex(1,1)
# Out: (1+1j)s

math.sqrt(complex(1,1))
# TypeError: can't convert complex to float

import cmath
cmath.sqrt(complex(1,1))
# Out: (1.09868411346781+0.45508986056222733j)

# you need numpy to get conjugate
np.conj(complex(0.5,0.5))
# Out: (0.5-0.5j)

# we can represent the basic meaning of the imaginary unit like so
assert cmath.sqrt(complex(-1, 0)) == complex(0,1)

As you can see, it uses j to denote the imaginary unit, instead of i.

The conjugate of a complex number is flipping the sign of the imaginary part.

If z is a python complex number, z.real gets the real part (exactly as an object attribute) and z.imag gets the imaginary part.

Just as complex numbers can be interpreted as a sort of wrapper around tuples of reals, a complex number data type wraps two floats. Numpy uses this to implement complex numbers of different sizes/precisions.

The syntax is close enough to cmath, but it comes with the power and convenience of numpy. Importantly, other numpy methods are better at casting to and from complex.

observe the following cube roots of unity

z1 = 0.5 * np.complex(-1, math.sqrt(3)) # Numpy's constructor is basically the same.  
z2 = np.conj(z1) # but numpy gives us a conjugation function, while the standard module does not. 

assert math.isclose(z1**3, z2**3)
# TypeError: can't convert complex to float

np.testing.assert_almost_equal(z1**3, z2**3)

Read on about numpy's complex numbers

dot & cross

The dot · and cross × symbols have different uses depending on context.

They might seem obvious, but it's important to understand the subtle differences before we continue into other sections.

scalar multiplication

Both symbols can represent simple multiplication of scalars. The following are equivalent:

dotcross1

In programming languages we tend to use asterisk for multiplication:

result = 5 * 4

Often, the multiplication sign is only used to avoid ambiguity (e.g. between two numbers). Here, we can omit it entirely:

dotcross2

If these variables represent scalars, the code would be:

result = 3 * k * j

vector multiplication

To denote multiplication of one vector with a scalar, or element-wise multiplication of a vector with another vector, we typically do not use the dot · or cross × symbols. These have different meanings in linear algebra, discussed shortly.

Let's take our earlier example but apply it to vectors. For element-wise vector multiplication, you might see an open dot to represent the Hadamard product.2

dotcross3

In other instances, the author might explicitly define a different notation, such as a circled dot or a filled circle .3

Here is how it would look in code, using arrays [x, y] to represent the 2D vectors.

s = 3
k = [1, 2]
j = [2, 3]

tmp = multiply(k, j)
result = multiply_scalar(tmp, s)
# Out: [6, 18]

Our multiply and multiply_scalar functions look like this:

def multiply(a, b):
  return [aa * bb for aa,bb in zip(a,b)


def multiply_scalar(scalar, a):
  return [scalar * aa for aa in a]

Similarly, matrix multiplication typically does not use the dot · or cross symbol ×.

Numpy's broadcasted syntax for scaling looks like this:

def multiply_scalar(scalar, a): 
  return scalar * np.array(a)

dot product

The dot symbol · can be used to denote the dot product of two vectors. Sometimes this is called the scalar product since it evaluates to a scalar.

dotcross4

It is a very common feature of linear algebra, and with a 3D vector it might look like this:

k = [0, 1, 0]
j = [1, 0, 0]

d = np.dot(k, j)
# Out: 0

The result 0 tells us our vectors are perpendicular. Here is a dot function for 3-component vectors:

def dot(a, b):
  return a[0] * b[0] + a[1] * b[1] + a[2] * b[2]

cross product

The cross symbol × can be used to denote the cross product of two vectors.

dotcross5

In code, it would look like this:

k = [0, 1, 0]
j = [1, 0, 0]

result = cross(k, j)
# Out: [ 0, 0, -1 ]

Here, we get [0, 0, -1], which is perpendicular to both k and j.

Our cross function:

def cross(a, b): 
  ''' take two 3D vectors and return their cross product. '''
  rx = a[1] * b[2] - a[2] * b[1]
  ry = a[2] * b[0] - a[0] * b[2]
  rz = a[0] * b[1] - a[1] * b[0]
  return rx, ry, rz

It's good to practice and grok these operations, but in real life you'll use Numpy.

sigma

The big Greek Σ (Sigma) is for Summation. In other words: summing up some numbers.

sigma

Here, i=1 says to start at 1 and end at the number above the Sigma, 100. These are the lower and upper bounds, respectively. The i to the right of the "E" tells us what we are summing. In code:

Hence, the big sigma is the for keyword.

sum([k for k in range(100)])
# Out: 5050

Tip: With whole numbers, this particular pattern can be optimized to the following (and try to grok the proof. The legend of how Gauss discovered I can only describe as "typical programmer antics"):

def sum_to_n(n):
  ''' return the sum of integers from 0 to n'''
  return 0.5 * n * (n + 1)

Here is another example where the i, or the "what to sum," is different:

sum2

In code:

sum([2*k + 1 for k in range(100)])
# Out: 10000

important: range in python has an inclusive lower bound and exclusive upper bound, meaning that ... for k in range(100) is equivalent to the sum of ... for k=0 to k=n.

If you're still not absolutely fluent in indexing for these applications, spend some time with Trev Tutor on youtube.

The notation can be nested, which is much like nesting a for loop. You should evaluate the right-most sigma first, unless the author has enclosed them in parentheses to alter the order. However, in the following case, since we are dealing with finite sums, the order does not matter.

sigma3

In code:

sum(
  [sum([3*i*j 
        for j 
        in range(4,7)]) 
   for i 
   in range(1,3)])
# Out: 135

capital Pi

The capital Pi or "Big Pi" is very similar to Sigma, except we are using multiplication to find the product of a sequence of values.

Take the following:

capitalPi

This was removed from vanilla python for python 3, but it's easy to recover with a generalization of the list accumulator.

def times(x, y): 
  ''' first, give a name to the multiplication operator '''
  return x * y

from functools import reduce

reduce(times, range(1,7))
# Out: 720

With reduce, you can actually repeatedly apply a binary function to items of a list and accumulate the value for any binary operator. Python gives and and or out of the box like sum, but keep reduce in mind if you encounter a less common binary operator out in the wild.

Note that in Numpy arrays, the syntax is different (and product is given out of the box)

import numpy as np

xs = np.array([2*k + 1 for k in range(100)])
ys = np.array(range(1,7))

xs.sum()
# Out: 10000

ys.prod()
# Out: 720

which is better on larger input, but you're always welcome to use functions for ordinary lists as you please.

pipes

Pipe symbols, known as bars, can mean different things depending on the context. Below are three common uses: absolute value, Euclidean norm, and determinant.

These three features all describe the length of an object.

absolute value

pipes1

For a number x, |x| means the absolute value of x. In code:

x = -5
abs(x)
# Out: 5

Euclidean norm

pipes4

For a vector v, ‖v‖ is the Euclidean norm of v. It is also referred to as the "magnitude" or "length" of a vector.

Often this is represented by double-bars to avoid ambiguity with the absolute value notation, but sometimes you may see it with single bars:

pipes2

Here is an example using an array [x, y, z] to represent a 3D vector.

v = [0, 4, -3]
length(v)
# Out: 5.0

The `length** function:

def length(vec):
  x = vec[0]
  y = vec[1]
  z = vec[2]
  return math.sqrt(x**2 + y**2 + z**2)

The implementation for arbitrary length'd vectors is left as an exercise for the reader.

In practice, you'll probably use the following numpy call

np.linalg.norm([0, 4, -3])
# Out: 5.0

Resources:

  • [numpy.linalg docs](get link to numpy.linalg docs)

determinant

pipes3

For a matrix A, |A| means the determinant of matrix A.

Here is an example computing the determinant of a 2x2 identity matrix

ident_2 = [[1, 0], 
           [0, 1]]

np.linalg.det(ident_2)
# Out: 1

You should watch 3blue1brown, but in short if a matrix (list of list of numbers) is interpreted as hitting a coordinate system with a squisher-stretcher-rotater, the determinant of that matrix is the measure of how much the unit area/volume of the coordinate system got squished-stretched-rotated.

np.linalg.det(np.identity(100)) # the determinant of the 100 x 100 identity matrix is still one, because the identity matrix doesn't squish, stretch, or rotate at all. 
# Out: 1.0

np.linalg.det(np.array([[0, -1], [1, 0]])) # 90 degree rotation. 
# Out: 1.0

The second matrix was the 2D rotation at 90 degrees.

hat

In geometry, the "hat" symbol above a character is used to represent a unit vector. For example, here is the unit vector of a:

hat

In Cartesian space, a unit vector is typically length 1. That means each part of the vector will be in the range of -1.0 to 1.0. Here we normalize a 3D vector into a unit vector:

a = [ 0, 4, -3 ]
normalize(a)
# Out: [ 0, 0.8, -0.6 ]

If a vector is that which has magnitude and direction, normalization of a vector is the operation that deletes magnitude and preserves direction.

Here is the normalize function, operating on 3D vectors:

def normalize(vec):
  x = vec[0]
  y = vec[1]
  z = vec[2]
  squaredLength = x * x + y * y + z * z

  if (squaredLength > 0):
    length = math.sqrt(squaredLength)
    vec[0] = x / length
    vec[1] = y / length
    vec[2] = z / length
  
  return vec

Which Numpy's broadcasting syntax sugar can do in fewer lines

You should try to generalize this to vectors of arbitrary length yourself, before reading this...

Go, I mean it!

def normalize(vec):
  ''' *sigh* if you don't do it yourself you'll never learn! '''
  vec = np.array(vec) # ensure that input is casted to numpy
  length = np.linalg.norm(vec)
  if length > 0:
    return vec / length

Notice that broadcasting here is just short for [x / length for x in vec]. But it's actually faster on large input, because arrays.

Read the Numpy docs. BE the Numpy docs

element

In set theory, the "element of" symbol and can be used to describe whether something is an element of a set. For example:

element1

Here we have a set of numbers A = { 3, 9, 14 } and we are saying 3 is an "element of" that set.

The in keyword plays the role of the elementhood function, giving a bool.

A = [ 3, 9, 14 ]

3 in A
# Out: True

Python also has set. You can wrap any iterable or generator with the set keyword to delete repeats.

set([3,3,3,2,4,3,3,3,1,2,4,5,3])
# Out: {1, 2, 3, 4, 5}

3 in set(range(1, 20, 4))
# Out: False

The backwards is the same, but the order changes:

element2

You can also use the "not an element of" symbols and like so:

element3

Which you know is represented by the convenient not keyword in python.

common number sets

You may see some some large Blackboard letters among equations. Often, these are used to describe sets.

For example, we might describe k to be an element of the set .

real

Listed below are a few common sets and their symbols.

real numbers

The large describes the set of real numbers. These include integers, as well as rational and irrational numbers.

Computers approximate with float.

You can use isinstance to check "k ∈ ℝ", where float and aren't really the same thing but the intuition is close enough.

isinstance(np.pi, float)
# Out: True

Again, you may elevate that bool to an assertion that makes-or-breaks the whole program with the assert keyword when you see fit.

Excellent resource on floats in python

rational numbers

Rational numbers are real numbers that can be expressed as a fraction, or ratio. Rational numbers cannot have zero as a denominator.

Imagine taking and removing radicals (like np.sqrt) and logarithms (in a family called transcendentals), that's basically what is, at least enough for a rough first approximation.

This also means that all integers are rational numbers, since the denominator can be expressed as 1.

An irrational number, on the other hand, is one that cannot be expressed as a ratio, like π (math.pi).

A reason a programmer might care about the difference between Q and R is in the design of unit tests--- fractions are terminating decimals, and sometimes when you're a 100% sure that a number will be a basic rational (like counting change, 0.25, 0.10, 0.05, etc.), you're allowed to use == in unit tests rather than isclose or assert_almost_equal. The point is that you know not to use exact equality == when anything like sqrt or log is involved!

You can work with rationals without dividing them into floatiness with the fractions standard module

integers

An integer is a whole number. Just imagine starting from zero and one and building out an inventory with addition and subtraction.

An integer has no division, no decimals.

assert isinstance(8/7, int), "GO DIRECTLY TO JAIL"

natural numbers

A natural number, a non-negative integer.

This is actually the only set invented by the flying spaghetti monster: as for the others, humans have themselves to blame.

Depending on the context and field of study, the set may or may not start with zero.

...ok but, between you and me, they 200% start with zero.

also happens to be the first inductive construction in the study of datatypes, consisting of a single axiom ("Zero is a ") and a single inference rule ("if n is a then n + 1 is also a ")

is not a datatype in python, we can't use typechecking to disambiguate int from non-negative int, but in a pinch you could easily write up something that combines x >= 0 judgments with isinstance(x, int).

complex numbers

As we saw earlier, the complex numbers are a particular wrapper around tuples of reals.

A complex number is a combination of a real and imaginary number, viewed as a co-ordinate in the 2D plane. For more info, see A Visual, Intuitive Guide to Imaginary Numbers.

We can say ℂ = {a + b*i | a,b ∈ ℝ}, which is a notation called

Set builder notation

Pythoners have a name for set builder notation; and the name is comprehension

  • { }: delimiter around iterable (curlybois for dict or set, [ for list)
  • a + b * i: an expression (for instance, earlier when we made a list of odd numbers this expression was 2*k + 1) to be evaluated for each item in source list.
  • |: for
  • a,b ∈ ℝ: this just shows that a,b are drawn from a particular place, in this case the real numbers.

So if you've been writing Python listcomps, that definition of the complex numbers wasn't so bad! Say it with me this time

`ℂ = {a + b*i | a,b ∈ ℝ}``

inhaaaaaaless unison "C IS THE SET OF a + b*i FOR REAL NUMBERS a AND b"

If you want, you can draw up a grainy picture of an interval of ℂ with zip and np.linspace, and of course list comprehension.

j = np.complex(0,1)

R = np.linspace(-2, 2, 100)

{a + b * j for a,b in zip(R, R)}
# too much to print but try it yourself. 

functions

Functions are fundamental features of mathematics, and the concept is fairly easy to translate into code.

A function transforms an input into an output value. For example, the following is a function:

function1

We can give this function a name. Commonly, we use ƒ to describe a function, but it could be named A or anything else.

function2

In code, we might name it square and write it like this:

def square(x): 
  return math.pow(x, 2)

Sometimes a function is not named, and instead the output is written.

function3

In the above example, x is the input, the transformation is squaring, and y is the output. We can express this as an equation because, conventionally, we think of x as input and y as output.

But we have a stronger idea called anonymous functions to generalize this.

Just as we can name strings x = "Alonzo" then call them with their names or we can just pass string literals, we also have function literals.

Math first, then python:

x ↦ x^2 is equivalent to the equational description above.

Nearly identical, but very different to the untrained eye, is λx.x^2, hence the python keyword

lambda x: x**2

Functions can also have multiple parameters, like in a programming language. These are known as arguments in mathematics, and the number of arguments a function takes is known as the arity of the function.

function4

dictionaries are functions

Sometimes mathematicians, like software developers, need to specify maps by enumerating each input-output pair* when there is no expression that computes output from input.

Note: formally, mathematicians require that functions not be ambiguous*, so when you have a function and you have an input, there can be no uncertainty as to what the output should be; you mustn't be confused about whether an apple is red or purple (in introductory algebra courses this is called the "vertical line test", but it applies to all maps). Notice that the implementation of hash maps already guarantees this in the case of dictionaries! Notice also that we make no such requirement on outputs, both an apple and a banana can land on purple! With caveats like these, we can study the properties of different kinds of functions into different kinds, important in compression and security engineering.

piecewise function

Some functions will use different relationships depending on the input value, x.

The following function ƒ chooses between two "sub functions" depending on the input value.

piecewise1

This is very similar to if / else in code. The right-side conditions are often written as "for x < 0" or "if x = 0". If the condition is true, the function to the left is used.

In piecewise functions, "otherwise" and "elsewhere" are analogous to the else statement in code.

def f(x):
  if (x >= 1):
    return (math.pow(x, 2) - x) / x
  else:
    return 0

common functions

There are some function names that are ubiquitous in mathematics. For a programmer, these might be analogous to functions "built-in" to the language (like parseInt in JavaScript).

One such example is the sgn function. This is the signum or sign function. Let's use piecewise function notation to describe it:

sgn

In code, it might look like this:

def signum(x):
  if (x < 0):
    return -1
  elif (x > 0):
    return 1
  else: 
    return 0

See signum for this function as a module.

Other examples of such functions: sin, cos, tan.

function notation

In some literature, functions may be defined with more explicit notation. For example, let's go back to the square function we mentioned earlier:

function2

It might also be written in the following form:

mapsto

The arrow here with a tail typically means "maps to," as in x maps to x2.

Sometimes, when it isn't obvious, the notation will also describe the domain and codomain of the function. A more formal definition of ƒ might be written as:

funcnot

A function's domain and codomain is a bit like its input and output types, respectively. Here's another example, using our earlier sgn function, which outputs an integer:

domain2

The arrow here (without a tail) is used to map one set to another.

In Python and other dynamically typed languages, you might use documentation and/or runtime checks to explain and validate a function's input/output. Example:

def square_ints(k): 
  ''' FEED ME INTEGER '''
  try: 
    assert isinstance(k, int), "I HUNGER FOR AN INTEGER! "
    return math.pow(k, 2)
  except AssertionError as e:
    raise e

The python of a more glorious future as described in pep484 proposes a static type checker for Python, but no one's proposed anything shrewd enough to prevent code with type errors from compiling for Python yet.

As we will see in the following sample of pep484 Python, the set interpretation of domain and codomain makes way for a types interpretation of domain and codomain

def square(x: float) -> float: 
  ''' a pep484 annotation isn't that different from if i declared in the docstring; 
  
  input/domain: a float
  output/codomain: another float 
  '''
  return x**2

Other languages, like Java, allow for true method overloading based on the static types of a function's input/output. This is closer to mathematics: two functions are not the same if they use a different domain. This is also called polymorphism and it explains why 'literally' + 'alonzo' concats two strings together but 1 + 1 is addition on numbers.

prime

The prime symbol () is often used in variable names to describe things which are similar, without giving it a different name altogether. It can describe the "next value" after some transformation.

For example, if we take a 2D point (x, y) and rotate it, you might name the result (x′, y′). Or, the transpose of matrix M might be named M′.

In code, we typically just assign the variable a more descriptive name, like transformedPosition.

For a mathematical function, the prime symbol often describes the derivative of that function. Derivatives will be explained in a future section. Let's take our earlier function:

function2

Its derivative could be written with a prime symbol:

prime1

In code:

def f(x): 
  return Math.pow(x, 2)

def f_prime(x):
  return 2 * x

Multiple prime symbols can be used to describe the second derivative ƒ′′ and third derivative ƒ′′′. After this, authors typically express higher orders with roman numerals ƒIV or superscript numbers ƒ(n).

floor & ceiling

The special brackets ⌊x⌋ and ⌈x⌉ represent the floor and ceil functions, respectively.

floor

ceil

In code:

math.floor(4.8)
math.ceil(3.1)
np.floor(4.9)
np.ceil(3.001)

Note: the Numpy version returns a float, in the above example 4.0, rather than the int 4

When the two symbols are mixed ⌊x⌉, it typically represents a function that rounds to the nearest integer:

round

Python automatically gives you a keyword round to call on a number.

arrows

Arrows are often used in function notation. Here are a few other areas you might see them.

material implication

Arrows like and are sometimes used in logic for material implication. That is, if A is true, then B is also true.

material1

Interpreting this as code might look like this:

def if_A_then_B:
  if A:
    assert B, "alas, not A!"
    return B

The arrows can go in either direction , or both . When A ⇒ B and B ⇒ A, they are said to be equivalent:

material-equiv

inequality

In math, the < > and are typically used in the same way we use them in code: less than, greater than, less than or equal to and greater than or equal to, respectively.

assert 50 > 2
assert 2 < 10
assert 3 <= 4
assert 4 >= 4

On rare occasions you might see a slash through these symbols, to describe not. As in, k is "not greater than" j.

ngt

The and are sometimes used to represent significant inequality. That is, k is an order of magnitude larger than j. Sometimes read "beats", when I say x^k ≫ log(x) what I'm really saying is that "polynomial functions grow an order of magnitude faster than logarithms; in a word, the polynomial beats the logarithm."

orderofmag

conjunction & disjunction

Another use of arrows in logic is conjunction and disjunction . They are analogous to a programmer's AND and OR operators, respectively.

The following shows conjunction , the logical AND.

and

In Python, we just say and. Assuming k is a natural number, the logic implies that k is 3:

lambda k: if (k > 2 and k < 4): assert k == 3, "Exercise: can this error ever be raised?"

Since both sides are equivalent , it also implies the following:

lambda k: if (k == 3): assert (k > 2 and k < 4), "I mean it, think through this exercise."

The down arrow is logical disjunction, like the OR operator.

logic-or

In Python, we have the or keyword. Like and, it is a function that will trade you one bool for two bools.

logical negation

Occasionally, the ¬, ~ and ! symbols are used to represent logical NOT. For example, ¬A is only true if A is false.

Here is a simple example using the not symbol:

negation

An example of how we might interpret this in code:

lambda x, y: if (x != y): assert not x == y, "arrr, buried treasure lost forever. "

Note: The tilde ~ has many different meanings depending on context. For example, row equivalence (matrix theory) or same order of magnitude (discussed in equality).

intervals

Sometimes a function deals with real numbers restricted to some range of values, such a constraint can be represented using an interval

For example we can represent the numbers between zero and one including/not including zero and/or one as:

  • Not including zero or one: interval-opened-left-opened-right
  • Including zero or but not one: interval-closed-left-opened-right
  • Not including zero but including one: interval-opened-left-closed-right
  • Including zero and one: interval-closed-left-closed-right

For example we to indicate that a point x is in the unit cube in 3D we say:

interval-unit-cube

In Python, we have to be sensitive about inclusive vs. exclusive boundaries in generators like range, but you already know that.

if you want to play with infinite lists in Python, learn more about generators

Intervals are used in conjunction with set operations:

  • intersection e.g. interval-intersection
  • union e.g. interval-union
  • difference e.g. interval-difference-1 and interval-difference-2

Integer versions in basic python look like this

# intersection of two int intervals
[x for x in range(3,5) if x in range(4, 6+1)]
# Out: [4]

# Union of two int intervals
[x for x in range(20) if x in range(3, 5) or x in range(4, 6+1)]
# Out: [3, 4, 5, 6]

# Set difference
[x for x in range(3, 5) if x not in range(4, 6+1)]
# Out: [3]

[x for x in range(4, 6+1) if x not in range(3, 5)]
# Out: [5, 6]

Using np.linspace, we can approximate what the real versions would look like.

R = np.linspace(-1, 9, 100)

# intersection of two float intervals
[x for x in R if 3 <= x < 5 and 4 <= x <= 6]

# Union of two float intervals
[x for x in R if 3 <= x < 5 or 4 <= x <= 6]

# set differences of two float intervals. 
[x for x in R if 3 <= x < 5 and not (4 <= x <= 6)]

[x for x in R if 4 <= x <= 6 and not (3 <= x < 5)]

You should definitely run these in repl and try to wrap your head around them.

more...

Like this guide? Suggest some more features or send us a Pull Request!

Contributing

For details on how to contribute, see CONTRIBUTING.md.

License

MIT, see LICENSE.md for details.