ellisonbg

def g(x):
  print x

b = 20

ellisonbg

def f(x):
  print x

a = 10

ellisonbg

ellisonbg

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Exploration of the Cal Poly Library API"
   ]
  },
  {

ellisonbg

"""Functions for reading binary files with the digits of pi stored
as 4-bit ints."""

def bin_file_to_digits(filename, the_type=str):
    raw_bytes = np.fromfile(filename, dtype=np.byte)
    for b in raw_bytes[0:10000]:
        digits = byte_to_digits(b)
        yield the_type(digits[0])
        yield the_type(digits[1])

ellisonbg

from sympy.physics.quantum import *
from sympy.physics.cartesian import *
from sympy.physics.piab import *
from sympy.physics.spin import *

# Basic manipulations with operators, bras, kets and commutators

psi, phi = Ket('psi'), Ket('phi')
psi, phi
A, B, C = Operator('A'), Operator('B'), Operator('C')

ellisonbg

"""Quantum mechanics on graphs.

Todo:

* Implement dirichlet boundary conditions.
* Use node attributes to label boundary points and impose dirichlet boundary
  conditions.
* Develop plotting/viz routines using matplotlib/mayavi2.
* Use node attributes to store the physics coordinates of the point.
* Write a bisect function to put additional nodes between existing ones.

ellisonbg

"""Solve the 1D diffusion equation using CN and finite differences."""
from time import sleep

import numpy as np
import matplotlib.pyplot as plt
import networkx as nx

# The total number of nodes
nnodes = 10
# The total number of times

ellisonbg

from numpy import *
from matplotlib.pyplot import *

from time import sleep

import numpy as np
import networkx as nx

dim = [3,2]

ellisonbg

"""Solve the 1D diffusion equation using CN and finite differences."""
from time import sleep

import numpy as np
import matplotlib.pyplot as plt
import networkx as nx

# The total number of nodes
numx = 5
numy = 5