ZGainsforth

# For Python files.
# Created 2017, Zack Gainsforth
import matplotlib
#matplotlib.use('Qt4Agg')
import matplotlib.pyplot as plt
import numpy as np
from rich.traceback import install; install()
import pdb
#pdb.set_trace()

ZGainsforth

# Created 2014, Zack Gainsforth
import matplotlib
matplotlib.use('Qt4Agg')
import matplotlib.pyplot as plt
import numpy as np
from scipy.interpolate import interp1d
from itertools import chain

def QuickPlot(x, y, linestyle='', xlabel=None, ylabel=None, title=None, legendstrs=None, legendloc=None, savefile=None, boldlevel=1, figax=None, plottype=None, xlim=None, ylim=None, **kwargs):
    """

ZGainsforth

from numba import jit
import numpy as np
from timeit import timeit

x = np.random.rand(100000)

@jit
def myfunc(x):
    for i in range(len(x)-1):
        x[i] **= x[i+1]

ZGainsforth

# Created 2014, Zack Gainsforth
import matplotlib
matplotlib.use('Qt4Agg')
import matplotlib.pyplot as plt
import numpy as np

def Beautify(xlabel=None, ylabel=None, title=None, legendstrs=None, savefile=None, boldlevel=1, figax=None):
    """
    :param xlabel: String (can include TeX) for the label of the x-axis.
    :param ylabel: String (can include TeX) for the label of the y-axis.

ZGainsforth

# coding: utf-8
# Created 2014, Zack Gainsforth

# Definitions:
# CLEN: CLinoENstatite, unit cell a = 9 Å
# OREN: ORthoENstatite, unit cell a = 18 Å
# HALF: Half unit cell, a = 4.5 Å
# PEN: ProtoENstatite, unit cell a = 18 Å
# An OREN -> CLEN transition always has to transform an entire 18 Å cell.  Hence always even numbers of CLEN cells.
# For PEN origin, you can get half-planes.  Don't know why.

ZGainsforth

import numpy as np
import Image
import scipy.ndimage as ndimage

def sobel(img):
    gx = ndimage.filters.convolve(img, np.array([[-1, 0, +1], [-2, 0, +2], [-1, 0, +1]]))
    gy = ndimage.filters.convolve(img, np.array([[+1, +2, +1], [0, 0, 0], [-1, -2, -1]]))
    
    imgout = np.sqrt(gx**2+gy**2)
    

ZGainsforth

View network usage by process.

Make sure to make the terminal window the width of the full screen (~200 chars wide is good). Press 'd' to switch to delta mode and then every second it reports what was sent in the last second. You can use right and left arrow to see all columns if you can't make the terminal window wide enough.

nettop

Find a file under the current directory matching a wildcard.

`

ZGainsforth

from numpy import *
from matplotlib.pyplot import *

# Read in the raw
x = fromfile('Stack 1 - 360s.raw', dtype='int8')

# reshape to x, y, E
x.shape=(300,300,2048)

# Make integrated image and show it.

ZGainsforth

First cell stores the path to a variable:

%%javascript
// This gets the name of the notebook so we can use 	relative paths.
var kernel = IPython.notebook.kernel;
var thename = window.document.getElementById("notebook_name").innerHTML;
var command = "NotebookPath = " + "'"+thename+"'";
kernel.execute(command);

ZGainsforth

%pylab inline
import matplotlib.pylab as pylab
import os
pylab.rcParams['figure.figsize'] = 8, 6  # that's default image size for this interactive session